550 Geowissenschaften
Refine
Year of publication
- 2019 (344) (remove)
Document Type
- Article (269)
- Doctoral Thesis (29)
- Postprint (21)
- Other (12)
- Review (7)
- Habilitation Thesis (4)
- Conference Proceeding (1)
- Master's Thesis (1)
Keywords
- climate (5)
- remote sensing (5)
- bryophytes (4)
- lichens (4)
- precipitation (4)
- Andes (3)
- Climate (3)
- Ecosystem services (3)
- Holocene (3)
- Permian (3)
Institute
- Institut für Geowissenschaften (253)
- Institut für Umweltwissenschaften und Geographie (66)
- Institut für Biochemie und Biologie (9)
- Mathematisch-Naturwissenschaftliche Fakultät (7)
- Institut für Physik und Astronomie (5)
- Extern (3)
- Institut für Mathematik (2)
- Zentrum für Umweltwissenschaften (2)
- Institut für Künste und Medien (1)
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
Geysers are hot springs whose frequency of water eruptions remain poorly understood. We set up a local broadband seismic network for 1 year at Strokkur geyser, Iceland, and developed an unprecedented catalog of 73,466 eruptions. We detected 50,135 single eruptions but find that the geyser is also characterized by sets of up to six eruptions in quick succession. The number of single to sextuple eruptions exponentially decreased, while the mean waiting time after an eruption linearly increased (3.7 to 16.4 min). While secondary eruptions within double to sextuple eruptions have a smaller mean seismic amplitude, the amplitude of the first eruption is comparable for all eruption types. We statistically model the eruption frequency assuming discharges proportional to the eruption multiplicity and a constant probability for subsequent events within a multituple eruption. The waiting time after an eruption is predictable but not the type or amplitude of the next one. <br /> Plain Language Summary Geysers are springs that often erupt in hot water fountains. They erupt more often than volcanoes but are quite similar. Nevertheless, it is poorly understood how often volcanoes and also geysers erupt. We created a list of 73,466 eruption times of Strokkur geyser, Iceland, from 1 year of seismic data. The geyser erupted one to six times in quick succession. We found 50,135 single eruptions but only 1 sextuple eruption, while the mean waiting time increased from 3.7 min after single eruptions to 16.4 min after sextuple eruptions. Mean amplitudes of each eruption type were higher for single eruptions, but all first eruptions in a succession were similar in height. Assuming a constant heat inflow at depth, we can predict the waiting time after an eruption but not the type or amplitude of the next one.
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.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
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.
Bank filtration (BF) is an established indirect water-treatment technology. The quality of water gained via BF depends on the subsurface capture zone, the mixing ratio (river water versus ambient groundwater), spatial and temporal distribution of subsurface travel times, and subsurface temperature patterns. Surface-water infiltration into the adjacent aquifer is determined by the local hydraulic gradient and riverbed permeability, which could be altered by natural clogging, scouring and artificial decolmation processes. The seasonal behaviour of a BF system in Germany, and its development during and about 6 months after decolmation (canal reconstruction), was observed with a long-term monitoring programme. To quantify the spatial and temporal variation in the BF system, a transient flow and heat transport model was implemented and two model scenarios, 'with' and 'without' canal reconstruction, were generated. Overall, the simulated water heads and temperatures matched those observed. Increased hydraulic connection between the canal and aquifer caused by the canal reconstruction led to an increase of similar to 23% in the already high share of BF water abstracted by the nearby waterworks. Subsurface travel-time distribution substantially shifted towards shorter travel times. Flow paths with travel times <200 days increased by similar to 10% and those with <300 days by 15%. Generally, the periodic temperature signal, and the summer and winter temperature extrema, increased and penetrated deeper into the aquifer. The joint hydrological and thermal effects caused by the canal reconstruction might increase the potential of biodegradable compounds to further penetrate into the aquifer, also by potentially affecting the redox zonation in the aquifer.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
The geochemical composition of oceanic basalts provides us with a window into the distribution of geochemical elements within the Earth’s mantle in space and time. In conjunction with a throughout knowledge on how the different elements behave e.g. during melt formation and evolution or on their partition behaviour between e.g. minerals and melts this information has been transformed into various models on how oceanic crust is formed along plume influenced or normal mid-ocean ridge segments, how oceanic crust evolves in response to seawater, on subduction recycling of oceanic crust and so forth. The work presented in this habilitation was aimed at refining existing models, putting further constraints on some of the major open questions in this field of research while at the same time trying to increase our knowledge on the behaviour of noble gases as a tracer for melt formation and evolution processes. In the line of this work the author and her co-workers were able to answer one of the major questions concerning the formation of oceanic crust along plume-influenced ridges – in which physical state does the plume material enter the ridge? Based on submarine volcanic glass He, Ne and Ar data, the author and her co-workers have shown that the interaction of mantle plumes with mid-ocean ridges occurs in the physical form of melts. In addition, the author and her co-workers have also put further constraints on one of the major questions concerning the formation of oceanic crust along normal mid-ocean ridges – namely how is the mid-ocean ridge system effectively cooled to form the lower oceanic crust? Based on Ne and Ar data in combination with Cl/K ratios of basaltic glass from the Mid-Atlantic ridge and estimates of crystallisation pressures they have shown, that seawater penetration reaches lower crustal levels close to the Moho, indicating that hydrothermal circulation might be an effective cooling mechanism even for the deep parts of the oceanic crust. Considering subduction recycling, the heterogeneity of the Earth’s mantle and mantle dynamic processes the key question is on which temporal and spatial scales is the Earth’s mantle geochemically heterogeneous? In the line of this work the author along with her co-workers have shown based on Cl/K ratios in conjunction with the Sr, Nd, and Pb isotopes of the OIBs representing the type localities for the different mantle endmembers that the quantity of Cl recycled into the mantle via subduction is not uniform and that neither the HIMU nor the EM1 and EM2 mantle components can be considered as distinct mantle endmembers. In addition, we have shown, based on He, Ne and Ar isotope and trace-element data from the Foundation hotspot that the near ridge seamounts of the Foundation seamount chain formed by the Foundation hotspot erupt lavas with a trace-element signature clearly characteristic of oceanic gabbro which indicates the existence of recycled, virtually unchanged lower oceanic crust in the plume source. This is a clear sign of the inefficiency of the stirring mechanism existing at mantle depth. Similar features are seen in other near-axis hotspot magmas around the world. Based on He, Sr, Nd, Pb and O isotopes and trace elements in primitive mafic dykes from the Etendeka flood basalts, NW Namibia the author along with her co-workers have shown that deep, less degassed mantle material carried up by a mantle plume contributed significantly to the flood basalt magmatism. The Etendeka flood basalts are part of the South Atlantic LIP, which is associated with the breakup of Gondwana, the formation of the Paraná-Etendeka flood basalts and the Walvis Ridge - Tristan da Cunha hotspot track. Thus reinforcing the lately often-challenged concept of mantle plumes and the role of mantle plumes in the formation of large igneous provinces. Studying the behaviour of noble gases during melt formation and evolution the author along with her co-workers has shown that He can be considerable more susceptible to changes during melt formation and evolution resulting not only in a complete decoupling of He isotopes from e.g. Ne or Pb isotopes but also in a complete loss of the primary mantle isotope signal. They have also shown that this decoupling occurs mainly during the melt formation processes requiring He to be more compatible during mantle melting than Ne. In addition, the author along with her co workers were able to show that incorporation of atmospheric noble gases into igneous rocks is in general a two-step process: (1) magma contamination by assimilation of altered oceanic crust results in the entrainment of air-equilibrated seawater noble gases; (2) atmospheric noble gases are adsorbed onto grain surfaces during sample preparation. This implies, considering the ubiquitous presence of the contamination signal, that magma contamination by assimilation of a seawater-sourced component is an integral part of mid-ocean ridge basalt evolution.
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.
A well-preserved sequence, by Archean standards, of mantle-derived metabasalts and metakomatiites forms large parts of the lower Onverwacht Group of the Barberton Greenstone Belt (South Africa). To elucidate the origin of mafic and ultramafic rocks from this 3.55 to 3.45 Ga sequence, we present a comprehensive geochemical dataset including major and trace elements as well as Lu-Hf and Sm-Nd isotope compositions for a variety of metavolcanic rocks. These include metabasalts of the amphibolite-facies Sandspruit and Theespruit Formations as well as metabasalts and metakomatiites of the lower greenschist-facies Komati Formation. Based on their incompatible trace element patterns, the basalts of the Sandspruit and Theespruit Formations can be subdivided into a light rare earth element (LREE) depleted group, a LREE-undepleted group, and a LREE-enriched group. Positive epsilon Hf-(t) and epsilon Nd-(t) values of ca. +3 to +4 and 0 to +2, respectively, together with depletions in Th and La-CN/Yb-CN indicate derivation of the LREE-depleted basalts from a depleted mantle source. However, chondritic epsilon Hf-(t) and epsilon Nd-(t) values combined with positive Th and La-CN/Yb-CN of the LREE-enriched samples indicate a contribution from older granitoid crust in the petrogenesis of these samples. Trace element patterns of komatiites and basalts of the Komati Formation are generally flat relative to primitive mantle with slight depletions in heavy rare earth elements and Th and overall positive epsilon Hf-(t) of + 2.5 +/- 3.5 (2 s.d.) and epsilon Nd-(t) of + 0.5 +/- 2.2 (2 s. d.). The coherence in trace element characteristics suggests a common magmatic origin for basalts and komatiites. This study reveals that the two lavas were derived from the same mantle plume, i. e. komatiites were formed by high degrees of melting of a depleted mantle source containing residual garnet and the basalts were formed by moderate degrees of partial melting in shallower regions of the mantle. Based on the current dataset, combined with published data, we propose a geodynamic model for the oldest units of the Barberton Greenstone Belt that describes the development from a submerged continental setting (for the Sandspruit and Theespruit Formations) to a submarine plateau setting (for the Komati Formation) as a consequence of continental rifting.
Climate science is highly interdisciplinary by nature, so understanding interactions between Earth processes inherently warrants the use of analytical software that can operate across the disciplines of Earth science. Toward this end, we present the Climate Data Toolbox for MATLAB, which contains more than 100 functions that span the major climate-related disciplines of Earth science. The toolbox enables streamlined, entirely scriptable workflows that are intuitive to write and easy to share. Included are functions to evaluate uncertainty, perform matrix operations, calculate climate indices, and generate common data displays. Documentation is presented pedagogically, with thorough explanations of how each function works and tutorials showing how the toolbox can be used to replicate results of published studies. As a well-tested, well-documented platform for interdisciplinary collaborations, the Climate Data Toolbox for MATLAB aims to reduce time spent writing low-level code, let researchers focus on physics rather than coding and encourage more efficacious code sharing. Plain Language Summary This article describes a collection of computer code that has recently been released to help scientists analyze many types of Earth science data. The code in this toolbox makes it easy to investigate things like global warming, El Nino, or other major climate-related processes such as how winds affect ocean circulation. Although the toolbox was designed to be used by expert climate scientists, its instruction manual is well written, and beginners may be able to learn a great deal about coding and Earth science, simply by following along with the provided examples. The toolbox is intended to help scientists save time, help them ensure their analysis is accurate, and make it easy for other scientists to repeat the results of previous studies.
The Cumbres Calchaquies Range forms part of the Famatinian metamorphic basement of the Eastern Sierras Pampeanas. The sedimentary protoliths of the metamorphic sequence were deposited in a marine basin alongside the western margin of Gondwana during the Neoproterozoic. New petrologic, geochemical and thermobarometric data give insight into the evolution of the sedimentary basin, its sediment source area, its later metamorphic overprint and its regional relationship to other parts of the Famatinian basement. The metamorphic series studied here consists of banded schists and gneisses and rare calcsilcate-rocks and migmatites that have been reworked by mid-to deep-crustal metamorphic and tectonic processes. The bulk rock compositions indicate shale, wacke, marl and litharenitic protoliths. The metamorphosed elastic sediments have major and trace element compositions indicating a continental granitoid-dominated source area with low sediment recycling. Low SiO2/Al2O3 ratios suggest a relatively low maturity of the sedimentary protoliths. Therefore, the Cumbres Calchaquies section represents a sequence of turbidity currents with progressive shallowing of the depositional environment, as indicated by quartz- and carbonate-rich sediments. The overall data are consistent with the geodynamic environment of a basin adjacent to a continental magmatic arc as the most probable scenario. Whereas the sedimentary protoliths of the metamorphic basement in the Sierra de Ancasti and Sierra de Aconquija, located ca 100-300 km south of the study area are interpreted as originating in an evolving back-arc basin, our results from the Cumbres Calchaquies region indicate a sedimentary source in a felsic continental arc with no significant influx of basic rocks. The Famatinian metamorphic evolution of the Cumbres Calchaquies rocks is of typical Barrow-type, culminating in partial melting of the metasediments. Conventional thermobarometry combined with thermodynamic models (pseudosections) reveal a prograde evolution reaching peak conditions of ca 665 degrees C/6.1 Kbar. This implies a geothermal gradient of ca 35 degrees C/km, which is slightly higher than the average for continental crust and suggests a period of crustal thinning, as known from back-arc basins, or additional heat supply by voluminous intrusions.
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).
Following the Middle Permian (Capitanian) mass extinction there was a global ‘reef eclipse’, and this event had an important role in the Paleozoic-Mesozoic transition of reef ecosystems. Furthermore, the recovery pattern of reef ecosystems in the Wuchiapingian of South China, before the radiation of Changhsingian reefs, is poorly understood. Here, we present a detailed sedimentological account of the Tieqiao section, South China, which records the only known Wuchiapingian reef setting from South China. Six reef growing phases were identified within six transgressive-regressive cycles. The cycles represent changes of deposition in a shallow basin to a subtidal outer platform setting, and the reefal build-ups are recorded in the shallowest part of the cycles above wave base in the euphotic zone. Our results show that the initial reef recovery started from the shallowing up part of the 1st cycle, within the Clarkina leveni conodont zone, which is two conodont zones earlier than previously recognized. In addition, even though metazoans, such as sponges, do become important in the development of the reef bodies, they are not a major component until later in the Wuchiapingian in the 5th and 6th transgressive-regressive cycles. This suggests a delayed recovery of metazoan reef ecosystems following the Middle Permian extinction. Furthermore, even though sponges do become abundant within the reefs, it is the presence and growth of the encrusters Archaeolithoporella and Tubiphytes and abundance of microbial micrites that play an important role in stabilizing the reef structures that form topographic highs.
Vermetid reefs and rocky shores are hot spots of biodiversity, often referred to as the subtropical equivalent of coral reefs. The development of the ecosystem depends on the activity of several reef builders, including red crustose coralline algae (CCA) such as Neogoniolithon brassica-florida. Despite its importance, little is known about Neogoniolithon sp. acclimation to rapid changes in light intensity and corresponding photosynthetic activity. To overcome the large spatial variability in the light field (due to location and the porous nature of the rocks) we grew Neogoniolithon sp. on glass slides and characterized its photosynthetic performance in response to various light intensities by following O-2 exchange and fluorescence parameters. This was also performed on rock-inhabiting thalli collected from the east Mediterranean basin. Generally, maximal photosynthetic rate was reached when Neogoniolithon sp. thalli grown under low illumination (such as in protected niches where the light intensity can be as low as 1% of surface illumination) were examined. When exposed to light intensities higher than those experienced during growth, Neogoniolithon sp. activates adaptive/protective mechanisms such as state transition and nonphotochemical fluorescence quenching and increases the dark respiration thereafter. We find that the Fv/Fm parameter (variable/maximal fluorescence) is not suitable to assess photosynthetic performance in Neogoniolithon sp. and propose using instead an alternative parameter recently developed. Our findings help to clarify why Neogoniolithon sp. is usually observed in shaded niches along the reef surfaces.
This study aims to identify the best-performing site characterization proxy alternative and complementary to the conventional 30 m average shear-wave velocity V-S30, as well as the optimal combination of proxies in characterizing linear site response. Investigated proxies include T-0 (site fundamental period obtained from earthquake horizontal-to-vertical spectral ratios), V-Sz (measured average shear-wave velocities to depth z, z = 5, 10, 20 and 30 m), Z(0.8) and Z(1.0) (measured site depths to layers having shear-wave velocity 0.8 and 1.0 km/s, respectively), as well as Z(x-infer) (inferred site depths from a regional velocity model, x = 0.8 and 1.0, 1.5 and 2.5 km/s). To evaluate the performance of a site proxy or a combination, a total of 1840 surface-borehole recordings is selected from KiK-net database. Site amplifications are derived using surface-to-borehole response-, Fourier- and cross-spectral ratio techniques and then are compared across approaches. Next, the efficacies of 7 single-proxies and 11 proxy-pairs are quantified based on the site-to-site standard deviation of amplification residuals of observation about prediction using the proxy or the pair. Our results show that T-0 is the best-performing single-proxy among T-0, Z(0.8), Z(1.0) and V-Sz. Meanwhile, T-0 is also the best-performing proxy among T-0, Z(0.8), Z(1.0) and Z(x-infer) complementary to V-S30 in accounting for the residual amplification after V-S30-correction. Besides, T-0 alone can capture most of the site effects and should be utilized as the primary site indicator. Though (T-0, V-S30) is the best-performing proxy pair among (V-S30, T-0), (V-S30, Z(0.8)), (V-S30, Z(1.0)), (V-S30, Z(x-infer)) and (T-0, V-Sz), it is only slightly better than (T-0, V-S20). Considering both efficacy and engineering utility, the combination of T-0 (primary) and V-S20 (secondary) is recommended. Further study is needed to test the performances of various proxies on sites in deep sedimentary basins.
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.
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.
During the earliest Triassic microbial mats flourished in the photic zones of marginal seas, generating widespread microbialites. It has been suggested that anoxic conditions in shallow marine environments, linked to the end-Permian mass extinction, limited mat-inhibiting metazoans allowing for this microbialite expansion. The presence of a diverse suite of proxies indicating oxygenated shallow sea-water conditions (metazoan fossils, biomarkers and redox proxies) from microbialite successions have, however, challenged the inference of anoxic conditions. Here, the distribution and faunal composition of Griesbachian microbialites from China, Iran, Turkey, Armenia, Slovenia and Hungary are investigated to determine the factors that allowed microbialite-forming microbial mats to flourish following the end-Permian crisis. The results presented here show that Neotethyan microbial buildups record a unique faunal association due to the presence of keratose sponges, while the Palaeotethyan buildups have a higher proportion of molluscs and the foraminifera Earlandia. The distribution of the faunal components within the microbial fabrics suggests that, except for the keratose sponges and some microconchids, most of the metazoans were transported into the microbial framework via wave currents. The presence of both microbialites and metazoan associations were limited to oxygenated settings, suggesting that a factor other than anoxia resulted in a relaxation of ecological constraints following the mass extinction event. It is inferred that the end-Permian mass extinction event decreased the diversity and abundance of metazoans to the point of significantly reducing competition, allowing photosynthesis-based microbial mats to flourish in shallow water settings and resulting in the formation of widespread microbialites.
The Raman spectra of pure N-2, CO2, and CH4 were analyzed over the range 10 to 500 bars and from -160 degrees C to 200 degrees C (N-2), 22 degrees C to 350 degrees C (CO2), and -100 degrees C to 450 degrees C (CH4). At constant temperature, Raman peak position, including the more intense CO2 peak (nu+), decreases (shifts to lower wave number) with increasing pressure for all three gases over the entire pressure and temperature (PT) range studied. At constant pressure, the peak position for CO2 and CH4 increases (shifts to higher wave number) with increasing temperature over the entire PT range studied. In contrast, N-2 first shows an increase in peak position with increasing temperature at constant pressure, followed by a decrease in peak position with increasing temperature. The inflection temperature at which the trend reverses for N-2 is located between 0 degrees C and 50 degrees C at pressures above similar to 50 bars and is pressure dependent. Below similar to 50 bars, the inflection temperature was observed as low as -120 degrees C. The shifts in Raman peak positions with PT are related to relative density changes, which reflect changes in intermolecular attraction and repulsion. A conceptual model relating the Raman spectral properties of N-2, CO2, and CH4 to relative density (volume) changes and attractive and repulsive forces is presented here. Additionally, reduced temperature-dependent densimeters and barometers are presented for each pure component over the respective PT ranges. The Raman spectral behavior of the pure gases as a function of temperature and pressure is assessed to provide a framework for understanding the behavior of each component in multicomponent N-2-CO2-CH4 gas systems in a future study.
The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressure-temperature (P-T) histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lu-Hf geochronological results for mafic rocks of the Halilbagi Complex (Anatolia) that evolved along different thermal gradients. Samples include a lawsonite-epidote blueschist, a lawsonite-epidote eclogite, and an epidote eclogite (all with counter-clockwise P-T paths), a prograde lawsonite blueschist with a "hairpin"-type P-T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at similar to 0.6-0.7 GPa and 800-850 degrees C, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420 degrees C to 2.6 GPa and 520 degrees C. Well-defined Lu-Hf isochrons were obtained for the epidote eclogite (92.38 +/- 0.22 Ma) and the lawsonite-epidote blueschist (90.19 +/- 0.54 Ma), suggesting rapid garnet growth. The lawsonite-epidote eclogite (87.30 +/- 0.39 Ma) and the prograde lawsonite blueschist (ca. 86 Ma) are younger, whereas the garnet amphibolite (104.5 +/- 3.5 Ma) is older. Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at similar to 104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lu-Hf garnet dates (between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims (as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbagi Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating. (C) 2018, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V.
Fire and grazing shape biodiversity in savannah landscapes. In land use management, knowing the effects of fire and grazing on biodiversity are important in order to ensure environmental sustainability. Beetles specifically are indicators of the biodiversity response to fire and grazing. A grazing exclusion and burning experiment in a split-plot design was used in order to investigate the interacting effects of fire and wildlife grazing on biomass, diversity, and species composition of darkling beetles (Coleoptera, Tenebrionidae) over time after fire. Darkling beetle species richness and diversity were responding in a three-way-interaction to fire, grazing, and time after fire, whereby biomass of darkling beetles remained unaffected and species compositional changes were attributed to seasonal changes of time only. Fire on ungrazed plots had a negative effect on species diversity and richness 2 weeks and 6 months post fire, whereas fire on grazed plots had no impact on species diversity and richness. Grazing only lowered species diversity and richness 6 months after fire treatments. Results suggest that grazing overrides the effects of fire and that the similar effects caused by fire and grazing are due to niche and assemblage simplification of the habitat.
In crop modeling and yield predictions, the heterogeneity of agricultural landscapes is usually not accounted for. This heterogeneity often arises from landscape elements like forests, hedges, or single trees and shrubs that cast shadows. Shading from forested areas or shrubs has effects on transpiration, temperature, and soil moisture, all of which affect the crop yield in the adjacent arable land. Transitional gradients of solar irradiance can be described as a function of the distance to the zero line (edge), the cardinal direction, and the height of trees. The magnitude of yield reduction in transition zones is highly influenced by solar irradiance-a factor that is not yet implemented in crop growth models on a landscape level. We present a spatially explicit model for shading caused by forested areas, in agricultural landscapes. With increasing distance to forest, solar irradiance and yield increase. Our model predicts that the shading effect from the forested areas occurs up to 15 m from the forest edge, for the simulated wheat yields, and up to 30 m, for simulated maize. Moreover, we estimated the spatial extent of transition zones, to calculate the regional yield reduction caused by shading of the forest edges, which amounted to 5% to 8% in an exemplary region.
Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben
(2019)
The Upper Rhine Graben (URG) is an active rift with a high geothermal potential. Despite being a well-studied area, the three-dimensional interaction of the main controlling factors of the thermal and hydraulic regime is still not fully understood. Therefore, we have used a data-based 3D structural model of the lithological configuration of the central URG for some conceptual numerical experiments of 3D coupled simulations of fluid and heat transport. To assess the influence of the main faults bordering the graben on the hydraulic and the deep thermal field, we carried out a sensitivity analysis on fault width and permeability. Depending on the assigned width and permeability of the main border faults, fluid velocity and temperatures are affected only in the direct proximity of the respective border faults. Hence, the hydraulic characteristics of these major faults do not significantly influence the graben-wide groundwater flow patterns. Instead, the different scenarios tested provide a consistent image of the main characteristics of fluid and heat transport as they have in common: (1) a topography-driven basin-wide fluid flow perpendicular to the rift axis from the graben shoulders to the rift center, (2) a N/NE-directed flow parallel to the rift axis in the center of the rift and, (3) a pronounced upflow of hot fluids along the rift central axis, where the streams from both sides of the rift merge. This upflow axis is predicted to occur predominantly in the center of the URG (northern and southern model area) and shifted towards the eastern boundary fault (central model area).
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.
Above and underground hydrological processes depend on soil moisture (SM) variability, driven by different environmental factors that seldom are well-monitored, leading to a misunderstanding of soil water temporal patterns. This study investigated the stability of the SM temporal dynamics to different monitoring temporal resolutions around the border between two soil types in a tropical watershed. Four locations were instrumented in a small-scale watershed (5.84 km(2)) within the tropical coast of Northeast Brazil, encompassing different soil types (Espodossolo Humiluvico or Carbic Podzol, and Argissolo Vermelho-Amarelo or Haplic Acrisol), land covers (Atlantic Forest, bush vegetation, and grassland) and topographies (flat and moderate slope). The SM was monitored at a temporal resolution of one hour along the 2013-2014 hydrological year and then resampled a resolutions of 6 h, 12 h, 1 day, 2 days, 4 days, 7 days, and 15 days. Descriptive statistics, temporal variability, time-stability ranking, and hierarchical clustering revealed uneven associations among SM time components. The results show that the time-invariant component ruled SM temporal variability over the time-varying parcel, either at high or low temporal resolutions. Time-steps longer than 2 days affected the mean statistical metrics of the SM time-variant parcel. Additionally, SM at downstream and upstream sites behaved differently, suggesting that the temporal mean was regulated by steady soil properties (slope, restrictive layer, and soil texture), whereas their temporal anomalies were driven by climate (rainfall) and hydrogeological (groundwater level) factors. Therefore, it is concluded that around the border between tropical soil types, the distinct behaviour of time-variant and time-invariant components of SM time series reflects different combinations of their soil properties.
Metamorphic geology
(2019)
Encroachment of shrubs into the unique pastoral grassland ecosystems of the Tibetan Plateau has significant impact on ecosystem services, especially forage production. We developed a process-based ecohydrological model to identify the relative importance of the main drivers of shrub encroachment for the alpine meadows within the Qinghai province. Specifically, we explored the effects of summer livestock grazing (intensity and type of livestock) together with the effects of climate warming, including interactions between herbaceous and woody vegetation and feedback loops between soil, water and vegetation. Under current climatic conditions and a traditional herd composition, an increasing grazing intensity above a threshold value of 0.32 +/- 0.10 large stock units (LSU) ha(-1) day(-1) changes the vegetation composition from herbaceous towards a woody and bare soil dominated system. Very high grazing intensity (above 0.8 LSU ha(-1) day(-1)) leads to a complete loss of any vegetation. Under warmer conditions, the vegetation showed a higher resilience against livestock farming. This resilience is enhanced when the herd has a higher browser : grazer ratio. A cooler climate has a shrub encroaching effect, whereas warmer conditions increase the cover of the herbaceous vegetation. This effect was primarily due to season length and an accompanied competitive loss of slower growing shrubs, rather than evaporative water loss leading to less soil water in deeper soil layers for deeper rooting shrubs. If climate warming is driving current shrub encroachment, we conclude it is only indirectly so. It would be manifest by an advancing shrubline and could be regarded as a climatic escape of specific shrub species such as Potentilla fruticosa. Under the recent high intensity of grazing, only herding by more browsing animals can potentially prevent both shrub encroachment and the complete loss of herbaceous vegetation.
The Schwarzenberg mining district in the western Erzgebirge hosts numerous skarn-hosted tin-polymetallic deposits, such as Breitenbrunn. The St. Christoph mine is located in the Breitenbrunn deposit and is the locus typicus of christophite, an iron-rich sphalerite variety, which can be associated with indium enrichment. This study presents a revision of the paragenetic scheme, a contribution to the indium behavior and potential, and discussion on the origin of the sulfur. This was achieved through reflected light microscopy, SEM-based MLA, EPMA, and bulk mineral sulfur isotope analysis on 37 sulfide-rich skarn samples from a mineral collection. The paragenetic scheme includes: a pre-mineralization stage of anhydrous calc-silicates and hydrous minerals; an oxide stage, dominated by magnetite; a sulfide stage of predominantly sphalerite, minor pyrite, chalcopyrite, arsenopyrite, and galena. Some sphalerite samples present elevated indium contents of up to 0.44 wt%. Elevated iron contents (4-10 wt%) in sphalerite can be tentatively linked to increased indium incorporation, but further analyses are required. Analyzed sulfides exhibit homogeneous delta S-34 values (-1 to +2 parts per thousand VCDT), assumed to be post-magmatic. They correlate with other Fe-Sn-Zn-Cu-In skarn deposits in the western Erzgebirge, and Permian vein-hosted associations throughout the Erzgebirge region.
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.
Water fluxes in highly impounded regions are heavily dependent on reservoir properties. However, for large and remote areas, this information is often unavailable. In this study, the geometry and volume of small surface reservoirs in the semi-arid region of Brazil were estimated using terrain and shape attributes extracted by remote sensing. Regression models and data classification were used to predict the volumes, at different water stages, of 312 reservoirs for which topographic information is available. The power function used to describe the reservoir shapes tends to overestimate the volumes; therefore, a modified shape equation was proposed. Among the methods tested, four were recommended based on performance and simplicity, for which the mean absolute percentage errors varied from 24 to 39%, in contrast to the 94% error achieved with the traditional method. Despite the challenge of precisely deriving the flooded areas of reservoirs, water management in highly reservoir-dense environments should benefit from volume prediction based on remote sensing.
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
Groundwater travel time distributions (TTDs) provide a robust description of the subsurface mixing behavior and hydrological response of a subsurface system. Lagrangian particle tracking is often used to derive the groundwater TTDs. The reliability of this approach is subjected to the uncertainty of external forcings, internal hydraulic properties, and the interplay between them. Here, we evaluate the uncertainty of catchment groundwater TTDs in an agricultural catchment using a 3-D groundwater model with an overall focus on revealing the relationship between external forcing, internal hydraulic properties, and TTD predictions. Eight recharge realizations are sampled from a high-resolution dataset of land surface fluxes and states. Calibration-constrained hydraulic conductivity fields (Ks fields) are stochastically generated using the null-space Monte Carlo (NSMC) method for each recharge realization. The random walk particle tracking (RWPT) method is used to track the pathways of particles and compute travel times. Moreover, an analytical model under the random sampling (RS) assumption is fit against the numerical solutions, serving as a reference for the mixing behavior of the model domain. The StorAge Selection (SAS) function is used to interpret the results in terms of quantifying the systematic preference for discharging young/old water. The simulation results reveal the primary effect of recharge on the predicted mean travel time (MTT). The different realizations of calibration-constrained Ks fields moderately magnify or attenuate the predicted MTTs. The analytical model does not properly replicate the numerical solution, and it underestimates the mean travel time. Simulated SAS functions indicate an overall preference for young water for all realizations. The spatial pattern of recharge controls the shape and breadth of simulated TTDs and SAS functions by changing the spatial distribution of particles' pathways. In conclusion, overlooking the spatial nonuniformity and uncertainty of input (forcing) will result in biased travel time predictions. We also highlight the worth of reliable observations in reducing predictive uncertainty and the good interpretability of SAS functions in terms of understanding catchment transport processes.
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.
Shear-waves are the most energetic body-waves radiated from an earthquake, and are responsible for the destruction of engineered structures. In both short-term emergency response and long-term risk forecasting of disaster-resilient built environment, it is critical to predict spatially accurate distribution of shear-wave amplitudes. Although decades’ old theory proposes a deterministic, highly anisotropic, four-lobed shear-wave radiation pattern, from lack of convincing evidence, most empirical ground-shaking prediction models settled for an oversimplified stochastic radiation pattern that is isotropic on average. Today, using the large datasets of uniformly processed seismograms from several strike, normal, reverse, and oblique-slip earthquakes across the globe, compiled specifically for engineering applications, we could reveal, quantify, and calibrate the frequency-, distance-, and style-of-faulting dependent transition of shear-wave radiation between a stochastic-isotropic and a deterministic-anisotropic phenomenon. Consequent recalibration of empirical ground-shaking models dramatically improved their predictions: with isodistant anisotropic variations of ±40%, and 8% reduction in uncertainty. The outcomes presented here can potentially trigger a reappraisal of several practical issues in engineering seismology, particularly in seismic ground-shaking studies and seismic hazard and risk assessment.
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.
Quantifying erosion rates, and how they compare to rock uplift rates, is fundamental for understanding landscape response to tectonics and associated sediment fluxes from upland areas. The erosional response to uplift is well-represented by river incision and the associated landslide activity. However, characterising the relationship between these processes remains a major challenge in tectonically active areas, in some cases because landslides can preclude obtaining reliable erosion rates from cosmogenic radionuclide (CRN) concentrations. Here, we quantify the control of tectonics and its coupled geomorphic response on the erosion rates of catchments in southern Italy that are experiencing a transient response to normal faulting. We analyse in-situ Be-10 concentrations for detrital sediment samples, collected along the strike of faults with excellent tectonic constraints and landslide inventories. We demonstrate that Be-10-derived erosion rates are controlled by fault throw rates and the extent of transient incision and associated landsliding in the catchments. We show that the low-relief sub-catchments above knickpoints erode at uniform background rates of similar to 0.10 mm/yr, while downstream of knickpoints, erosion removes similar to 50% of the rock uplifted by the faults, at rates of 0.10-0.64 mm/yr. Despite widespread landsliding, CRN samples provide relatively consistent and accurate erosion rates, most likely because landslides are frequent, small, and shallow, and represent the integrated record of landsliding over several seismic cycles. Consequently, we combine these validated Be-10 erosion rates and data from a geomorphological landslide inventory in a published numerical model, to gain further insight into the long-term landslide rates and sediment mixing, highlighting the potential of CRN data to study landslide dynamics. (C) 2018 Elsevier B.V. All rights reserved.
Several hydraulic fracturing tests were performed in boreholes located in central Hungary in order to determine the in-situ stress for a geological site investigation. At a depth of about 540m, the observed pressure versus time curves in mica schist with low dip angle foliation shows atypical pressure versus time results. After each pressurization cycle, the fracture breakdown pressure in the first fracturing cycle is lower than the refracturing or reopening pressure in the subsequent pressurizations. It is assumed that the viscosity of the drilling mud and observed foliation of the mica schist have a significant influence on the pressure values. In order to study this problem, numerical modeling was performed using the distinct element code particle flow code, which has been proven to be a valuable tool to investigate rock engineering problems such as hydraulic fracturing. The two-dimensional version of the code applied in this study can simulate hydro-mechanically coupled fluid flow in crystalline rock with low porosity and pre-existing fractures. In this study, the effect of foliation angle and fluid viscosity on the peak pressure is tested. The atypical characteristics of the pressure behaviour are interpreted so that mud with higher viscosity penetrates the sub-horizontal foliation plane, blocks the plane of weakness and makes the partly opened fracture tight and increase the pore pressure which decreases slowly with time. We see this viscous blocking effect as one explanation for the observed increase in fracture reopening pressure in subsequent pressurization cycles.
The occurrence of mounds dominated by siliceous sponges and microbialites is often related to distal, deep settings of middle ramps and shelves. This paper presents evidence for Bajocian (Garanliana garantiana Zone) microbial-siliceous sponge mounds formed in open marine but relatively shallow settings of a ramp from the Iberian Basin of eastern Spain. Marked differences in mound spacing, morphology, and composition of the related intermound facies are observed from distal to more proximal settings. The distal (below storm wave base) settings are characterized by alternating tabular-bedded marls and limestones rich in pelagic fossils (ammonites, belemnites), open-marine thin-shelled bivalves (Bositra-like), as well as peloids, which include widely or randomly spaced isolated, small (up to 0.4 m high) and larger (up to 2.5 m high) mounds with upward accretion. The intermediate (near to above storm wave base) settings show tabular, thickened beds of peloidal and/or intraclastic limestones with closely spaced mounds (similar to 1 m high), which often coalesce laterally, forming extensive lenticular structures (up to 10 m wide). The proximal (above storm wave base) depositional settings consist of tabular to irregular beds of intraclastic limestones with widely spaced small (up to 0.4 m high) mounds with mainly tabular geometries. The mound framework contains variable proportions of microbialites (dense to clotted peloidal thrombolitic fabrics) and siliceous sponges (hexactinellids and lithistids in similar proportion) ranging from planar to conic shapes. These morphological and compositional changes allow characterizing three shallowing-upward sequences (sequences 1-3) developed in the overall regressive trend of a basin-wide, upper Bajocian T-R cycle. Episodic wave reworking of the early-cemented mounds resulted in the formation of peloids, small rounded intraclasts, and large, rounded or subangular intraclasts. These nonskeletal micritic grains show internal fabrics related to those of the mound and/or microbialites. A progressive textural gradation towards greater size and lesser roundness of the nonskeletal grains in the areas in the vicinity of the main mound factory is documented (i.e., from large, subangular intraclasts in the areas close to the main mound factory to peloids in the areas that are far from it). We discuss the alternative model of internal waves (instead of storm-induced waves) as the hydrodynamic agent providing the high-energy events needed to explain the origin of the peloidal-intraclastic intermound facies and, likely, also the nutrients needed by the microbialites and siliceous sponges to grow.
Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity
(2019)
The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity.
The northward indentation of the Pamir salient into the Tarim basin at the western syntaxis of the India-Asia collision zone is the focus of controversial models linking lithospheric to surface and atmospheric processes. Here we report on tectonic events recorded in the most complete and best-dated sedimentary sequences from the western Tarim basin flanking the eastern Pamir (the Aertashi section), based on sedimentologic, provenance, and magnetostratigraphic analyses. Increased tectonic subsidence and a shift from marine to continental fluvio-deltaic deposition at 41Ma indicate that far-field deformation from the south started to affect the Tarim region. A sediment accumulation hiatus from 24.3 to 21.6Ma followed by deposition of proximal conglomerates is linked to fault propagation into the Tarim basin. From 21.6 to 15.0Ma, increasing accumulation rates of fining upward clastics is interpreted as the expression of a major dextral transtensional system linking the Kunlun to the Tian Shan ahead of the northward Pamir indentation. At 15.0Ma, the appearance of North Pamir-sourced conglomerates followed at 11Ma by Central Pamir-sourced volcanics coincides with a shift to E-W compression, clockwise vertical-axis rotations and the onset of growth strata associated with the activation of the local east vergent Qimugen thrust wedge. Together, this enables us to interpret that Pamir indentation into Tarim had started by 24.3Ma, reached the study location by 15.0Ma and had passed it by 11Ma, providing kinematic constraints on proposed tectonic models involving intracontinental subduction and delamination.
The improvement of process representations in hydrological models is often only driven by the modelers' knowledge and data availability. We present a comprehensive comparison between two hydrological models of different complexity that is developed to support (1) the understanding of the differences between model structures and (2) the identification of the observations needed for model assessment and improvement. The comparison is conducted on both space and time and by aggregating the outputs at different spatiotemporal scales. In the present study, mHM, a process‐based hydrological model, and ParFlow‐CLM, an integrated subsurface‐surface hydrological model, are used. The models are applied in a mesoscale catchment in Germany. Both models agree in the simulated river discharge at the outlet and the surface soil moisture dynamics, lending their supports for some model applications (drought monitoring). Different model sensitivities are, however, found when comparing evapotranspiration and soil moisture at different soil depths. The analysis supports the need of observations within the catchment for model assessment, but it indicates that different strategies should be considered for the different variables. Evapotranspiration measurements are needed at daily resolution across several locations, while highly resolved spatially distributed observations with lower temporal frequency are required for soil moisture. Finally, the results show the impact of the shallow groundwater system simulated by ParFlow‐CLM and the need to account for the related soil moisture redistribution. Our comparison strategy can be applied to other models types and environmental conditions to strengthen the dialog between modelers and experimentalists for improving process representations in Earth system models.
The Value of Empirical Data for Estimating the Parameters of a Sociohydrological Flood Risk Model
(2019)
In this paper, empirical data are used to estimate the parameters of a sociohydrological flood risk model. The proposed model, which describes the interactions between floods, settlement density, awareness, preparedness, and flood loss, is based on the literature. Data for the case study of Dresden, Germany, over a period of 200years, are used to estimate the model parameters through Bayesian inference. The credibility bounds of their estimates are small, even though the data are rather uncertain. A sensitivity analysis is performed to examine the value of the different data sources in estimating the model parameters. In general, the estimated parameters are less biased when using data at the end of the modeled period. Data about flood awareness are the most important to correctly estimate the parameters of this model and to correctly model the system dynamics. Using more data for other variables cannot compensate for the absence of awareness data. More generally, the absence of data mostly affects the estimation of the parameters that are directly related to the variable for which data are missing. This paper demonstrates that combining sociohydrological modeling and empirical data gives additional insights into the sociohydrological system, such as quantifying the forgetfulness of the society, which would otherwise not be easily achieved by sociohydrological models without data or by standard statistical analysis of empirical data.
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.
Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.
Terrestrial gravimetry is increasingly used to monitor mass transport processes in geophysics boosted by the ongoing technological development of instruments. Resolving a particular phenomenon of interest, however, requires a set of gravity corrections of which the uncertainties have not been addressed up to now. In this study, we quantify the time domain uncertainty of tide, global atmospheric, large-scale hydrological, and nontidal ocean loading corrections. The uncertainty is assessed by comparing the majority of available global models for a suite of sites worldwide. The average uncertainty expressed as root-mean-square error equals 5.1nm/s(2), discounting local hydrology or air pressure. The correction-induced uncertainty of gravity changes over various time periods of interest ranges from 0.6nm/s(2) for hours up to a maximum of 6.7nm/s(2) for 6months. The corrections are shown to be significant and should be applied for most geophysical applications of terrestrial gravimetry. From a statistical point of view, however, resolving subtle gravity effects in the order of few nanometers per square second is challenged by the uncertainty of the corrections. Plain Language Summary Many scientists are exploring ways to benefit from gravity measurements in fields of high societal relevance such as monitoring of volcanoes or measuring the amount of water in underground. Any application of such new methods, however, requires careful preparation of the gravity measurements. The intention of the preparation process is to ensure that the measurements do not contain information about processes that are not of interest. For that reason, the influence of atmosphere, ocean, tides, and hydrology needs to be reduced from the gravity. In this study, we investigate how this reduction process influences the quality of the measurement. We found that the precision degrades especially owing to the hydrology. The ocean plays an important role at sites close to the coast and the atmosphere at sites located in mountains. The overall errors of the reductions may complicate a reliable use of gravity measurements in certain studies focusing on very small signals. Nevertheless, the precision of gravity reductions alone does not obstruct a meaningful use of gravity measurements in most research fields. Details specifying the reduction precision are provided in this study allowing scientist dealing with gravity measurements to decide if their signal of interest can be reliably resolved.
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.
Numerous preflight investigations were necessary prior to the exposure experiment BIOMEX on the International Space Station to test the basic potential of selected microorganisms to resist or even to be active under Mars-like conditions. In this study, methanogenic archaea, which are anaerobic chemolithotrophic microorganisms whose lifestyle would allow metabolism under the conditions on early and recent Mars, were analyzed. Some strains from Siberian permafrost environments have shown a particular resistance. In this investigation, we analyzed the response of three permafrost strains (Methanosarcina soligelidi SMA-21, Candidatus Methanosarcina SMA-17, Candidatus Methanobacterium SMA-27) and two related strains from non-permafrost environments (Methanosarcina mazei, Methanosarcina barkeri) to desiccation conditions (-80 degrees C for 315 days, martian regolith analog simulants S-MRS and P-MRS, a 128-day period of simulated Mars-like atmosphere). Exposure of the different methanogenic strains to increasing concentrations of magnesium perchlorate allowed for the study of their metabolic shutdown in a Mars-relevant perchlorate environment. Survival and metabolic recovery were analyzed by quantitative PCR, gas chromatography, and a new DNA-extraction method from viable cells embedded in S-MRS and P-MRS. All strains survived the two Mars-like desiccating scenarios and recovered to different extents. The permafrost strain SMA-27 showed an increased methanogenic activity by at least 10-fold after deep-freezing conditions. The methanogenic rates of all strains did not decrease significantly after 128 days S-MRS exposure, except for SMA-27, which decreased 10-fold. The activity of strains SMA-17 and SMA-27 decreased after 16 and 60 days P-MRS exposure. Non-permafrost strains showed constant survival and methane production when exposed to both desiccating scenarios. All strains showed unaltered methane production when exposed to the perchlorate concentration reported at the Phoenix landing site (2.4 mM) or even higher concentrations. We conclude that methanogens from (non-)permafrost environments are suitable candidates for potential life in the martian subsurface and therefore are worthy of study after space exposure experiments that approach Mars-like surface conditions.
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.
Interannual variations in ecosystem primary productivity are dominated by water availability. Until recently, characterizing the photosynthetic response of different ecosystems to soil moisture anomalies was hampered by observational limitations. Here, we use a number of satellite-based proxies for productivity, including spectral indices, sun-induced chlorophyll fluorescence, and data-driven estimates of gross primary production, to reevaluate the relationship between terrestrial photosynthesis and water. In contrast to nonwoody vegetation, we find a resilience of forested ecosystems to reduced soil moisture. Sun-induced chlorophyll fluorescence and data-driven gross primary production indicate an increase in photosynthesis as a result of the accompanying higher amounts of light and temperature despite lowered light-use-efficiency. Conversely, remote sensing indicators of greenness reach their detection limit and largely remain stable. Our study thus highlights the differential responses of ecosystems along a tree cover gradient and illustrates the importance of differentiating photosynthesis indicators from those of greenness for the monitoring and understanding of ecosystems. Plain Language Summary The capacity of vegetation to thrive and to sequester carbon depends on how much water they can have access to. In this work, we evaluate how different types of satellite observations can describe the response of vegetation to changes in soil moisture over the entire planet. The first source of observation measures only the greenness of the land surface, the second measures light that is emitted by pigments in plants which are photosynthetically active (chlorophyll fluorescence), and the third are simulations of gross carbon uptake derived from machine learning techniques. For periods of water shortage all three indicate a reduction of growth in ecosystems with few trees. However, in cold boreal forests, when soil moisture is particularly low, we still detect an increase in photosynthesis due to higher light and temperature conditions, but this is not reflected in the greenness indicator. This work illustrates how lack of water is not necessarily harmful for catching carbon through photosynthesis, but to monitor this effect, we need remote sensing indicators that measure more than just how green the plants are, and fluorescence is likely a good candidate.
Breakup Without Borders
(2019)
Relative plate motions during continental rifting result from the interplay of local with far-field forces. Here we study the dynamics of rifting and breakup using large-scale numerical simulations of mantle convection with self-consistent evolution of plate boundaries. We show that continental separation follows a characteristic evolution with four distinctive phases: (1) an initial slow rifting phase with low divergence velocities and maximum tensional stresses, (2) a synrift speed-up phase featuring an abrupt increase of extension rate with a simultaneous drop of tensional stress, (3) the breakup phase with inception of fast sea-floor spreading, and (4) a deceleration phase occurring in most but not all models where extensional velocities decrease. We find that the speed-up during rifting is compensated by subduction acceleration or subduction initiation even in distant localities. Our study illustrates new links between local rift dynamics, plate motions, and subduction kinematics during times of continental separation.
Changes in the steepness of river profiles or abrupt vertical steps (i.e. waterfalls) are thought to be indicative of changes in erosion rates, lithology or other factors that affect landscape evolution. These changes are referred to as knickpoints or knickzones and are pervasive in bedrock river systems. Such features are thought to reveal information about landscape evolution and patterns of erosion, and therefore their locations are often reported in the geomorphic literature. It is imperative that studies reporting knickpoints and knickzones use a reproducible method of quantifying their locations, as their number and spatial distribution play an important role in interpreting tectonically active landscapes. In this contribution we introduce a reproducible knickpoint and knickzone extraction algorithm that uses river profiles transformed by integrating drainage area along channel length (the so-called integral or chi method). The profile is then statistically segmented and the differing slopes and step changes in the elevations of these segments are used to identify knickpoints, knickzones and their relative magnitudes. The output locations of identified knickpoints and knickzones compare favourably with human mapping: we test the method on Santa Cruz Island, CA, using previously reported knickzones and also test the method against a new dataset from the Quadrilatero Ferrifero in Brazil. The algorithm allows for the extraction of varying knickpoint morphologies, including stepped, positive slope-break (concave upward) and negative slope-break knickpoints. We identify parameters that most affect the resulting knickpoint and knickzone locations and provide guidance for both usage and outputs of the method to produce reproducible knickpoint datasets.
Climate and weather conditions in the mid-latitudes are strongly driven by the large-scale atmosphere circulation. Observational data indicate that important components of the large-scale circulation have changed in recent decades, including the strength and the width of the Hadley cell, jets, storm tracks and planetary waves. Here, we use a new statistical-dynamical atmosphere model (SDAM) to test the individual sensitivities of the large-scale atmospheric circulation to changes in the zonal temperature gradient, meridional temperature gradient and global-mean temperature. We analyze the Northern Hemisphere Hadley circulation, jet streams, storm tracks and planetary waves by systematically altering the zonal temperature asymmetry, the meridional temperature gradient and the global-mean temperature. Our results show that the strength of the Hadley cell, storm tracks and jet streams depend, in terms of relative changes, almost linearly on both the global-mean temperature and the meridional temperature gradient, whereas the zonal temperature asymmetry has little or no influence. The magnitude of planetary waves is affected by all three temperature components, as expected from theoretical dynamical considerations. The width of the Hadley cell behaves nonlinearly with respect to all three temperature components in the SDAM. Moreover, some of these observed large-scale atmospheric changes are expected from dynamical equations and are therefore an important part of model validation.
Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.
Higher resilience to climatic disturbances in tropical vegetation exposed to more variable rainfall
(2019)
With ongoing global warming, the amount and frequency of precipitation in the tropics is projected to change substantially. While it has been shown that tropical forests and savannahs are sustained within the same intermediate mean annual precipitation range, the mechanisms that lead to the resilience of these ecosystems are still not fully understood. In particular, the long-term impact of rainfall variability on resilience is as yet unclear. Here we present observational evidence that both tropical forest and savannah exposed to a higher rainfall variability-in particular on interannual scales-during their long-term past are overall more resilient against climatic disturbances. Based on precipitation and tree cover data in the Brazilian Amazon basin, we constructed potential landscapes that enable us to systematically measure the resilience of the different ecosystems. Additionally, we infer that shifts from forest to savannah due to decreasing precipitation in the future are more likely to occur in regions with a precursory lower rainfall variability. Long-term rainfall variability thus needs to be taken into account in resilience analyses and projections of vegetation response to climate change.
Remote Sensing technologies allow to map biophysical, biochemical, and earth surface parameters of the land surface. Of especial interest for various applications in environmental and urban sciences is the combination of spectral and 3D elevation information. However, those two data streams are provided separately by different instruments, namely airborne laser scanner (ALS) for elevation and a hyperspectral imager (HSI) for high spectral resolution data. The fusion of ALS and HSI data can thus lead to a single data entity consistently featuring rich structural and spectral information. In this study, we present the application of fusing the first pulse return information from ALS data at a sub-decimeter spatial resolution with the lower-spatial resolution hyperspectral information available from the HSI into a hyperspectral point cloud (HSPC). During the processing, a plausible hyperspectral spectrum is assigned to every first-return ALS point. We show that the complementary implementation of spectral and 3D information at the point-cloud scale improves object-based classification and information extraction schemes. This improvements have great potential for numerous land cover mapping and environmental applications.
Private precaution is an important component in contemporary flood risk management and climate adaptation. However, quantitative knowledge about vulnerability reduction via private precautionary measures is scarce and their effects are hardly considered in loss modeling and risk assessments. However, this is a prerequisite to enable temporally dynamic flood damage and risk modeling, and thus the evaluation of risk management and adaptation strategies. To quantify the average reduction in vulnerability of residential buildings via private precaution empirical vulnerability data (n = 948) is used. Households with and without precautionary measures undertaken before the flood event are classified into treatment and nontreatment groups and matched. Postmatching regression is used to quantify the treatment effect. Additionally, we test state-of-the-art flood loss models regarding their capability to capture this difference in vulnerability. The estimated average treatment effect of implementing private precaution is between 11 and 15 thousand EUR per household, confirming the significant effectiveness of private precautionary measures in reducing flood vulnerability. From all tested flood loss models, the expert Bayesian network-based model BN-FLEMOps and the rule-based loss model FLEMOps perform best in capturing the difference in vulnerability due to private precaution. Thus, the use of such loss models is suggested for flood risk assessments to effectively support evaluations and decision making for adaptable flood risk management.
Using optical absorption and Raman spectroscopic measurements, in conjunction with the first-principles calculations, a pressure-induced high-spin (HS)-to-low-spin (LS) state electronic transition of Fe2+ (M2-octahedral site) was resolved around 76-80GPa in a natural triphylite-lithiophilite sample with chemical composition (LiFe0.7082+Mn0.292PO4)-Li-M1-Fe-M2 (theoretical composition (LiFe0.52+Mn0.5PO4)-Li-M1-Fe-M2). The optical absorption spectra at ambient conditions consist of a broad doublet band with two constituents (1) (similar to 9330cm(-1)) and (2) (similar to 7110cm(-1)), resulting from the electronic spin-allowed transition (T2gEg)-T-5-E-5 of octahedral (HSFe2+)-Fe-M2. Both (1) and (2) bands shift non-linearly with pressure to higher energies up to similar to 55GPa. In the optical absorption spectrum measured at similar to 81GPa, the aforementioned HS-related bands disappear, whereas a new broadband with an intensity maximum close to 16,360cm(-1) appears, superimposed on the tail of the high-energy ligand-to-metal O2-Fe2+ charge-transfer absorption edge. We assign this new band to the electronic spin-allowed dd-transition (1)A(1g)(1)T(1g) of LS Fe2+ in octahedral coordination. The high-pressure Raman spectra evidence the Fe2+ HS-to-LS transition mainly from the abrupt shift of the P-O symmetric stretching modes to lower frequencies at similar to 76GPa, the highest pressure achieved in the Raman spectroscopic experiments. Calculations indicated that the presence of Mn-M2(2+) simply shifts the isostructural HS-to-LS transition to higher pressures compared to the triphylite Fe-M2(2+) end-member, in qualitative agreement with our experimental observations.
Ecosystem services inherently involve people, whose values help define the benefits of nature's services. It is thus important for researchers to involve stakeholders in ecosystem services research. However, a simple and practicable framework to guide such engagement, and in particular to help researchers anticipate and consider key issues and challenges, has not been well explored. Here, we use experience from the 12 case studies in the European Operational Potential of Ecosystem Research Applications (OPERAs) project to propose a stakeholder engagement framework comprising three key elements: creating space, aligning motivations, and building trust. We argue that involving stakeholders in research demands thoughtful reflection from the researchers about what kind of space they want to create, including if and how they want to bring different interests together, how much space they want to allow for critical discussion, and whether there is a role for particular stakeholders to serve as conduits between others. In addition, understanding their own motivations—including values, knowledge, goals, and desired benefits—will help researchers decide when and how to involve stakeholders, identify areas of common ground and potential disagreement, frame the project appropriately, set expectations, and ensure each party is able to see benefits of engaging with each other. Finally, building relationships with stakeholders can be difficult but considering the roles of existing relationships, time, approach, reputation, and belonging can help build mutual trust. Although the three key elements and the paths between them can play out differently depending on the particular research project, we suggest that a research design that considers how to create the space in which researchers and stakeholders will meet, align motivations between researchers and stakeholders, and build mutual trust will help foster productive researcher–stakeholder relationships.
Growing attention to phytoplankton mixotrophy as a trophic strategy has led to significant revisions of traditional pelagic food web models and ecosystem functioning. Although some empirical estimates of mixotrophy do exist, a much broader set of in situ measurements are required to (i) identify which organisms are acting as mixotrophs in real time and to (ii) assess the contribution of their heterotrophy to biogeochemical cycling. Estimates are needed through time and across space to evaluate which environmental conditions or habitats favour mixotrophy: conditions still largely unknown. We review methodologies currently available to plankton ecologists to undertake estimates of plankton mixotrophy, in particular nanophytoplankton phago-mixotrophy. Methods are based largely on fluorescent or isotopic tracers, but also take advantage of genomics to identify phylotypes and function. We also suggest novel methods on the cusp of use for phago-mixotrophy assessment, including single-cell measurements improving our capacity to estimate mixotrophic activity and rates in wild plankton communities down to the single-cell level. Future methods will benefit from advances in nanotechnology, micromanipulation and microscopy combined with stable isotope and genomic methodologies. Improved estimates of mixotrophy will enable more reliable models to predict changes in food web structure and biogeochemical flows in a rapidly changing world.
We studied bacterial abundance and community structure of five soil cores using high-throughput sequencing of the 16S rRNA gene. Shifts in the soil bacterial composition were more pronounced within a vertical profile than across the landscape. Soil organic carbon (SOC) and nitrogen (N) concentrations decreased exponentially with soil depth and revealed a buried carbon-rich horizon between 0.8 and 1.3 m across all soil cores. This buried horizon was phylogenetically similar to its surrounding subsoils supporting the idea that the type of carbon, not necessarily the amount of carbon was driving the apparent similarities. In contrast to other studies, Nitrospirae was one of our major phyla with relatively high abundances throughout the soil profile except for the surface soil. Although depth is the major driver shaping soil bacterial community structure, positive correlations with SOC and N concentrations, however, were revealed with the bacterial abundance of Acidobacteria, one of the major, and Gemmatimonadetes, one of the minor phyla in our study. Our study showed that bacterial diversity in soils below 2.0 m can be still as high if not higher than in the above laying subsurface soil suggesting that various bacteria throughout the soil profile influence major biogeochemical processes in floodplain soils.
The early Cretaceous Koegel Fontein intrusive complex is situated near the Atlantic coast in South Africa, about 350 km northwest of Cape Town. The complex comprises felsic units of granite and syenite with compositionally related dykes, and a single intrusive plug of diorite. Existing zircon U-Pb ages of 144 +/- 2 Ma for the syenite and 133.9 +/- 1.3 Ma for the granite suggest that the emplacement of the complex took place over a period of about 10 My. This study provides additional and independent ages of the Koegel Fontein complex by Ar-40/Ar-39 dating to confirm the onset and duration of magmatism and better define the sequence of igneous units that comprise it. New laser step-heating Ar-40/Ar-3(9) ages on plagioclase and biotite from the main intrusive units in the complex are presented here, including samples previously dated by U-Pb dating. The Ar-40/Ar-39 ages for the granite and syenite units (131.1 +/- 0.9 Ma and 143.3 +/- 0.9, respectively) are in good agreement with the zircon U-Pb ages. Other units not previously dated include the Rooivleitjie alkaline granite (150.7 +/- 0.6 Ma), two quartz-porphyry dykes (143.0 +/- 0.9 and 139.4 +/- 1.7 Ma) and the Zout Rivier diorite plug (133.0 +/- 1.0 Ma). The new results confirm an early onset of magmatism at Koegel Fontein relative to that of the Etendeka Province some 1000 km to the north, which is consistent with the regional south-to-north propagation of South Atlantic rifting. The youngest Ar-40/Ar-3(9) ages at Koegel Fontein (134 to 131 Ma, Rietpoort Granite and 133 Ma, Zout Rivier diorite) correspond to the age of the first magnetic seafloor-spreading anomaly offshore, and we suggest that the longevity of Koegel Fontein magmatism relates to a superposition of pre-drift magmatism onshore and spreading-related magmatism as continental separation began.
The occurrence of deep low-frequency (DLF) microearthquakes beneath volcanoes is commonly attributed to mass transport in the volcanic plumbing system and used to infer feeding channels from and into magma reservoirs. The key question is how magmas migrate from depth to the shallow crust and whether magma reservoirs are currently being recharged. For the first time since the improvement of the local seismic networks in the East Eifel region (Rhineland-Palatinate, Germany), we detect and locate recurrent DLF earthquakes in the lower crust and upper mantle beneath the Laacher See Volcano (LSV), using a joint data set of permanent sensors and a temporary deployment. So far, eight DLF earthquake sequences were observed in four distinct clusters between 10 and 40 km depth. These clusters of weak events (M-L< 2) align along an approximately 80. southeast dipping line south of the LSV. Moment tensor solutions of these events have large shear components, and the irregular dispersion and long coda of body waves indicate interaction processes between shear cracks and fluids. We find a rotation of P-axes orientation for shallow tectonic earthquakes compared to DLF events, indicating that the stress field in the depth interval of DLF events might favour a vertical migration of magma or magmatic fluids. The caldera of the LSV was formed by the last major eruption of the East Eifel Volcanic Field only 12.9 kyr ago, fed by a shallow magma chamber at 5-8 km depth and erupting a total magma volume of 6.7 km(3). The observed DLF earthquake activity and continuous volcanic gas emissions around the LSV indicate an active magmatic system, possibly connected with an upper mantle melt zone.
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.
Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (> 40 degrees N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.
The link between streamflow extremes and climatology has been widely studied in recent decades. However, a study investigating the effect of large-scale circulation variations on the distribution of seasonal discharge extremes at the European level is missing. Here we fit a climate-informed generalized extreme value (GEV) distribution to about 600 streamflow records in Europe for each of the standard seasons, i.e., to winter, spring, summer and autumn maxima, and compare it with the classical GEV distribution with parameters invariant in time. The study adopts a Bayesian framework and covers the period 1950 to 2016. Five indices with proven influence on the European climate are examined independently as covariates, namely the North Atlantic Oscillation (NAO), the east Atlantic pattern (EA), the east Atlantic-western Russian pattern (EA/WR), the Scandinavia pattern (SCA) and the polar-Eurasian pattern (POL). It is found that for a high percentage of stations the climate-informed model is preferred to the classical model. Particularly for NAO during winter, a strong influence on streamflow extremes is detected for large parts of Europe (preferred to the classical GEV distribution for 46% of the stations). Climate-informed fits are characterized by spatial coherence and form patterns that resemble relations between the climate indices and seasonal precipitation, suggesting a prominent role of the considered circulation modes for flood generation. For certain regions, such as northwestern Scandinavia and the British Isles, yearly variations of the mean seasonal climate indices result in considerably different extreme value distributions and thus in highly different flood estimates for individual years that can also persist for longer time periods.
We analyze a large transient strainmeter signal recorded at 62.5 m depth along the southern shore of the eastern Sea of Marmara region in northwestern Turkey. This region represents a passage of stress transfer from the Izmit rupture to the Marmara seismic gap. The strain signal was recorded at the Esenkoy site by one of the ICDP-GONAF (International Continental Drilling Programme - Geophysical Observatory at the North Anatolian Fault) strainmeters on the Armutlu peninsula with a maximum amplitude of 5 microstrain and lasting about 50 days. The onset of the strain signal coincided with the origin time of a M-w 4.4 earthquake offshore Yalova, which occurred as part of a seismic sequence including eight M-w >= 3.5 earthquakes. The Mw 4.4 event occurred at a distance of about 30 km from Esenkoy on June 25th 2016 representing the largest earthquake in this region since 2008. Before the event, the maximum horizontal strain was subparallel to the regional maximum horizontal stress derived from stress inversion of local seismicity. During the strain transient, we observe a clockwise rotation in the local horizontal strain field of about 20 degrees. The strain signal does not correlate with known environmental parameters such as annual changes of sea level, rainfall or temperature. The strain signal could indicate local slow slip on the Cinarcik fault and thus a transfer of stress to the eastern Marmara seismic gap.
The establishment and evolution of the Asian monsoons and arid interior have been linked to uplift of the Tibetan Plateau, retreat of the inland proto-Paratethys Sea and global cooling during the Cenozoic. However, the respective role of these driving mechanisms remains poorly constrained. This is partly due to a lack of continental records covering the key Eocene epoch marked by the onset of Tibetan Plateau uplift, proto-Paratethys Sea incursions and long-term global cooling. In this study, we reconstruct paleoenvironments in the Xining Basin, NE Tibet, to show a long-term drying of the Asian continental interior from the early Eocene to the Oligocene. Superimposed on this trend are three alternations between arid mudflat and wetter saline lake intervals, which are interpreted to reflect atmospheric moisture fluctuations in the basin. We date these fluctuations using magnetostratigraphy and the radiometric age of an intercalated tuff layer. The first saline lake interval is tentatively constrained to the late Paleocene-early Eocene. The other two are firmly dated between similar to 46 Ma (top magnetochron C21n) and similar to 41 Ma (base C18r) and between similar to 40 Ma (base C18n) and similar to 37 Ma (top C17n). Remarkably, these phases correlate in time with highstands of the proto-Paratethys Sea. This strongly suggests that these sea incursions enhanced westerly moisture supply as far inland as the Xining Basin. We conclude that the proto-Paratethys Sea constituted a key driver of Asian climate and should be considered in model and proxy interpretations. (C) 2019 Elsevier B.V. All rights reserved.
Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.
A strong temperature increase in the Arctic is expected to lead to latitudinal treeline shift. This tundra-taiga turnover would cause a positive vegetation-climate feedback due to albedo decrease. However, reliable estimates of tree migration rates are currently lacking due to the complex processes involved in forest establishment, which depend strongly on seed dispersal. We aim to fill this gap using LAVESI, an individual-based and spatially explicit Larix vegetation simulator. LAVESI was designed to simulate plots within homogeneous forests. Here, we improve the implementation of the seed dispersal function via field-based investigations. We inferred the effective seed dispersal distances of a typical open-forest stand on the southern Taymyr Peninsula (northern central Siberia) from genetic parentage analysis using eight nuclear microsatellite markers. The parentage analysis gives effective seed dispersal distances (median similar to 10 m) close to the seed parents. A comparison between simulated and observed effective seed dispersal distances reveals an overestimation of recruits close to the releasing tree and a shorter dispersal distance generally. We thus adapted our model and used the newly parameterised version to simulate south-to-north transects; a slow-moving treeline front was revealed. The colonisation of the tundra areas was assisted by occasional long-distance seed dispersal events beyond the treeline area. The treeline (similar to 1 tree ha(-1)) advanced by similar to 1.6 m yr(-1), whereas the forest line (similar to 100 trees ha(-1)) advanced by only similar to 0.6 m yr(-1). We conclude that the treeline in northern central Siberia currently lags behind the current strong warming and will continue to lag in the near future.
Measuring (dis)similarity between ecosystem states is a key theme in ecology. Much of community and ecosystem ecology is devoted to searching for patterns in ecosystem similarity from an external observer's viewpoint, using variables such as species abundances, measures of diversity and complexity. However, from the point of view of organisms in the ecosystem, proportional population growth rates are the only relevant aspect of ecosystem state, because natural selection acts on groups of organisms with different proportional population growth rates. We therefore argue that two ecosystem states are equivalent if and only if, for each species they contain, the proportional population growth rate does not differ between the states. Based on this result, we develop species-level and aggregated summary measures of ecosystem state and discuss their ecological meaning. We illustrate our approach using a long-term dataset on the plankton community from the Central European Lake Constance. We show that the first three principal components of proportional population growth rates describe most of the variation in ecosystem state in Lake Constance. We strongly recommend using proportional population growth rates and the derived equivalence classes for comparative ecosystem studies. This opens up new perspectives on important existing topics such as alternative stable ecosystem states, community assembly, and the processes generating regularities in ecosystems.
Proxy records from climate archives provide evidence about past climate changes, but the recorded signal is affected by non-climate-related effects as well as time uncertainty. As proxy-based climate reconstructions are frequently used to test climate models and to quantitatively infer past climate, we need to improve our understanding of the proxy record signal content as well as the uncertainties involved. In this study, we empirically estimate signal-to-noise ratios (SNRs) of temperature proxy records used in global compilations of the middle to late Holocene (last 6000 years). This is achieved through a comparison of the correlation of proxy time series from nearby sites of three compilations and model time series extracted at the proxy sites from two transient climate model simulations: a Holocene simulation of the ECHAM5/MPI-OM model and the Holocene part of the TraCE-21ka simulation. In all comparisons, we found the mean correlations of the proxy time series on centennial to millennial timescales to be low (R < 0.2), even for nearby sites, which resulted in low SNR estimates. The estimated SNRs depend on the assumed time uncertainty of the proxy records, the timescale analysed, and the model simulation used. Using the spatial correlation structure of the ECHAM5/MPI-OM simulation, the estimated SNRs on centennial timescales ranged from 0.05 - assuming no time uncertainty - to 0.5 for a time uncertainty of 400 years. On millennial timescales, the estimated SNRs were generally higher. Use of the TraCE-21ka correlation structure generally resulted in lower SNR estimates than for ECHAM5/MPI-OM. As the number of available high-resolution proxy records continues to grow, a more detailed analysis of the signal content of specific proxy types should become feasible in the near future. The estimated low signal content of Holocene temperature compilations should caution against over-interpretation of these multi-proxy and multisite syntheses until further studies are able to facilitate a better characterisation of the signal content in paleoclimate records.
Southern Patagonia is a prime example of ongoing oceanic ridge collision and slab-window formation sustained over several million years. The impact of these phenomena on the thermal structure and exhumation of the crust have been mainly assessed with low-temperature thermochronology of bedrock samples. Here, we infer thermal histories from new and existing thermochronological data from the region of most recent ridge collision. In particular, we evaluate the potential far-reaching thermal effects of the evolving slab window, which have previously been considered responsible for patterns of late Miocene reheating associated with back-arc alkaline volcanism. Our model results define protracted cooling since similar to 15 Ma and stepwise exhumation since the late Miocene. The pattern of stepwise exhumation closely matches the onset of Patagonian glaciation at 7 Ma and the successive pulse of glacial incision coeval with neotectonic activity since 3-4 Ma that are also documented by independent geological and geomorphological evidence in the region. Importantly, our findings challenge the recently suggested lack of glacial erosion and incision since 5 Ma in this region. Furthermore, in contrast to previous modelling studies, we find that the available data do not evidence a previously proposed northward-propagating heating event associated with alkaline volcanism. We hypothesize that the anomalous alkaline volcanism in the Patagonian back-arc might be related to trench-orthogonal tears aligned with transform faults in the subducting plate. The substantial differences from the previous modelling procedure on some of the same samples is demonstrated to result from an important lack of convergence in model runs. (C) 2019 Elsevier B.V. All rights reserved.
Stochastic modeling is a common practice for modeling uncertainty in hydrogeology. In stochastic modeling, aquifer properties are characterized by their probability density functions (PDFs). The Bayesian approach for inverse modeling is often used to assimilate information from field measurements collected at a site into properties’ posterior PDFs. This necessitates the definition of a prior PDF, characterizing the knowledge of hydrological properties before undertaking any investigation at the site, and usually coming from previous studies at similar sites. In this paper, we introduce a Bayesian hierarchical algorithm capable of assimilating various information–like point measurements, bounds and moments–into a single, informative PDF that we call ex-situ prior. This informative PDF summarizes the ex-situ information available about a hydrogeological parameter at a site of interest, which can then be used as a prior PDF in future studies at the site. We demonstrate the behavior of the algorithm on several synthetic case studies, compare it to other methods described in the literature, and illustrate the approach by applying it to a public open-access hydrogeological dataset.
The intermontane Humahuaca Basin in the Eastern Cordillera of the northwest Argentine Andes lies leeward of an orographic barrier to easterly derived moisture. An average of >2000 mm/yr of rainfall along the eastern flanks of the barrier contrasts with <200 mm/yr in the orogen interior. Paleoenvironmental reconstructions suggest that the basin became disconnected from the foreland during the Miocene-Pliocene by the growth of fault-bounded mountain ranges. Fossil records, sedimentology, and stable isotope data imply that rerouting of the fluvial network by 4.2 Ma and reduced rainfall by ca. 3 Ma were consequences of that range uplift. Here, we present cosmogenic nuclide-derived (Be-10) paleodenudation rates from 6 to 2 Ma fluvial deposits collected from the Humahuaca Basin. Despite increased tectonic activity, our Be-10 data show a tenfold decrease in denudation rates at ca. 3 Ma, documenting a link between uplift-induced semiarid conditions and decreasing hillslope denudation rates. This new data set thus demonstrates the influence of hydrological change on spatiotemporal denudation patterns in tectonically active mountain areas.
The application of electrical resistivity tomography to peatlands supports conventional coring by providing data on the current condition of peatlands, including data on stratigraphy, peat properties and thickness of organic deposits. Data on the current condition of drained peatlands are particularly required to improve estimates of carbon storage as well as losses and emissions from agriculturally used peatlands. However, most of the studies focusing on electrical resistivity tomography surveys have been conducted on natural peatlands with higher groundwater levels. Peatlands drained for agriculture have not often been studied using geophysical techniques. Drained sites are characterized by low groundwater levels and high groundwater fluctuations during the year, which lead to varying levels of water saturation. To validate better electrical resistivity tomography surveys of drained peatlands, the aim of this laboratory study is to investigate the influence of varying water saturation levels on electrical conductivity (reciprocal of resistivity) for a variety of peat and gyttja types, as well as for different degrees of peat decomposition. Results show that different levels of water saturation strongly influence bulk electrical conductivity. Distinct differences in this relationship exist between peat and gyttja substrates and between different degrees of peat decomposition. Peat shows an exponential relationship for all degrees of decomposition, whereas gyttja, in particular organic-rich gyttja, is characterized by a rather unimodal relationship. The slopes for the relationship between electrical conductivity and water content are steeper at high degrees of decomposition than for peat of low degrees of decomposition. These results have direct implications for field electrical resistivity tomography surveys. In drained peatlands that are strongly susceptible to drying, electrical resistivity tomography surveys have a high potential to monitor the actual field water content. In addition, at comparable water saturations, high or low degrees of decomposition can be inferred from electrical conductivity.
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.
Soil degradation by water is a serious environmental problem worldwide, with specific climatic factors being the major causes. We investigated the relationships between synoptic atmospheric patterns (i.e. weather types, WTs) and runoff, erosion and sediment yield throughout the Mediterranean basin by analyzing a large database of natural rainfall events at 68 research sites in 9 countries. Principal Component Analysis (PCA) was used to identify spatial relationships of the different WTs including three hydro-sedimentary variables: rainfall, runoff, and sediment yield (SY, used to refer to both soil erosion measured at plot scale and sediment yield registered at catchment scale). The results indicated 4 spatial classes of rainfall and runoff: (a) northern sites dependent on North (N) and North West (NW) flows; (b) eastern sites dependent on E and NE flows; (c) southern sites dependent on S and SE flows; and, finally, (d) western sites dependent on W and SW flows. Conversely, three spatial classes are identified for SY characterized by: (a) N and NE flows in northern sites (b) E flows in eastern sites, and (c) W and SW flows in western sites. Most of the rainfall, runoff and SY occurred during a small number of daily events, and just a few WTs accounted for large percentages of the total. Our results confirm that characterization by WT improves understanding of the general conditions under which runoff and SY occur, and provides useful information for understanding the spatial variability of runoff, and SY throughout the Mediterranean basin. The approach used here could be useful to aid of the design of regional water management and soil conservation measures.
Water stable isotope signatures can provide valuable insights into the catchment internal runoff processes. However, the ability of the water isotope data to constrain the internal apportionments of runoff components in hydrological models for glacierized basins is not well understood. This study developed an approach to simultaneously model the water stable isotopic compositions and runoff processes in a glacierized basin in Central Asia. The fractionation and mixing processes of water stable isotopes in and from the various water sources were integrated into a glacio-hydrological model. The model parameters were calibrated on discharge, snow cover and glacier mass balance data, and additionally isotopic composition of streamflow. We investigated the value of water isotopic compositions for the calibration of model parameters, in comparison to calibration methods without using such measurements. Results indicate that: (1) The proposed isotope-hydrological integrated modeling approach was able to reproduce the isotopic composition of streamflow, and improved the model performance in the evaluation period; (2) Involving water isotopic composition for model calibration reduced the model parameter uncertainty, and helped to reduce the uncertainty in the quantification of runoff components; (3) The isotope-hydrological integrated modeling approach quantified the contributions of runoff components comparably to a three-component tracer-based end-member mixing analysis method for summer peak flows, and required less water tracer data. Our findings demonstrate the value of water isotopic compositions to improve the quantification of runoff components using hydrological models in glacierized basins.
Speleothems, secondary cave carbonates, are important tools for climate reconstruction, especially as they often can be very precisely dated with the UTh method. If the U-Th method fails, dating becomes difficult, and often results in abandonment of samples and study sites. Radiocarbon dating is the only other radiometric dating technique applicable to the last similar to 50 ka, but presents complexities related to temporal variability of the reservoir effect in speleothems. Thus, radiocarbon dating of speleothems is not straightforward, and there are currently no publicly available tools to define proper age-depth relationships with this method. Here, we present an improved version of a previously published radiocarbon based age-depth modelling approach (star, Lechleitner et al., 2016b), which is now made freely available. The software is easy to use and provides the possibility to obtain reliable age-depth relationships, without prior knowledge of reservoir effects and their variability. In addition, star is able to detect and handle growth stops and phases with different growth rates. We test star on artificially constructed data sets and illustrate steps to improve the model performance. Furthermore, we apply the new approach to published radiocarbon data of U-Th dated stalagmites. This offers the possibility to investigate the strengths and weaknesses of the new approach with respect to potentially significant long term trends in the radiocarbon reservoir effect, which might otherwise remain undetected. In summary, we have produced a valuable software, which easily enables to construct age-depth relationships on the basis of reservoir effect disturbed radiocarbon measurements.
Models of ring current electron dynamics unavoidably contain uncertainties in boundary conditions, electric and magnetic fields, electron scattering rates, and plasmapause location. Model errors can accumulate with time and result in significant deviations of model predictions from observations. Data assimilation offers useful tools which can combine physics-based models and measurements to improve model predictions. In this study, we systematically analyze performance of the Kalman filter applied to a log-transformed convection model of ring current electrons and Van Allen Probe data. We consider long-term dynamics of mu = 2.3 MeV/G and K = 0.3 G(1/2) R-E electrons from 1 February 2013 to 16 June 2013. By using synthetic data, we show that the Kalman filter is capable of correcting errors in model predictions associated with uncertainties in electron lifetimes, boundary conditions, and convection electric fields. We demonstrate that reanalysis retains features which cannot be fully reproduced by the convection model such as storm-time earthward propagation of the electrons down to 2.5 R-E. The Kalman filter can adjust model predictions to satellite measurements even in regions where data are not available. We show that the Kalman filter can adjust model predictions in accordance with observations for mu = 0.1, 2.3, and 9.9 MeV/G and constant K = 0.3 G(1/2) R-E electrons. The results of this study demonstrate that data assimilation can improve performance of ring current models, better quantify model uncertainties, and help deeper understand the physics of the ring current particles.
Thawing of subsea permafrost can impact offshore infrastructure, affect coastal erosion, and release permafrost organic matter. Thawing is usually modeled as the result of heat transfer, although salt diffusion may play an important role in marine settings. To better quantify nearshore subsea permafrost thawing, we applied the CryoGRID2 heat diffusion model and coupled it to a salt diffusion model. We simulated coastline retreat and subsea permafrost evolution as it develops through successive stages of a thawing sequence at the Bykovsky Peninsula, Siberia. Sensitivity analyses for seawater salinity were performed to compare the results for the Bykovsky Peninsula with those of typical Arctic seawater. For the Bykovsky Peninsula, the modeled ice-bearing permafrost table (IBPT) for ice-rich sand and an erosion rate of 0.25m/year was 16.7 m below the seabed 350m offshore. The model outputs were compared to the IBPT depth estimated from coastline retreat and electrical resistivity surveys perpendicular to and crossing the shoreline of the Bykovsky Peninsula. The interpreted geoelectric data suggest that the IBPT dipped to 15-20m below the seabed at 350m offshore. Both results suggest that cold saline water forms beneath grounded ice and floating sea ice in shallow water, causing cryotic benthic temperatures. The freezing point depression produced by salt diffusion can delay or prevent ice formation in the sediment and enhance the IBPT degradation rate. Therefore, salt diffusion may facilitate the release of greenhouse gasses to the atmosphere and considerably affect the design of offshore and coastal infrastructure in subsea permafrost areas.
An essential, respected, and critical aspect of the modern practice of science and scientific publishing is peer review. The process of peer review facilitates best practices in scientific conduct and communication, ensuring that manuscripts published are as accurate, valuable, and clearly communicated. The over 216 papers published in Tectonics in 2018 benefit from the time, effort, and expertise of our reviewers who have provided thoughtfully considered advice on each manuscript. This role is critical to advancing our understanding of the evolution of the continents and their margins, as these reviews lead to even clearer and higher-quality papers. In 2018, the over 443 papers submitted to Tectonics were the beneficiaries of more than 1,010 reviews provided by 668 members of the tectonics community and related disciplines. To everyone who has volunteered their time and intellect to peer reviewing, thank you for helping Tectonics and all other AGU Publications provide the best science possible.
Pluvial flood risk is mostly excluded in urban flood risk assessment. However, the risk of pluvial flooding is a growing challenge with a projected increase of extreme rainstorms compounding with an ongoing global urbanization. Considered as a flood type with minimal impacts when rainfall rates exceed the capacity of urban drainage systems, the aftermath of rainfall-triggered flooding during Hurricane Harvey and other events show the urgent need to assess the risk of pluvial flooding. Due to the local extent and small-scale variations, the quantification of pluvial flood risk requires risk assessments on high spatial resolutions. While flood hazard and exposure information is becoming increasingly accurate, the estimation of losses is still a poorly understood component of pluvial flood risk quantification. We use a new probabilistic multivariable modeling approach to estimate pluvial flood losses of individual buildings, explicitly accounting for the associated uncertainties. Except for the water depth as the common most important predictor, we identified the drivers for having loss or not and for the degree of loss to be different. Applying this approach to estimate and validate building structure losses during Hurricane Harvey using a property level data set, we find that the reliability and dispersion of predictive loss distributions vary widely depending on the model and aggregation level of property level loss estimates. Our results show that the use of multivariable zero-inflated beta models reduce the 90% prediction intervalsfor Hurricane Harvey building structure loss estimates on average by 78% (totalling U.S.$3.8 billion) compared to commonly used models.
Lineament mapping, which is an important part of any structural geological investigation, is made more efficient and easier by the availability of optical as well as radar remote sensing data, such as Landsat and Sentinel with medium and high spatial resolutions. However, the results from these multi-resolution data vary due to their difference in spatial resolution and sensitivity to soil occupation. The accuracy and quality of extracted lineaments depend strongly on the spatial resolution of the imagery. Therefore, the aim of this study was to compare the optical Landsat-8, Sentinel-2A, and radar Sentinel-1A satellite data for automatic lineament extraction. The framework of automatic approach includes defining the optimal parameters for automatic lineament extraction with a combination of edge detection and line-linking algorithms and determining suitable bands from optical data suited for lineament mapping in the study area. For the result validation, the extracted lineaments are compared against the manually obtained lineaments through the application of directional filtering and edge enhancement as well as to the lineaments digitized from the existing geological maps of the study area. In addition, a digital elevation model (DEM) has been utilized for an accuracy assessment followed by the field verification. The obtained results show that the best correlation between automatically extracted lineaments, manual interpretation, and the preexisting lineament map is achieved from the radar Sentinel-1A images. The tests indicate that the radar data used in this study, with 5872 and 5865 lineaments extracted from VH and VV polarizations respectively, is more efficient for structural lineament mapping than the Landsat-8 and Sentinel-2A optical imagery, from which 2338 and 4745 lineaments were extracted respectively.
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.
We present a model of the electrical resistivity structure of the lithosphere in the Central Andes between 20 degrees and 24 degrees S from 3-D inversion of 56 long-period magnetotelluric sites. Our model shows a complex resistivity structure with significant variability parallel and perpendicular to the trench direction. The continental forearc is characterized mainly by high electrical resistivity (>1,000m), suggesting overall low volumes of fluids. However, low resistivity zones (LRZs, <5m) were found in the continental forearc below areas where major trench-parallel faults systems intersect NW-SE transverse faults. Forearc LRZs indicate circulation and accumulation of fluids in highly permeable fault zones. The continental crust along the arc shows three distinctive resistivity domains, which coincide with segmentation in the distribution of volcanoes. The northern domain (20 degrees-20.5 degrees S) is characterized by resistivities >1,000m and the absence of active volcanism, suggesting the presence of a low-permeability block in the continental crust. The central domain (20.5 degrees-23 degrees S) exhibits a number of LRZs at varying depths, indicating different levels of a magmatic plumbing system. The southern domain (23 degrees-24 degrees S) is characterized by resistivities >1,000m, suggesting the absence of large magma reservoirs below the volcanic chain at crustal depths. Magma reservoirs located below the base of the crust or in the backarc may fed active volcanism in the southern domain. In the subcontinental mantle, the model exhibits LRZs in the forearc mantle wedge and above clusters of intermediate-depth seismicity, likely related to fluids produced by serpentinization of the mantle and eclogitization of the slab, respectively.
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.
Study region: Tisza and Prut catchments, originating on the slopes of the Carpathian mountains. Study focus: The study reported here investigates (i) climate change impacts on flood risk in the region, and (ii) uncertainty related to hydrological modelling, downscaling techniques and climate projections. The climate projections used in the study were derived from five GCMs, downscaled either dynamically with RCMs or with the statistical downscaling model XDS. The resulting climate change scenarios were applied to drive the eco-hydrological model SWIM, which was calibrated and validated for the catchments in advance using observed climate and hydrological data. The changes in the 30-year flood hazards and 98 and 95 percentiles of discharge were evaluated for the far future period (2071-2100) in comparison with the reference period (1981-2010). New hydrological insights for the region: The majority of model outputs under RCP 4.5 show a small to strong increase of the 30-year flood level in the Tisza ranging from 4.5% to 62%, and moderate increase in the Prut ranging from 11% to 22%. The impact results under RCP 8.5 are more uncertain with changes in both directions due to high uncertainties in GCM-RCM climate projections, downscaling methods and the low density of available climate stations.
Design flood estimation is an essential part of flood risk assessment. Commonly applied are flood frequency analyses and design storm approaches, while the derived flood frequency using continuous simulation has been getting more attention recently. In this study, a continuous hydrological modelling approach on an hourly time scale, driven by a multi-site weather generator in combination with a -nearest neighbour resampling procedure, based on the method of fragments, is applied. The derived 100-year flood estimates in 16 catchments in Vorarlberg (Austria) are compared to (a) the flood frequency analysis based on observed discharges, and (b) a design storm approach. Besides the peak flows, the corresponding runoff volumes are analysed. The spatial dependence structure of the synthetically generated flood peaks is validated against observations. It can be demonstrated that the continuous modelling approach can achieve plausible results and shows a large variability in runoff volume across the flood events.
The propagation of a seismic rupture on a fault introduces spatial variations in the seismic wave field surrounding the fault. This directivity effect results in larger shaking amplitudes in the rupture propagation direction. Its seismic radiation pattern also causes amplitude variations between the strike-normal and strike-parallel components of horizontal ground motion. We investigated the landslide response to these effects during the 2016 Kumamoto earthquake (M-w 7.1) in central Kyushu (Japan). Although the distribution of some 1500 earthquake-triggered landslides as a function of rupture distance is consistent with the observed Arias intensity, the landslides were more concentrated to the northeast of the southwest-northeast striking rupture. We examined several landslide susceptibility factors: hillslope inclination, the median amplification factor (MAF) of ground shaking, lithology, land cover, and topographic wetness. None of these factors sufficiently explains the landslide distribution or orientation (aspect), although the landslide head scarps have an elevated hillslope inclination and MAF. We propose a new physics-based ground-motion model (GMM) that accounts for the seismic rupture effects, and we demonstrate that the low-frequency seismic radiation pattern is consistent with the overall landslide distribution. Its spatial pattern is influenced by the rupture directivity effect, whereas landslide aspect is influenced by amplitude variations between the fault-normal and fault-parallel motion at frequencies < 2 Hz. This azimuth dependence implies that comparable landslide concentrations can occur at different distances from the rupture. This quantitative link between the prevalent landslide aspect and the low-frequency seismic radiation pattern can improve coseismic landslide hazard assessment.
Core Ideas
3D MRI relaxation time maps reflect water mobility in root, rhizosphere, and soil.
3D NCT water content maps of the same plant complement relaxation time maps.
The relaxation time T1 decreases from soil to root, whereas water content increases.
Parameters together indicate modification of rhizosphere pore space by gel phase.
The zone of reduced T1 corresponds to the zone remaining dry after rewetting.
In situ investigations of the rhizosphere require high‐resolution imaging techniques, which allow a look into the optically opaque soil compartment. We present the novel combination of magnetic resonance imaging (MRI) and neutron computed tomography (NCT) to achieve synergistic information such as water mobility in terms of three‐dimensional (3D) relaxation time maps and total water content maps. Besides a stationary MRI scanner for relaxation time mapping, we used a transportable MRI system on site in the NCT facility to capture rhizosphere properties before desiccation and after subsequent rewetting. First, we addressed two questions using water‐filled test capillaries between 0.1 and 5 mm: which root diameters can still be detected by both methods, and to what extent are defined interfaces blurred by these imaging techniques? Going to real root system architecture, we demonstrated the sensitivity of the transportable MRI device by co‐registration with NCT and additional validation using X‐ray computed tomography. Under saturated conditions, we observed for the rhizosphere in situ a zone with shorter T1 relaxation time across a distance of about 1 mm that was not caused by reduced water content, as proven by successive NCT measurements. We conclude that the effective pore size in the pore network had changed, induced by a gel phase. After rewetting, NCT images showed a dry zone persisting while the MRI intensity inside the root increased considerably, indicating water uptake from the surrounding bulk soil through the still hydrophobic rhizosphere. Overall, combining NCT and MRI allows a more detailed analysis of the rhizosphere's functioning.
Human-driven fragmentation of landscapes leads to the formation of transition zones between ecosystems that are characterised by fluxes of matter, energy and information. These transition zones may offer rather inhospitable habitats that could jeopardise biodiversity. On the other hand, transition zones are also reported to be hotspots for biodiversity and even evolutionary processes. The general mechanisms and influence of processes in transition zones are poorly understood. Although heterogeneity and diversity of land use of fragments and the transition zones between them play an important role, most studies only refer to forested transition zones. Often, only an extrapolation of measurements in the different fragments themselves is reported to determine gradients in transition zones. This paper contributes to a quantitative understanding of agricultural landscapes beyond individual ecotopes, and towards connected ecosystem mosaics that may be beneficial for the provision of ecosystem services.
Spinning up large-scale coupled surface-subsurface numerical models can be a time and resource consuming task. If an uninformed initial condition is chosen, the spin-up can easily require 20 years of repeated simulations on high-performance computing machines. In this paper we compare the classical approach of starting from a fixed shallow depth to groundwater (here 3 m) with three more informed approaches for the definition of initial conditions in the spin up. In the first of these three approaches, we start from a known-steady state groundwater table, calculated with a 2-D groundwater model and the yearly net recharge, and combine it with an unsaturated zone that assumes hydrostatic conditions. In the second approach, we start from the same groundwater table combined with vertical profiles in the unsaturated zone with uniform vertical flow identical to the groundwater recharge. In the third approach we calculate a dynamic steady state from a simplified subsurface model combining a transient 2-D groundwater model with a limited number of 1-D transient unsaturated zone columns on top. Results for spinning-up a 3-D Parflow-CLM model using the different initial conditions show that large gains can be made by considering states in groundwater and the vadose zone that are consistent, i.e. where groundwater recharge and the vertical flux in the vadose zone agree. By this, the spin-up time was reduced from about 10 years to about 3 years of simulated time. In the light of seasonal fluctuations of net recharge, using the transient approach showed more stable results.
As a consequence of the rapid growing worldwide seismic data set, a huge variety of automatized data-processing methods have been developed. To perform automatized waveform-based seismological studies aiming for magnitudes or source process inversion, it is crucial to identify network stations with erroneous transfer functions, gain factors, or component orientations. We developed a new tool dedicated to automated station quality control of dense seismic networks and arrays. The python-based AutoStatsQ toolbox uses the pyrocko seismic data-processing environment. The toolbox automatically downloads data and metadata for selected teleseismic events and performs different tests. As a result, relative gain factors, sensor orientation corrections, and reliable frequency bands are computed for all stations in a chosen time period. Relative gain factors are calculated for all stations and events in a time domain based on maximum P-phase amplitudes. A Rayleigh-wave polarization analysis is used to identify deviating sensor orientations. The power spectra of all stations in a given frequency range are compared with synthetic ones, accessing Global Centroid Moment Tensor (CMT) solutions. Frequency ranges of coinciding synthetic and recorded power spectral densities (PSDs) may serve as guidelines for choosing band-pass filters for moment tensor (MT) inversion and help confirm the corner frequency of the instrument. The toolbox was applied to the permanent and temporary AlpArray networks as well as to the denser SWATH-D network, a total of over 750 stations. Stations with significantly deviating gain factors were identified, as well as stations with inverse polarity and misorientations of the horizontal components. The tool can be used to quickly access network quality and to omit or correct stations before MT inversion. Electronic Supplement: List of teleseismic events and tables of median, mean, and standard deviation of relative gain factors, and figures of relative gain factors of all event-station pairs, waveform example showing inverse polarity of horizontal components on ZS.D125, histograms of median, mean, and standard deviation of the correction angles, examples of synthetic and recorded frequency spectra of ZS.D046 and NI.VINO.