Deterministic sensitivity analysis for a model for flow in porous media

A deterministic method for sensitivity analysis is developed and applied to a mathematical model for the simulation of flow in porous media. The method is based on the singular value decomposition (SVD) of the Jacobian matrix of the model. It is a local approach to sensitivity analysis providing a hierarchical classification of the directions in both the input space and of those in the output space reflecting the degree of sensitiveness of the latter to the former. Its low computational cost, in comparison with that of statistical approaches, allows the study of the variability of the results of the sensitivity analysis due to the variations of the input parameters of the model, and thus it can provide a quality criterion for the validity of more classical probabilistic global approaches. For the example treated here, however, this variability is weak, and deterministic and statistical methods yield similar sensitivity results.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Balancing deformation in NW Borneo: Quantifying plate-scale vs. gravitational tectonics in a delta and deepwater fold-thrust belt system

Recent GPS measurements demonstrate that NW Borneo undergoes 4–6 mm of plate-scale shortening a year, which is not accommodated by plate-scale structures. The only geological structure in NW Borneo described to accommodate on-going shortening is the Baram Delta System located on the outer shelf to basin floor. Delta toe fold-thrust belts are commonly thought to be caused by margin-normal compressional stresses generated by margin-parallel upslope gravitational extension.     

Wave climate variability in the North-East Atlantic Ocean over the last six decades

Ocean surface gravity waves play a major role in many engineering and environmental problems, both in the open ocean and in coastal zones. Therefore, it is essential to improve our knowledge on spatial and temporal variability of wave climate. This study aims at investigating this variability in the North-East Atlantic Ocean (25°W–0°W and 30°N–60° N), using a 57-year hindcast (1953–2009) obtained with a spectral wave model forced with reanalysis wind fields. The hindcast analysis reveals firstly strong seasonal fluctuations of wave climate, with winters characterized by large and long-period waves of mean direction spreading from south-west to north-west, and summers characterized by smaller and shorter-period waves originating from norther directions. From northern (55°N) to southern (35°N) latitudes, the significant wave height (Hs) decreases by roughly 40%, the mean wave direction (Mwd) rotates clockwise by about 25% while the peak period (Tp) only grows by 5%. These three parameters also exhibit a strong inter-annual variability, particularly when winter-means (from 1st of December to 1st of April) are considered. Linear trend analysis over the studied period shows spatially variable long-term trends, with a significant increase of Hs (up to 0.02 m yr^?1) and a counterclockwise shift of Mwd (up to ?0.1° yr^?1) at northern latitude, contrasting with a fairly constant trend for Hs and a clockwise shift of Mwd (up to +0.15° yr^?1) at southern latitudes. Long-term trends of Tp are less significant, with still a slight increase in the north-eastern part of the study area (up to +0.01 s yr^?1). Eventually, a comparison between the inter-annual variability of the winter-means of the three selected wave parameters and the North Atlantic Oscillation (NAO) reveals: (1) a strong positive correlation between Hs and the NAO index at northern latitudes (correlation coefficient up to R = 0.91) and a significant negative correlation at southern latitudes (up to R = ?0.6); (2) no significant correlation for Mwd north of 40°N and a clear positive correlation southward of 40°N (up to R = 0.8) and (3) a northward increasing correlation for Tp (up to R = 0.8). Long-term trends for Hs, Mwd and Tp are finally explained by a significant increase in the NAO index over the studied period.     

Paleogene evolution and demise of the proto‐Paratethys Sea in Central Asia (Tarim and Tajik basins): Role of intensified tectonic activity at ca. 41 Ma

The proto‐Paratethys Sea covered a vast area extending from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous and early Paleogene. Climate modelling and proxy studies suggest that Asian aridification has been governed by westerly moisture modulated by fluctuations of the proto‐Paratethys Sea. Transgressive and regressive episodes of the proto‐Paratethys Sea have been previously recognized but their timing, extent and depositional environments remain poorly constrained. This hampers understanding of their driving mechanisms (tectonic and/or eustatic) and their contribution to Asian aridification. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed palaeoenvironmental analysis for the Paleogene proto‐Paratethys Sea incursions in the Tajik and Tarim basins. This enables us to identify the major drivers of marine fluctuations and their potential consequences on Asian aridification. A major regional restriction event, marked by the exceptionally thick (≤ 400 m) shelf evaporites is assigned a Danian‐Selandian age (ca. 63–59 Ma) in the Aertashi Formation. This is followed by the largest recorded proto‐Paratethys Sea incursion with a transgression estimated as early Thanetian (ca. 59–57 Ma) and a regression within the Ypresian (ca. 53–52 Ma), both within the Qimugen Formation. The transgression of the next incursion in the Kalatar and Wulagen formations is now constrained as early Lutetian (ca. 47–46 Ma), whereas its regression in the Bashibulake Formation is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian–Priabonian (ca. 39.7–36.7 Ma). We interpret the long‐term westward retreat of the proto‐Paratethys Sea starting at ca. 41 Ma to be associated with far‐field tectonic effects of the Indo‐Asia collision and Pamir/Tibetan plateau uplift. Short‐term eustatic sea level transgressions are superimposed on this long‐term regression and seem coeval with the transgression events in the other northern Peri‐Tethyan sedimentary provinces for the 1st and 2nd sea incursions. However, the 3rd sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto‐Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto‐Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Palaeomagnetism of the Permian Volcanic Rocks of Moissey (French Jura): Implications For the Palaeofield and Tectonic Evolution

Rocks from the Massif de la Serre in the French Jura (latitude: 47.3°N longitude: 5.6°E) belonging to an ignimbritic assemblage dominated by vitrophyric rhyolites, and whose age of formation is probably Permian (Autunian to Saxonian) have been studied by applying thermal and alternating field demagnetization. the characteristic magnetization has a mean direction derived from 89 samples of D= 170°, I = - 16°, k = 26.2°, α₉₅= 3° and a corresponding north palaeopole at 41°N, 172°E, A ₉₅= 5°. the pole, which is very close to the Permian European poles, can thus be considered as a new contribution. Some samples are found to carry a unique normal polarity magnetization, others carry both normal and reverse polarities. It therefore seems that, similar to Permian series in the USSR, these west European rocks have registered a normal event in the Kiaman interval. From a structural point of view, we may conclude that during the Alpine tectonic phases the Massif de la Serre has not been subjected to substantial rotation.     

Co-use of existing scenario sets to extend and quantify the shared socioeconomic pathways

More often than not, assessments of future climate risks are based on future climatic conditions superimposed on current socioeconomic conditions only. The new IPCC-guided alternative global development trends, the shared socioeconomic pathways (SSPs), have the potential to enhance the integration of future socioeconomic conditions—in the form of socioeconomic scenarios—within assessments of future climate risks. Being global development pathways, the SSPs lack regional and sectoral details. To increase their suitability in sectoral and/or regional studies and their relevance for local stakeholders, the SSPs have to be extended. We propose here a new method to extend the SSPs that makes use of existing scenario studies, the (re)use of which has been underestimated so far. Our approach lies in a systematic matching of multiple scenario sets that facilitates enrichment of the global SSPs with regional and sectoral information, in terms of both storylines and quantitative projections. We apply this method to develop extended SSPs of human vulnerability in Europe and to quantify them for a number of key indicators at the sub-national level up to 2050, based on the co-use of the matched scenarios’ quantitative outputs. Results show that such a method leads to internally consistent extended SSPs with detailed and highly quantified narratives that are tightly linked to global contexts. This method also provides multiple entry points where the relevance of scenarios to local stakeholders can be tested and strengthened. The extended SSPs can be readily employed to explore future populations’ vulnerability to climate hazards under varying levels of socioeconomic development.     

Anharmonicity and high-temperature heat capacity of crystals: the examples of Ca2GeO4, Mg2GeO4 and CaMgGeO4 olivines

Drop-calorimetry determinations of the isobaric heat capacity (C_P) of Mg₂GeO₄, Ca₂GeO₄ and CaMgGeO₄ have been made up to 1700 K. The thermal expansion coefficient (α) of these olivine germanates has been determined from high-temperature X-ray measurements up to 1500 K. From these measurements and available compressibility data, one calculates that the isochoric heat capacity (C_V) exceeds the harmonic limit of Dulong and Petit above 1000–1200 K. Such an intrinsic anharmonic behaviour can be accounted for by introducing anharmonic parameters a_i=(? ln v _i/?T)_V in vibrational modelling of C_V. These parameters are calculated from pressure and temperature shifts of the vibrational frequencies as measured by Raman spectroscopy up to 10 GPa at room temperature and up to 1300 K at 1 bar. A comparison of the Raman spectra of the three germanates with those of natural olivines justifies once again the use of germanates as silicate analogues. Extensive Ca,Mg disordering likely takes place in CaMgGeO₄, beginning at about 1100 K and leading to unusually high increases of the heat capacity and thermal expansion coefficient.