Refine
Has Fulltext
- no (240) (remove)
Year of publication
- 2019 (240) (remove)
Document Type
- Article (216)
- Other (9)
- Review (7)
- Doctoral Thesis (6)
- Habilitation Thesis (2)
Is part of the Bibliography
- yes (240)
Keywords
- climate (4)
- Holocene (3)
- paleoclimate (3)
- remote sensing (3)
- AMOC (2)
- Asia (2)
- Central Andes (2)
- Climate (2)
- Climate change (2)
- Climate change adaptation (2)
Institute
- Institut für Geowissenschaften (240) (remove)
Stable carbon isotopes of sediment organic matter (delta C-13(OM)) are widely applied in paleoenvironmental studies. Interpretations of delta C-13(OM), however, remain challenging and factors that influence delta C-13(OM) may not apply across all lakes. Common explanations for stratigraphic shifts in delta C-13(OM) include changes in lake productivity or changes in inputs of allochthonous OM. We investigated the influence of different oxygen conditions (oxic versus anoxic) on the delta C-13(OM) values in the sediments of Lake Tiefer See. We analysed (1) a long sediment core from the deepest part of the lake, (2) two short, sediment-water interface cores from shallower water depths, and (3) OM in the water column, i.e. from sediment traps. Fresh OM throughout the entire water column showed a relatively constant delta C-13(OM) value of approximately -30.5 parts per thousand. Similar values, about -31 parts per thousand, were obtained for well-varved sediments in both the long and short, sediment-water interface cores. In contrast, delta C-13(OM) values from non-varved sediments in all cores were significantly less negative (-29 parts per thousand). The delta C-13(OM) values in the sediment-water interface cores from different water depths differ for sediments of the same age, if oxygen conditions at the time of deposition were different at these sites, as suggested by the state of varve preservation. Sediments deposited from AD 1924 to 1980 at 62m water depth are varved and exhibit delta C-13(OM) values around -31 parts per thousand, whereas sediments of the same age in the core from 35m water depth are not varved and show less negative delta C-13(OM) values of about -29 parts per thousand. The relation between varve occurrence and delta C-13(OM) values suggests that delta C-13(OM) is associated with oxygen conditions because varve preservation depends on hypolimnetic anoxia. A mechanism that likely influences delta C-13(OM) is selective degradation of OM under oxic conditions, such that organic components with more negative delta C-13(OM) are preferably decomposed, leading to less negative delta C-13(OM) values in the remaining, undegraded OM pool. Greater decomposition of OM in non-varved sediments is supported by lower TOC concentrations in these deposits (similar to 5%) compared to well-varved sediments (similar to 15%). Even in lakes that display small variations in productivity and terrestrial OM input through time, large spatial and temporal differences in hypolimnetic oxygen concentrations may be an important factor controlling sediment delta C-13(OM).
Late Quaternary landscapes of unglaciated Beringia were largely shaped by ice-wedge polygon tundra. Ice Complex (IC) strata preserve such ancient polygon formations. Here we report on the Yukagir IC from Bol'shoy Lyakhovsky Island in northeastern Siberia and suggest that new radioisotope disequilibria (230Th/U) dates of the Yukagir IC peat confirm its formation during the Marine Oxygen Isotope Stage (MIS) 7a–c interglacial period. The preservation of the ice-rich Yukagir IC proves its resilience to last interglacial and late glacial–Holocene warming. This study compares the Yukagir IC to IC strata of MIS 5, MIS 3, and MIS 2 ages exposed on Bol'shoy Lyakhovsky Island. Besides high intrasedimental ice content and syngenetic ice wedges intersecting silts, sandy silts, the Yukagir IC is characterized by high organic matter (OM) accumulation and low OM decomposition of a distinctive Drepanocladus moss-peat. The Yukagir IC pollen data reveal grass-shrub-moss tundra indicating rather wet summer conditions similar to modern ones. The stable isotope composition of Yukagir IC wedge ice is similar to those of the MIS 5 and MIS 3 ICs pointing to similar atmospheric moisture generation and transport patterns in winter. IC data from glacial and interglacial periods provide insights into permafrost and climate dynamics since about 200 ka.
Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments, studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for 3 to 9 years and analysed for stable isotopes (delta O-18 and delta H-2), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two-component lumped parameter models (LPMs) that are solved using delta O-18 records. The study illustrates that two-component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low-permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration, short mTTs (<40 weeks) were determined for sites closer to the river (<200 m), whereas for the precipitation infiltration, the mTTs were longer (>80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface-groundwater interaction can be conceptualized by exploiting two-component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.
Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M-w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.
We derive a set of regional ground-motion prediction equations (GMPEs) in the Fourier amplitude spectra (FAS-GMPE) and in the spectral acceleration (SA-GMPE) domains for the purpose of interpreting the between-event residuals in terms of source parameter variability. We analyze a dataset of about 65,000 recordings generated by 1400 earthquakes (moment magnitude 2: 5 <= M-w <= 6: 5, hypocentral distance R-hypo <= 150 km) that occurred in central Italy between January 2008 and October 2017. In a companion article (Bindi, Spallarossa, et al., 2018), the nonparametric acceleration source spectra were interpreted in terms of omega-square models modified to account for deviations from a high-frequency flat plateau through a parameter named k(source). Here, the GMPEs are derived considering the moment (M-w), the local (M-L), and the energy (M-e) magnitude scales, and the between-event residuals are computed as random effects. We show that the between-event residuals for the FAS-GMPE implementing M-w are correlated with stress drop, with correlation coefficients increasing with increasing frequency up to about 10 Hz. Contrariwise, the correlation is weak for the FAS-GMPEs implementing M-L and M-e, in particular between 2 and 5 Hz, where most of the corner frequencies lie. At higher frequencies, all models show a strong correlation with k(source). The correlation with the source parameters reflects in a different behavior of the standard deviation tau of the between-event residuals with frequency. Although tau is smaller for the FAS-GMPE using M-w below 1.5 Hz, at higher frequencies, the model implementing either M-L or M-e shows smaller values, with a reduction of about 30% at 3 Hz (i.e., from 0.3 for M-w to 0.1 for M-L). We conclude that considering magnitude scales informative for the stress-drop variability allows to reduce the between-event variability with a significant impact on the hazard assessment, in particular for studies in which the ergodic assumption on site is removed.
Flood disasters severely impact human subjective well-being (SWB). Nevertheless, few studies have examined the influence of flood events on individual well-being and how such impacts may be limited by flood protection measures. This study estimates the long term impacts on individual subjective well-being of flood experiences, individual subjective flood risk perceptions, and household flood preparedness decisions. These effects are monetised and placed in context through a comparison with impacts of other adverse events on well-being. We collected data from households in flood-prone areas in France. The results indicate that experiencing a flood has a large negative impact on subjective well-being that is incompletely attenuated over time. Moreover, individuals do not need to be directly affected by floods to suffer SWB losses since subjective well-being is lower for those who expect their flood risk to increase or who have seen a neighbour being flooded. Floodplain inhabitants who prepared for flooding by elevating their home have a higher subjective well-being. A monetisation of the aforementioned well-being impacts shows that a flood requires Euro150,000 in immediate compensation to attenuate SWB losses. The decomposition of the monetised impacts of flood experience into tangible losses and intangible effects on SWB shows that intangible effects are about twice as large as the tangible direct monetary flood losses. Investments in flood protection infrastructure may be under funded if the intangible SWB benefits of flood protection are not taken into account.
Increased landslide activity on forested hillslopes following two recent volcanic eruptions in Chile
(2019)
Large explosive eruptions can bury landscapes beneath thick layers of tephra. Rivers subsequently overloaded with excess pyroclastic sediments have some of the highest reported specific sediment yields. Much less is known about how hillslopes respond to tephra loads. Here, we report a pulsed and distinctly delayed increase in landslide activity following the eruptions of the Chaiten (2008) and Puyehue-Cordon Caulle (2011) volcanoes in southern Chile. Remote-sensing data reveal that land-slides clustered in densely forested hillslopes mostly two to six years after being covered by tephra. This lagged instability is consistent with a gradual loss of shear strength of decaying tree roots in areas of high tephra loads. Surrounding areas with comparable topography, forest cover, rainfall and lithology maintained landslide rates roughly ten times lower. The landslides eroded the landscape by up to 4.8 mm on average within 30 km of both volcanoes, mobilizing up to 1.6 MtC at rates of about 265 tC km(-2) yr(-1). We suggest that these yields may reinforce the elevated river loads of sediment and organic carbon in the decade after the eruptions. We recommend that studies of post-eruptive mass fluxes and hazards include lagged landslide responses of tephra-covered forested hillslopes, to avoid substantial underestimates.
Caldera unrest can lead to major volcanic eruptions. Analysis of subtle subsidence or inflation at calderas helps understanding of their subsurface volcanic processes and related hazards. Several subsiding calderas have shown similar patterns of ground deformation composed of broad subsidence affecting the entire volcanic edifice and stronger localized subsidence focused inside the caldera. Physical models of internal deformation sources used to explain these observations typically consist of two magma reservoirs at different depths in an elastic half-space. However, such models ignore important subsurface structures, such as ring faults, that may influence the deformation pattern. Here we use both analog subsidence experiments and boundary element modeling to study the three-dimensional geometry and kinematics of caldera subsidence processes, evolving from an initial downsag to a later collapse stage. We propose that broad subsidence is mainly caused by volume decrease within a single magma reservoir, whereas buried ring-fault activity localizes the deformation within the caldera. Omitting ring faulting in physical models of subsiding calderas and using multiple point/sill-like sources instead can result in erroneous estimates of magma reservoir depths and volume changes. (C) 2019 Elsevier B.V. All rights reserved.
Millennial-scale Dansgaard Oeschger (DO) variability at northern high latitudes has influenced climatic and environmental conditions in the Mediterranean during the last glacial period. There is evidence that the hemispheric transmission of the DO variability occurred at the end of DO event 25; however, the exact timing and the trigger that activated the environmental response in the Mediterranean remains incompletely understood. Here, we provide evidence that the clear millennial-scale teleconnection between Greenland and the Mediterranean started at similar to 111.4 ka BP and was initiated by a sub-millennial scale cooling in Greenland (GI-25b). High-resolution sediment proxies and the pollen record of Lago Grande di Monticchio (MON), Italy, reflect climatic instability during the last millennium of the last interglacial, which was characterised by a first and short cooling episode (MON 1) at 111.44 +/- 0.69 ka BP, coinciding with the Greenland cold sub-event GI-25b in duration and timing (within dating uncertainties). MON and Greenland (NorthGRIP ice core) also agree in recording a subsequent warm rebound phase that abruptly culminated in the stadial MON 2/GS-25, marking the transition into the last glacial period. Our results show that the GI-25b triggered an early environmental response at MON to centennial-scale climate change in Greenland as a prelude to the millennial-scale teleconnection that was maintained during the glacial period.
The 2015 Paris Agreement (PA) has been widely hailed as a diplomatic triumph and a breakthrough in global climate cooperation. However, it is commonly accepted that the PA's collective goal—keeping global warming “well below” 2°C above preindustrial levels—remains ambitious. Making matters even more challenging, in 2017, global CO2 emissions resumed growth after 3 years of near standstill. In 2018, this growth accelerated. It is therefore extremely important that the PA's institutional architecture meet expectations concerning its ability to induce member countries to promise and deliver emissions reductions. This study offers a review of the rapidly growing literature on the PA, to assess its strengths and weaknesses, its significance, and its prospects. We focus on evaluations of its institutional structure and its ability to induce member countries to implement policies. We frame the issues as a trilemma: the challenge of simultaneously satisfying all three main conditions for effectiveness—broad participation, deep commitments, and satisfactory compliance rates. Based on our review, we conclude that the key challenge for the PA will likely be to facilitate sufficiently fast ratcheting‐up of nationally determined contributions, while keeping compliance rates high.
A comprehensive description of electromagnetic processes related to equatorial plasma depletions (EPDs) is essential for understanding their evolution and day-to-day variability. Recently, field-aligned currents (FACs) flowing at both western and eastern edges of EPDs were observed to be interhemispheric rather than anti-parallel about the dip equator, as suggested by previous theoretical studies. In this paper, we investigate the spatial and temporal variability of the FACs orientation using simultaneous measurements of electron density and magnetic field gathered by ESA’s Swarm constellation mission. By using empirical models, we assess the role of the Pedersen conductance in the preference of the FACs to close either in the northern or southern magnetic hemisphere. Here we show that the closure of the FACs agrees with an electrostatic regime determined by a hemispherical asymmetry of the Pedersen conductance. That is, the EPD-related FACs close at lower altitudes in the hemisphere with the highest conductivity. The evidence of this conclusion stands on the general agreement between the longitudinal and seasonal variability of both the conductivity and the FACs orientation.
Introducing PebbleCounts
(2019)
Grain-size distributions are a key geomorphic metric of gravel-bed rivers. Traditional measurement methods include manual counting or photo sieving, but these are achievable only at the 1–10 ㎡ scale. With the advent of drones and increasingly high-resolution cameras, we can now generate orthoimagery over hectares at millimeter to centimeter resolution. These scales, along with the complexity of high-mountain rivers, necessitate different approaches for photo sieving. As opposed to other image segmentation methods that use a watershed approach, our open-source algorithm, PebbleCounts, relies on k-means clustering in the spatial and spectral domain and rapid manual selection of well-delineated grains. This improves grain-size estimates for complex riverbed imagery, without post-processing. We also develop a fully automated method, PebbleCountsAuto, that relies on edge detection and filtering suspect grains, without the k-means clustering or manual selection steps. The algorithms are tested in controlled indoor conditions on three arrays of pebbles and then applied to 12 × 1 ㎡ orthomosaic clips of high-energy mountain rivers collected with a camera-on-mast setup (akin to a low-flying drone). A 20-pixel b-axis length lower truncation is necessary for attaining accurate grain-size distributions. For the k-means PebbleCounts approach, average percentile bias and precision are 0.03 and 0.09 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, and 0.07 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. The automatic approach has higher bias and precision of 0.13 and 0.15 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, but similar values of −0.06 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. For the automatic approach, only at best 70 % of the grains are correct identifications, and typically around 50 %. PebbleCounts operates most effectively at the 1 ㎡ patch scale, where it can be applied in ∼5–10 min on many patches to acquire accurate grain-size data over 10–100 ㎡ areas. These data can be used to validate PebbleCountsAuto, which may be applied at the scale of entire survey sites (102–104 ㎡ ). We synthesize results and recommend best practices for image collection, orthomosaic generation, and grain-size measurement using both algorithms.
In this paper, we examine the influence of the 27 October 2012, M-w 7.8 earthquake on landslide occurrence in the southern half of Haida Gwaii (formerly Queen Charlotte Islands), British Columbia, Canada. Our 1350 km(2) study area is undisturbed, primarily forested terrain that has not experienced road building or timber harvesting. Our inventory of landslide polygons is based on optical airborne and spaceborne images acquired between 2007 and 2018, from which we extracted and mapped 446 individual landslides (an average of 33 landslides per 100 km(2)). The landslide rate in years without major earthquakes averages 19.4 per year, or 1.4/100 km(2)/year, and the annual average area covered by non-seismically triggered landslides is 35 ha/year. The number of landslides identified in imagery closely following the 2012 earthquake, and probably triggered by it, is 244 or an average of about 18 landslides per 100 km(2). These landslides cover a total area of 461 ha. In the following years-2013-2016 and 2016-2018-the number of landslides fell, respectively, to 26 and 13.5 landslides per year. In non-earthquake years, most landslides happen on south-facing slopes, facing the prevailing winds. In contrast, during or immediately after the earthquake, up to 32% of the landslides occurred on north and northwest-facing slopes. Although we could not find imagery from the day after the earthquake, overview reconnaissance flights 10 and 16 days later showed that most of the landslides were recent, suggesting they were co-seismic.
Soil degradation is a severe and growing threat to ecosystem services globally. Soil loss is often nonlinear, involving a rapid deterioration from a stable eco-geomorphic state once a tipping point is reached. Soil loss thresholds have been studied at plot scale, but for landscapes, quantitative constraints on the necessary and sufficient conditions for tipping points are rare. Here, we document a landscape-wide eco-geomorphic tipping point at the edge of the Tibetan Plateau and quantify its drivers and erosional consequences. We show that in the upper Kali Gandaki valley, Nepal, soil formation prevailed under wetter conditions during much of the Holocene. Our data suggest that after a period of human pressure and declining vegetation cover, a 20% reduction of relative humidity and precipitation below 200 mm/year halted soil formation after 1.6 ka and promoted widespread gullying and rapid soil loss, with irreversible consequences for ecosystem services.
We analyzed the Lomas de Carabajal area in the intermontane Lerma valley of the Cordillera Oriental to assess the level of neotectonic activity in a densely populated region of northwestern Argentina. In this region, Plio-Pleistocene synorogenic conglomerates are deformed, locally associated with high-angle faults, and NNW-SSE oriented en-echelon folds characterized by wavelengths of < 1 km. The deformed Quaternary units follow the same pattern of deformation as observed in the underlying Neogene deposits; growth-strata geometries are observed near faults. This configuration is compatible with local left-lateral transpressional tectonism driven by ENE-WSW buttressing against the NW-oriented border of a Cretaceous extensional basin (Alemania sub-basin). Optically Stimulated Luminescence analysis of sandy-silty layers interbedded within the folded late Pleistocene conglomeratic sequence helps to determine uplift rates of 0.83-0.87 mm/a during the last 30-40 ka. Nearby the Lomas de Carabajal, a WNW-striking, 3-m-high fault scarp disrupts radiocarbon dated, 10-ka-old loessic deposits providing a Holocene mean uplift rate of 0.30 mm/a. Our data unambiguously show that shallow crustal deformation in the intermontane Lerma valley is ongoing; some of this deformation may be associated with seismicity. Our findings support the notion of temporally and spatially disparate deformation processes in the broken foreland of the northwestern Argentinean Andes.
Issue Despite their rather similar climatic conditions, eastern Eurasia and northern North America are largely covered by different plant functional types (deciduous or evergreen boreal forest) composed of larch or pine, spruce and fir, respectively. I propose that these deciduous and evergreen boreal forests represent alternative quasi-stable states, triggered by their different northern tree refugia that reflect the different environmental conditions experienced during the Last Glacial. Evidence This view is supported by palaeoecological and environmental evidence. Once established, Asian larch forests are likely to have stabilized through a complex vegetation-fire-permafrost soil-climate feedback system. Conclusion With respect to future forest developments, this implies that Asian larch forests are likely to be governed by long-term trajectories and are therefore largely resistant to natural climate variability on time-scales shorter than millennia. The effects of regional human impact and anthropogenic global warming might, however, cause certain stability thresholds to be crossed, meaning that irreversible transitions occur and resulting in marked consequences for ecosystem services on these human-relevant time-scales.
The marine sedimentary record contains unique information about the history of erosion, uplift and climate of the adjacent continent. Inverting this record has been the purpose of many numerical studies. However, limited attention has been given to linking continental erosion to marine sediment transport and deposition in large-scale surface process evolution models. Here we present a new numerical method for marine sediment transport and deposition that is directly coupled to a landscape evolution algorithm solving for the continental fluvial and hillslope erosion equations using implicit and O(N) algorithms. The new method takes into account the sorting of grain sizes (e.g., silt and sand) in the marine domain using a non-linear multiple grain-size diffusion equation and assumes that the sediment flux exported from the continental domain is proportional to the bathymetric slope. Specific transport coefficients and compaction factors are assumed for the two different grain sizes to simulate the stratigraphic architecture. The resulting set of equations is solved using an efficient (O(N) and implicit) algorithm. It can thus be used to invert stratigraphic geometries using a Bayesian approach that requires a large number of simulations. This new method is used to invert the sedimentary geometry of a natural example, the Ogooue Delta (Gabon), over the last similar to 5 Myr. The objective is to unravel the set of erosional histories of the adjacent continental domain compatible with the observed geometry of the offshore delta. For this, we use a Bayesian inversion scheme in which the misfit function is constructed by comparing four geometrical parameters between the natural and the simulated delta: the volume of sediments stored in the delta, the surface slope, the initial and the final shelf lengths. We find that the best-fit values of the transport coefficients for silt in the marine domain are in the range of 300 - 500 m(2)/yr, in agreement with previous studies on offshore diffusion. We also show that, in order to fit the sedimentary geometry, erosion rate on the continental domain must have increased by a factor of 6 to 8 since 5.3 Ma. (C) 2019 Elsevier B.V. All rights reserved.
Geomorphic strain markers accumulating the effects of many earthquake cycles help to constrain the mechanical behaviour of continental rift systems as well as the related seismic hazards. In the Corinth Rift (Greece), the unique record of onshore and offshore markers of Pleistocene similar to 100-ka climate cycles provides an outstanding possibility to constrain rift mechanics over a range of timescales. Here we use high-resolution topography to analyse the 3D geometry of a sequence of Pleistocene emerged marine terraces associated with flexural rift-flank uplift. We integrate this onshore dataset with offshore seismic data to provide a synoptic view of the flexural deformation across the rift. This allows us to derive an average slip rate of 4.5-9.0 mm.yr(-1) on the master fault over the past similar to 610 ka and an uplift/ subsidence ratio of 1:1.1-2.4. We reproduce the observed flexure patterns, using 3 and 5-layered lithospheric scale finite element models. Modelling results imply that the observed elastic flexure is produced by coseismic slip along 40-60 degrees planar normal faults in the elastic upper crust, followed by postseismic viscous relaxation occurring within the basal lower crust or upper mantle. We suggest that such a mechanism may typify rapid localised extension of continental lithosphere.
Isolated extreme habitats are ideally suited to investigate pivotal ecological processes such as niche use, local adaptation and dispersal. Extremophilic animals living in isolated habitats face the problem that dispersal is limited through the absence of suitable dispersal corridors, which in turn facilitates local adaptation. We used five rotifer isolates from extremely acidic mining lakes with a pH of below 3 as model organisms to test whether these isolates are acidotolerant or acidophilic, whether they survive and reproduce at their niche edges (here pH 2 and circum-neutral pH) and whether local adaptation has evolved. To evaluate potential dispersal limitation, we tested whether animals and their parthenogenetic eggs survive and remain reproductive or viable at unfavourable pH-conditions. All five isolates were acidophilic with a pH-optimum in the range of 4-6, which is well above the pH (< 3) of their lakes of origin. At unfavourable high pH, in four out of the five isolates parthenogenetic females produced a high number of non-viable eggs. Females and eggs produced at favourable pH (4) remained vital at an otherwise unfavourable pH of 7, indicating that for dispersal no acidic dispersal corridors are necessary. Common garden experiments revealed no clear evidence for local adaptation in any of the five isolates. Despite their acidophilic nature, all five isolates can potentially disperse via circum-neutral water bodies as long as their residence time is short, suggesting a broader dispersal niche than their realized niche. Local adaptation might have been hampered by the low population sizes of the rotifers in their isolated habitat and the short time span the mining lakes have existed.
During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.
Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.
In active mountain belts with steep terrain, bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high-magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analysing earthquake- and monsoon-induced landslide inventories across different timescales. With time series of 5 m satellite images over four main valleys in central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured similar to 5-8 times more landsliding than the pre-earthquake mean rate. The longterm size-frequency distribution of monsoon-induced landsliding (MIL) was derived from these inventories and from an inventory of landslides larger than similar to 0.1 km(2) that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size-frequency distribution for earthquakeinduced landsliding (EQIL). These two distributions are dominated by infrequent, large and giant landslides but under-predict an estimated Holocene frequency of giant landslides (> 1 km(3)) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquakes of variable magnitude and when considering that a shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of similar to 2 +/- 0.75 mm yr(-1) in agreement with most thermo-chronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size-frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling timescale required to adequately capture the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic Be-10 methods. This observation presents a strong caveat when interpreting spatial or temporal variability in erosion rates from this method. Thus, in areas where a very large, rare landslide contributes heavily to long-term erosion (as the Himalayas), we recommend Be-10 sample in catchments with source areas > 10 000 km(2) to reduce the method mean bias to below similar to 20 % of the long-term erosion.
The Arctic is directly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems, the subsistence economy of the local population, and the climate because of the transformation of organic matter into greenhouse gases. Yet, the patterns of sediment dispersal in the nearshore zone are not well known, because ships do not often reach shallow waters and satellite remote sensing is traditionally focused on less dynamic environments. The goal of this study is to use the extensive Landsat archive to investigate sediment dispersal patterns specifically on an exemplary Arctic nearshore environment, where field measurements are often scarce. Multiple Landsat scenes were combined to calculate means of sediment dispersal and sea surface temperature under changing seasonal wind conditions in the nearshore zone of Herschel Island Qikiqtaruk in the western Canadian Arctic since 1982. We use observations in the Landsat red and thermal wavebands, as well as a recently published water turbidity algorithm to relate archive wind data to turbidity and sea surface temperature. We map the spatial patterns of turbidity and water temperature at high spatial resolution in order to resolve transport pathways of water and sediment at the water surface. Our results show that these pathways are clearly related to the prevailing wind conditions, being ESE and NW. During easterly wind conditions, both turbidity and water temperature are significantly higher in the nearshore area. The extent of the Mackenzie River plume and coastal erosion are the main explanatory variables for sediment dispersal and sea surface temperature distributions in the study area. During northwesterly wind conditions, the influence of the Mackenzie River plume is negligible. Our results highlight the potential of high spatial resolution Landsat imagery to detect small-scale hydrodynamic processes, but also show the need to specifically tune optical models for Arctic nearshore environments.
Ground and space-based geomagnetic data were used in the investigation of the longitudinal, seasonal and lunar phase dependence of the equatorial counter electrojet (CEJ) occurrence in the Peruvian, Brazilian, African, Indian and Philippine sectors during geomagnetically quiet days from the solar cycle 24 (2008 to 2018). We found that CEJ events occur more frequently during the morning (MCEJ) than in the afternoon (ACEJ). The highest MCEJ and ACEJ occurrence rates were observed for the Brazilian sector. Distinct seasonal dependence was found for each longitudinal sector under investigation. The lunar phase dependence was determined for the first time for the Philippine sector (longitude 125 degrees E), and it was shown to be less pronounced than in the Peruvian, Brazilian and African sectors. We demonstrate that differences in CEJ rates derived from ground-based and satellite data can arise from the longitudinal separation between low-latitude and equatorial stations that are used to determine the signal and its consequent time delay in their sunrise/sunset times at ionospheric heights.
SiO(2 )is the main component of silicate melts and thus controls their network structure and physical properties. The compressibility and viscosities of melts at depth are governed by their short range atomic and electronic structure. We measured the O K-edge and the Si L-2,L-3-edge in silica up to 110 GPa using X-ray Raman scattering spectroscopy, and found a striking match to calculated spectra based on structures from molecular dynamic simulations. Between 20 and 27 GPa, Si-[4] species are converted into a mixture of Si-[5] and Si-[6] species and between 60 and 70 GPa, Si-[6] becomes dominant at the expense of Si-[5] with no further increase up to at least 110 GPa. Coordination higher than 6 is only reached beyond 140 GPa, corroborating results from Brillouin scattering. Network modifying elements in silicate melts may shift this change in coordination to lower pressures and thus magmas could be denser than residual solids at the depth of the core-mantle boundary.
This study presents the first suite of apatite fission-track (AFT) ages from the SE part of the Western Sudetes. AFT cooling ages from the Orlica-snie(z) over dotnik Dome and the Upper Nysa Klodzka Graben range from Late Cretaceous (84 Ma) to Early Palaeocene-Middle Eocene (64-45 Ma). The first stage of basin evolution (similar to 100-90 Ma) was marked by the formation of a local extensional depocentre and disruption of the Mesozoic planation surface. Subsequent far-field convergence of European microplates resulted in Coniacian-Santonian (similar to 89-83 Ma) thrust faulting. AFT data from both metamorphic basement and Mesozoic sedimentary cover indicate homogenous Late Cretaceous burial of the entire Western Sudetes. Thermal history modeling suggests that the onset of cooling could be constrained between 89 and 63 Ma with a climax during the Palaeocene-Middle Eocene basin inversion phase.
Metamorphic geology
(2019)
Secondary mica minerals collected from the Santa Helena (W- (Cu) mineralization) and Venise (W-Mo mineralization) endogenic breccia structures were Ar-40/Ar-39 dated. The muscovite Ar-40/Ar-39 data yielded 286.8 +/- 1.2 (+/- 1 sigma) Ma (samples 6Ha and 11Ha) which reflect the age of secondary muscovite formation probably from magmatic biotite or feldspar alteration. Sericite Ar-40/Ar-39 data yielded 280.9 +/- 1.2 (+/- 1 sigma) Ma to 279.0 +/- 1.1 (+/- 1 sigma) Ma (samples 6Hb and 11Hb) reflecting the age of greisen alteration (T similar to 300 degrees C) where the W- disseminated mineralization occurs. The muscovite 40Ar/39Ar data of 277.3 +/- 1.3 (+/- 1 sigma) Ma and 281.3 +/- 1.2 (+/- 1 sigma) Ma (samples 5 and 6) also reflect the age of muscovite (selvage) crystallized adjacent to molybdenite veins within the Venise breccia. Geochronological data obtained confirmed that the W mineralization at Santa Helena breccia is older than Mo-mineralization at Venise breccia. Also, the timing of hydrothermal circulation and the cooling history for the W-stage deposition was no longer than 7 Ma and 4 Ma for Mo-deposition.
Microbial community composition and abundance after millennia of submarine permafrost warming
(2019)
Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 degrees C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore-offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 degrees C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes delta O-18 and delta D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 degrees C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.
We present new U-Pb LA-ICP-MS data from the Central Andean foreland basins combined with new and published stratigraphic information in order to reconstruct the Miocene fragmentation of the Andean foreland between 26 and 28 degrees S. The disruption of this foreland basin and the subsequent development of elevated intermountain basins have been the focus of several studies. However, the absence of temporal constraints in the Miocene to Pliocene sedimentary record of the low elevation Choromoro and Tucuman foreland basins has presented an obstacle for precise paleogeographic reconstructions. We describe 11 discontinuous stratigraphic sections and use the U-Pb LA-ICP-MS method to date 10 pyroclastic-bearing sediments in order to reconstruct the stratigraphic evolution of the Choromoro and Tucuman basins. We combine our results with published strati graphic and thermochronologic data from adjacent basins to present a refined Miocene paleogeographic model. In a first stage, a continuous Early Miocene foreland lacustrine basin developed, filling up the preexisting Paleogene topography. The second stage is characterized by basin unroofing around similar to 12 Ma; the easily eroded sedimentary cover was removed, leading to the uplift of the underlying basement rocks and the segmentation of the lacustrine system. In the third stage, relief increase took place after similar to 6 Ma due to the low erodibility of the basement blocks; as a result, stable fluvial systems developed. Progressive relief development caused pronounced unconformities in the basins and the development of proximal fluvial-gravitational depositional systems after 3 Ma. This model emphasizes on the relations between tectonics, climate, and erodibility, and their control on the evolution of the depositional systems and relief.
The Andean Plateau of NW Argentina is a prominent example of a high-elevation orogenic plateau characterized by internal drainage, arid to hyper-arid climatic conditions and a compressional basin-and-range morphology comprising thick sedimentary basins. However, the development of the plateau as a geomorphic entity is not well understood. Enhanced orographic rainout along the eastern, windward plateau flank causes reduced fluvial run-off and thus subdued surface-process rates in the arid hinterland. Despite this, many Puna basins document a complex history of fluvial processes that have transformed the landscape from aggrading basins with coalescing alluvial fans to the formation of multiple fluvial terraces that are now abandoned. Here, we present data from the San Antonio de los Cobres (SAC) area, a sub-catchment of the Salinas Grandes Basin located on the eastern Puna Plateau bordering the externally drained Eastern Cordillera. Our data include: (a) new radiometric U-Pb zircon data from intercalated volcanic ash layers and detrital zircons from sedimentary key horizons; (b) sedimentary and geochemical provenance indicators; (c) river profile analysis; and (d) palaeo-landscape reconstruction to assess aggradation, incision and basin connectivity. Our results suggest that the eastern Puna margin evolved from a structurally controlled intermontane basin during the Middle Miocene, similar to intermontane basins in the Mio-Pliocene Eastern Cordillera and the broken Andean foreland. Our refined basin stratigraphy implies that sedimentation continued during the Late Mio-Pliocene and the Quaternary, after which the SAC area was subjected to basin incision and excavation of the sedimentary fill. Because this incision is unrelated to baselevel changes and tectonic processes, and is similar in timing to the onset of basin fill and excavation cycles of intermontane basins in the adjacent Eastern Cordillera, we suspect a regional climatic driver, triggered by the Mid-Pleistocene Climate Transition, caused the present-day morphology. Our observations suggest that lateral orogenic growth, aridification of orogenic interiors, and protracted plateau sedimentation are all part of a complex process chain necessary to establish and maintain geomorphic characteristics of orogenic plateaus in tectonically active mountain belts.
Environmental factors shape the spatial distribution and dynamics of populations. Understanding how these factors interact with movement behavior is critical for efficient conservation, in particular for migratory species. Adult female green sea turtles, Chelonia mydas, migrate between foraging and nesting sites that are generally separated by thousands of kilometers. As an emblematic endangered species, green turtles have been intensively studied, with a focus on nesting, migration, and foraging. Nevertheless, few attempts integrated these behaviors and their trade‐offs by considering the spatial configurations of foraging and nesting grounds as well as environmental heterogeneity like oceanic currents and food distribution. We developed an individual‐based model to investigate the impact of local environmental conditions on emerging migratory corridors and reproductive output and to thereby identify conservation priority sites. The model integrates movement, nesting, and foraging behavior. Despite being largely conceptual, the model captured realistic movement patterns which confirm field studies. The spatial distribution of migratory corridors and foraging hot spots was mostly constrained by features of the regional landscape, such as nesting site locations, distribution of feeding patches, and oceanic currents. These constraints also explained the mixing patterns in regional forager communities. By implementing alternative decision strategies of the turtles, we found that foraging site fidelity and nesting investment, two characteristics of green turtles' biology, are favorable strategies under unpredictable environmental conditions affecting their habitats. Based on our results, we propose specific guidelines for the regional conservation of green turtles as well as future research suggestions advancing spatial ecology of sea turtles. Being implemented in an easy to learn open‐source software, our model can coevolve with the collection and analysis of new data on energy budget and movement into a generic tool for sea turtle research and conservation. Our modeling approach could also be useful for supporting the conservation of other migratory marine animals.
Germany and the United Kingdom have domestic shale gas reserves which they may exploit in the future to complement their national energy strategies. However gas production releases volatile organic compounds (VOC) and nitrogen oxides (NOx), which through photochemical reaction form ground-level ozone, an air pollutant that can trigger adverse health effects e.g. on the respiratory system. This study explores the range of impacts of a potential shale gas industry in these two countries on local and regional ambient ozone. To this end, comprehensive emission scenarios are used as the basis for input to an online-coupled regional chemistry transport model (WRF-Chem). Here we simulate shale gas scenarios over summer (June, July, August) 2011, exploring the effects of varying VOC emissions, gas speciation, and concentration of NOx emissions over space and time, on ozone formation. An evaluation of the model setup is performed, which exhibited the model’s ability to predict surface meteorological and chemical variables well compared with observations, and consistent with other studies. When different shale gas scenarios were employed, the results show a peak increase in maximum daily 8-hour average ozone from 3.7 to 28.3 μg m–3. In addition, we find that shale gas emissions can force ozone exceedances at a considerable percentage of regulatory measurement stations locally (up to 21% in Germany and 35% in the United Kingdom) and in distant countries through long-range transport, and increase the cumulative health-related metric SOMO35 (maximum percent increase of ~28%) throughout the region. Findings indicate that VOC emissions are important for ozone enhancement, and to a lesser extent NOx, meaning that VOC regulation for a future European shale gas industry will be of especial importance to mitigate unfavorable health outcomes. Overall our findings demonstrate that shale gas production in Europe can worsen ozone air quality on both the local and regional scales.
The ‘bomb-pulse’ method is a chronological approach to further constrain the age of speleothems that grew between 1950 CE – present. Establishing dependable chronological constraints is crucial for modern calibration studies of speleothems to instrumental climate records, which provides the basis for paleoclimate interpretations. However, a large unknown is how 14C is transferred from the atmosphere to any individual speleothem owing to the site-specific residence times of organic matter above cave systems. Here, we employ the bomb-pulse method to build chronologies from 14C measurements in combination with a new unsaturated zone C model which considers C decomposition as a continuum, to better understand unsaturated zone 14C dynamics. The bomb-pulse curves of eight speleothems from southern Australia in three contrasting climatic regions; the semi-arid Wellington Caves site, the mediterranean Golgotha Cave site and the montane Yarrangobilly Caves site, are investigated. Overall, the modelled 14C bomb-pulse curves produce excellent fits with measured 14C speleothem data (r2 = 0.82–0.99). The C modelling reveals that unsaturated zone C is predominately young at the semi-arid site, with a weighted-mean residence time of 32 years and that tree root respiration is likely an important source of vadose CO2. At the montane site, ∼39% of C is young (<1 years), but the weighted-mean C ages are older (145–220 years). The mediterranean site has very little contribution from young C (<12%: 0–1 years), with weighted-mean ages between 157 and 245 years, likely due to greater adsorption of organic matter in the upper vadose zone during matrix flow, and remobilisation of C from young syngenetic karst. New end members for low speleothem Dead Carbon Proportion (DCP) are identified (2.19% and 1.65%, respectively) for Australian montane and semi-arid zone speleothems, where oversupply of modern CO2 in the vadose zone leads to lower DCP. It is also demonstrated that DCP can be quite variable over small time scales, that processes may be difficult to untangle and a constant DCP assumption is likely invalid. DCP variability over time is mainly controlled by the changes vadose zone CO2, where vegetation regeneration, wild-fires and karst hydrology play an important role.
In this paper we report a rare and fortunate event of fast magnetosonic (MS, also called equatorial noise) waves modulated by compressional ultralow frequency (ULF) waves measured by Van Allen Probes. The characteristics of MS waves, ULF waves, proton distribution, and their potential correlations are analyzed. The results show that ULF waves can modulate the energetic ring proton distribution and in turn modulate the MS generation. Furthermore, the variation of MS intensities is attributed to not only ULF wave activities but also the variation of background parameters, for example, number density. The results confirm the opinion that MS waves are generated by proton ring distribution and propose a new modulation phenomenon.
Due to increasing demands for irrigation using groundwater as a source there is an urgent need for efficient methods that shed light on the resulting anthropogenic impacts on the connected aquifers. Thus an innovative approach is introduced, that aims to identify predominant geochemical changes in the groundwater system. The approach involves a principal component analysis as a promising tool to disentangle the effects of different impacts and even to give a quantitative assessment of the respective effect strength at each site. The study was applied in an irrigation region of the Nuthe River Basin, State Brandenburg, Northeast Germany. The results identify the negative impacts on the groundwater quality in the aquifer used for irrigation. A decrease of shallow groundwater quality under irrigation due to contamination with fertilizers (NO3, Cl, K, Na) and a slight shift in the redox system is indicated. Beside this direct impact on the shallow groundwater a long-term impact on a deeper groundwater resource could be identified. There is clear evidence, that the contamination is not restricted to the shallow groundwater but that extraction from deeper wells increasingly includes deeper, uncontaminated groundwater resources into the local irrigation cycle. The approach can be used as a basic tool for the adaptation of sustainable agricultural irrigation management strategies.
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
Adjustment of median ground motion prediction equations (GMPEs) from one region to another region is one of the major challenges within the current practice of seismic hazard analysis. In our approach of generating response spectra, we derive two separate empirical models for a) Fourier amplitude spectrum (FAS) and b) duration of ground motion. To calculate response spectra, the two models are combined within the random vibration theory (RVT) framework. The models are calibrated on recordings obtained from shallow crustal earthquakes in active tectonic regions. We use a subset of NGA-West2 database with M3.2-7.9 earthquakes at distances 0-300 km. The NGA-West2 database expanded over a wide magnitude and distance range facilitates a better constraint over derived models. A frequency-dependent duration model is derived to obtain adjustable response spectral ordinates. Excellent comparison of our approach with other NGA-West2 models implies that it can also be used as a stand-alone model.
Ongoing and past biome transitions are generally assigned to climate and atmospheric changes (e.g. temperature, precipitation, CO2), but the major regional factors or factor combinations that drive vegetation change often remain unknown. Modelling studies applying ensemble runs can help to partition the effects of the different drivers. Such studies require careful validation with observational data. In this study, fossil pollen records from 741 sites in Europe, 728 sites in North America, and 418 sites in Asia (extracted from terrestrial archives including lake sediments) are used to reconstruct biomes at selected time slices between 40 cal ka BP (calibrated thousand years before present) and today. These results are used to validate Northern Hemisphere biome distributions (>30 degrees N) simulated by the biome model BIOME4 that has been forced with climate data simulated by a General Circulation model. Quantitative comparisons between pollen- and model-based results show a generally good fit at a broad spatial scale. Mismatches occur in central-arid Asia with a broader extent of grassland throughout the last 40 ka (likely due to the over-representation of Artemisia and Chenopodiaceae pollen) and in Europe with over-estimation of tundra at 0 cal ka BP (likely due to human impacts to some extent). Sensitivity analysis reveals that broad-scale biome changes follow the global signal of major postglacial temperature change, although the climatic variables vary in their regional and temporal importance. Temperature is the dominant variable in Europe and other rather maritime areas for biome changes between 21 and 14 ka, while precipitation is highly important in the arid inland regions of Asia and North America. The ecophysiological effect of changes in the atmospheric CO2-concentration has the highest impact during this transition than in other intervals. With respect to modern vegetation in the course of global warming, our findings imply that vegetation change in the Northern Hemisphere may be strongly limited by effective moisture changes, i.e. the combined effect of temperature and precipitation, particularly in inland areas. (C) 2019 Elsevier Ltd. All rights reserved.
The evolution of the radiation belts in L-shell (L), energy (E), and equatorial pitch angle (alpha(0)) is analyzed during the calm 11-day interval (4-15 March) following the 1 March 2013 storm. Magnetic Electron and Ion Spectrometer (MagEIS) observations from Van Allen Probes are interpreted alongside 1D and 3D Fokker-Planck simulations combined with consistent event-driven scattering modeling from whistler mode hiss waves. Three (L, E, alpha(0)) regions persist through 11 days of hiss wave scattering; the pitch angle-dependent inner belt core (L similar to <2.2 and E < 700 keV), pitch angle homogeneous outer belt low-energy core (L > similar to 5 and E similar to < 100 keV), and a distinct pocket of electrons (L similar to [4.5, 5.5] and E similar to [0.7, 2] MeV). The pitch angle homogeneous outer belt is explained by the diffusion coefficients that are roughly constant for alpha(0) similar to <60 degrees, E > 100 keV, 3.5 < L < L-pp similar to 6. Thus, observed unidirectional flux decays can be used to estimate local pitch angle diffusion rates in that region. Top-hat distributions are computed and observed at L similar to 3-3.5 and E = 100-300 keV.
In this study, we detect high percentile rainfall events in the eastern central Andes, based on Tropical Rainfall Measuring Mission (TRMM) with a spatial resolution of 0.25 × 0.25°, a temporal resolution of 3 h, and for the duration from 2001 to 2018. We identify three areas with high mean accumulated rainfall and analyze their atmospheric behaviour and rainfall characteristics with specific focus on extreme events. Extreme events are defined by events above the 95th percentile of their daily mean accumulated rainfall. Austral summer (DJF) is the period of the year presenting the most frequent extreme events over these three regions. Daily statistics show that the spatial maxima, as well as their associated extreme events, are produced during the night. For the considered period, ERA-Interim reanalysis data, provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) with 0.75° x0.75° spatial and 6-hourly temporal resolutions, were used for the analysis of the meso- and synoptic-scale atmospheric patterns. Night- and day-time differences indicate a nocturnal overload of northerly and northeasterly low-level humidity flows arriving from tropical South America. Under these conditions, cooling descending air from the mountains may find unstable air at the surface, giving place to the development of strong local convection. Another possible mechanism is presented here: a forced ascent of the low-level flow due to the mountains, disrupting the atmospheric stratification and generating vertical displacement of air trajectories. A Principal Component Analysis (PCA) in T-mode is applied to day- and night-time data during the maximum and extreme events. The results show strong correlation areas over each subregion under study during night-time, whereas during day-time no defined patterns are found. This confirms the observed nocturnal behavior of rainfall within these three hotspots.
OpenForecast
(2019)
The development and deployment of new operational runoff forecasting systems are a strong focus of the scientific community due to the crucial importance of reliable and timely runoff predictions for early warnings of floods and flashfloods for local businesses and communities. OpenForecast, the first operational runoff forecasting system in Russia, open for public use, is presented in this study. We developed OpenForecast based only on open-source software and data-GR4J hydrological model, ERA-Interim meteorological reanalysis, and ICON deterministic short-range meteorological forecasts. Daily forecasts were generated for two basins in the European part of Russia. Simulation results showed a limited efficiency in reproducing the spring flood of 2019. Although the simulations managed to capture the timing of flood peaks, they failed in estimating flood volume. However, further implementation of the parsimonious data assimilation technique significantly alleviates simulation errors. The revealed limitations of the proposed operational runoff forecasting system provided a foundation to outline its further development and improvement.
Thermokarst lake landscapes are permafrost regions, which are prone to rapid (on seasonal to decadal time scales) changes, affecting carbon and nitrogen cycles. However, there is a high degree of uncertainty related to the balance between carbon and nitrogen cycling and storage. We collected 12 permafrost soil cores from six drained thermokarst lake basins (DTLBs) along a chronosequence north of Teshekpuk Lake in northern Alaska and analyzed them for carbon and nitrogen contents. For comparison, we included three lacustrine cores from an adjacent thermokarst lake and one soil core from a non thermokarst affected remnant upland. This allowed to calculate the carbon and nitrogen stocks of the three primary landscape units (DTLB, lake, and upland), to reconstruct the landscape history, and to analyze the effect of thermokarst lake formation and drainage on carbon and nitrogen stocks. We show that carbon and nitrogen contents and the carbon-nitrogen ratio are considerably lower in sediments of extant lakes than in the DTLB or upland cores indicating degradation of carbon during thermokarst lake formation. However, we found similar amounts of total carbon and nitrogen stocks due to the higher density of lacustrine sediments caused by the lack of ground ice compared to DTLB sediments. In addition, the radiocarbon-based landscape chronology for the past 7,000years reveals five successive lake stages of partially, spatially overlapping DTLBs in the study region, reflecting the dynamic nature of ice-rich permafrost deposits. With this study, we highlight the importance to include these dynamic landscapes in future permafrost carbon feedback models. Plain Language Summary When permanently frozen soils (permafrost) contain ice-rich sediments, the thawing of this permafrost causes the surface to sink, which may result in lake formation. This process, the thaw of ice-rich permafrost and melting of ground ice leads to characteristic landforms-known as thermokarst. Once such a thaw process is initiated in ice-rich sediments, a thaw lake forms and grows by shoreline erosion, eventually expanding until a drainage pathway is encountered and the lake eventually drains, resulting in a drained thermokarst lake basin. In our study, we show that such a thermokarst-affected landscape north of Teshekpuk Lake in northern Alaska is shaped by repeated thaw lake formation and lake drainage events during the past 7,000years, highlighting the dynamic nature of these landscapes. These landscape-scale processes have a big effect on the carbon and nitrogen stored in permafrost soils. We show that large amounts of carbon (>45kg C/m(2)) and nitrogen (>2.6kg N/m(2)) are stored in unfrozen lake sediments and in frozen soil sediments. The findings are important when considering the potential effect that permafrost thaw has for the global climate through releasing carbon and nitrogen, which was frozen and therefore locked away for millennia, from the active carbon cycle.
With the advent of the two Sentinel-1 (S1) satellites, Synthetic Aperture Radar (SAR) data with high temporal and spatial resolution are freely available. This provides a promising framework to facilitate detailed investigations of surface instabilities and movements on large scales with high temporal resolution, but also poses substantial processing challenges because of storage and computation requirements. Methods are needed to efficiently detect short term changes in dynamic environments. Approaches considering pair-wise processing of a series of consecutive scenes to retain maximum temporal resolution in conjunction with time series analyses are required. Here we present OSARIS, the “Open Source SAR Investigation System,” as a framework to process large stacks of S1 data on high-performance computing clusters. Based on Generic Mapping Tools SAR, shell scripts, and the workload manager Slurm, OSARIS provides an open and modular framework combining parallelization of high-performance C programs, flexible processing schemes, convenient configuration, and generation of geocoded stacks of analysis-ready base data, including amplitude, phase, coherence, and unwrapped interferograms. Time series analyses can be conducted by applying automated modules to the data stacks. The capabilities of OSARIS are demonstrated in a case study from the northwestern Tien Shan, Central Asia. After merging of slices, a total of 80 scene pairs were processed from 174 total input scenes. The coherence time series exhibits pronounced seasonal variability, with relatively high coherence values prevailing during the summer months in the nival zone. As an example of a time series analysis module, we present OSARIS' “Unstable Coherence Metric” which identifies pixels affected by significant drops from high to low coherence values. Measurements of motion provided by LOSD measurements require careful evaluation because interferometric phase unwrapping is prone to errors. Here, OSARIS provides a series of modules to detect and mask unwrapping errors, correct for atmospheric disturbances, and remove large-scale trends. Wall clock processing time for the case study (area ~9,000 km2) was ~12 h 4 min on a machine with 400 cores and 2 TB RAM. In total, ~12 d 10 h 44 min (~96%) were saved through parallelization. A comparison of selected OSARIS datasets to results from two state-of-the-art SAR processing suites, ISCE and SNAP, shows that OSARIS provides products of competitive quality despite its high level of automatization. OSARIS thus facilitates efficient S1-based region-wide investigations of surface movement events over multiple years.
The progress of science is tied to the standardization of measurements, instruments, and data. This is especially true in the Big Data age, where analyzing large data volumes critically hinges on the data being standardized. Accordingly, the lack of community-sanctioned data standards in paleoclimatology has largely precluded the benefits of Big Data advances in the field. Building upon recent efforts to standardize the format and terminology of paleoclimate data, this article describes the Paleoclimate Community reporTing Standard (PaCTS), a crowdsourced reporting standard for such data. PaCTS captures which information should be included when reporting paleoclimate data, with the goal of maximizing the reuse value of paleoclimate data sets, particularly for synthesis work and comparison to climate model simulations. Initiated by the LinkedEarth project, the process to elicit a reporting standard involved an international workshop in 2016, various forms of digital community engagement over the next few years, and grassroots working groups. Participants in this process identified important properties across paleoclimate archives, in addition to the reporting of uncertainties and chronologies; they also identified archive-specific properties and distinguished reporting standards for new versus legacy data sets. This work shows that at least 135 respondents overwhelmingly support a drastic increase in the amount of metadata accompanying paleoclimate data sets. Since such goals are at odds with present practices, we discuss a transparent path toward implementing or revising these recommendations in the near future, using both bottom-up and top-down approaches.
A number of sedimentary provenance studies have been undertaken in order to determine whether the palaeo-Red River was once a river of continental proportions into which the upper reaches of the Yangtze, Salween, Mekong, Irrawaddy, and Yarlung drained. We have assessed the evidence that the Yarlung originally flowed into the palaeo-Red river, and then sequentially into the Irrawaddy and Brahmaputra, connecting to the latter first via the Lohit and then the Siang. For this river system, we have integrated our new data from the Paleogene-Recent Irrawaddy drainage basin (detrital zircon U-Pb with Hf and fission track, rutile U-Pb, mica Ar-Ar, bulk rock Sr-Nd, and petrography) with previously published data, to produce a palaeodrainage model that is consistent with all datasets. In our model, the Yarlung never flowed into the Irrawaddy drainage: during the Paleogene, the Yarlung suture zone was an internally drained basin, and from Neogene times onwards the Yarlung drained into the Brahmaputra in the Bengal Basin. The Central Myanmar Basin, through which the Irrawaddy River flows today, received predominantly locally-derived detritus until the Middle Eocene, the Irrawaddy initiated as a through-going river draining the Mogok Metamorphic Belt and Bomi-Chayu granites to the north sometime in the Late Eocene to Early Oligocene, and the river was dominated by a stable MMB-dominated drainage throughout the Neogene to present day. Existing evidence does not support any connection between the Yarlung and the Red River in the past, but there is a paucity of suitable palaeo-Red River deposits with which to make a robust comparison. We argue that this limitation also precludes a robust assessment of a palaeo-connection between the Yangtze/ Salween/Mekong and the Red River; it is difficult to unequivocally interpret the recorded provenance changes as the result of specific drainage reorganisations. We highlight the palaeo-Red River deposits of the Hanoi Basin as a potential location for future research focus in view of the near-complete Cenozoic record of palaeo-Red River deposits at this location. A majority of previous studies consider that if a major continental-scale drainage ever existed at all, it fragmented early in the Cenozoic. Such a viewpoint would agree with the growing body of evidence from palaeoaltitude studies that large parts of SE Tibet were uplifted by this period. This then leads towards the intriguing question as to the mechanisms which caused the major period of river incision in the Miocene in this region.
Paleogene evolution and demise of the proto-Paratethys Sea in Central Asia (Tarim and Tajik basins)
(2019)
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.
Rapid local adaptation frequently occurs during the spread of invading species. It remains unclear, however, how consistent, and therefore potentially predictable, such patterns of local adaptation are. One approach to this question is to measure patterns of local differentiation in functional traits and plasticity levels in invasive species in multiple regions. Finding consistent patterns of local differentiation in replicate regions suggests that these patterns are adaptive. Further, this outcome indicates that the invading species likely responds predictably to selection along environmental gradients, even though standing genetic variation is likely to have been reduced during introduction. We studied local differentiation in the invasive annual plant Erodium cicutarium in two invaded regions, California and Chile. We collected seeds from across strong gradients in precipitation and temperature in Mediterranean-climate parts of the two regions (10 populations per region). We grew seeds from maternal families from these populations through two generations and exposed the second generation to contrasting levels of water and nutrient availability. We measured growth, flowering time and leaf functional traits across these treatments to obtain trait means and plasticity measures. We found strong differentiation among populations in all traits. Plants from drier environments flowered earlier, were less plastic in flowering time and reached greater size in all treatments. Correlations among traits within regions suggested a coordinated evolutionary response along environmental gradients associated with growing season length. There was little divergence in traits and trait intercorrelations between regions, but strongly parallel divergence in traits within regions. Similar, statistically consistent patterns of local trait differentiation across two regions suggest that local adaptation to environmental gradients has aided the spread of this invasive species, and that the formation of ecotypes in newly invaded environments has been relatively consistent and predictable.