Botswana gravity reference net

Gravity reference stations for the National Gravity Survey of Botswana have been established at twenty-three sites throughout the country in a net linked to existing bases in South Africa, Kenya and Zambia with an internal accuracy of better than 0.5 gravity units (one gravity unit, gu, equals an acceleration of 10⁻⁶ m.s⁻²). The field procedure and reduction of data are explained and a list is given of the gravity values.     

Absolute gravity measurements in Switzerland: Definition of a base network for geodynamic investigations and for the Swiss fundamental gravity net

Summary Results of two absolute gravity surveys performed in Switzerland between 1978 and 1979 are presented and discussed in the framework of the uplift history of the Swiss Alps. Five absolute stations have been established as a contribution to the Swiss fundamental gravity net as well as to geodynamic investigations on the Alpine uplift. Two sites (Interlaken—Jungfraujoch) form the end points of a calibration line for field gravimeters. The gravity range of this line amounts to 605×10⁻⁵ ms⁻² (=605 mgal). It can be traversed in a relatively short time interval of less than 3 hours. Two other sites (Brig and Chur) are located in the area of the most negative gravity anomalies and highest uplift rates encountered in Switzerland. They serve as reference stations for a more extended gravity net for studying non—periodic secular gravity variations associated with the Alpine uplift. Institut für Geod?sie und Photogrammetrie, ETH-Zürich, Separata No. 13. Institut für Geophysik, ETH-Zürich, Contribution No. 333.     

Combination of temporal gravity variations resulting from superconducting gravimeter (SG) recordings, GRACE satellite observations and global hydrology models

Gravity recovery and climate experiment (GRACE)-derived temporal gravity variations can be resolved within the μgal (10^?8 m/s ²) range, if we restrict the spatial resolution to a half-wavelength of about 1,500 km and the temporal resolution to 1 month. For independent validations, a comparison with ground gravity measurements is of fundamental interest. For this purpose, data from selected superconducting gravimeter (SG) stations forming the Global Geodynamics Project (GGP) network are used. For comparison, GRACE and SG data sets are reduced for the same known gravity effects due to Earth and ocean tides, pole tide and atmosphere. In contrast to GRACE, the SG also measures gravity changes due to load-induced height variations, whereas the satellite-derived models do not contain this effect. For a solid spherical harmonic decomposition of the gravity field, this load effect can be modelled using degree-dependent load Love numbers, and this effect is added to the satellite-derived models. After reduction of the known gravity effects from both data sets, the remaining part can mainly be assumed to represent mass changes in terrestrial water storage. Therefore, gravity variations derived from global hydrological models are applied to verify the SG and GRACE results. Conversely, the hydrology models can be checked by gravity variations determined from GRACE and SG observations. Such a comparison shows quite a good agreement between gravity variation derived from SG, GRACE and hydrology models, which lie within their estimated error limits for most of the studied SG locations. It is shown that the SG gravity variations (point measurements) are representative for a large area within the accuracy, if local gravity effects are removed. The individual discrepancies between SG, GRACE and hydrology models may give hints for further investigations of each data series.     

A model for adjustment of differential gravity measurements with simultaneous gravimeter calibration

 A mathematical model is proposed for adjustment of differential or relative gravity measurements, involving simultaneously instrumental readings, coefficients of the calibration function, and gravity values of selected base stations. Tests were performed with LaCoste and Romberg model G gravimeter measurements for a set of base stations located along a north–south line with 1750 mGal gravity range. This line was linked to nine control stations, where absolute gravity values had been determined by the free-fall method, with an accuracy better than 10 μGal. The model shows good consistence and stability. Results show the possibility of improving the calibration functions of gravimeters, as well as a better estimation of the gravity values, due to the flexibility admitted to the values of the calibration coefficients. Received: 15 November 1999 / Accepted: 31 October 2000     

On computing ellipsoidal harmonics using Jekeli’s renormalization

Gravity data observed on or reduced to the ellipsoid are preferably represented using ellipsoidal harmonics instead of spherical harmonics. Ellipsoidal harmonics, however, are difficult to use in practice because the computation of the associated Legendre functions of the second kind that occur in the ellipsoidal harmonic expansions is not straightforward. Jekeli’s renormalization simplifies the computation of the associated Legendre functions. We extended the direct computation of these functions—as well as that of their ratio—up to the second derivatives and minimized the number of required recurrences by a suitable hypergeometric transformation. Compared with the original Jekeli’s renormalization the associated Legendre differential equation is fulfilled up to much higher degrees and orders for our optimized recurrences. The derived functions were tested by comparing functionals of the gravitational potential computed with both ellipsoidal and spherical harmonic syntheses. As an input, the high resolution global gravity field model EGM2008 was used. The relative agreement we found between the results of ellipsoidal and spherical syntheses is 10^?14, 10^?12 and 10^?8 for the potential and its first and second derivatives, respectively. Using the original renormalization, this agreement is 10^?12, 10^?8 and 10^?5, respectively. In addition, our optimized recurrences require less computation time as the number of required terms for the hypergeometric functions is less.     

First results with the transportable absolute gravity meter JILAG-3

A transportable absolute gravity meter has been built at the Joint Institute for Laboratory Astrophysics (JILA) of the University of Colorado and delivered in January 1986 to the Institut für Erdmessung (IFE), Universität Hannover. Instrumental investigations, software improvements and absolute gravity determinations on 15 stations performed byIFE in 1986 are reported here. The main conclusion from one year experience with the instrument is that absolute gravity can be observed on a station within one day with a precision of a few μgal and an accuracy of about 10 μgal.     

The Gauss–Listing geopotential value W 0 and its rate from altimetric mean sea level and GRACE

In geopotential space, the fundamental geodetic parameter W ₀ defines the Gauss–Listing geoid which can be used to best represent the Earth’s mean sea level (MSL) and hence specifies a conventional zero height level to unify vertical datums employed by mapping agencies throughout the world. Further, W ₀ cannot be considered invariant as the parameter varies temporally as a direct response to sea level change and mass redistributions. This study determines W ₀ and its rate, dW ₀/dt, by utilizing altimetric MSL models and an independent mean dynamic topography (MDT) model to define points on the geoid. W ₀ and dW ₀/dt are estimated by two approaches: (i) by means of a global gravity field model (GGM) and (ii) within normal gravity field space as the geopotential value of the best fitting reference ellipsoid. The study shows that uncertainty in W ₀ is mainly influenced by MDT while the choice of methodology, GGM and MSL data coverage are not significant within reason. Our estimate W ₀ =?62636854.2 ± 0.2 m²?s^?2 at epoch 2005.0 differs by 1.8?m²s^?2 from the International Astronomical Union reference value. This study shows that, at a sub-decadal time scale, the time variation dW ₀/dt stems mainly from sea level change with negligible effect from gravity field variations. dW ₀/dt =?(?2.70 ± 0.03)?×?10^?2?m²?s^?2?year^?1, corresponding to a MSL rise of 2.9?mm?year^?1, is evaluated from sea level change based on 16?years of TOPEX and Jason-1 data.     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

Spatially-integrated estimates of net ecosystem exchange and methane fluxes from Canadian peatlands

Background

Peatlands are an important component of Canada’s landscape, however there is little information on their national-scale net emissions of carbon dioxide [Net Ecosystem Exchange (NEE)] and methane (CH₄). This study compiled results for peatland NEE and CH₄ emissions from chamber and eddy covariance studies across Canada. The data were summarized by bog, poor fen and rich-intermediate fen categories for the seven major peatland containing terrestrial ecozones (Atlantic Maritime, Mixedwood Plains, Boreal Shield, Boreal Plains, Hudson Plains, Taiga Shield, Taiga Plains) that comprise >?96% of all peatlands nationally. Reports of multiple years of data from a single site were averaged and different microforms (e.g., hummock or hollow) within these peatland types were kept separate. A new peatlands map was created from forest composition and structure information that distinguishes bog from rich and poor fen. National Forest Inventory k-NN forest structure maps, bioclimatic variables (mean diurnal range and seasonality of temperatures) and ground surface slope were used to construct the new map. The Earth Observation for Sustainable Development map of wetlands was used to identify open peatlands with minor tree cover.

Results

The new map was combined with averages of observed NEE and CH₄ emissions to estimate a growing season integrated NEE (±?SE) at ??108.8 (±?41.3) Mt CO₂ season^?1 and CH₄ emission at 4.1 (±?1.5) Mt CH₄ season^?1 for the seven ecozones. Converting CH₄ to CO₂ equivalent (CO₂e; Global Warming Potential of 25 over 100 years) resulted in a total net sink of ??7.0 (±?77.6) Mt CO₂e season^?1 for Canada. Boreal Plains peatlands contributed most to the NEE sink due to high CO₂ uptake rates and large peatland areas, while Boreal Shield peatlands contributed most to CH₄ emissions due to moderate emission rates and large peatland areas. Assuming a winter CO₂ emission of 0.9 g CO₂ m^?2 day^?1 creates an annual CO₂ source (24.2 Mt CO₂ year^?1) and assuming a winter CH₄ emission of 7 mg CH₄ m^?2 day^?1 inflates the total net source to 151.8 Mt CO₂e year^?1.

Conclusions

This analysis improves upon previous basic, aspatial estimates and discusses the potential sources of the high uncertainty in spatially integrated fluxes, indicating a need for continued monitoring and refined maps of peatland distribution for national carbon and greenhouse gas flux estimation.

     

Applications of an orbiting gravity gradiometer

Considering present attempts to develop a gradiometer with an accuracy between 10⁻³ E and 10⁻⁴ E, two applications for such a device have been studied: (a) mapping the gravitational field of the Earth, and (b) estimating the geocentric distance of a satellite carrying the instrument. Given a certain power spectrum for the signal and 10⁻⁴ E (rms) of white measurement noise, the results of an error analysis indicate that a six-month mission in polar orbit at a height of 200 km, with samples taken every three seconds, should provide data for estimating the spherical harmonic potential coefficients up to degree and order 300 with less than 50% error, and improve the coefficients through degree 30 by up to four orders of magnitude compared to existing models. A simulation study based on numerical orbit integrations suggests that a simple adjustment of the initial conditions based on gradiometer data could produce orbits where the geocentric distance is accurate to 10 cm or better, provided the orbits are 2000 km high and some improvement in the gravity field up to degree 30 is first achieved. In this sense, the gravity-mapping capability of the gradiometer complements its use in orbit refinement. This idea can be of use in determining orbits for satellite altimetry. Furthermore, by tracking the gradiometer-carrying spacecraft when it passes nearly above a terrestrial station, the geocentric distance of this station can also be estimated to about one decimeter accuracy. This principle could be used in combination with VLBI and other modern methods to set up a world-wide 3-D network of high accuracy.     

Unification of vertical datums by GPS and gravimetric geoid models with application to Fennoscandia

The second Baltic Sea Level (BSL) GPS campaign was run for one week in June 1993. Data from 35 tide gauge sites and five fiducial stations were analysed, for three fiducial stations (Onsala, Mets?hovi and Wettzell) fixed at the ITRF93 system. On a time-scale of 5 days, precision was several parts in 10⁹ for the horizontal and vertical components. Accuracies were about 1 cm in comparison with the International GPS Geodynamical Service (IGS) coordinates in three directions. To connect the Swedish and the Finnish height systems, our numerical application utilises three approaches: a rigorous approach, a bias fit and a three-parameter fit. The results between the Swedish RH70 and the Finnish N 60 systems are estimated to −19.3 ± 6.5, −17 ± 6 and −15 ± 6 cm, respectively, by the three approaches. The results of the three indirect methods are in an agreement with those of a direct approach from levelling and gravity measurements. Received: 3 April 1996 / Accepted: 4 August 1997     

The deflection of the vertical in the Western and Central Mediterranean area

The deflection of the vertical at 23 stations was determined by gravimetrical methods and compared to the difference of astronomic observations and geodetic values of the European adjustment. The results show a fair agreement between the two independent sets of values, indicating that doubts about the relative reliability of submarine gravity observations are not justified. A systematic discrepancy in the meridian components of the deflection may reflect on theHayford ellipsoid as used in the European adjustment.

---

Zusammenfassung Auf 23 Stationen wurde die Lotabweichung gravimetrisch bestimmt und den astro-geod?tischen Werten der europ?ischen Ausgleichung gegenübergestellt. Die beiden Ergebnisse stimmen sehr gut überein. Damit werden die Bedenken hinf?llig, die sich gegen den Wert von Messungen in Unterseebooten richteten. Systematische Abweichungen in den meridionalen Komponenten beruhen m?glicherweise auf der art der Anwendung des Hayfordschen Ellipsoids in der Ausgleichung des Zentraleurop?ischen Netzes.

---

Resumen La desviación de la vertical en 23 estaciones ha sido determinada por el método gravimétrico y compara con los resultados astrogeodéticos de la compensación europea. Existe muy buen acuerdo entre los dos conjuntos de datos, que demuestra que las dudas emitidas en cuanto al valor de las observaciones en submarino son injustificadas. Una discordancia sistemática en las componentes meridianas podria ser debida al empleo del elipsoide deHayford en la compensación de la red europea.

---

Résumé La déviation de la verticale en 23 stations a été déterminée par la méthode gravimétrique et comparée aux résultats astrogéodésiques de la compensation européenne. Il y a très bon accord entre les 2 ensembles de données, ce qui montre que les doutes émis quant à la valeur des observations en sous-marins sont injustifiés. Une discordance systématique dans les composantes méridiennes pourrait être due à la fa?on dont on a employé l'ellipso?de deHayford dans la compensation du réseau central européen.

---

Sommario E'stata determinata col metodo gravimetrico la deviazione della verticale in 23 stazioni e comparata ai risultati astro-geodetici della compensazione europea. Vi é molto buon accordo tra i due insiemi di dati, ciò che dimostra che i dubbi emessi circa il valore delle osservazioni in sottomarino sono ingiustificati. Una discordanza sistematica tra le componenti maridiane potrebbe essere dovuta all'impiego nella compensazione della rete central europea dell'ellisoide diHayford.

     

Estimating vertical plant area density profile and growth parameters of a wheat canopy at different growth stages using three-dimensional portable lidar imaging

Vertical plant area density profiles of wheat (Triticum aestivum L.) canopy at different growth stages (tillering, stem elongation, flowering, and ripening stages) were estimated using high-resolution portable scanning lidar based on the voxel-based canopy profiling method. The canopy was scanned three-dimensionally by laser beams emitted from several measuring points surrounding the canopy. At the ripening stage, the central azimuth angle was inclined about 23^° to the row direction to avoid obstruction of the beam into the lower canopy by the upper part. Plant area density profiles were estimated, with root mean square errors of 0.28–0.79 m² m^?3 at each growth stage and of 0.45 m² m^?3 across all growth stages. Plant area index was also estimated, with absolute errors of 4.7%–7.7% at each growth stage and of 6.1% across all growth stages. Based on lidar-derived plant area density, the area of each type of organ (stem, leaves, ears) per unit ground area was related to the actual dry weight of each organ type, and regression equations were obtained. The standard errors of the equations were 4.1 g m^?2 for ears and 26.6 g m^?2 for stems and leaves. Based on these equations, the estimated total dry weight was from 63.3 to 279.4 g m^?2 for ears and from 35.8 to 375.3 g m^?2 for stems and leaves across the growth stages. Based on the estimated dry weight at ripening and the ratio of carbon to dry weight in wheat plants, the carbon stocks were 76.3 g C m^?2 for grain, 225.0 g C m^?2 for aboveground residue, and 301.3 g C m^?2 for all aboveground organs.     

Correcting superconducting gravity time-series using rainfall modelling at the Vienna and Membach stations and application to Earth tide analysis

We demonstrate the possibility to improve the signal-to-noise ratio of superconducting gravity time-series by correcting for local hydrological effects. Short-term atmospheric events associated with heavy rain induce step-like gravity signals that deteriorate the frequency spectrum estimates. Based on 4D modeling constrained by high temporal resolution rain gauge data, rainfall admittances for the Vienna and Membach superconducting gravity stations are calculated. This allows routine correction for Newtonian rain water effects, which reduces the standard deviation of residuals after tidal parameter adjustment by 10%. It also improves the correction of steps of instrumental origin when they coincide with step-like water mass signals.     

我国大地测量学的进展和展望

回顾了我国大地测量工作的进展。面向 2 1世纪前期的我国经济和国防建设及科技和社会发展 ,展望了我国新世纪的大地测量 ,提出应逐步进入精确、动态、实时的现代化体系 ,即完善国家三维空间大地网 ;建立 GPS综合服务体系 ;提供导航和定位服务 ;测定地壳运动、电离层参数、大气中可降水份等信息 ;精化中国地区重力场参数 ;建立新的国家重力基准网 ;完成分米级精度的中国似大地水准面的推算 ;积极开展海洋和空间大地测量 ,为资源、环境的管理以及防灾监测做出应有的贡献。     

Gravity field approximation in the area of Greece with emphasis on local characteristics

The accuracy of the gravity field approximation depends on the amount of the available data and their distribution as well as on the variation of the gravity field. The variation of the gravity field in the Greek mainland, which is the test area in this study, is very high (the variance of point free air gravity anomalies is 3191.5mgal ²). Among well known reductions used to smooth the gravity field, the complete isostatic reduction causes the best possible smoothing, however remain strong local anomalies which disturb the homogeneity of the gravity field in this area. The prediction of free air gravity anomalies using least squares collocation and regional covariance function is obtained within a ±4 ... ±19mgal accuracy depending on the local peculiarities of the free air gravity field. By taking into account the topography and its isostatic compensation with the usual remove-restore technique, the accuracy of the prediction mentioned obove was increased by about a factor of 4 and the prediction results become quite insensitive to the covariance function used (local or regional). But when predicting geoidal heights, in spite of using the smoothed field, the prediction results remain still depend on the covariance function used in such a way that differences up to about 50cm/100km result between relative geoidal heights computed with regional or local covariance functions.     

A preliminary error analysis of the gravity field recovery from a lunar Satellite-to-Satellite mission

Summary A low cost lunar Satellite-to-Satellite radio tracking mission in a low-low configuration could considerably improve the existing knowledge about the lunar gravity field. The impact of various mission parameters that may contribute to the recovery of the gravity field, such as satellite altitude, satellite separation, mission duration, measurement precision and sampling interval were quantified using the Jekeli-Rapp algorithm. Preliminary results indicate that the gravity field resolution up to harmonic degree 40 to 80 is feasible depending on various mission configurations. Radio tracking data from a six-month mission with a precision of 1 mm s^–1 every 10 s and 300 km satellite separation at 150 km altitude will permit the determination of 5^o×5^o mean gravity anomalies with an error of approximately 15 mgals. Consideration of other unaccounted error sources of instrumental, operational as well as environmental nature may lower this resolution.     

Accuracy analysis for the NAD 83 geodetic reference system

The North American Datum of 1983 (NAD 83) provides horizontal coordinates for more than 250,000 geodetic stations. These coordinates were derived by a least squares adjustment of existing terrestrial and space-based geodetic data. For pairs of first order stations with interstation distances between 10km and 100km, therms discrepancy between distances derived fromNAD 83 coordinates and distances derived from independentGPS data may be suitably approximated by the empirical rulee=0.008 K^0.7 where e denotes therms discrepancy in meters and K denotes interstation distance in kilometers. For the same station pairs, therms discrepancy in azimuth may be approximated by the empirical rule e=0.020 K^0.5. Similar formulas characterize therms discrepancies for pairs involving second and third order stations. Distance and orientation accuracies, moreover, are well within adopted standards. While these expressions indicate that the magnitudes of relative positional accuracies depend on station order, absolute positional accuracies are similar in magnitude for first, second, and third order stations. Adjustment residuals reveal a few local problems with theNAD 83 coordinates and with the weights assigned to certain classes of observations.     

Estimation of Partial Pressure of Carbon Dioxide and Air-Sea Fluxes in Hooghly Estuary Based on In Situ and Satellite Observations

An empirical model is developed and used with remotely sensed predictors: sea surface temperature (SST) and chlorophyll-a concentration (Chl-a), to compute surface water partial pressure of carbon dioxide (pCO_2w) and air-sea fluxes of CO₂ in the Hooghly estuary and its adjacent coastal oceans. In situ observations used here were based on measurements carried out in this region during winter and summer periods in 2008. The estimated pCO_2w compares well with the in situ observations at root mean square error ±18 μatm. In winter, estimated pCO_2w ranges between 320 and 500 μatm with large values (>400 μatm) on the south-western and south-eastern flanks of the coastal domain and lower values (340–375 μatm) on the main-channel. In summer, it remained spatially uniform at 450 μatm. Extrapolation of the results over the study region based on the Moderate Imaging Specroradiometer (MODIS) measured SST and Chl-a suggests that the region is a strong source of atmospheric CO₂ during the summer with net release of 0.095 Tg C year^?1 (equivalent to mean flux of 90 molC m^?2 year^?1) and is a weak source during the winter with net release of 0.006 Tg C yr^?1 (0.5 molC m^?2 year^?1) from the geographical extent of 6000 Km² area.     

应用抗差估计理论分析2000国家重力基本网

应用抗差估计理论，结合2000国家重力基本网的实际数据。通过大量的试算对2000国家重力基本网的数据进行分析比较，提出了符合2000国家重力基本网平差计算的抗差估计方法。