A new integrated geophysical-petrological global three dimensional model of upper-mantle electrical conductivity validated by the Swarm M2 tidal magnetic field

Summary A new global model of the present-day thermochemical state of the lithosphere and upper mantle based on global waveform inversion, satellite gravity and gradiometry measurements, surface elevation and heat flow data has been recently presented: WINTERC-G (Fullea et al., 2021). WINTERC-G is built within an integrated geophysical-petrological framework where the mantle seismic velocity and density fields are computed in a thermodynamically self-consistent framework, allowing for a direct parameterisation in terms of the temperature, pressure and composition of the subsurface rocks. In this paper, we combine WINTERC-G thermal and compositional fields along with laboratory experiments constraining the electrical conductivity of mantle minerals, melt and water, and derive a set of new global three dimensional electrical conductivity models of the upper mantle. The new conductivity models, WINTERC-e, consist of two end-members corresponding to minimum and maximum conductivity of the in-situ rock aggregate accounting for mantle melting, mineral water content and the individual conductivities of the main stable mantle mineral phases. The end-member models are validated over oceans by simulating the magnetic field induced by the ocean M2 tidal currents and comparing the predicted fields with the M2 tidal magnetic field estimated from six-year Swarm satellite observations. Our new conductivity model, derived independently from any surface or satellite magnetic data sets, is however able to predict tidal magnetic fields that are in good agreement with the Swarm M2 tidal magnetic field models estimated by Sabaka et al. (2018, 2020) and Grayver & Olsen (2019). Our predicted M2 tidal magnetic fields differ in amplitudes by about 5-20% from the Swarm M2 tidal magnetic field, with the high conductivity WINTERC-e end-member model accounting for mantle melt and water content capturing the structure of Swarm data better than the low conductivity end-member model. Spherically symmetric conductivity models derived by averaging new WINTERC-e conductivities over oceanic areas are slightly more conductive than the recent global conductivity models AA17 by Grayver et al. (2017) derived from Swarm and CHAMP satellite data in the 60-140 km depth range, while they are less conductive deeper in the mantle. The conductivities in WINTERC-e are about 3-4 times smaller than the AA17 conductivities at a depth of 400 km. Despite the differences in electrical conductivity, our spherically symmetric high conductivity end-member model WINTERC-e captures the structure of Swarm M2 tidal magnetic field almost the same as a state of the art 1D conductivity models derived entirely from magnetic data (AA17, (Grayver et al., 2017). Moreover, we show that realistic lateral electrical conductivity inhomogeneities of the oceanic upper mantle derived from the temperature, melt and water distributions in WINTERC-e contribute to the M2 tidal magnetic field up to ±0.3 nT at 430 km altitude.

Numerical modelling of landscape and sediment flux response to precipitation rate change

Laboratory-scale experiments of erosion have demonstrated that landscapes have a natural (or intrinsic) response time to a change in precipitation rate. In the last few decades there has been growth in the development of numerical models that attempt to capture landscape evolution over long timescales. However, there is still an uncertainty regarding the validity of the basic assumptions of mass transport that are made in deriving these models. In this contribution we therefore return to a principal assumption of sediment transport within the mass balance for surface processes; we explore the sensitivity of the classic end-member landscape evolution models and the sediment fluxes they produce to a change in precipitation rates. One end-member model takes the mathematical form of a kinetic wave equation and is known as the stream power model, in which sediment is assumed to be transported immediately out of the model domain. The second end-member model is the transport model and it takes the form of a diffusion equation, assuming that the sediment flux is a function of the water flux and slope. We find that both of these end-member models have a response time that has a proportionality to the precipitation rate that follows a negative power law. However, for the stream power model the exponent on the water flux term must be less than one, and for the transport model the exponent must be greater than one, in order to match the observed concavity of natural systems. This difference in exponent means that the transport model generally responds more rapidly to an increase in precipitation rates, on the order of 105 years for post-perturbation sediment fluxes to return to within 50 % of their initial values, for theoretical landscapes with a scale of 100×100 km. Additionally from the same starting conditions, the amplitude of the sediment flux perturbation in the transport model is greater, with much larger sensitivity to catchment size. An important finding is that both models respond more quickly to a wetting event than a drying event, and we argue that this asymmetry in response time has significant implications for depositional stratigraphies. Finally, we evaluate the extent to which these constraints on response times and sediment fluxes from simple models help us understand the geological record of landscape response to rapid environmental changes in the past, such as the Paleocene–Eocene thermal maximum (PETM). In the Spanish Pyrenees, for instance, a relatively rapid (10 to 50 kyr) duration of the deposition of gravel is observed for a climatic shift that is thought to be towards increased precipitation rates. We suggest that the rapid response observed is more easily explained through a diffusive transport model because (1) the model has a faster response time, which is consistent with the documented stratigraphic data, (2) there is a high-amplitude spike in sediment flux, and (3) the assumption of instantaneous transport is difficult to justify for the transport of large grain sizes as an alluvial bedload. Consequently, while these end-member models do not reproduce all the complexity of processes seen in real landscapes, we argue that variations in long-term erosional dynamics within source catchments can fundamentally control when, how, and where sedimentary archives can record past environmental change.

Numerical modelling landscape and sediment flux response to precipitation rate change

Laboratory-scale experiments of erosion have demonstrated that landscapes have a natural (or intrinsic) response time to a change in precipitation rate. In the last few decades there has been a growth in the development of numerical models that attempt to capture landscape evolution over long time-scales. Recently, a sub-set of these numerical models have been used to invert river profiles for past tectonic conditions even during variable climatic conditions. However, there is still an uncertainty over validity of the basic assumption of mass transport that are made in deriving these models. In this contribution we therefore return to a principle assumption of sediment transport within the mass balance for surface processes, and explore the sensitivity of the classic end-member landscape evolution models to change in precipitation rates. One end-member model takes the mathematical form of a kinetic wave equation and is known as the stream power model, where sediment is assumed to be transported immediately out of the model domain. The second end-member model takes the form of a diffusion equation, and assumes that the sediment flux is a function of the water flux and slope. We find that both of these end-member models have a response time that has a proportionality to the precipitation rate that follows a negative power law. For the stream power model the exponent on the water flux term must be less than one, and for the sediment transport model the exponent must be greater than one in order to match the observed concavity of natural systems. This difference in exponent means that sediment transport model responds more rapidly to an increase in precipitation rates, on the order of 105 years for a landscape with a scale of 105 m. In nature, landscape response times to a rapid environmental change have been estimated for events such as the Paleocene-Eocene thermal maximum (PETM). In the Spanish Pyrenees, a relatively rapid, 20 to 100 kyr, duration of deposition of gravel during the PETM is observed for a climatic shift that is thought to be towards increased precipitation rates. We suggest the rapid response observed is more easily explained through a diffusive sediment transport model, as (1) this model has a faster response time, consistent with the documented stratigraphic data, and (2) the assumption of instantaneous transport is difficult to justify for the transport of large grain sizes as an alluvial bed-load.

The role of log jams and exceptional flood events in mobilizing coarse particulate organic matter in a steep headwater stream

Coarse particulate organic matter (CPOM) fulfills important functions in the physical and ecological system of a stream. CPOM delivery to and export from the stream has implications for the stream's morphology and sediment transport capacity as well as the energy budget and food availability. Export rates of CPOM from mountain catchments have been observed to strongly increase with rising discharge, but the mechanism leading to this strong relationship is unclear. Here, we show that log jams in the Erlenbach, a steep headwater stream in the Swiss Prealps, are an effective barrier for the transport of CPOM pieces, and thus become sites of storage of large quantities of material over time. Exceptional discharge events with return periods exceeding 20 years play a dual role in CPOM transport in the Erlenbach. First, they appear to destroy existing log jams, releasing the stored material (wood and sediment). Second, they intensify channel–hillslope coupling, thereby recruiting new logs to the channel, around which new jams can form. This allows for the formulation of a new, fully episodic end-member in a four-end-member model of CPOM dynamics of steep mountain streams based on wood delivery and export.