Repeat station crustal biases and accuracy determined from regional field models

1. Repeat station crustal biases and accuracy determined from regional field models M. Korte, E. Thébault* and M. Mandea, GeoForschungsZentrum Potsdam(*now at Institut de Physique du Globe, Paris)

2. Outline • The German repeat station network and data processing • R-SCHA regional magnetic anomaly modelling technique • Preparation of data- satellite data- aeromagnetic data- repeat station data • Repeat station uncertainties • Repeat station crustal biases • Conclusions blue: this talk black: details in talk by E. Thébault

3. The German repeat station network • Ca. 45 repeat stations, subset of denserground vector surveys established in thepast. • Since 2000 two classes of stations: • Variometer stations Portable LEMI variometer installed nearby, 3 to 7 days • Other stations surveyed while variometer is recording at nearest variometer station

4. Data processing I • Determine “local baseline” of nearest variometer from absolute measurements (like standard observatory practice)=> Series of full field vector variations for a few days for each station • Use quiet night times with constant differences between recordings at the station and observatory for reduction to “annual means”. C(xi,tmean) = C(xi,ti) – C(O,ti) + C(O, tmean) Repeat station“annual mean”of component C Repeat stationmeasurementvalue at time ti Observatoryannual meanof component C Observatoryrecording attime ti This difference determined robustlyfrom quiet night time values

5. Data processing II Differences between recordings of the LEMI variometer at station EBH to NGK observatory recordings for 3 days Quiet hours used for final datareduction to annual means

6. R-SCHA modellingRevised spherical cap harmonic analysis • Method developed by E. Thébault, based on SCHA (spherical cap harmonic analysis), a “regional spherical harmonic analysis” • Data from different altitudes can be taken into account because the vector field is modelled inside a cone • Well-suited for joint inversion of ground, aeromagnetic and satellite data • Ideally suited for regional models of the vector crustal magnetic anomaly field (first example: France by Thébault et al.) • Details given in presentation by Erwan Thébault This work: study crustal field influence at repeat station locations (“crustal biases”) and use compatibility of different data types in modelling for estimations of data uncertainty

7. Data for R-SCHA anomaly models Crustal anomaly data, different data types: • Satellite data+ full vector information+ long wavelength information- low spatial resolution (no better than ~300km) due to altitude of satellite • Aeromagnetic data+ high spatial resolution (up to km scale)- only intensity data- long wavelength information missing • Ground data, repeat station data+ full vector information- very localized information, even ground vector surveys mostly notdense enough for detailed crustal field mapping => Only combination of all data types give full information

8. Data preparation • Satellite data- selection of quiet night time data - corrections for external fields- subtraction of core field model up to SH degree 14 • Aeromagnetic data- existing anomaly compilation for Germany (grid) used(assumption: external variation corrections, upward-/downward continuation for individual surveys to common altitude, reductions of individual surveys to common date have been done properly before)- addition of the originally subtracted DGRF 1980 model- subtraction of core field model up to SH degree 14 • Repeat station data- Results from one year would be enough, but use as many data as possible for higher confidence and to get an idea of uncertainties

9. Repeat station data used in this work • All available ground survey andrepeat station data within the time span for which a continuous core field model is available (CM4 model by Sabaka et al., valid 1960 – 2002) • Only stations, where measurements had been taken at least three times in this time interval (for statistics)

10. Repeat station data preparation • Subtract core field model up to SH degree 14 • Afterwards, observatory annual means show systematic variations, assumed to be external field influences which did not average out. • Correction for these influences based on a simple empirical template, using the homogeneity of these variations in a small area • Average residual is crustal field contribution (assumed to be constant over 40 years) • Standard deviation about average is good estimate of uncertainty • In most cases, standard deviation of average is reduced by the empirical external field correction

11. Empirical external field correction for (repeat station) annual means Left: Observatory annual meansafter subtraction of core field model, black line is best fit linear trend. Right: Residual variations in annualmeans after subtraction of linear trends, black line is average of the three= template for external field correction of all annual means Red: NGK Green: WNG Blue: FUR

12. Repeat station data preparation • Subtract core field model up to SH degree 14 • Afterwards, observatory annual means show systematic variations, assumed to be external field influences which did not average out. • Correction for these influences based on a simple empirical template, using the homogeneity of these variations in a small area • Average residual is crustal field contribution (assumed to be constant over 40 years) • Standard deviation about average is good estimate of uncertainty • In most cases, standard deviation of average is reduced by the empirical external field correction

13. Repeat station crustal biases and uncertainties Avg. 4.7 8.3 4.5 These uncertainties also represent repeat station data uncertainties

14. Repeat station crustal biases and uncertainties Uncertainties are smaller at timesand locations, where a local variometer was used. These uncertainties might containinfluences from insufficient repre-sentation of core field and secularvariation by the subtracted model.Systematic deviationssuggesting a significant influence of such an effect could not be found. However, the data mostly come from measurements without localvariometers, a study of small-scalesecular variation requires highest accuracy data!

15. The R-SCHA model Y-Anomaly X-Anomaly Z-Anomaly

16. Model fit to the data Average rms: Model 1 X 33.8nT Y 22.8nT Z 45.1nTAerom. 13.2nT Model 2 X 23.4nTY 17.0nT Z 26.0nT Aerom. 13.6nT Rms misfit larger than uncertainty estimates,probably due to resolution limit of the model (~35 km)

17. Conclusions • We have obtained estimates of repeat station data uncertainties:on average 4.7 nT for X 8.3 nT for Y 4.5 nT for ZThe misfit of a regional R-SCHA model is significantly larger, but that is probably due to the limited model resolution (~35 km) • The crustal field influence at the German repeat stations reaches up to +/-150 nT in the different components, which has to be considered when the data are used for core field studies(see example on Poster, F in 2004 and 2006) • Repeat station data are useful for joint inversions of ground, satellite and aeromagnetic data for magnetic anomaly mapping.