Search Results
Local stratopause temperature variabilities and their embedding in the global context
The stratopause is by definition the transition between the stratosphere and mesosphere. During winter the circulation at mid-latitudes and high latitudes in the stratosphere is mainly driven by quasi-stationary planetary waves (PWs), while the circulation in the mesosphere is mainly driven by gravity waves (GWs). The question arises of whether PWs or GWs dominate the variability of the stratopause. The most famous and dramatic variability of the middle atmosphere is a sudden stratospheric warming (SSW) generated by PWs interacting with the polar vortex. A similar phenomenon but smaller in magnitude and more regional is stratopause temperature enhancements (STEs) initially observed by local measurements and generated by breaking PWs. Thus it seems that PWs dominate the variability of the stratopause. In this study we want to quantify to which extent quasi-stationary PWs contribute to the stratopause variability. To do that we combine local lidar observations at Kühlungsborn (54∘ N, 11∘ E) and Andenes (69∘ N, 16∘ E) with global MERRA-2 reanalysis data bringing the local variability of the stratopause into the global context. Therefore we compare the temperature time series at Kühlungsborn and Andenes at 2 hPa, the altitude where STEs maximize, with characteristics (amplitude and phase) of PWs with wave numbers 1, 2 and 3. We found that for Kühlungsborn and Andenes 98 % of the local day-to-day variability of the stratopause can be explained by the variability of PWs with wave number 1, 2 and 3. Thus, the winter stratopause day-to-day variability is highly dominated by the variability of PWs.
Earth system data cubes unravel global multivariate dynamics
Understanding Earth system dynamics in light of ongoing human intervention and dependency remains a major scientific challenge. The unprecedented availability of data streams describing different facets of the Earth now offers fundamentally new avenues to address this quest. However, several practical hurdles, especially the lack of data interoperability, limit the joint potential of these data streams. Today, many initiatives within and beyond the Earth system sciences are exploring new approaches to overcome these hurdles and meet the growing interdisciplinary need for data-intensive research; using data cubes is one promising avenue. Here, we introduce the concept of Earth system data cubes and how to operate on them in a formal way. The idea is that treating multiple data dimensions, such as spatial, temporal, variable, frequency, and other grids alike, allows effective application of user-defined functions to co-interpret Earth observations and/or model-data integration. An implementation of this concept combines analysis-ready data cubes with a suitable analytic interface. In three case studies, we demonstrate how the concept and its implementation facilitate the execution of complex workflows for research across multiple variables, and spatial and temporal scales: (1) summary statistics for ecosystem and climate dynamics; (2) intrinsic dimensionality analysis on multiple timescales; and (3) model-data integration. We discuss the emerging perspectives for investigating global interacting and coupled phenomena in observed or simulated data. In particular, we see many emerging perspectives of this approach for interpreting large-scale model ensembles. The latest developments in machine learning, causal inference, and model-data integration can be seamlessly implemented in the proposed framework, supporting rapid progress in data-intensive research across disciplinary boundaries. © 2020 Institute of Electrical and Electronics Engineers Inc.. All rights reserved.
Earth system modeling with endogenous and dynamic human societies: the copan:CORE open World–Earth modeling framework
Analysis of Earth system dynamics in the Anthropocene requires explicitly taking into account the increasing magnitude of processes operating in human societies, their cultures, economies and technosphere and their growing feedback entanglement with those in the physical, chemical and biological systems of the planet. However, current state-of-the-art Earth system models do not represent dynamic human societies and their feedback interactions with the biogeophysical Earth system and macroeconomic integrated assessment models typically do so only with limited scope. This paper (i) proposes design principles for constructing world-Earth models (WEMs) for Earth system analysis of the Anthropocene, i.e., models of social (world)-ecological (Earth) coevolution on up to planetary scales, and (ii) presents the copan:CORE open simulation modeling framework for developing, composing and analyzing such WEMs based on the proposed principles. The framework provides a modular structure to flexibly construct and study WEMs. These can contain biophysical (e.g., carbon cycle dynamics), socio-metabolic or economic (e.g., economic growth or energy system changes), and sociocultural processes (e.g., voting on climate policies or changing social norms) and their feedback interactions, and they are based on elementary entity types, e.g., grid cells and social systems. Thereby, copan:CORE enables the epistemic flexibility needed for contributions towards Earth system analysis of the Anthropocene given the large diversity of competing theories and methodologies used for describing socio-metabolic or economic and sociocultural processes in the Earth system by various fields and schools of thought. To illustrate the capabilities of the framework, we present an exemplary and highly stylized WEM implemented in copan:CORE that illustrates how endogenizing sociocultural processes and feedbacks such as voting on climate policies based on socially learned environmental awareness could fundamentally change macroscopic model outcomes. © Author(s) 2020.
Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)
The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4mm of global sea level rise, with a standard deviation of 3.7mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 °C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles.We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade. © Author(s) 2020.
A multi-model analysis of teleconnected crop yield variability in a range of cropping systems
Climate oscillations are periodically fluctuating oceanic and atmospheric phenomena, which are related to variations in weather patterns and crop yields worldwide. In terms of crop production, the most widespread impacts have been observed for the El Niño-Southern Oscillation (ENSO), which has been found to impact crop yields on all continents that produce crops, while two other climate oscillations - the Indian Ocean Dipole (IOD) and the North Atlantic Oscillation (NAO) - have been shown to especially impact crop production in Australia and Europe, respectively. In this study, we analyse the impacts of ENSO, IOD, and NAO on the growing conditions of maize, rice, soybean, and wheat at the global scale by utilising crop yield data from an ensemble of global gridded crop models simulated for a range of crop management scenarios. Our results show that, while accounting for their potential co-variation, climate oscillations are correlated with simulated crop yield variability to a wide extent (half of all maize and wheat harvested areas for ENSO) and in several important crop-producing areas, e.g. in North America (ENSO, wheat), Australia (IOD and ENSO, wheat), and northern South America (ENSO, soybean). Further, our analyses show that higher sensitivity to these oscillations can be observed for rainfed and fully fertilised scenarios, while the sensitivity tends to be lower if crops were to be fully irrigated. Since the development of ENSO, IOD, and NAO can potentially be forecasted well in advance, a better understanding about the relationship between crop production and these climate oscillations can improve the resilience of the global food system to climate-related shocks. © 2020 American Institute of Physics Inc.. All rights reserved.
Multiphase, decoupled faulting in the southern German Molasse Basin – evidence from 3-D seismic data
We use three-dimensional seismic reflection data from the southern German Molasse Basin to investigate the structural style and evolution of a geometrically decoupled fault network in close proximity to the Alpine deformation front. We recognise two fault arrays that are vertically separated by a clay-rich layer – lower normal faults and upper normal and reverse faults. A frontal thrust fault partially overprints the upper fault array. Analysis of seismic stratigraphy, syn-kinematic strata, throw distribution, and spatial relationships between faults suggest a multiphase fault evolution: (1) initiation of the lower normal faults in the Upper Jurassic carbonate platform during the early Oligocene, (2) development of the upper normal faults in the Cenozoic sediments during the late Oligocene, and (3) reverse reactivation of the upper normal faults and thrusting during the mid-Miocene. These distinct phases document the evolution of the stress field as the Alpine orogen propagated across the foreland. We postulate that interplay between the horizontal compression and vertical stresses due to the syn-sedimentary loading resulted in the intermittent normal faulting. The vertical stress gradients within the flexed foredeep defined the independent development of the upper faults above the lower faults, whereas mechanical behaviour of the clay-rich layer precluded the subsequent linkage of the fault arrays. The thrust fault must have been facilitated by the reverse reactivation of the upper normal faults, as its maximum displacement and extent correlate with the occurrence of these faults. We conclude that the evolving tectonic stresses were the primary mechanism of fault activation, whereas the mechanical stratigraphy and pre-existing structures locally governed the structural style.
Tropical and mid-latitude teleconnections interacting with the Indian summer monsoon rainfall: a theory-guided causal effect network approach
The alternation of active and break phases in Indian summer monsoon (ISM) rainfall at intraseasonal timescales characterizes each ISM season. Both tropical and mid-latitude drivers influence this intraseasonal ISM variability. The circumglobal teleconnection observed in boreal summer drives intraseasonal variability across the mid-latitudes, and a two-way interaction between the ISM and the circumglobal teleconnection pattern has been hypothesized. We use causal discovery algorithms to test the ISM circumglobal teleconnection hypothesis in a causal framework. A robust causal link from the circumglobal teleconnection pattern and the North Atlantic region to ISM rainfall is identified, and we estimate the normalized causal effect (CE) of this link to be about 0.2 (a 1 standard deviation shift in the circumglobal teleconnection causes a 0.2 standard deviation shift in the ISM rainfall 1 week later). The ISM rainfall feeds back on the circumglobal teleconnection pattern, however weakly. Moreover, we identify a negative feedback between strong updraft located over India and the Bay of Bengal and the ISM rainfall acting at a biweekly timescale, with enhanced ISM rainfall following strong updraft by 1 week. This mechanism is possibly related to the boreal summer intraseasonal oscillation. The updraft has the strongest CE of 0.5, while the Madden–Julian oscillation variability has a CE of 0.2–0.3. Our results show that most of the ISM variability on weekly timescales comes from these tropical drivers, though the mid-latitude teleconnection also exerts a substantial influence. Identifying these local and remote drivers paves the way for improved subseasonal forecasts.
Incremental improvements of 2030 targets insufficient to achieve the Paris Agreement goals
Current global mitigation ambition up to 2030 under the Paris Agreement, reflected in the National Determined Contributions (NDCs), is insufficient to achieve the agreement's 1.5 °C long-term temperature limit. As governments are preparing new and updated NDCs for 2020, the question as to how much collective improvement is achieved is a pivotal one for the credibility of the international climate regime. The recent Special Report on Global Warming of 1.5 °C by the Intergovernmental Panel on Climate Change has assessed a wide range of scenarios that achieve the 1.5 °C limit. Those pathways are characterised by a substantial increase in near-term action and total greenhouse gas (GHG) emission levels about 50 % lower than what is implied by current NDCs. Here we assess the outcomes of different scenarios of NDC updating that fall short of achieving this 1.5 °C benchmark. We find that incremental improvements in reduction targets, even if achieved globally, are insufficient to align collective ambition with the goals of the Paris Agreement. We provide estimates for global mean temperature increase by 2100 for different incremental NDC update scenarios and illustrate climate impacts under those median scenarios for extreme temperature, long-term sea-level rise and economic damages for the most vulnerable countries. Under the assumption of maintaining ambition as reflected in current NDCs up to 2100 and beyond, we project a reduction in the gross domestic product (GDP) in tropical countries of around 60 % compared to a no-climate-change scenario and median long-term sea-level rise of close to 2 m in 2300. About half of these impacts can be avoided by limiting warming to 1.5 °C or below. Scenarios of more incremental NDC improvements do not lead to comparable reductions in climate impacts. An increase in aggregated NDC ambition of big emitters by 33 % in 2030 does not reduce presented climate impacts by more than about half compared to limiting warming to 1.5 °C. Our results underscore that a transformational increase in 2030 ambition is required to achieve the goals of the Paris Agreement and avoid the worst impacts of climate change. © 2020 SPIE. All rights reserved.
Bilinear pressure diffusion and termination of bilinear flow in a vertically fractured well injecting at constant pressure
This work studies intensively the flow in fractures with finite hydraulic conductivity intersected by a well injecting or producing at constant pressure, either during an injection or production well test or the operation of a production well. Previous investigations showed that for a certain time the reciprocal of flow rate is proportional to the fourth root of time, which is characteristic of the flow regime known as bilinear flow. Using a 2D numerical model, we demonstrated that during the bilinear flow regime the transient propagation of isobars along the fracture is proportional to the fourth root of time. Moreover, we present relations to calculate the termination time of bilinear flow under constant injection or production well pressure as well as an expression for the bilinear hydraulic diffusivity of fractures with finite hydraulic conductivity. To determine the termination of bilinear flow regime, two different methods were used: (a) numerically measuring the transient flow rate in the well and (b) analyzing the propagation of isobars along the fracture. Numerical results show that for low dimensionless fracture conductivities the transition from bilinear flow to another flow regime (e.g., pseudo-radial flow) occurs before the pressure front reaches the fracture tip, and for high dimensionless fracture conductivities it occurs when the pressure front arrives at the fracture tip. Hence, this work complements and advances previous research on the interpretation and evaluation of well test analysis under different reservoir conditions. Our results aim to improve the understanding of the hydraulic diffusion in fractured geologic media, and as a result they can be utilized for the interpretation of hydraulic tests, for example to estimate the fracture length.