SUMMARY
We present a new global thermochemical model of the lithosphere and underlying upper mantle constrained by state of the art seismic waveform inversion, satellite gravity (geoid and gravity ...anomalies and gradiometric measurements from ESA's GOCE mission), surface elevation and heat flow data: WINTERC-G. The model is based upon an integrated geophysical–petrological approach where seismic velocities and density in the mantle are computed within a thermodynamically self-consistent framework, allowing for a direct parametrization in terms of the temperature and composition variables. The complementary sensitivities of the data sets allow us to constrain the geometry of the lithosphere–asthenosphere boundary, to separate thermal and compositional anomalies in the mantle, and to obtain a proxy for dynamic surface topography. At long spatial wavelengths, our model is generally consistent with previous seismic (or seismically derived) global models and earlier integrated studies incorporating surface wave data at lower lateral resolution. At finer scales, the temperature, composition and density distributions in WINTERC-G offer a new state of the art image at a high resolution globally (225 km average interknot spacing). Our model shows that the deepest lithosphere–asthenosphere boundary is associated with cratons and, also, some tectonically active areas (Andes, Persian Gulf). Among cratons we identify considerable differences in temperature and composition. The North American and Siberian Cratons are thick (>260 km) and compositionally refractory, whereas the Sino-Korean, Aldan and Tanzanian Cratons have a thinner, fertile lithosphere, similar to younger continental lithosphere elsewhere. WINTERC-G shows progressive thickening of oceanic lithosphere with age, but with significant regional differences: the lithospheric mantle beneath the Atlantic and Indian Oceans is, on average, colder, more fertile and denser than that beneath the Pacific Ocean. Our results suggest that the composition, temperature and density of the oceanic mantle lithosphere are related to the spreading rate for the rates up to 50–60 mm yr–1: the lower spreading rate, the higher the mantle fertility and density, and the lower the temperature. At greater spreading rates, the relationship disappears. The 1-D radial average of WINTERC-G displays a mantle geothermal gradient of 0.55–0.6 K km–1 and a potential temperature of 1300–1320 °C for depths >200 km. At the top of the mantle transition zone the amplitude of the maximum lateral temperature variations (cratons versus hotspots) is about 120 K. The isostatic residual topography values, a proxy for dynamic topography, are large (>1 km) mostly in active subduction settings. The residual isostatic bathymetry from WINTERC-G is remarkably similar to the pattern independently determined based on oceanic crustal data compilations. The amplitude of the continental residual topography is relatively large and positive (>600 m) in the East European Craton, Greenland, and the Andes and Himalayas. By contrast, central Asia, most of Antarctica, southern South America and, to a lesser extent, central Africa are characterized by negative residual topography values (>–400 m). Our results show that a substantial part of the topography signal previously identified as residual (or dynamic) is accounted for, isostatically, by lithospheric density variations.
SUMMARY
The interaction of the oceanic tidal flow with the Earth’s main magnetic field provides a powerful natural source of electromagnetic (EM) energy suitable for suboceanic upper-mantle ...electrical conductivity sounding. In this paper, we have developed and tested a new frequency-domain, spherical harmonic-finite element approach to the inverse problem of global EM induction. It is set up for an effective inversion of satellite-observed tidally induced magnetic field in terms of 3-D structure of the electrical conductivity in the suboceanic upper mantle. Before proceeding to the inversion of Swarm-derived models of tidal magnetic signatures, we have performed a series of parametric studies, using the 3-D conductivity model WINTERC-e as a testbed. The WINTERC-e model has been derived using state-of-the-art laboratory conductivity measurements of mantle minerals, and thermal and compositional model of the lithosphere and upper mantle WINTERC-G. The latter model is based on the inversion of global surface waveforms, satellite gravity and gradiometry measurements, surface elevation and heat flow data in a thermodynamically self-consistent framework. Therefore, the WINTERC-e model, independent of any EM data, represents an ideal target for synthetic tests of the 3-D EM inversion. We tested the impact of the truncation degree of the spherical-harmonic expansion of the M2 tidal signal, the effect of random noise in synthetic data and inclusion of the N2 and O1 tidal constituents on the ability to recover the suboceanic upper-mantle conductivity structure. We demonstrate that with suitable regularization we can successfully reconstruct the 3-D upper-mantle conductivity beneath world oceans. In the ideal noise-free case, the correlation coefficient between the target and recovered conductivity is greater than 0.8 in the 150–270 km depth range.
Display omitted
•Local Green's functions to estimate hydrological loading displacements.•Consideration of heterogeneities in Earth's crust.•Scatter plots show dependency between loading response and ...local crustal structure.•Crustal stratification influences near-field displacements.•Impact of crustal structure is relevant for hydrological loading displacements.
The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500km2, as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.
We derive a set of local Green's functions defined globally on a 1°×1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.
We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 ...to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 plus or minus 18 Gt yr super(-1). Therefore, our AIS mass balance of -114 plus or minus 23 Gt yr super(-1) is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (-26 plus or minus 3 Gt yr super(-1) and -127 plus or minus 7 Gt yr super(-1), respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 plus or minus 13 Gt yr super(-1)), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 plus or minus 7 Gt yr super(-1) (ca. 0.4 mm yr super(-1)) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.
Interannual ice-mass variations along the Antarctic Peninsula (AP) and in the Amundsen Sea Sector (AS) are obtained for the years 2002 until 2009 using satellite data of the Gravity Recovery and ...Climate Experiment, that correlate well (r
≈
0.7) with accumulation variations based on the net precipitation from the European Centre for Medium Range Weather Forecasts. Moreover, mass signals for AP and AS are anti-correlated in time (r
≈
−
0.4) and contain El Niño Southern Oscillation signatures related to the strength of the Amundsen Sea Low pressure system, that has a dominant influence on West Antarctic atmospheric moisture transports. The GRACE interannual mass variations exhibit root-mean squared amplitudes of: 16.4
±
4.1
Gt (AP) and 28.6
±
10.5
Gt (AS), which are significant compared to the mean annual mass loss of −
110.2
±
6.7
Gt/a in coastal West Antarctica.
In this study, the memory of the Greenland Ice Sheet (GIS) with respect to its past states is analyzed. According to ice core reconstructions, the present‐day GIS reflects former climatic conditions ...dating back to at least 250 thousand years before the present (kyr BP). This fact must be considered when initializing an ice sheet model. The common initialization techniques are paleoclimatic simulations driven by atmospheric forcing inferred from ice core records and steady state simulations driven by the present‐day or past climatic conditions. When paleoclimatic simulations are used, the information about the past climatic conditions is partly reflected in the resulting present‐day state of the GIS. However, there are several important questions that need to be clarified. First, for how long does the model remember its initial state? Second, it is generally acknowledged that, prior to 100 kyr BP, the longest Greenland ice core record (GRIP) is distorted by ice‐flow irregularities. The question arises as to what extent do the uncertainties inherent in the GRIP‐based forcing influence the resulting GIS? Finally, how is the modeled thermodynamic state affected by the choice of initialization technique (paleo or steady state)? To answer these questions, a series of paleoclimatic and steady state simulations is carried out. We conclude that (1) the choice of an ice‐covered initial configuration shortens the initialization simulation time to 100 kyr, (2) the uncertainties in the GRIP‐based forcing affect present‐day modeled ice‐surface topographies and temperatures only slightly, and (3) the GIS forced by present‐day climatic conditions is overall warmer than that resulting from a paleoclimatic simulation.
Global navigation satellite systems (GNSSs) have revealed that a mega-thrust earthquake that occurs in an island-arc trench system causes post-seismic crustal deformation. Such crustal deformation ...data have been interpreted by combining three mechanisms: afterslip, poroelastic rebound and viscoelastic relaxation. It is seismologically important to determine the contribution of each mechanism because it provides frictional properties between the plate boundaries and viscosity estimates in the asthenosphere which are necessary to evaluate the stress behaviour during earthquake cycles. However, the observation sites of GNSS are mostly deployed over land and can detect only a small part of the large-scale deformation, which precludes a clear separation of the mechanisms. To extend the spatial coverage of the deformation area, recent studies started to use satellite gravity data that can detect long-wavelength deformations over the ocean. To date, compared with theoretical models for calculating the post-seismic crustal deformation, a few models have been proposed to interpret the corresponding gravity variations. Previous approaches have adopted approximations for the effects of compressibility, sphericity and self-gravitation when computing gravity changes. In this study, a new spectral-finite element approach is presented to consider the effects of material compressibility for Burgers viscoelastic earth model with a laterally heterogeneous viscosity distribution. After the basic principles are explained, it is applied to the 2004 Sumatra–Andaman earthquake. For this event, post-seismic deformation mechanisms are still a controversial topic. Using the developed approach, it is shown that the spatial patterns of gravity change generated by the above three mechanisms clearly differ from one another. A comparison of the theoretical simulation results with the satellite gravity data obtained from the Gravity Recovery and Climate Experiment reveals that both afterslip and viscoelastic relaxation are occurring. Considering the spatial patterns in satellite gravity fields is an effective method for investigating post-seismic deformation mechanisms.