An accuracy assessment of the fan wavelet coherence method for elastic thickness estimation

Abstract

[1] A variety of methods exist to estimate the elastic thickness (T e ) of the lithosphere. In this contribution, we attempt to provide an indication of how well the fan wavelet coherence method recovers T e , through synthetic modeling. The procedure involves simulating initial topographic and subsurface loads and emplacing them on a thin elastic plate of known T e , generating the postloading topography and gravity. We then attempt to recover that T e distribution from the gravity and topography through the wavelet method, hence discovering where its strengths and weaknesses lie. The T e distributions we use here have elliptical and fractal geometries, while the initial loads are fractal. Importantly, we have found that this widely used synthetic loading calibration method will tend to result in underestimates of T e no matter which recovery method is used. This is due to random correlations between the initial loads which, on average, serve to increase their coherence at all wavelengths and spatial locations. For the fan wavelet method, the degree of underestimation from this ''background'' source is approximately 10% of the true T e . In addition, the fan wavelet coherence method will provide underestimates of (1) the true T e when the study area size is of the order of the highest flexural wavelength or less, (2) relative T e differences when the T e anomaly is narrow compared to its flexural wavelength, and (3) steep T e gradients. Significantly, we find that the recovery is not greatly affected by the assumption of uniform T e in the inversion of the coherence. We also find that T e recovery from the coherence is only weakly dependent upon the initial subsurface-to-surface loading ratio ( f ). In contrast to the coherence, T e recovery from the admittance is highly ''noisy,'' with discontinuities and overestimates of T e frequently arising. This is most likely due to the high sensitivity of the admittance to f and is likely to apply to real data as well.

Estimates ofTefor continental regions using GOCE gravity

“…The same is true of those obtained from surface data by M c Kenzie and Fairhead (1997). However, Kirby and Swain (2008) did not find this behaviour when they analysed synthetic data sets. The principal reason why we believe those from the admittance method are likely to be the more accurate than that from the Bouguer coherence concerns the ratio of T e to T s , the seismogenic thickness.…”

Assessment of flexural analysis applied to the Sumatra-Java subduction zone

“…Recently, Kirby and Swain (2007) demonstrated that T e was better recovered in all settings (i.e., land only, sea only and mixed land sea) when using the scaling (equivalent topography) and landloading equations. Hence we converted bathymetry to equivalent topography and free-air anomalies to complete Bouguer anomalies whose details are given here.…”

Spatial variations of effective elastic thickness of the lithosphere in Central America and surrounding regions

“…The Ikol value used in this study is 2.668, which give a space-domain wavelet whose first sidelobes is 1/16 of the magnitude of the central am plitude (Kirby and Swain, 2011 deconvolution equations to the entire data set. Kirby and Swain (2008) showed that although this approach may bias Te in ocean areas, the bias is small. Values of the densities are given in Table 1.…”