Abstract: Spherical cap harmonic analysis (SCHA) has become a common tool for the regional modeling of potential fields since its introduction by Haines (1985). The fact that SCHA satisfies Laplace equation and the possibility of representing high-frequency fields with a small number of coefficients (compared to the global spherical harmonic analysis) made SCHA the preferred choice for the development, for example, of magnetic field models at national scale. However, Thébault et al. (2006a) demonstrated that the traditional SCHA presented some deficiencies, in particular related to the inversion of multilevel data sets. The authors presented the R-SCHA technique as an alternative method in which the introduction of a new set of basis functions and boundary conditions solved this issue. In this paper we present some numerical comparisons between the SCHA and R-SCHA techniques applied with different synthetic vector data sets, from near-surface main field, main difference, and crustal field data simulating a World Digital Magnetic Anomaly Map subset. Other analyses are carried out with synthetic vector data set that mimics the expected data distribution from a multisatellite mission like the forthcoming European Swarm mission. No regularization, weighting, or ad hoc procedures are applied to the synthetic vector data, and a cap of 7° aperture is considered. The numerical analyses show that SCHA is a satisfying approximation in a band-limited spectral region that depends on the cap's size. It does not work correctly either for main field or for the short-scale crustal field modeling. These aspects are supported by equations illustrating why SCHA may fail. On the contrary, R-SCHA converges more slowly than SCHA but is valid in all cases. It gives a consistent set of regional coefficients and fits the radial variation of the field in a realistic way. At last, the special case of data incompatibility shows that R-SCHA does not fit incompatible data while SCHA assimilates most of them. These results should help the scientific community to evaluate the level of approximation needed for the development of regional magnetic field models in the era of the European Space Agency Swarm mission.

Abstract: Located both on the Mid-Atlantic Ridge and above a mantle plume, Iceland is subject to horizontal and vertical motions. Many studies described these deformations in terms of rifting episodes that have combined both extensional tectonics and magmatism. However, few studies have described the glacio-isostatic response induced by the retreat of the Weichselian ice cap. The melting of this ice cap induced a postglacial rebound for the whole of Iceland that may be controlled by the geodynamic setting and the rheological layering of the lithosphere. This study is devoted to (1) understanding the Holocene rebound on the southwestern coast and (2) estimating the asthenosphere viscosity and depth beneath Iceland. Two stages of holocene evolution were determined by means of GPS profiles, morphological observations, and data compilation. The first stage corresponds to a vertical uplift of 67.5 to 157.5 m. It started at 10,000 years BP and ended at 8500 years BP implying uplift rates between 4.5 and 10.5 cm/a. It was a quick isostatic response to the fast ice retreat. The second stage had vertical motion of tens of meters with a probable tectonic origin and started at 8500 years BP. The uplift rate is 1 to 2 orders of magnitude slower than the one during the first stage. Uplift partitioning during the first stage was controlled by the thermal state of the lithosphere, the highest geothermal flux inducing the maximum uplift rates. The relaxation time for uplift provides a viscosity estimate of 5.4–5.8 × 1019 Pa s for the asthenosphere. This value is similar to those determined for glacial areas in different continental contexts. However, the flexural wavelength indicates a shallower asthenosphere than that occurring in continental domains. Therefore this study highlights a coupling between the thermal structure of the Icelandic asthenosphere and the glacial rebound.

Abstract: A large data set of fundamental mode Rayleigh wave amplitudes is analyzed to derive a new global three-dimensional model of shear wave attenuation in the upper mantle. The amplitude observations span a range of periods between 50 and 250 s and are derived from earthquakes with M W > 6.0 that occurred between 1993 and 2005. Four separate factors may influence an amplitude anomaly: intrinsic attenuation along the raypath, elastic focusing effects along the raypath, uncertainties in the strength of excitation, and uncertainties in the response at the station. In an earlier paper (Dalton and Ekstrm, 2006a), dependence of the retrieved attenuation structure on these terms was shown to be significant and an approach was developed to invert the amplitudes simultaneously for each term. The new three-dimensional attenuation model QRFSI12, which is the subject of this paper, is derived using this method. The model contains large lateral variations in upper-mantle attenuation, 60% to 100%, and exhibits strong agreement with surface tectonic features at depths shallower than 200 km. At greater depth, QRFSI12 is dominated by high attenuation in the southeastern Pacific and eastern Africa and low attenuation along many subduction zones in the western Pacific. Resolution tests confirm that the change in pattern of attenuation above and below 200-km depth can be determined with confidence using the fundamental mode data set. The new model is highly correlated with global models of shear wave velocity, particularly in the uppermost mantle, suggesting that the same factors may control both seismic attenuation and velocity in this depth range. However, forcing the lateral perturbations in attenuation to match those found in global velocity models decreases the data variance reduction, which suggests that subtle differences between patterns of attenuation and velocity are robust.

Abstract: The relationship between intraplate volcanism and continental tectonics has been investigated for North and East Africa using a high resolution three-dimensional anisotropic tomographic model derived from seismic data of a French experiment “Horn of Africa” and existing broadband data. The joint inversion for seismic velocity and anisotropy of the upper 400 km of the mantle, and geochemical data reveals a complex interaction between mantle upwellings, and lithosphere. Two kinds of mantle upwellings can be distinguished: The first one, the Afar “plume” originates from deeper than 400 km and is characterized by enrichment in primordial 3He and 3He/4He ratios higher than those along mid-ocean ridges (MOR). The second one, associated with other Cenozoic volcanic provinces (Darfur, Tibesti, Hoggar, Cameroon), with 3He/4He ratios similar to, or lower than MOR, is a consequence of shallower upwelling. The presumed asthenospheric convective instabilities are oriented in an east-west direction, resulting from interaction between south-north asthenospheric mantle flow, main plume head and topography on the base of lithosphere.