Impact Estimation and Filtering of Disturbances in FG 5 Absolute Gravimeter Observations

Instrumental and environmental disturbances do affect FG5 absolute gravimeter observations and the estimated gravity values, sometimes to the degree that entire measurement campaigns are discarded. We propose a method which moves towards the re-assessment of previously discarded observations. Once an estimate of the frequency and amplitude of a disturbance in a FG5 data set exists, the proposed method can estimate its impact on the estimated gravity value. This is performed through a Gaussian Bell Summation approach of the functional relationship between disturbance frequency and standard deviation of gravity. The filtering of the identified disturbance is realized through a modification of the functional model of the equation of motion in the least squares adjustment of FG5 observations. The results reveal that the Gaussian Bell Summation approximates the frequency—gravity impact relationship sufficiently well with negligible uncertainties, while the accuracy of the detected disturbance frequency defines a limiting factor for the gravity impact estimation. A realistic disturbance of 15 Hz with an amplitude of 1.5 nm had an impact of ≈48 [μGal] on the gravity estimate. The proposed filter approach reduced the impact to ≈12 [μGal], with the remaining effect being almost entirely associated to the uncertainty in disturbance frequency detection.


Introduction
In recent decades, a large amount of FG5 absolute gravimeter data have been accumulated in North America by agencies such as the National Geodetic Survey (USA) and the Geological Survey of Canada.A part of these observations has not been fully exploited yet for different reasons.For example, the FG5 requires an intensive and regular maintenance program and, sometimes, measurements performed at the end of a maintenance cycle are of lower quality.In addition, disturbances (defined as a signal plus noise) during a measurement campaign contaminate absolute gravimeter observations, e.g., signals/ noise caused by construction, microseismicity, or instrumental disturbances.Many studies have addressed individual components of FG5 instrumental or environmental disturbances, cf., [1][2][3][4].Since FG5 observations are very time consuming and costly, the re-assessment of contaminated data through improved de-noising and analysis methods is promising.In addition to the economical aspects, the re-assessment is also beneficial for the gravimetry and geoid communities, which rely on these fundamental observations for calibration or validation of satellite or airborne gravity observations or for the definition of a gravimetric vertical datum.Further, existing time series could be extended and the determination of gravity could be improved [4,5].It is often the length of the gravity time series which determines its usefulness as temporal changes in gravity are often caused by slow geodynamic processes such as glacial isostatic adjustment [6][7][8].The re-assessment of previously discarded FG5 data is the only means to extend a gravity time series back in time.The benefits of being able to filter out detected disturbances are twofold, namely, 1) the determination of gravity is more accurate; and 2) noisy time series of the FG5 observations, which have not been used for analysis due to their high noise level or contamination level could be revitalized.
In this study, we will demonstrate how a detected disturbance in a FG5 time series impacts the estimated gravity value.This approach makes the assumption that a disturbance has been identified in the time series, e.g. using methods such as Lomb-Scargle periodogram analysis or wavelet analysis as demonstrated in [9,10].A very simple way of filtering out a signal (not considering noise as in a disturbance) from FG5 observations can be achieved by modifying the equation of motion used in the least squares fit, cf., [4,10].[4,11] proposed a technique which includes a damping factor and a summation over several signals.Herein, no damping factor is included and only one signal is considered at a time.The reason for neglecting a damping factor is because the uncertainty in the required prior frequency detection presents the limiting factor and dominates the damping effects.The analysed time-distance measurements are composed by adding a synthetic signal to a real FG5 data set.This ensures that any real signals originating from the instrument or environment as well as their noise is contained in the data and does not need to be simulated.Only the added signal is perfectly known.This combination allows us to consider the analysed data set as representative for the FG5 measurement process and the following adjustment to estimate gravity.
A detected disturbance frequency, D f , can be considered by adding a sinusoidal term Herein, 0 , 0 z g , and 0 are the unknown initial parameters and and are the laser modulation parameters (for details on Equation 1 refer to [10,12,13]).The detected frequency is treated as a 'true' parameter (indicated by the tilde), i.e., v B C D f  is not introduced as an unknown parameter which is estimated by the fitting algorithm, cf., the usage of "true" parameter in [14].However, if Equation ( 1) is used as a functional model, the following vectors and matrices must be modified within the least squares fit algorithm, cf.[10], and,

X
and the design matrix becomes .
For a signal amplitude of guess Before starting a time consuming reprocessing with a modified observation equation for large amounts of data, it is beneficial to investigate what level of improvement can be achieved.In other words, if a signal in a data set was detected, e.g., by spectral analysis tools such as Lomb-Scargle or wavelets [10], an estimation of its impact on the derived gravity value would enable the quantification of potential improvements through filtering, an initial step in the process of re-assessment.In this study, we present a numerical tool for the estimation of the disturbance signal impact on gravity.There are multiple ways of filtering a detected disturbance in order to decrease its impact on gravity estimates.Herein, only one filtering approach is presented, as the focus of this study lies on the impact assessment.This includes the important consideration that the FG5 observations at hand first go through the transfer function of the FG5 measurement process, then through an adjustment process, before estimating a gravity value.In order to keep most of these processes in the impact assessment, we used real observations together with known synthetic disturbances in this study.
The FG5 data used herein were provided by Natural Resources Canada.The measurements used in the text were collected July 8, 2008 at the Pacific Geoscience Centre in Sidney, British Columbia, Canada.A data subset of 100 drops was used for most computations (FG5#236: file 2008a0708 ({*}.raw,{*}.project),Sets 37-40, 25 Drops/Set, total 100 Drops).

Signal Impact Analysis
For the following synthetic analysis a real FG5 data set of 100 drops with all time-distance measurements, i and i , was used.The 100 drops were analyzed before by Lomb-Scargle Periodogram analysis to be sure that no significant disturbance exists in the data set [9].A sinusoidal signal was added to the distance measurement sin 2π where A is the amplitude in   nm and f the signal . Further, in Equation ( 2) the index denotes the used fringes, .Furthermore, from drop to drop the phase was chanced to cover the range of within the 100 drops.Hence, the phase is defined by , , , j  with .A full parameter space study of how sensitive the FG5 measurement and adjustment process is for varying disturbance frequency, amplitude and phase was performed in [9].We used the results of that study for defining relevant frequencies, amplitudes and phase.A least squares fit was performed for the 100 drops with both types of distances (without an added signal and with an added signal ).This can be expressed by , , , , , , , , , .
and the standard deviation of the set of differences was calculated.This analysis was performed for the set of frequencies   0.2, ,99.8,100  given in Hz .Figure 1 shows the results of a frequency analysis up to 100 [Hz] with three different amplitudes.It also validates the results of [15] and provides an analysis tool similar to the results of [16], who have conducted a similar impact analysis for geodetic velocity estimation.The most important conclusions from analyzing Figure 1 are: 1) the lower signal frequencies have the largest impact on gravity; and 2) small changes in frequency can have significantly different impacts on gravity.These facts are caused intrinsically by the way the FG5 transfers a disturbance to time-distance measurements and finally to the gravity estimate.It is worth to note that changing the number of selected fringes will lead to a different drop frequency, cf equation in next page, which in turn will change the impact on gravity.Hence, the curves in Fig- ure 1 must be re-evaluated every time, the fundamental parameters drop frequency, number of drops, selected start and end fringes, amplitude and phase change.

Estimation of Signal Impact on Gravity
In order to perform the task of data re-assessment in an effective manner, a tool for the estimation of the signal impact on gravity is very useful.Figure 1 shows the frequency sensitivity of the FG5 and its impact on gravity.Since the standard deviation is based on a set of gravity The differences between each pair of affected and unaffected gravity estimates were computed for each drop, i.e., differences, not on absolute values, the plot is of a general character.Figure 1 can be used to estimate the impact on gravity for a detected frequency, D f , and for a certain amplitude, which has to be estimated.Since the curves in Figure 1 have a distinct shape and are continuously differentiable, it should be possible to find an appropriate approximation using analytical functions.It is obvious that this desired function depends on the disturbance amplitude and drop period/frequency.Furthermore, it seems possible that the curves shown in Figure 1 can be approximated by a summation of several Gaussian bells.In general, the Gaussian bell summation is given by In order to determine a relationship between the Gaussian parameters i i i and the amplitude of the disturbance A (in [nm]) and the drop period and ef refer to the start fringe and end fringe), several curves as in Figure 1 with different amplitudes and drop periods were analyzed.A parameter space of nm was used.For each constellation of amplitude and frequency the phase was evenly spread between, i.e.,

 
2π 1 where is the drop index and the number of used drops.With this parameter space, the curves were generated and a least squares fit for each curve was performed to obtain the three parameters for each Gaussian bell.The results for each parameter were analyzed.The heuristically identified relationship between a b is described by with and is valid for frequencies up to 100 Hz 14 n  a b , which corresponds to a summation of Gaussian bells.Parameter i corresponds to the amplitude of the i -th Gaussian bell, and is controlled by the amplitude of the synthetic signal and the drop period used.Parameter, i , describes the position of the center of each Gaussian bell.A closer look reveals that the first peak is not exactly at the drop frequency 1 1 but has a shift of 0.5149 Hz   on average.It is very likely that this is caused by the intrinsic data window of the FG5 observations.A periodic behavior of this shift was detected between  0.39 -0.65 Hz 1 peak st , and it depends on the drop period used.The differences over the fringe space of

DF
of this length is taken to define each bell width.In Figure 3, the approximations are shown for two different amplitudes and for two different .In general, the approximation at the beginning is worse than at the end; this can be observed in the second maximum in each subplot.However, this discrepancy is acceptable, since the goal is to estimate the impact on gravity, and also the disturbance amplitude.Pseudo-Voigt functions representing sums of Gaussian and Lorentzian functions [17,18] were tested instead of Gaussian bells, which provided a better approximation.However, this introduced more parameters, which made it much more difficult to find a relationship between the parameters, and the signal amplitude.In Table 1, a numerical example for a single drop is given.The generated synthetic signal of   .The remaining difference is caused by the fact that Lomb-Scargle is not able to detect the exact frequency of  

15
, cf. [9].The results of the esti-Hz   4) are acceptable and the misfit described in Figure 3 has a negligible impact on gravity compared to the limitations in frequency detection using Lomb-Scargle.Because of the strong slope of the function at this point, a difference of 0. .
The numerical example in Table 1 demonstrates that the filtering efficiency depends mostly on the detection accuracy.The guessed amplitude has a minor impact on g  and can be considered as very robust.

Conclusion
A Gaussian Bell Summation   GBS T 

GBS
was developed to approximate the frequency dependent impact of disturbances on gravity estimates (Figure 1 and 3).The approximation requires only two parameters, the drop period , and a rough estimate of the disturbance amplitude, guess .In addition to the analytical term, a numerical comparison of estimation and filtering is made for a single drop to demonstrate the performance.The filter is realized by a modification of the FG5 least-squares observation equation and uses the detected frequency as a "true observation" in the adjustment algorithm.The modified least-squares fit algorithm estimates the unknown amplitude of the signal quite well and does not depend on the guessed initial amplitude value.The developed can be used to obtain estimates of the impact of known or detected disturbances on gravity.It is a useful tool for the re-assessment of previously discarded FG5 observations.It can be used to improve the estimated gravity values markedly.The limiting factor of this approach is the uncertainty in disturbance frequency detection, e.g. by using Lomb-Scargle periodogram analysis, which prevends the proposed filter to remove all of the disturbance impact.Initial tests using a disturbance frequency detection algorithm based on wavelets [10] shows improved accuracy and would further reduce the gravity impact of detectable disturbances.

Figure 2 .
For the sake of convenience, the average of keep Equation (4) simple.Parameter i controls the width of each Gaussian bell.Since local maxima are around the multiples of , the space between two local maxima can be well -Scargle and filtering via Equation (1) the impact can be reduced considerably to   12.39 μGal 

Figure 2 .Figure 3 .
Figure 2. Δf is the difference of frequency at the first peak (in the plot of standard deviation analysis of g  , cf., max peak in Figure 1) and DF depending on the start and end fringes used.Data used: FG5#236, file 2008a0708, sets 37-40, (25 drops/set, total 100 drops).