Correction of Discretization Errors Simulated at Supply Wells

Many hydrogeology problems require predictions of hydraulic heads in a supply well. In most cases, the regional hydraulic response to groundwater withdrawal is best approximated using a numerical model; however, simulated hydraulic heads at supply wells are subject to errors associated with model discretization and well loss. An approach for correcting the simulated head at a pumping node is described here. The approach corrects for errors associated with model discretization and can incorporate the user's knowledge of well loss. The approach is model independent, can be applied to finite difference or finite element models, and allows the numerical model to remain somewhat coarsely discretized and therefore numerically efficient. Because the correction is implemented external to the numerical model, one important benefit of this approach is that a response matrix, reduced model approach can be supported even when nonlinear well loss is considered.     

VLBI-derived troposphere parameters during CONT08

Time-series of zenith wet and total troposphere delays as well as north and east gradients are compared, and zenith total delays (ZTD) are combined on the level of parameter estimates. Input data sets are provided by ten Analysis Centers (ACs) of the International VLBI Service for Geodesy and Astrometry (IVS) for the CONT08 campaign (12?C26 August 2008). The inconsistent usage of meteorological data and models, such as mapping functions, causes systematics among the ACs, and differing parameterizations and constraints add noise to the troposphere parameter estimates. The empirical standard deviation of ZTD among the ACs with regard to an unweighted mean is 4.6?mm. The ratio of the analysis noise to the observation noise assessed by the operator/software impact (OSI) model is about 2.5. These and other effects have to be accounted for to improve the intra-technique combination of VLBI-derived troposphere parameters. While the largest systematics caused by inconsistent usage of meteorological data can be avoided and the application of different mapping functions can be considered by applying empirical corrections, the noise has to be modeled in the stochastic model of intra-technique combination. The application of different stochastic models shows no significant effects on the combined parameters but results in different mean formal errors: the mean formal errors of the combined ZTD are 2.3?mm (unweighted), 4.4?mm (diagonal), 8.6?mm [variance component (VC) estimation], and 8.6?mm (operator/software impact, OSI). On the one hand, the OSI model, i.e. the inclusion of off-diagonal elements in the cofactor-matrix, considers the reapplication of observations yielding a factor of about two for mean formal errors as compared to the diagonal approach. On the other hand, the combination based on VC estimation shows large differences among the VCs and exhibits a comparable scaling of formal errors. Thus, for the combination of troposphere parameters a combination of the two extensions of the stochastic model is recommended.     

EOP and scale from continuous VLBI observing: CONT campaigns to future VGOS networks

Continuous (CONT) VLBI campaigns have been carried out about every 3 years since 2002. The basic idea of these campaigns is to acquire state-of-the-art VLBI data over a continuous time period of about 2 weeks to demonstrate the highest accuracy of which the current VLBI system is capable. In addition, these campaigns support scientific studies such as investigations of high-resolution Earth rotation, reference frame stability, and daily to sub-daily site motions. The size of the CONT networks and the observing data rate have increased steadily since 1994. Performance of these networks based on reference frame scale precision and polar motion/LOD comparison with global navigation satellite system (GNSS) earth orientation parameters (EOP) has been substantially better than the weekly operational R1 and R4 series. The precisions of CONT EOP and scale have improved by more than a factor of two since 2002. Polar motion precision based on the WRMS difference between VLBI and GNSS for the most recent CONT campaigns is at the 30 \(\upmu \)as level, which is comparable to that of GNSS. The CONT campaigns are a natural precursor to the planned future VLBI observing networks, which are expected to observe continuously. We compare the performance of the most recent CONT campaigns in 2011 and 2014 with the expected performance of the future VLBI global observing system network using simulations. These simulations indicate that the expected future precision of scale and EOP will be at least 3 times better than the current CONT precision.     

Monitoring the TOPEX and Jason-1 Microwave Radiometers with GPS and VLBI Wet Zenith Path Delays

Monitoring of altimeter microwave radiometer measurements is necessary in order to identify radiometer drifts or offsets that if uncorrected will introduce systematic errors into ocean height measurements. To examine TOPEX Microwave Radiometer (TMR) and Jason-1 Microwave Radiometer (JMR) behavior, we have used coincident wet zenith delay estimates from Very Long Baseline Interferometry (VLBI) and Global Positioning System (GPS) geodetic sites near altimeter ground tracks. We derived a TMR path delay drift rate of -1.1 ± 0.1 mm/yr using GPS data for the period from 1993.0-1999.0 and -1.2 ± 0.5 mm/yr using VLBI data. Thereafter, the drift appears to have leveled off. Already after 2.3 years (82 cycles) of the Jason-1 mission, it is clear that there have been significant systematic errors in the JMR path delay measurements. From comparison with GPS wet delays, there is an offset of -5.2 ± 0.6 mm at about cycle 30 and a more abrupt offset of -11.5 ± 0.8 mm at cycle 69. If we look at the behavior of the JMR coldest brightness temperatures, we see that the offsets near cycle 30 and cycle 69 are mainly caused by corresponding offsets in the 23.8 GHz channel of -0.49 ± 0.12 K and -1.18 ± 0.13 K, although there is a small 34.0 GHz offset at cycle 69 of 0.75 ± 0.22 K. Drifts in the 18.0 and 34.0 GHz channels produce a small path delay drift of 0.3 ± 0.5 mm/yr.     

Precise geodesy with the Very Long Baseline Array

We report on a program of geodetic measurements between 1994 and 2007 which used the Very Long Baseline Array (VLBA) and up to ten globally distributed antennas. One of the goals of this program was to monitor positions of the array at a 1 mm level of accuracy and to tie the VLBA into the international terrestrial reference frame. We describe the analysis of these data and report several interesting geophysical results including measured station displacements due to crustal motion, earthquakes, and antenna tilt. In terms of both formal errors and observed scatter, these sessions are among the very best geodetic very long baseline interferometry experiments.     

Combination of long time-series of troposphere zenith delays observed by VLBI

Within the International Very Long Baseline Interferometry (VLBI) Service for Geodesy and Astrometry (IVS), long time-series of zenith wet and total troposphere delays have been combined at the level of parameter estimates. The data sets were submitted by eight IVS Analysis Centers (ACs) and cover January 1984 to December 2004. In this paper, the combination method is presented and the time-series submitted by the eight IVS ACs are compared with each other. The combined zenith delays are compared with time-series provided by the International Global Navigation Satellite System (GNSS) Service (IGS), and with zenith delays derived from the European Centre for Medium-Range Weather Forecasts (ECMWF). Before the combination, outliers are eliminated from the individual time-series using the robust BIBER (bounded influence by standardized residuals) estimator. For each station and AC, relative weight factors are obtained by variance component estimation. The mean bias of the IVS ACs’ time-series with respect to the IVS combined time-series is 0.89 mm and the mean root mean square is 7.67 mm. Small differences between stations and ACs can be found, which are due to the inhomogeneous analysis options, different parameterizations, and different treatment of missing in-situ pressure records. Compared to the IGS zenith total delays, the combined IVS series show small positive mean biases and different long-term trends. Zenith wet delays from the ECMWF are used to validate the IVS combined series. Inconsistencies, e.g., long-term inhomogeneity of the in-situ pressure data used for the determination of VLBI zenith delays, are identified.     

Monte Carlo simulations of the impact of troposphere,clock and measurement errors on the repeatability of VLBI positions

Within the International VLBI Service for Geodesy and Astrometry (IVS) Monte Carlo simulations have been carried out to design the next generation VLBI system (“VLBI2010”). Simulated VLBI observables were generated taking into account the three most important stochastic error sources in VLBI, i.e. wet troposphere delay, station clock, and measurement error. Based on realistic physical properties of the troposphere and clocks we ran simulations to investigate the influence of the troposphere on VLBI analyses, and to gain information about the role of clock performance and measurement errors of the receiving system in the process of reaching VLBI2010’s goal of mm position accuracy on a global scale. Our simulations confirm that the wet troposphere delay is the most important of these three error sources. We did not observe significant improvement of geodetic parameters if the clocks were simulated with an Allan standard deviation better than 1 × 10⁻¹⁴ at 50 min and found the impact of measurement errors to be relatively small compared with the impact of the troposphere. Along with simulations to test different network sizes, scheduling strategies, and antenna slew rates these studies were used as a basis for the definition and specification of VLBI2010 antennas and recording system and might also be an example for other space geodetic techniques.     

Protecting window‐size female American lobster,Homarus americanus,to increase egg production

A new conservation measure aimed at protecting the larger broodstock (females with a carapace length between 115 and 129 mm, i.e., window‐size females, WSF) to enhance egg production was implemented in 2003 for the American lobster (Homarus americanus) fishery in the southern Gulf of St. Lawrence, Canada. To assess the effectiveness of this new measure, a tagging project was carried out between 2004 and 2007 in the Hillsborough Bay area with the collaboration of local harvesters. Based on a single mark‐recapture Petersen model, the estimated catch rates of WSF not carrying eggs under the abdomen in the study area varied between 40% and 53%. The catch‐per‐unit‐effort (animals per 100 traps) of both berried (egg‐bearing) and non‐berried WSF increased from 0.2 to 0.8 and 1.0 to 1.6, respectively, between 2004 and 2007. Based on the estimated number of non‐berried WSF in the population observed during the spring fishery between 2004 and 2007, the potential total egg production for WSF was estimated at 1.1 billion in the study area. From this total, 515 million eggs would have been produced by WSF directly protected by the conservation measure, i.e., females that would have been harvested if not for the new regulation. Finally, the contribution of berried WSF to the annual egg production in the study area increased from 4% to 23% between 2003 and 2007 based on the at‐sea sampling programme and the fecundity curve. WSF that represented 10% of the total berried female population in 2007 contributed 23% of the annual egg production. Hence, the protection of WSF seems to be an effective conservation measure to increase egg production based on the catch‐per‐unit‐effort and egg production.