Boedecker, G.

Francis, Olivier

Kenyeres, A.

Open Repository and Bibliography - Luxembourg

Unified European Gravity Reference Network 2002 (UEGN02) - Status 2004 G. Boedecker Bavarian Academy of Sciences and Humanities / BEK, Miinchen, Germany. Fax +49 8923031100 boe@bek.badw-muenchen.de O. Francis . European Centre for Geodynamics and Seismology / University of Luxembourg,Luxembourg, Fax +352 33 14 87 88, olivierfglecgs.lu A. Kenyeres FOMI Satellite Geodetic Observatory, Penc, Hungary; Fax +36 27 374982, kenyeres@sgo.fomi.hu Abstract. The last previous realisation of the Unified European Gravity Reference Network was completed in 1994 (UEGN94) covering 11 West European countries. Since that time significant work has been done e.g. in the frame of the UNIGRACE project that focused on the establishment of absolute gravity stations in Central and East European countries. The past European Subcommission of the IGGC recommended to continue with the unification of the gravity reference networks at the continental level. During the IAG2001 Scientific Assembly in Budapest, it was decided to establish a new continental gravimetric reference network. The total UEGN2002 network area covers 25 countries, the number of stations expected was around 1000, but has already exeeded this number. The data processing starts from raw absolute and relative gravity observations. The tidal corrections are computed by the ECGS in Luxembourg in a uniform manner and based on state-of-the-art models. At the current status (2004) the work is focused on 12 countries with 404 absolute and 33000 relative gravity observations at some 1500 stations. Main problems are data bugs. Hence, appropriate tools had to be developed. The paper presents details of the adjustment, earth tide corrections as well as analyses of the data and adjusted gravity values. After this partial UEGN02 adjustment, the remaining countries that also showed their interest to participate are again invited to submit data. Keywords. Reference network, gravity reference 1 Introduction 1.1 Role of Gravity Reference Networks Because bulk field gravimetry still is carried out by relative spring gravimeters, it is necessary to link these to an absolute reference. Decades ago, only a few absolute observations worldwide provided the basis for extended relative gravity reference net-works. Today, as absolute gravity meters become more easily available, they become more fieldwor-thy and portable, also their number is increasing. Consequently, the number of absolute stations is increasing and the contribution of relative meters to reference networks is decreasing. In each case where a new national reference or an extended field campaign is envisaged, the benefit-cost-relation will lead to a new optimum design appropriate to the specific situation, see e.g. Boedecker (2002). Rela-tive observations are not only important for refer-ence networks, their use in combined absolute / relative networks will also benefit in their calibra-tion. 1.2 European cooperation / Preceding works After the old global gravity reference network IGSN71 (Morelli et al. 1971) proved to be not a sufficiently accurate basis for new gravity observa-tions any more, many countries designed and ob-served national networks. In 1994, an attempt was made on behalf of the IGC Subcommission Western Europe to unify a number of European national networks (Boedecker, Marson, Wenzel 1994) and named UEGN94 ('Unified European Gravity Net-works 1994'). The network covered 11 countries 286 comprising 499 stations with 123 absolute and 14532 relative gravity observations; see fig. 1. From 1998 to 2002 the project UNIGRACE ('Uni-Figure 1: UEGN94 fication of Gravity Systems in Central and Eastern Europe') was carried out, coordinated by the Ger-man BKG ('Bundesamt fiir Kartographie und Geodasie') with contractors from all participating countries, see fig. 2. See Reinhart, Richter, Wilmes 1998. In the framework of that project, absolute observations at 19 stations in 12 countries in central and eastern Europe have been observed and were made available to the national survey agencies and to the UEGN project. Figure 2: UNIGRACE These works for a unification of European gravity reference networks were continued by the UEGN02. A first meeting to agree on the guide-lines took place on May 13-14, 2001, in Vienna. The work was delayed because of problems with details of the ties between absolute and relative observation - a key problem for the correct gravity level of the network. Some further delay is caused by incorrect data: As a matter of fact, this is a bigger problem than in other much larger data sets, because in this case we are dealing with (currently) 33000 observations each of which has its own evolution and is not the result of one unique process as is the case in other type of geodetc observations like GPS observation series or similar. For this reason, some software tools had to be developed in order to make the observation series and procedures more transparent. The current status of the UEGN02 activities is illustrated in figure 3: Figure 3: UEGN02: Activity status 08/04 dark: adjustment computation, data cleaning medium: some data available light: data envisaged 2 Data collection The data were collected in four types of formatted files: station data, absolute observations, relative observations. Some national agencies do not agree to have their data published. For this reason, no original data set will be transferred without permis-sion of the original owner. 3 Tidal reduction All the raw observations of the relative gravity meters have been corrected for the body Earth tides and ocean loading and attraction effect using state-ofthe art models. The tidal parameters for the body Earth tides come from the Dehant-Defraigne-Wahr (1999) model for an inelastic non-hydrostatic Earth (Dehant et al., 1999) including a delta factor of 1.16 for the long periods tides. The body tides prediction 287 is computed with the etgtab software written by Wenzel (http://www-geod.kugi.kvoto-u.ac.ip/iagetc/etcdat/etgtab/etgtab .txt On the other hand, the ocean loading and attraction parameters were calculated using the ocean tides model CSR3.0 and the Green's function of the PREM with the Load89 software (Francis O. and P. Mazzega, 1993). The oceanic loading tidal predic-tion is then calculated using the ocean loading pa-rameters in a separate program. The tidal correction is the sum of the contributions: body Earth tides and oceanic loading and attraction gravity effect. The data from the absolute gravimeters were pro-vided already corrected for tides. It means that the tides correction for these data could be slightly different. It was outside the scope of this work to reprocess all the absolute gravity data from the raw observations.). 4 Adjustment model The adjustment model is identical to the adjustment of UEGN94. Absolute gravity observations are immediate observations of station gravity. Absolute gravity observations : ^i=gj-gi Vj residual of observation i g j unknown station gravity gj observed gravity incl. tidal reduction Relative gravity meter obervations are taken as a series of consecutive observations within one day by one gravity meter. A series is modeled by station gravity values, an offset (per series), linear drift of the readings with time (per series) and a linear calibration factor for the entire set of observations of the meter. Relative observations are the raw observations - corrected for e.g.periodic errors, if necessary or available, converted by the manu-facturers calibration table to care for nonlinear com-ponents, and reduced by tidal effects (see below). Relative gravity observation : ^i =gj -^Oi+Zifg+tjCii -ri Sj oi g residual of observation i unknown gravity of station j orientation of series 1 raw gravimeter reading i linear calibration factor of gravimeter g time of observation i linear drift for series 1 gravimeter reading i incl. tides & pre - calib. The intention was to use the original observations at definition height of the absolute observation and to connect to all other stations including the immediate ground marker by relative observations. This would have avoided introducing prior 'gradient corrections' etc. However, frequently the original relative gravity observations were not carried out at the height of the absolute observations; hence, the ties between absolute observations and the network is a mix of immediate relative observations and intermediate use of gradients. To account for the variations in accuracy of the absolute observations, the internal error variances estimations were introduced, as communicated by the agencies providing the data. The weights of the relative observations were assumed constant within one series and the series weights were estimated from an iterative variance component estimation. Initial weights : Or, zero variance O, variance of series 1 Weight iteration: w / + e m zero variance ^ 0 [vvJi squared sum of residuals of series 1 ni number of obsns. of series 1 £ small number 5 Execution of adjustment The adjustment programme was developed in MATLAB. The challenges of this work are in the data organization etc. E.g., relative observation series have to be identified from the contiguous observation set, station data checked; series of only two observations do not provide network information and are removed. The data were processed to a certain point, then metadatafiles were transferred to ECS Luxembourg (O. Francis) for tidal reduction. After this, the processing continued. Variants of the adjustment are based on the same data set, hence the tidal 288 reduction computation was necessary only once so far. Because of the restrictions of this paper format, it is not possible to provide a full and detailed image of the screen for the practical work on the network. Rather, a sequence of cutout figures tries to provide at least a rough impression, see appendix. The network figures (see appendix) show e.g. normal stations, absolute stations, stations with observations by parties from other countries, station names, station gravity standard deviation estimates, absolute and relative residuals above some threshold (arrows), observation series with residuals above some threshold (red dotted line). The graphical user interface is not just a nice toy, but it is essential because otherwise a clarification of bugs in the large data set of individual observations would not be possible. A standard programme run including all original data file input and output, graphics, parameter esti-mation for about 10.000 parameters, iterative co-variance component estimates in six iterations, but without tidal reduction and inversion of the normals takes less than 1 minute on a normal PC under Windows. The inversion for computing the standard deviations of the parameters takes between half a minute to several minutes depending on the quality of the normal matrix. 1.4 1.2 1 •E 0.8 .5* '(U 5 0.6 0.4 0.2 Q ^ ^_ ^^ a a Series weight estimation iteration mtSStt . , „ «. ^w-- <_»• "*" gas _ « , s s ^ = - -^ _ ^„^ ' init ' it=1 it=2 . it=3 ' it=6 . _ ——————— ***" """^  « B B S = 2.015 2.02 2.025 2.03 2.035 #obsn 2.04 x 1 0 ' Figure 4: Iterative weight estimation Further explanation see text The weight estimation for relative observations starts with unit weight (black in the sample figure 4), except the series is downweighted to zero as a result of a preceding bug analysis (see figure). The iteration subsequently yields the weights depicted above. As can be recognized, the convergence is quite stable and fast. The key numbers of the sample run are given be-low. The histogram of residuals shows that 50% of the residuals are less than 0.007 mGal. Absolute obsns: 405 Relative obsns: 33092 Rel. gravimeters 44 Series 3989 Stations: 1446 Paranneters: 9445 14000 12000 10000 «, 8000 O =«= 6000 4000 2000 0 Histogram of residuals Hl^ _ --' ---0.1 -0.05 0 mGal 0.05 0.1 Figure 5: Histogram of residuals 6 Conclusions and outlook • A subnet of the UEGN02 has successfully been processed • The national network structures are very di-verse • Observation data bugs are a major problem • This project completes the UNIGRACE objec-tives • The accuracy is satisfactory • The benefits are very manifold • All European countries are called to participate • Some of the countries contributing already are called to deliberate to submit more complete data in order to have greater benefit 289 • It is envisaged to complete the project in 2005 • The further procedure, e.g. publication (also stations ?) has to be agreed by the participating countries 7 Acknowledgements This work has been made possible through the con-tributions of many colleagues in institutes and agencies in the participating countries, Austria, Czech Republic, Croatia, Estonia, Finland, France, Germany, Greece, Hungary, Italy, Latvia, Lithua-nia, Netherlands, Poland, Slovakia, Slovenia, Spain, Switzerland, Ukraine. Others offered there active participation in computations such as Volgyesi and Csapo from Hungary, Simek from Czech Republic and Klobusiak from Slovakia; these offers are ap-preciated. References Boedecker, G.: World gravity standards - present status and friture challenges. Metrologia 2002, 39 no.5, 429-434 Boedecker, G.; Marson, I.; Wenzel, H.-G. : The Adjustment of the Unified European Grav-ity Network 1994 (UEGN 94). In: Siinkel, H. / Marson, I. (eds.): Gravity and Geoid. Joint IGC/ICG Symposium Graz 1994. lAG Sympo-sium No. 113, Springer 1995, 82-91 Dehant, V., P. Defraigne and J. M. Wahr, Tides for a convective Erath, Journal of Geophysical Research, 104, 1035-1058,1999. Francis O. and Mazzega P., Global charts of ocean tide loading effects. Journal of Geophysical Re-search, 95(C7), 11411-11424, 1990. Morelli, C , Gantor, C, Honkasalo, T., McConnell, R.K., Tanner, J.G., Szabo, B., Votila, V. & Whalen, C.T.-1971- The International Gravity Standardisation Net. Pub. Spec. no. 4 du Bulletin Geodesique., 1974 Reinhart, E.; Richter, B.; Wilmes, H.: "UNIGRACE - Ein Projekt zur Vereinheitlichung der Schwere-Referenzsysteme in Zentral-Europa" Mitteilungen des Bundesamtes fur Kartographie und Geodasie, Bandl, Frankfurt a. M. 1998; see also: INCO: International Scientific Cooperation Projects (1998-2002) Contract Number: ERBIC15CT970805 FINAL REPORT ngrJorn DeMirnEir iiiiiii^mH g^T)/ ItaK/ % f * W T + % , . 2 -Appendix: Impres-sions of Network-ing: 290 UEGN results 2004 8 19 19 12 Stations files read: D:\UEGN02\at\stations_at_03.txt4 D:\UEGN02\cz\stations cz 02.txt4 Appendix: Sample output Absolute observations files read: D:\UEGN02\at\absobs_at_03.txt 5 D:\UEGN02\cz\absobs_,cz__02.txt 5 Relative observations files read: D:\UEGN02\at\relobs_at_4b.txt 5 D:\UEGN02\cz\relobs cz 2a.txt 5 Stations with identical coordinates 3D atO-059-20 atA-059-20 CZ107.10 hu107.10 CZ107.10 sk107.10 czl 08.00 hu108 ; - - -Parameters not solved for: Station Gravity: i 794 gr299IRAFI 1 797gr309KYMH { Observations per instrument: # instr, # obsns 1 -900.002 1008 2 -625.001 110 3 -90.002 178 4 3.004 302 5 4.001 26 6 51.002 1934 Parameters : names (seq: i atO-021-00 atO-021-01 i atO-050-00 atO-058-20 180.041 191.001 192.001 192.05 193.001 . . . grav, lineal orient, drift), adjusted values, 980867.380 0.002 980868.928 0.003 980741.267 0.002 980842.055 0.003 -0.089 6.083 -3.727 2.749 -23.301 i Observations per stations: / # rel obsns, # abs obsns ] atO-021-00 65 0 ; atO-021-01 17 0 ! Iiu87 7 0 hu88 0 2 hu88.0 54 2 iiu89 0 1 : Absolute outliers: 0.000 0.001 0.001 0.006 0.001 1 \ hu88 930826 980765.513000 0.048 AXIS-FG5-107 0.002 Gschwind hu88 971213 980765.610000 -0.049 JILAG-6 0.002 Ruess mse: Bad series: Bad series: mse, maxr: 0.057230 0.060223 (StatName, date, time, instr, raw, obs(inc cal/tid), tides, resid, status) atll-SCHLb 010929 1837 -900.002 119.817 134.397 -0.056 +0.0310 at2-171-20 010929 1946 -900.002 17.798 19.951 -0.032-0.013 0 at2-171-21 010929 2012 -900.002 42.309 47.471 -0.023-0.060 0 atll-SCHLb 010929 2114 -900.002 119.777 134.402 -0.006 +0.042 0 291 

Unified European Gravity Reference Network 2002 (UEGN02)—Status 2004

http://orbilu.uni.lu/bitstream/10993/2518/1/10.1007_3-540-26932-0_50.pdf

Unified European Gravity Reference Network 2002 (UEGN02)—Status 2004

Abstract

Similar works

Full text

Available Versions

Open Repository and Bibliography - Luxembourg