2021, Wolf, Frederik, Voigt, Aiko, Donner, Reik V.

The intertropical convergence zone (ITCZ) is an important component of the tropical rain belt. Climate models continue to struggle to adequately represent the ITCZ and differ substantially in its simulated response to climate change. Here we employ complex network approaches, which extract spatiotemporal variability patterns from climate data, to better understand differences in the dynamics of the ITCZ in state-of-the-art global circulation models (GCMs). For this purpose, we study simulations with 14 GCMs in an idealized slab-ocean aquaplanet setup from TRACMIP – the Tropical Rain belts with an Annual cycle and a Continent Model Intercomparison Project. We construct network representations based on the spatial correlation patterns of monthly surface temperature anomalies and study the zonal-mean patterns of different topological and spatial network characteristics. Specifically, we cluster the GCMs by means of the distributions of their zonal network measures utilizing hierarchical clustering. We find that in the control simulation, the distributions of the zonal network measures are able to pick up model differences in the tropical sea surface temperature (SST) contrast, the ITCZ position, and the strength of the Southern Hemisphere Hadley cell. Although we do not find evidence for consistent modifications in the network structure tracing the response of the ITCZ to global warming in the considered model ensemble, our analysis demonstrates that coherent variations of the global SST field are linked to ITCZ dynamics. This suggests that climate networks can provide a new perspective on ITCZ dynamics and model differences therein.

2021, Haruzi, Peleg, Katsman, Regina, Halisch, Matthias, Waldmann, Nicolas, Spiro, Baruch

This paper presents a detailed description and evaluation of a multi-methodological petrophysical approach for the comprehensive multi-scale characterization of reservoir sandstones. The suggested methodology enables the identification of links between Darcy-scale permeability and an extensive set of geometrical, textural and topological rock descriptors quantified at the pore scale. This approach is applied to the study of samples from three consecutive sandstone layers of Lower Cretaceous age in northern Israel. These layers differ in features observed at the outcrop, hand specimen, petrographic microscope and micro-CT scales. Specifically, laboratory porosity and permeability measurements of several centimetre-sized samples show low variability in the quartz arenite (top and bottom) layers but high variability in the quartz wacke (middle) layer. The magnitudes of this variability are also confirmed by representative volume sizes and by anisotropy evaluations conducted on micro-CT-imaged 3-D pore geometries. Two scales of directional porosity variability are revealed in quartz arenite sandstone of the top layer: the pore size scale of ∼0.1 mm in all directions and ∼3.5 mm scale related to the occurrence of high- and low-porosity horizontal bands occluded by Fe oxide cementation. This millimetre-scale variability controls the laboratory-measured macroscopic rock permeability. More heterogeneous pore structures were revealed in the quartz wacke sandstone of the intermediate layer, which shows high inverse correlation between porosity and clay matrix in the vertical direction attributed to depositional processes and comprises an internal spatial irregularity. Quartz arenite sandstone of the bottom layer is homogenous and isotropic in the investigated domain, revealing porosity variability at a ∼0.1 mm scale, which is associated with the average pore size. Good agreement between the permeability upscaled from the pore-scale modelling and the estimates based on laboratory measurements is shown for the quartz arenite layers. The proposed multi-methodological approach leads to an accurate petrophysical characterization of reservoir sandstones with broad ranges of textural, topological and mineralogical characteristics and is particularly applicable for describing anisotropy and heterogeneity of sandstones on various rock scales. The results of this study also contribute to the geological interpretation of the studied stratigraphic units.

2021, Alberti, Tommaso, Donner, Reik V., Vannitsem, Stéphane

Atmosphere and ocean dynamics display many complex features and are characterized by a wide variety of processes and couplings across different timescales. Here we demonstrate the application of multivariate empirical mode decomposition (MEMD) to investigate the multivariate and multiscale properties of a reduced order model of the ocean–atmosphere coupled dynamics. MEMD provides a decomposition of the original multivariate time series into a series of oscillating patterns with time-dependent amplitude and phase by exploiting the local features of the data and without any a priori assumptions on the decomposition basis. Moreover, each oscillating pattern, usually named multivariate intrinsic mode function (MIMF), represents a local source of information that can be used to explore the behavior of fractal features at different scales by defining a sort of multiscale and multivariate generalized fractal dimensions. With these two complementary approaches, we show that the ocean–atmosphere dynamics presents a rich variety of features, with different multifractal properties for the ocean and the atmosphere at different timescales. For weak ocean–atmosphere coupling, the resulting dimensions of the two model components are very different, while for strong coupling for which coupled modes develop, the scaling properties are more similar especially at longer timescales. The latter result reflects the presence of a coherent coupled dynamics. Finally, we also compare our model results with those obtained from reanalysis data demonstrating that the latter exhibit a similar qualitative behavior in terms of multiscale dimensions and the existence of a scale dependency of the statistics of the phase-space density of points for different regions, which is related to the different drivers and processes occurring at different timescales in the coupled atmosphere–ocean system. Our approach can therefore be used to diagnose the strength of coupling in real applications.

2021, Armstrong McKay, David I., Cornell, Sarah E., Richardson, Katherine, Rockström, Johan

The Earth's oceans are one of the largest sinks in the Earth system for anthropogenic CO2 emissions, acting as a negative feedback on climate change. Earth system models project that climate change will lead to a weakening ocean carbon uptake rate as warm water holds less dissolved CO2 and as biological productivity declines. However, most Earth system models do not incorporate the impact of warming on bacterial remineralisation and rely on simplified representations of plankton ecology that do not resolve the potential impact of climate change on ecosystem structure or elemental stoichiometry. Here, we use a recently developed extension of the cGEnIE (carbon-centric Grid Enabled Integrated Earth system model), ecoGEnIE, featuring a trait-based scheme for plankton ecology (ECOGEM), and also incorporate cGEnIE's temperature-dependent remineralisation (TDR) scheme. This enables evaluation of the impact of both ecological dynamics and temperature-dependent remineralisation on particulate organic carbon (POC) export in response to climate change. We find that including TDR increases cumulative POC export relative to default runs due to increased nutrient recycling (+∼1.3 %), whereas ECOGEM decreases cumulative POC export by enabling a shift to smaller plankton classes (−∼0.9 %). However, interactions with carbonate chemistry cause opposite sign responses for the carbon sink in both cases: TDR leads to a smaller sink relative to default runs (−∼1.0 %), whereas ECOGEM leads to a larger sink (+∼0.2 %). Combining TDR and ECOGEM results in a net strengthening of POC export (+∼0.1 %) and a net reduction in carbon sink (−∼0.7 %) relative to default. These results illustrate the degree to which ecological dynamics and biodiversity modulate the strength of the biological pump, and demonstrate that Earth system models need to incorporate ecological complexity in order to resolve non-linear climate–biosphere feedbacks.