The processes responsible for land surface subsidence in the Mississippi Delta (MD) have been vigorously debated. Numerous studies have postulated a dominant role for isostatic subsidence associated with sediment loading. Previous computational modeling of present-day vertical land motion has been carried out in order to understand geodetic data. While the magnitudes of these measured rates have been reproduced, the model parameter values required have often been extreme and, in some cases, unrealistic. In contrast, subsidence rates in the MD on the 10 3 year timescale due to delta loading estimated from relative sea level reconstructions are an order of magnitude lower. In an attempt to resolve this conflict, a sensitivity analysis was carried out using a spherically symmetric viscoelastic solid Earth deformation model with sediment, ice, and ocean load histories. The model results were compared with geologic and geodetic observations that provide a record of vertical land motion over three distinctly different timescales (past 80 kyr, past 7 kyr, and past~15 years). It was found that glacial isostatic adjustment is likely to be the dominant contributor to vertical motion of the Pleistocene and underlying basement. Present-day basement subsidence rates solely due to sediment loading are found to be less than~0.5 mm yr À1 . The analysis supports previous suggestions in the literature that Earth rheology parameters are time dependent. Specifically, the effective elastic thickness of the lithosphere may be <50 km on a 10 5 year timescale, but closer to 100 km over 10 3 to 10 4 year timescales.

The measurement of ongoing ice-mass loss and associated melt water contribution to sea-level change from regions such as West Antarctica is dependent on a combination of remote sensing methods. A key method, the measurement of changes in Earth's gravity via the GRACE satellite mission, requires a potentially large correction to account for the isostatic response of the solid Earth to ice-load changes since the Last Glacial Maximum. In this study, we combine glacial isostatic adjustment modelling with a new GPS dataset of solid Earth deformation for the southern Antarctic Peninsula to test the current understanding of ice history in this region. A sufficiently complete history of past ice-load change is required for glacial isostatic adjustment models to accurately predict the spatial variation of ongoing solid Earth deformation, once the independently-constrained effects of present-day ice mass loss have been accounted for. Comparisons between the GPS data and glacial isostatic adjustment model predictions reveal a substantial misfit. The misfit is localized on the southwestern Weddell Sea, where current ice models under-predict uplift rates by approximately 2 mm yr −1. This under-prediction suggests that either the retreat of the ice sheet grounding line in this region occurred significantly later in the Holocene than currently assumed, or that the region previously hosted more ice than currently assumed. This finding demonstrates the need for further fieldwork to obtain direct constraints on the timing of Holocene grounding line retreat in the southwestern Weddell Sea and that GRACE estimates of ice sheet mass balance will be unreliable in this region until this is resolved.

No abstract

Sea level rise presents a hazard for coastal populations, and the Mississippi Delta (MD) is a region particularly at risk due to the high rates of land subsidence. We apply a gravitationally self‐consistent model of glacial and sediment isostatic adjustment (SIA) along with a realistic sediment load reconstruction in this region for the first time to determine isostatic contributions to relative sea level (RSL) and land motion. We determine optimal model parameters (Earth rheology and ice history) using a new high‐quality compaction‐free sea level indicator database. Using the optimal model parameters, we show that SIA can lower predicted RSL in the MD area by several meters over the Holocene and so should be taken into account when modeling these data. We compare modeled contemporary rates of vertical land motion with those inferred using GPS. This comparison indicates that isostatic processes can explain the majority of the observed vertical land motion north of latitude 30.7°N, where subsidence rates average about 1 mm/yr; however, subsidence south of this latitude shows large data‐model discrepancies of greater than 3 mm/yr, indicating the importance of nonisostatic processes. This discrepancy extends to contemporary RSL change, where we find that the SIA contribution in the Delta is on the order of 10−1 mm/yr. We provide estimates of the isostatic contributions to 20th and 21st century sea level rates at Gulf Coast Permanent Service for Mean Sea Level tide gauge locations as well as vertical and horizontal land motion at GPS station locations near the MD.

No abstract

Abstract. We investigate the influence on mantle convection of the negative Clapeyron slope ringwoodite to perovskite and ferro-periclase mantle phase transition, which is correlated with the seismic discontinuity at 660 km depth. In particular, we focus on understanding the influence of the magnitude of the Clapeyron slope (as measured by the Phase Buoyancy parameter, P ) and the vigour of convection (as measured by the Rayleigh number, Ra) on mantle convection. We have undertaken 76 simulations of isoviscous mantle convection in spherical geometry, varying Ra and P . Three domains of behaviour were found: layered convection for high Ra and more negative P , whole mantle convection for low Ra and less negative P , and transitional behaviour in an intervening domain. The boundary between the layered and transitional domain was fit by a curve P = αRa β where α = −1.05, and β = −0.1, and the fit for the boundary between the transitional and whole mantle convection domain was α = −4.8, and β = −0.25. These two curves converge at Ra ≈2.5×10 4 (well below Earth mantle vigour) and P ≈ −0.38. Extrapolating to high Ra, which is likely earlier in Earth history, this work suggests a large transitional domain. It is therefore likely that convection in the Archean would have been influenced by this phase change, with Earth being at least in the transitional domain, if not the layered domain.

The evolution of the planetary interior during plate tectonics is controlled by slow convection within the mantle. Global-scale geochemical differences across the upper mantle are known, but how they are preserved during convection has not been adequately explained. We demonstrate that the geographic patterns of chemical variations around the Earth’s mantle endure as a direct result of whole-mantle convection within largely isolated cells defined by subducting plates. New 3D spherical numerical models embedded with the latest geological paleo-tectonic reconstructions and ground-truthed with new Hf-Nd isotope data, suggest that uppermost mantle at one location (e.g. under Indian Ocean) circulates down to the core-mantle boundary (CMB), but returns within ≥100 Myrs via large-scale convection to its approximate starting location. Modelled tracers pool at the CMB but do not disperse ubiquitously around it. Similarly, mantle beneath the Pacific does not spread to surrounding regions of the planet. The models fit global patterns of isotope data and may explain features such as the DUPAL anomaly and long-standing differences between Indian and Pacific Ocean crust. Indeed, the geochemical data suggests this mode of convection could have influenced the evolution of mantle composition since 550 Ma and potentially since the onset of plate tectonics.

We investigate the influence on mantle convection of the negative Clapeyron slope ringwoodite to perovskite and ferro-periclase mantle phase transition, which is correlated with the seismic discontinuity at 660 km depth. In particular, we focus on understanding the influence of the magnitude of the Clapeyron slope (as measured by the Phase Buoyancy parameter, <i>P</i>) and the vigour of convection (as measured by the Rayleigh number, <i>Ra</i>) on mantle convection. We have undertaken 76 simulations of isoviscous mantle convection in spherical geometry varying <i>Ra</i> and <i>P</i>. Three domains of behaviour were found: layered convection for high <i>Ra</i> and more negative <i>P</i>, whole mantle convection for low <i>Ra</i> and less negative <i>P</i> and transitional behaviour in an intervening domain. The boundary between the layered and transitional domain was fit by a curve <i>P</i> = &alpha;<i>Ra</i>^&beta; where α = −1.05, and β = −0.1, and the fit for the boundary between the transitional and whole mantle convection domain was α = −4.8, and β = −0.25. These two curves converge at <i>Ra</i>≈2.5×10⁴ and <i>P</i>&approx;&minus;0.38. Extrapolating to high <i>Ra</i>, which is likely earlier in Earth history, this work suggests a large transitional domain. It is therefore likely that convection in the Archean would have been influenced by this phase change, with Earth being at least in the transitional domain, if not the layered domain