  Astrometric accuracy of complex modern VLBI arrays cannot be calculated
analytically. We study the astrometric accuracy of phase-referenced VLBI
observations for the VLBA, EVN and global VLBI array by simulating VLBI data
for targets at declinations -25$^\circ$, 0$^\circ$, 25$^\circ$, 50$^\circ$,
75$^\circ$ and 85$^\circ$. The systematic error components considered in this
study are calibrator position, station coordinate, Earth orientation and
troposphere parameter uncertainties. We provide complete tables of the
astrometric accuracies of these arrays for a source separation of 1$^\circ$
either along the right ascension axis or along the declination axis.
Astrometric accuracy is 50microas at mid declination and is 300microas at low
(-25$^\circ$) and high (85$^\circ$) declinations for the VLBA and EVN. In
extending our simulations to source separations of 0.5$^\circ$ and 2$^\circ$,
we establish the formula for the astrometric accuracy of the VLBA: Delta =
(Delta_1$^\circ$-14)*d+ 14 (microas) where Delta_1$^\circ$ is the astrometric
accuracy for a separation d=1$^\circ$ provided in our tables for various
declinations and conditions of the wet troposphere. We argue that this formula
is also valid for the astrometric accuracy of the EVN and global VLBI array.
