Refine
Year of publication
Document Type
- Article (2818) (remove)
Keywords
- Holocene (35)
- climate change (30)
- Earthquake source observations (21)
- Climate change (18)
- Pollen (18)
- erosion (18)
- permafrost (18)
- Himalaya (17)
- Tibetan Plateau (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2818) (remove)
The region of the Andean back-arc of northwestern Argentina has been struck by several magnitude >= 6 crustal earthquakes since the first historically recorded event in 1692. One of these events corresponds to the Anta earthquake on 25 August 1948, with epicenter in the Santa Barbara System causing three deaths and severe damage in Salta and Jujuy provinces with maximum Modified Mercalli seismic intensities (MMI) of IX. We collected and digitized analog seismograms of this earthquake from worldwide seismic observatories in order to perform first-motion analysis and modeling of long-period teleseismic P-waveforms. Our results indicate a simple seismic source of M0 = 2.85 x 1019 N m consistent with a moment magnitude Mw = 6.9. We have also tested for the focal depth determining a shallow source at 8 km with a reverse focal mechanism solution with a minor dextral strike-slip component (strike 20 degrees, dip 30 degrees, rake 120 degrees) from the best fit of waveforms. Using magnitude size empirical relationships, the comparison of the obtained Mw 6.9 magnitude value and the ca. 10,000 km2 area of MMI >= IX from our seismic intensity map, which was obtained from newspaper and many historical reports, indicates a rupture length of 42 +/- 8 km for the Anta earthquake. We show our results in a 3D geological model around the epicentral area, which integrates modern seismicity, geological data, and information of a previously studied east-west cross section located a few kilometers south of the 1948 epicenter. The integration of all available information provides evidence of the re-activation of the Pie de la Sierra del Gallo fault during the 1948 Mw 6.9 shallow earthquake; this thrust fault bounds the Santa Barbara System along its western foothill.
Since the beginning of the Anthropocene, lacustrine biodiversity has been influenced by climate change and human activities. These factors advance the spread of harmful cyanobacteria in lakes around the world, which affects water quality and impairs the aquatic food chain. In this study, we assessed changes in cyanobacterial community dynamics via sedimentary DNA (sedaDNA) from well-dated lake sediments of Lake Tiefer See, which is part of the Klocksin Lake Chain spanning the last 350 years. Our diversity and community analysis revealed that cyanobacterial communities form clusters according to the presence or absence of varves. Based on distance-based redundancy and variation partitioning analyses (dbRDA and VPA) we identified that intensified lake circulation inferred from vegetation openness reconstructions, delta C-13 data (a proxy for varve preservation) and total nitrogen content were abiotic factors that significantly explained the variation in the reconstructed cyanobacterial community from Lake Tiefer See sediments. Operational taxonomic units (OTUs) assigned to Microcystis sp. and Aphanizomenon sp. were identified as potential eutrophication-driven taxa of growing importance since circa common era (ca. CE) 1920 till present. This result is corroborated by a cyanobacteria lipid biomarker analysis. Furthermore, we suggest that stronger lake circulation as indicated by non-varved sediments favoured the deposition of the non-photosynthetic cyanobacteria sister clade Sericytochromatia, whereas lake bottom anoxia as indicated by subrecent- and recent varves favoured the Melainabacteria in sediments. Our findings highlight the potential of high-resolution amplicon sequencing in investigating the dynamics of past cyanobacterial communities in lake sediments and show that lake circulation, anoxic conditions, and human-induced eutrophication are main factors explaining variations in the cyanobacteria community in Lake Tiefer See during the last 350 years.
In-depth understanding of the reorganization of the hydrological cycle in response to global climate change is crucial in highly sensitive regions like the eastern Mediterranean, where water availability is a major factor for socioeconomic and political development.
The sediments of Lake Lisan provide a unique record of hydroclimatic change during the last glacial to Holocene transition (ca. 24-11 ka) with its tremendous water level drop of similar to 240 m that finally led to its transition into the present hypersaline water body-the Dead Sea.
Here we utilize high-resolution sedimentological analyses from the marginal terraces and deep lake to reconstruct an unprecedented seasonal record of the last millennia of Lake Lisan. Aragonite varve formation in intercalated intervals of our record demonstrates that a stepwise long-term lake level decline was interrupted by almost one millennium of rising or stable water level.
Even periods of pronounced water level drops indicated by gypsum deposition were interrupted by decades of positive water budgets.
Our results thus highlight that even during major climate change at the end of the last glacial, decadal to millennial periods of relatively stable or positive moisture supply occurred which could have been an important premise for human sedentism.
Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering.
However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood.
In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects.
We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes.
We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations.
Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect.
We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.
The Arctic is rich in aquatic systems and experiences rapid warming due to climate change. The accelerated warming causes permafrost thaw and the mobilization of organic carbon.
When dissolved organic carbon is mobilized, this DOC can be transported to aquatic systems and degraded in the water bodies and further downstream. Here, we analyze the influence of different landscape components on DOC concentrations and export in a small (6.45 km(2)) stream catchment in the Lena River Delta.
The catchment includes lakes and ponds, with the flow path from Pleistocene yedoma deposits across Holocene non-yedoma deposits to the river outlet. In addition to DOC concentrations, we use radiocarbon dating of DOC as well as stable oxygen and hydrogen isotopes (delta O-18 and delta D) to assess the origin of DOC.
We find significantly higher DOC concentrations in the Pleistocene yedoma area of the catchment compared to the Holocene non-yedoma area with medians of 5 and 4.5 mg L-1 (p < 0.05), respectively. When yedoma thaw streams with high DOC concentration reach a large yedoma thermokarst lake, we observe an abrupt decrease in DOC concentration, which we attribute to dilution and lake processes such as mineralization. The DOC ages in the large thermokarst lake (between 3,428 and 3,637 C-14 y BP) can be attributed to a mixing of mobilized old yedoma and Holocene carbon.
Further downstream after the large thermokarst lake, we find progressively younger DOC ages in the stream water to its mouth, paired with decreasing DOC concentrations. This process could result from dilution with leaching water from Holocene deposits and/or emission of ancient yedoma carbon to the atmosphere. Our study shows that thermokarst lakes and ponds may act as DOC filters, predominantly by diluting incoming waters of higher DOC concentrations or by re-mineralizing DOC to CO2 and CH4.
Nevertheless, our results also confirm that the small catchment still contributes DOC on the order of 1.2 kg km(-2) per day from a permafrost landscape with ice-rich yedoma deposits to the Lena River.
Mechanical behaviors of granite after thermal treatment under loading and unloading conditions
(2021)
Understanding the mechanical behaviors of granite after thermal treatment under loading and unloading conditions is of utmost relevance to deep geothermal energy recovery. In the present study, a series of loading and unloading triaxial compression tests (20, 40 and 60 MPa) on granite specimens after exposure to different temperatures (20, 200, 300, 400, 500 and 600 degrees C) was carried out to quantify the combined effects of thermal treatment and loading/unloading stress conditions on granite strength and deformation. Changes in the microstructure of granite exposed to high temperatures were revealed by optical microscopy. The experimental results indicate that both, thermal treatment and loading/unloading stress conditions, degrade the mechanical behaviors and further decrease the carrying capacity of granite. The gradual degradation of the mechanical characteristics of granite after thermal treatment is mainly associated with the evolution of thermal micro-cracks based on optical microscopy observations. The unloading stress state induces the extension of tension cracks parallel to the axial direction, and thus, the mechanical properties are degraded. Temperatures above 400 degrees C have a more significant influence on the mechanical characteristics of granite than the unloading treatment, whereby 400 degrees C can be treated as a threshold temperature for the delineation of significant deterioration. This study is expected to support feasibility and risk assessments by means of providing data for analytical calculations and numerical simulations on granite exposed to high temperatures during geothermal energy extraction.
It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics(1-4). However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process(5,6) and mantle-scale tectonic model(7). End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state(8), strength controlled (Bm approximate to 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.
In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).
Dryland xeric conditions exert a deterministic effect on microbial communities, forcing life into refuge niches. Deposited rocks can form a lithic niche for microorganisms in desert regions. Mineral weathering is a key process in soil formation and the importance of microbial-driven mineral weathering for nutrient extraction is increasingly accepted. Advances in geobiology provide insight into the interactions between microorganisms and minerals that play an important role in weathering processes. In this study, we present the examination of the microbial diversity in dryland rocks from the Tsauchab River banks in Namibia. We paired culture-independent 16S rRNA gene amplicon sequencing with culture-dependent (isolation of bacteria) techniques to assess the community structure and diversity patterns. Bacteria isolated from dryland rocks are typical of xeric environments and are described as being involved in rock weathering processes. For the first time, we extracted extra- and intracellular DNA from rocks to enhance our understanding of potentially rock-weathering microorganisms. We compared the microbial community structure in different rock types (limestone, quartz-rich sandstone and quartz-rich shale) with adjacent soils below the rocks. Our results indicate differences in the living lithic and sublithic microbial communities.
In some catchments, the distribution of annual maximum streamflow shows heavy tail behavior, meaning the occurrence probability of extreme events is higher than if the upper tail decayed exponentially. Neglecting heavy tail behavior can lead to an underestimation of the likelihood of extreme floods and the associated risk. Partly contradictory results regarding the controls of heavy tail behavior exist in the literature and the knowledge is still very dispersed and limited. To better understand the drivers, we analyze the upper tail behavior and its controls for 480 catchments in Germany and Austria over a period of more than 50 years. The catchments span from quickly reacting mountain catchments to large lowland catchments, allowing for general conclusions. We compile a wide range of event and catchment characteristics and investigate their association with an indicator of the tail heaviness of flood distributions, namely the shape parameter of the GEV distribution. Following univariate analyses of these characteristics, along with an evaluation of different aggregations of event characteristics, multiple linear regression models, as well as random forests, are constructed. A novel slope indicator, which represents the relation between the return period of flood peaks and event characteristics, captures the controls of heavy tails best. Variables describing the catchment response are found to dominate the heavy tail behavior, followed by event precipitation, flood seasonality, and catchment size. The pre-event moisture state in a catchment has no relevant impact on the tail heaviness even though it does influence flood magnitudes.
The metastable paragenesis of corundum and quartz is rare in nature but common in laboratory experiments where according to thermodynamic predictions aluminum-silicate polymorphs should form. We demonstrate here that the existence of a hydrous, silicon-bearing, nanometer-thick layer (called "HSNL") on the corundum surface can explain this metastability in experimental studies without invoking unspecific kinetic inhibition. We investigated experimentally formed corundum reaction products synthesized during hydrothermal and piston-cylinder experiments at 500-800 degrees C and 0.25-1.8 GPa and found that this HSNL formed inside and on the corundum crystals, thereby controlling the growth behavior of its host. The HSNL represents a substitution of Al with Si and H along the basal plane of corundum. Along the interface of corundum and quartz, the HSNL effectively isolates the bulk phases corundum and quartz from each other, thus apparently preventing their reaction to the stable aluminum silicate. High temperatures and prolonged experimental duration lead to recrystallization of corundum including the HSNL and to the formation of quartz + fluid inclusions inside the host crystal. This process reduces the phase boundary area between the bulk phases, thereby providing further opportunity to expand their coexistence. In addition to its small size, its transient nature makes it difficult to detect the HSNL in experiments and even more so in natural samples. Our findings emphasize the potential impact of nanometer-sized phases on geochemical reaction pathways and kinetics under metamorphic conditions in one of the most important chemical systems of the Earth's crust.
Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.
The release of greenhouse gases from the large organic carbon stock in permafrost deposits in the circumarctic regions may accelerate global warming upon thaw. The extent of this positive climate feedback is thought to be largely controlled by the microbial degradability of the organic matter preserved in these sediments. In addition, weathering and oxidation processes may release inorganic carbon preserved in permafrost sediments as CO2, which is generally not accounted for. We used C-13 and C-14 analysis and isotopic mass balances to differentiate and quantify organic and inorganic carbon released as CO2 in the field from an active retrogressive thaw slump of Pleistocene-age Yedoma and during a 1.5-years incubation experiment. The results reveal that the dominant source of the CO2 released from freshly thawed Yedoma exposed as thaw mound is Pleistocene-age organic matter (48-80%) and to a lesser extent modern organic substrate (3-34%). A significant portion of the CO2 originated from inorganic carbon in the Yedoma (17-26%). The mixing of young, active layer material with Yedoma at a site on the slump floor led to the preferential mineralization of this young organic carbon source. Admixtures of younger organic substrates in the Yedoma thaw mound were small and thus rapidly consumed as shown by lower contributions to the CO2 produced during few weeks of aerobic incubation at 4 degrees C corresponding to approximately one thaw season. Future CO2 fluxes from the freshly thawed Yedoma will contain higher proportions of ancient inorganic (22%) and organic carbon (61-78%) as suggested by the results at the end, after 1.5 years of incubation. The increasing contribution of inorganic carbon during the incubation is favored by the accumulation of organic acids from microbial organic matter degradation resulting in lower pH values and, in consequence, in inorganic carbon dissolution. Because part of the inorganic carbon pool is assumed to be of pedogenic origin, these emissions would ultimately not alter carbon budgets. The results of this study highlight the preferential degradation of younger organic substrates in freshly thawed Yedoma, if available, and a substantial release of CO2 from inorganic sources.
The 100 km wide Merida Andes extend from the Colombian/Venezuelan border to the Coastal Cordillera. The mountain chain and its associated major strike-slip fault systems in western Venezuela formed due to oblique convergence of the Caribbean with the South American Plates and the north-eastwards expulsion of the North Andean Block. Due to the limited knowledge of lithospheric structures related to the formation of the Merida Andes research projects have been developed to illuminate this zone with deep geophysical data. In this study, we present three-dimensional inversion of broadband magnetotelluric data, collected along a 240 km long profile crossing the Merida Andes and the Maracaibo and Barinas-Apure foreland basins. The distribution of the stations limits resolution of the model to off-profile features. Combining 3D inversion of synthetic data sets derived from 3D modelling with 3D inversion of measured data, we could derive a 10 to 15 km wide corridor with good lateral resolution to develop hypotheses about the origin of deep-reaching anomalies of high electrical conductivity. The Merida Andes appear generally as electrically resistive structures, separated by anomalies associated with the most important fault systems of the region, the Bocono and Valera faults. Sensitivity tests suggest that the Valera Fault reaches to depths of up to 12 km and the Bocono Fault to more than 35 km depth. Both structures are connected to a sizeable conductor located east of the profile at 12-15 km depth. We propose that the high conductivity associated with this off-profile conductor may be related to the detachment of the Trujillo Block. We also identified a conductive zone that correlates spatially with the location of a gravity low, possibly representing a SE tilt of the Maracaibo Triangular Block under the mountain chain to great depths (>30 km). The relevance of these tectonic blocks in our models at crustal depths seems to be consistent with proposed theories that describe the geodynamics of western Venezuela as dominated by floating blocks or orogens. Our results stress the importance of the Trujillo Block for the current tectonic evolution of western Venezuela and confirm the relevance of the Bocono Fault carrying deformation to the lower crust and upper mantle. The Barinas-Apure and the Maracaibo sedimentary basins are imaged as electrically conductive with depths of 4 to 5 km and 5 to 10 km, respectively. The Barinas-Apure basin is imaged as a simple 1D structure, in contrast to the Maracaibo Basin, where a series of conductive and resistive bodies could be related to active deformation causing the juxtaposition of older geological formations and younger basin sediments.
The Altiplano-Puna plateau, in Central Andes, is the second-largest continental plateau on Earth, extending between 22 degrees and 27 degrees S at an average altitude of 4400 m. The Puna plateau has been formed in consequence of the subduction of the oceanic Nazca Plate beneath the continental South American plate, which has an average crustal thickness of 50 km at this location. A large seismicity cluster, the Jujuy cluster, is observed at depth of 150-250 km beneath the central region of the Puna plateau. The cluster is seismically very active, with hundreds of earthquakes reported and a peak magnitude MW 6.6 on 25th August 2006. The cluster is situated in one of three band of intermediate-depth focus seismicity, which extend parallel to the trench roughly North to South. It has been hypothesized that the Jujuy cluster could be a seismic nest, a compact seismogenic region characterized by a high stationary activity relative to its surroundings. In this study, we collected more than 40 years of data from different catalogs and proof that the cluster meets the three conditions of a seismic nest. Compared to other known intermediate depth nests at Hindu Kush (Afganisthan) or Bucaramanga (Colombia), the Jujuy nest presents an outstanding seismicity rate, with more than 100 M4+ earthquakes per year. We additionally performed a detailed analysis of the rupture process of some of the largest earthquakes in the nest, by means of moment tensor inversion and directivity analysis. We focused on the time period 2017-2018, where the seismic monitoring was the most extended. Our results show that earthquakes in the nest take place within the eastward subducting oceanic plate, but rupture along sub-horizontal planes dipping westward. We suggest that seismicity at Jujuy nest is controlled by dehydration processes, which are also responsible for the generation of fluids ascending to the crust beneath the Puna volcanic region. We use the rupture plane and nest geometry to provide a constraint to maximal expected magnitude, which we estimate as MW -6.7.
Purpose:
Soil erosion by water yields sediment to surface reservoirs, reducing their storage capacities, changing their geometry, and degrading water quality. Sediment reuse, i.e., fertilization of agricultural soils with the nutrient-enriched sediment from reservoirs, has been proposed as a recovery strategy. However, the sediment needs to meet certain criteria. In this study, we characterize sediments from the densely dammed semiarid Northeast Brazil by VNIR-SWIR spectroscopy and assess the effect of spectral resolution and spatial scale on the accuracy of N, P, K, C, electrical conductivity, and clay prediction models.
Methods
Sediment was collected in 10 empty reservoirs, and physical and chemical laboratory analyses as well as spectral measurements were performed. The spectra, initially measured at 1 nm spectral resolution, were resampled to 5 and 10 nm, and samples were analysed for both high and low spectral resolution at three spatial scales, namely (1) reservoir, (2) catchment, and (3) regional scale.
Results
Partial least square regressions performed from good to very good in the prediction of clay and electrical conductivity from reservoir (<40 km(2)) to regional (82,500 km(2)) scales. Models for C and N performed satisfactorily at the reservoir scale, but degraded to unsatisfactory at the other scales. Models for P and K were more unstable and performed from unsatisfactorily to satisfactorily at all scales. Coarsening spectral resolution by up to 10 nm only slightly degrades the models' performance, indicating the potential of characterizing sediment from spectral data captured at lower resolutions, such as by hyperspectral satellite sensors.
Conclusion:
By reducing the costly and time-consuming laboratory analyses, the method helps to promote the sediment reuse as a practice of soil and water conservation.
Urban surface runoff management via best management practices (BMP) and low impact development (LID) has earned significant recognition owing to positive environmental and ecological impacts. However, due to the complexity of the parameters involved, the estimation of LID efficiency in attenuating the urban surface runoff at the watershed scale is challenging. A planning analysis of employing Green Roofs and Infiltration Trenches as BMPs/LIDs practices for urban surface runoff control is presented in this study. A multi-objective optimization decision-making framework is established by coupling SWMM (Storm Water Management Model) with NSGA-II models to check the performance of BMPs/LIDs concerning the cost-benefit analysis of LID at the watershed scale. Two urbanized areas belonging to Central Delhi in India were used as case studies. The results showed that the SWMM model is useful in simulating optimization problems for managing urban surface runoff. The optimum scenarios efficiently minimized the urban runoff volume while maintaining the BMPs/LIDs implementation costs and size. With BMPs/LIDs implementation, the reduction in runoff volume increases as expenses increase initially; however, there is no noticeable reduction in flood volume after a certain threshold. Contrasted with the haphazard arrangement of BMPs/LIDs, the proposed approach demonstrates 22%-24% runoff reductions for the same expenditures in watershed 1 and 23%-26% in watershed 2. The result of the study provides insights into planning and management of the urban surface runoff control with LID practices. The proposed framework assists the hydrologists in optimum selection and placements of BMPs/LIDs practices to acquire the most extreme ecological advantages with the least expenses.
Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.
We present two new empirical models of radiation belt electron flux at geostationary orbit. GOES-15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES-15 flux values and an upper boundary condition scaling factor (BF). The GOES-15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES-15 flux and therefore were used as the exogenous inputs. The NARX-derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July-November 1990, independent of in-situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite-era reconstruction showcases the potential to use these neural network-derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17-24 (Saikin et al., 2021, ), for which the results featured in this paper were used.
Regional flood frequency analysis (RFFA) is a powerful method for interrogating hydrological series since it combines observational time series from several sites within a region to estimate risk-relevant statistical parameters with higher accuracy than from single-site series. Since RFFA extreme value estimates depend on the shape of the selected distribution of the data-generating stochastic process, there is need for a suitable goodness-of-distributional-fit measure in order to optimally utilize given data. Here we present a novel, least-squares-based measure to select the optimal fit from a set of five distributions, namely Generalized Extreme Value (GEV), Generalized Logistic, Gumbel, Log-Normal Type III and Log-Pearson Type III. The fit metric is applied to annual maximum discharge series from six hydrological stations along the Sava River in South-eastern Europe, spanning the years 1961 to 2020. Results reveal that (1) the Sava River basin can be assessed as hydrologically homogeneous and (2) the GEV distribution provides typically the best fit. We offer hydrological-meteorological insights into the differences among the six stations. For the period studied, almost all stations exhibit statistically insignificant trends, which renders the conclusions about flood risk as relevant for hydrological sciences and the design of regional flood protection infrastructure.
Climate Benefits of Cleaner Energy Transitions in East and South Asia Through Black Carbon Reduction
(2022)
The state of air pollution has historically been tightly linked to how we produce and use energy. Air pollutant emissions over Asia are now changing rapidly due to cleaner energy transitions; however, magnitudes of benefits for climate and air quality remain poorly quantified. The associated risks involve adverse health impacts, reduced agricultural yields, reduced freshwater availability, contributions to climate change, and economic costs. We focus particularly on climate benefits of energy transitions by making first-time use of two decades of high quality observations of atmospheric loading of light-absorbing black carbon (BC) over Kanpur (South Asia) and Beijing (East Asia) and relating these observations to changing energy, emissions, and economic trends in India and China. Our analysis reveals that absorption aerosol optical depth (AAOD) due to BC has decreased substantially, by 40% over Kanpur and 60% over Beijing between 2001 and 2017, and thus became decoupled from regional economic growth. Furthermore, the resultant decrease in BC emissions and BC AAOD over Asia is regionally coherent and occurs primarily due to transitions into cleaner energies (both renewables and fossil fuels) and not due to the decrease in primary energy supply or decrease in use of fossil use and biofuels and waste. Model simulations show that BC aerosols alone contribute about half of the surface temperature change (warming) of the total forcing due to greenhouse gases, natural and internal variability, and aerosols, thus clearly revealing the climate benefits due to a reduction in BC emissions, which would significantly reduce global warming. However, this modeling study excludes responses from natural variability, circulation, and sea ice responses, which cause relatively strong temperature fluctuations that may mask signals from BC aerosols. Our findings show additional benefits for climate (beyond benefits of CO2 reduction) and for several other issues of sustainability over South and East Asia, provide motivation for ongoing cleaner energy production, and consumption transitions, especially when they are associated with reduced emissions of air pollutants. Such an analysis connecting the trends in energy transitions and aerosol absorption loading, unavailable so far, is crucial for simulating the aerosol climate impacts over Asia which is quite uncertain.
In this study, 3-D models of P-wave velocity (Vp) and P-wave and S-wave ratio (Vp/Vs) of the crust and upper mantle in the Eastern and eastern Southern Alps (northern Italy and southern Austria) were calculated using local earthquake tomography (LET). The data set includes high-quality arrival times from well-constrained hypocenters observed by the dense, temporary seismic networks of the AlpArray AASN and SWATH-D. The resolution of the LET was checked by synthetic tests and analysis of the model resolution matrix. The small inter-station spacing (average of similar to 15 km within the SWATH-D network) allowed us to image crustal structure at unprecedented resolution across a key part of the Alps. The derived P velocity model revealed a highly heterogeneous crustal structure in the target area. One of the main findings is that the lower crust is thickened, forming a bulge at 30-50 km depth just south of and beneath the Periadriatic Fault and the Tauern Window. This indicates that the lower crust decoupled both from its mantle substratum as well as from its upper crust. The Moho, taken to be the iso-velocity contour of Vp = 7.25 km/s, agrees with the Moho depth from previous studies in the European and Adriatic forelands. It is shallower on the Adriatic side than on the European side. This is interpreted to indicate that the European Plate is subducted beneath the Adriatic Plate in the Eastern and eastern Southern Alps.
Phase transitions in molecular crystals are often determined by intermolecular interactions. The cage complex of [Co(C12H30N8)](3+) . 3 NO3- is reported to undergo a disorder-order phase transition at T-c1 approximate to 133 K upon cooling. Temperature-dependent neutron and synchrotron diffraction experiments revealed satellite reflections in addition to main reflections in the diffraction patterns below T-c1. The modulation wave vector varies as function of temperature and locks in at T-c3 approximate to 98 K. Here, we demonstrate that the crystal symmetry lowers from hexagonal to monoclinic in the incommensurately modulated phases in T-c1<T<T-c3. Distinctive levels of competitions: trade-off between longer N-H...O and shorter C-H...O hydrogen bonds; steric constraints to dense C-H...O bonds give rise to pronounced modulation of the basic structure. Severely frustrated crystal packing in the incommensurate phase is precursor to optimal balance of intermolecular interactions in the lock-in phase.
Volcano-seismic signals such as long-period events and tremor are important indicators for volcanic activity and unrest. However, their wavefield is complex and characterization and location using traditional seismological instrumentation is often difficult.
In 2019 we recorded the full seismic wavefield using a newly developed 3C rotational sensor co-located with a 3C traditional seismometer on Etna, Italy. We compare the performance of the rotational sensor, the seismometer and the Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE) seismic network with respect to the analysis of complex volcano-seismic signals. We create event catalogs for volcano-tectonic (VT) and long-period (LP) events combining a STA/LTA algorithm and cross-correlations.
The event detection based on the rotational sensor is as reliable as the seismometer-based detection. The LP events are dominated by SH-type waves. Derived SH phase velocities range from 500 to 1,000 m/s for LP events and 300-400 m/s for volcanic tremor. SH-waves compose the tremor during weak volcanic activity and SH- and SV-waves during sustained strombolian activity.
We derive back azimuths using (a) horizontal rotational components and (b) vertical rotation rate and transverse acceleration. The estimated back azimuths are consistent with the INGV-OE event location for (a) VT events with an epicentral distance larger than 3 km and some closer events, (b) LP events and tremor in the main crater area. Measuring the full wavefield we can reliably analyze the back azimuths, phase velocities and wavefield composition for VT, LP events and tremor in regions that are difficult to access such as volcanoes.
Extracting information about past tectonic or climatic environmental changes from sedimentary records is a key objective of provenance research. Interpreting the imprint of such changes remains challenging as signals might be altered in the sediment-routing system.
We investigate the sedimentary provenance of the Oligocene/Miocene Upper Austrian Northern Alpine Foreland Basin and its response to the tectonically driven exhumation of the Tauern Window metamorphic dome (28 +/- 1 Ma) in the Eastern European Alps by using the unprecedented combination of Nd isotopic composition of bulk-rock clay-sized samples and partly previously published multi-proxy (Nd isotopic composition, trace-element geochemistry, U-Pb dating) sand-sized apatite single-grain analysis.
The basin offers an excellent opportunity to investigate environmental signal propagation into the sedimentary record because comprehensive stratigraphic and seismic datasets can be combined with present research results. The bulk-rock clay-sized fraction epsilon Nd values of well-cutting samples from one well on the northern basin slope remained stable at similar to-9.7 from 27 to 19 Ma but increased after 19 Ma to similar to-9.1. In contrast, apatite single-grain distributions, which were extracted from 22 drill-core samples, changed significantly around 23.3 Ma from apatites dominantly from low-grade (<upper amphibolite-facies) metamorphic sources with Permo-Mesozoic and late Variscan U-Pb ages and epsilon Nd values of -4.4 to dominantly high-grade metamorphic apatites with late Variscan U-Pb ages and epsilon Nd values of -2.2.
The change in apatite single-grain distributions at 23.3 Ma is interpreted to result from the exposure of a new Upper Austroalpine source nappe with less negative epsilon Nd values triggered by the ongoing Tauern Window exhumation. Combining these data with the clay-sized bulk-rock epsilon Nd values reveals that the provenance changed 4-5 Myrs later at 19 Ma in the clay-sized fraction.
Reasons for the delayed provenance-change recording are rooted in the characteristics of the applied methods.
Whereas single-grain distributions of orogen-wide sediment-routing systems can be dominated by geographically small areas with high erosion and mineral fertility rates, bulk-rock methods integrate over the entire drainage basin, thus diminishing extreme values. Hence, by combining these two methods, spatial information are uncovered, enabling a previously unattained understanding of the underlying environmental change.
White mica and tourmaline are the dominant hydrothermal alteration minerals at the world-class Panasqueira W-Sn-Cu deposit in Portugal. Thus, understanding the controls on their chemical composition helps to constrain ore formation processes at this deposit and determine their usefulness as pathfinder minerals for mineralization in general. We combine whole-rock geochemistry of altered and unaltered metasedimentary host rocks with in situ LA-ICP-MS measurements of tourmaline and white mica from the alteration halo. Principal component analysis (PCA) is used to better identify geochemical patterns and trends of hydrothermal alteration in the datasets. The hydrothermally altered metasediments are enriched in As, Sn, Cs, Li, W, F, Cu, Rb, Zn, Tl, and Pb relative to unaltered samples. In situ mineral analyses show that most of these elements preferentially partition into white mica over tourmaline (Li, Rb, Cs, Tl, W, and Sn), whereas Zn is enriched in tourmaline. White mica has distinct compositions in different settings within the deposit (greisen, vein selvages, wall rock alteration zone, late fault zone), indicating a compositional evolution with time. In contrast, tourmaline from different settings overlaps in composition, which is ascribed to a stronger dependence on host rock composition and also to the effects of chemical zoning and microinclusions affecting the LA-ICP-MS analyses. Hence, in this deposit, white mica is the better recorder of the fluid composition. The calculated trace-element contents of the Panasqueira mineralizing fluid based on the mica data and estimates of mica-fluid partition coefficients are in good agreement with previous fluid-inclusion analyses. A compilation of mica and tourmaline trace-element compositions from Panasqueira and other W-Sn deposits shows that white mica has good potential as a pathfinder mineral, with characteristically high Li, Cs, Rb, Sn, and W contents. The trace-element contents of hydrothermal tourmaline are more variable. Nevertheless, the compiled data suggest that high Sn and Li contents are distinctive for tourmaline from W-Sn deposits.
The Devonian Las Chacras-Potrerillos batholith comprises six nested monzonitic to granitic intrusions with metaluminous to weakly peraluminous composition and a Sr-Nd isotopic signature indicating a dominantly juvenile mantle-derived source. The chemically most evolved units in the southern batholith contain a large number of intra-granitic, pod-shaped tourmaline-bearing pegmatites. This study uses in situ chemical and boron isotopic analyses of tourmaline from nine of these pegmatites to discuss their relationship to the respective host intrusions and the implications of their B-isotope composition for the source and evolution of the magmas. The tourmalines reveal a diversity in element composition (e.g., FeO, MgO, TiO2, CaO, MnO, F) which distinguishes individual pegmatites from one another. However, all have a narrow 5 11 B range of -13.7 to -10.5%0 (n = 100) which indicates a relatively uniform magmatic system and similar temperature conditions during tourmaline crystallization. The average delta(11) B value of -11.7%0 is typical for S-type granites and is within the range reported for peraluminous granites. pegmatites, and metamorphic units of the Ordovician basement into which the Las Chacras-Potrerillos batholith intruded. The B-isotope evidence argues for a crustal boron source like that of the Ordovician basement, in contrast to the metaluminous to weakly peraluminous composition and juvenile initial Sr and Nd isotope ratios of the Las Chacras-Potrerillos batholith magmas. We propose that the boron was not derived from the magma source region but was incorporated from dehydration melting of elastic metasedimentary rocks higher up in the crustal column.
The Altiplano-Puna Plateau holds several shallow lakes, which are very sensitive to climate changes. This work is focused on a high-altitude lake system called Lagunas de Vilama (LVS), located in a complex climatic transition area with scarcity of continuous and homogeneous instrumental records. The objective of this study is to determine the regional spatial-temporal variability of precipitation and evaluate the seasonal and interannual lake responses. We use a lake-surfaces record derived from Landsat images to investigate links with regional precipitations and different climatic forcings. The results reveal that austral summer and autumn precipitations control the variability of the annual lake-surfaces. Also, we found intra-annual and interannual lags in the lake responses to precipitations, and identified several wet and dry stages. Our results show negative trends in precipitations and lake-surfaces, whose were strengthened by a shift to a warm phase of the Atlantic Multidecadal Oscillation in the 1990s. The El Nino Southern Oscillation, Pacific Decadal Oscillation, and Southern Annular Mode also exert a strong influence in the region. This study demonstrates that the variability of LVS lakes is strongly related to the South American Monsoon System dynamics and large-scale climate fordngs from the Pacific and Atlantic Oceans. This work provides novel indices which demonstrated to be good indicators of regional hydroclimatological variability for this region of South America.
Long- and short-term monitoring of a dam in response to seasonal changes and ground motion loading
(2021)
An experimental multi-parameter structural monitoring system has been installed on the Kurpsai dam, western Kyrgyz Republic. This system consists of equipment for seismic and strain measurements for making longer- (days, weeks, months) and shorter- (minutes, hours) term observations, dealing with, for example seasonal (longer) effects or the response of the dam to ground motion from noise or seismic events. Fibre-optic strain sensors allow the seasonal and daily opening and closing of the spaces between the dam's segments to be tracked. For the seismic data, both amplitude (in terms of using differences in amplitudes in the Fourier spectra for mapping the modes of vibration of the dam) and their time-frequency distribution for a set of small to moderate seismic events are investigated and the corresponding phase variabilities (in terms of lagged coherency) are evaluated. Even for moderate levels of seismic-induced ground motion, some influence on the structural response can be detected, which then sees the dam quickly return to its original state. A seasonal component was identified in the strain measurements, while levels of noise arising from the operation of the dam's generators and associated water flow have been provisionally identified.
Centroid moment tensor (CMT) parameters can be estimated from seismic waveforms. Since these data indirectly observe the deformation process, CMTs are inferred as solutions to inverse problems which are generally underdetermined and require significant assumptions, including assumptions about data noise. Broadly speaking, we consider noise to include both theory and measurement errors, where theory errors are due to assumptions in the inverse problem and measurement errors are caused by the measurement process. While data errors are routinely included in parameter estimation for full CMTs, less attention has been paid to theory errors related to velocity-model uncertainties and how these affect the resulting moment-tensor (MT) uncertainties. Therefore, rigorous uncertainty quantification for CMTs may require theory-error estimation which becomes a problem of specifying noise models. Various noise models have been proposed, and these rely on several assumptions. All approaches quantify theory errors by estimating the covariance matrix of data residuals. However, this estimation can be based on explicit modelling, empirical estimation and/or ignore or include covariances. We quantitatively compare several approaches by presenting parameter and uncertainty estimates in nonlinear full CMT estimation for several simulated data sets and regional field data of the M-1 4.4, 2015 June 13 Fox Creek, Canada, event. While our main focus is at regional distances, the tested approaches are general and implemented for arbitrary source model choice. These include known or unknown centroid locations, full MTs, deviatoric MTs and double-couple MTs. We demonstrate that velocity-model uncertainties can profoundly affect parameter estimation and that their inclusion leads to more realistic parameter uncertainty quantification. However, not all approaches perform equally well. Including theory errors by estimating non-stationary (non-Toeplitz) error covariance matrices via iterative schemes during Monte Carlo sampling performs best and is computationally most efficient. In general, including velocity-model uncertainties is most important in cases where velocity structure is poorly known.
Ice-rich permafrost has been subject to abrupt thaw and thermokarst formation in the past and is vulnerable to current global warming. The ice-rich permafrost domain includes Yedoma sediments that have never thawed since deposition during the late Pleistocene and Alas sediments that were formed by previous thermokarst processes during the Lateglacial and Holocene warming. Permafrost thaw unlocks organic carbon (OC) and minerals from these deposits and exposes OC to mineralization. A portion of the OC can be associated with iron (Fe), a redox-sensitive element acting as a trap for OC. Post-depositional thaw processes may have induced changes in redox conditions in these deposits and thereby affected Fe distribution and interactions between OC and Fe, with knock-on effects on the role that Fe plays in mediating present day OC mineralization. To test this hypothesis, we measured Fe concentrations and proportion of Fe oxides and Fe complexed with OC in unthawed Yedoma and previously thawed Alas deposits. Total Fe concentrations were determined on 1,292 sediment samples from the Yedoma domain using portable X-ray fluorescence; these concentrations were corrected for trueness using a calibration based on a subset of 144 samples measured by inductively coupled plasma optical emission spectrometry after alkaline fusion (R (2) = 0.95). The total Fe concentration is stable with depth in Yedoma deposits, but we observe a depletion or accumulation of total Fe in Alas deposits, which experienced previous thaw and/or flooding events. Selective Fe extractions targeting reactive forms of Fe on unthawed and previously thawed deposits highlight that about 25% of the total Fe is present as reactive species, either as crystalline or amorphous oxides, or complexed with OC, with no significant difference in proportions of reactive Fe between Yedoma and Alas deposits. These results suggest that redox driven processes during past thermokarst formation impact the present-day distribution of total Fe, and thereby the total amount of reactive Fe in Alas versus Yedoma deposits. This study highlights that ongoing thermokarst lake formation and drainage dynamics in the Arctic influences reactive Fe distribution and thereby interactions between Fe and OC, OC mineralization rates, and greenhouse gas emissions.
How biased are our models?
(2021)
Geophysical process simulations play a crucial role in the understanding of the subsurface. This understanding is required to provide, for instance, clean energy sources such as geothermal energy. However, the calibration and validation of the physical models heavily rely on state measurements such as temperature. In this work, we demonstrate that focusing analyses purely on measurements introduces a high bias. This is illustrated through global sensitivity studies. The extensive exploration of the parameter space becomes feasible through the construction of suitable surrogate models via the reduced basis method, where the bias is found to result from very unequal data distribution. We propose schemes to compensate for parts of this bias. However, the bias cannot be entirely compensated. Therefore, we demonstrate the consequences of this bias with the example of a model calibration.
This study is trying to understand the pre-eruptive magma storage and crystallization conditions of the Middle Miocene aged, silica-saturated trachytic rocks of the Afyon Volcanic Complex (AVC) in Western Anatolia, Turkey. Those rocks can be divided by their high K2O, K2O/Na2O ratio and Mg# into two groups, namely the intermediate-potassic (IPG) and the ultrapotassic (UPG). Here we are comparing calculated pressure (P) - temperature (T) conditions derived from geothermobarometric calculations of natural samples with results of high-pressure, high-temperature phase equilibria experiments. IPG samples are richer in silica (57-64 wt% SiO2), whereas UPG samples show intermediate SiO2 contents of 56-58 wt%. UPG are having high K2O contents ((>)9 wt %), K2O/Na2O ratios ((>)10 wt%) and Mg# values (75-77). IPG phenocrysts comprise plagioclase + biotite + amphibole + clinopyroxene +/- orthopyroxene +/- sanidine +/- phlogopite and oxides, while UPG mineralogical assemblage consists of amphibole + phlogopite + clinopyroxene + olivine + sanidine and oxides. IPG and UPG are enriched in Large-Ion Lithophile Elements (LILE), and both have negative anomalies in Nb, Sr, Zr and Ti elements. Additionally, IPG shows positive anomalies in Pb. Both IPG and UPG display enrichment in Light Rare Earth Elements (LREE), while IPG shows a more significant negative anomaly in Eu when compared to UPG. Plagioclase fractionation may play a role in magma generation. In IPG samples Ni and Cr values range between (3.3-18.8 ppm) and (2.6-27.8 ppm), respectively; whereas UPG samples have (119.1-120.7 ppm) Ni and (212.1-219.9 ppm) Cr. Dy/Yb ratios of IPG and UPG are higher than 2 and may indicate that garnet was present in the source. Geothermobarometric calculations for natural IPG clinopyroxene-melt pairs imply higher PT conditions (Dogan-Kulahci et al., 2015), while in this study high-pressure/high-temperature (HP/HT) phase equilibria experiments recreated the natural mineral assemblage at 2-4 kbar, 6-9 km and c. 900 degrees C. New plagioclase-melt calculations have confirmed lower mean magma storage temperatures, which are closer to the experimental results but still slightly elevated. Thus, trace element results of the natural rocks and experimental data may imply that a deep garnet-bearing magma source mixed with shallower magmas (IPG) was feeding the volcanic eruption.
Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.
The Kohat fold and thrust belt in Pakistan shows a significantly different structural style due to the structural evolution on the double décollement compared to the rest of the Subhimalaya. In order to better understand the spatio-temporal structural evolution of the Kohat fold and thrust belt, we combine balanced cross sections with apatite (U?Th-Sm)/He (AHe) and apatite fission track (AFT) dating. The AHe and AFT ages appear to be totally reset, allowing us to date exhumation above structural ramps. The results suggest that deformation began on the frontal Surghar thrust at-15 Ma, predating or coeval with the development of the Main Boundary thrust at-12 Ma. Deformation propagated southward from the Main Boundary thrust on double de?collements between 10 Ma and 2 Ma, resulting in a disharmonic structural style inside the Kohat fold and thrust belt. Thermal modeling of the thermochronologic data suggest that samples inside Kohat fold and thrust belt experienced cooling due to formation of the duplexes; this deformation facilitated tectonic thickening of the wedge and erosion of the Miocene to Pliocene foreland strata. The spatial distribution of AHe and AFT ages in combination with the structural forward model suggest that, in the Kohat fold and thrust belt, the wedge deformed in-sequence as a supercritical wedge (-15-12 Ma), then readjusted by out-sequence deformation (-12-0 Ma) within the Kohat fold and thrust belt into a sub-critical wedge.
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
The correct orientation of seismic sensors is critical for studies such as full moment tensor inversion, receiver function analysis, and shear-wave splitting. Therefore, the orientation of horizontal components needs to be checked and verified systematically. This study relies on two different waveform-based approaches, to assess the sensor orientations of the broadband network of the Kandilli Observatory and Earthquake Research Institute (KOERI). The network is an important backbone for seismological research in the Eastern Mediterranean Region and provides a comprehensive seismic data set for the North Anatolian fault. In recent years, this region became a worldwide field laboratory for continental transform faults. A systematic survey of the sensor orientations of the entire network, as presented here, facilitates related seismic studies. We apply two independent orientation tests, based on the polarization of P waves and Rayleigh waves to 123 broadband seismic stations, covering a period of 15 yr (2004-2018). For 114 stations, we obtain stable results with both methods. Approximately, 80% of the results agree with each other within 10 degrees. Both methods indicate that about 40% of the stations are misoriented by more than 10 degrees. Among these, 20 stations are misoriented by more than 20 degrees. We observe temporal changes of sensor orientation that coincide with maintenance work or instrument replacement. We provide time-dependent sensor misorientation correction values for the KOERI network in the supplemental material.
Water bodies are a highly abundant feature of Arctic permafrost ecosystems and strongly influence their hydrology, ecology and biogeochemical cycling. While very high resolution satellite images enable detailed mapping of these water bodies, the increasing availability and abundance of this imagery calls for fast, reliable and automatized monitoring. This technical work presents a largely automated and scalable workflow that removes image noise, detects water bodies, removes potential misclassifications from infrastructural features, derives lake shoreline geometries and retrieves their movement rate and direction on the basis of ortho-ready very high resolution satellite imagery from Arctic permafrost lowlands. We applied this workflow to typical Arctic lake areas on the Alaska North Slope and achieved a successful and fast detection of water bodies. We derived representative values for shoreline movement rates ranging from 0.40-0.56 m yr(-1) for lake sizes of 0.10 ha-23.04 ha. The approach also gives an insight into seasonal water level changes. Based on an extensive quantification of error sources, we discuss how the results of the automated workflow can be further enhanced by incorporating additional information on weather conditions and image metadata and by improving the input database. The workflow is suitable for the seasonal to annual monitoring of lake changes on a sub-meter scale in the study areas in northern Alaska and can readily be scaled for application across larger regions within certain accuracy limitations.
This study deals with the East Beni Suef Basin (Eastern Desert, Egypt) and aims to evaluate the source-generative potential, reconstruct the burial and thermal history, examine the most influential parameters on thermal maturity modeling, and improve on the models already published for the West Beni Suef to ultimately formulate a complete picture of the whole basin evolution.
Source rock evaluation was carried out based on TOC, Rock-Eval pyrolysis, and visual kerogen petrography analyses. Three kerogen types (II, II/III, and III) are distinguished in the East Beni Suef Basin, where the Abu Roash "F" Member acts as the main source rock with good to excellent source potential, oil-prone mainly type II kerogen, and immature to marginal maturity levels.
The burial history shows four depositional and erosional phases linked with the tectonic evolution of the basin. A hiatus (due to erosion or non-deposition) has occurred during the Late Eocene-Oligocene in the East Beni Suef Basin, while the West Beni Suef Basin has continued subsiding.
Sedimentation began later (Middle to Late Albian) with lower rates in the East Beni Suef Basin compared with the West Beni Suef Basin (Early Albian). The Abu Roash "F" source rock exists in the early oil window with a present-day transformation ratio of about 19% and 21% in the East and West Beni Suef Basin, respectively, while the Lower Kharita source rock, which is only recorded in the West Beni Suef Basin, has reached the late oil window with a present-day transformation ratio of about 70%.
The magnitude of erosion and heat flow have proportional and mutual effects on thermal maturity.
We present three possible scenarios of basin modeling in the East Beni Suef Basin concerning the erosion from the Apollonia and Dabaa formations.
Results of this work can serve as a basis for subsequent 2D and/or 3D basin modeling, which are highly recommended to further investigate the petroleum system evolution of the Beni Suef Basin.
The subsurface is a temporally dynamic and spatially heterogeneous compartment of the Earth's critical zone, and biogeochemical transformations taking place in this compartment are crucial for the cycling of nutrients.
The impact of spatial heterogeneity on such microbially mediated nutrient cycling is not well known, which imposes a severe challenge in the prediction of in situ biogeochemical transformation rates and further of nutrient loading contributed by the groundwater to the surface water bodies.
Therefore, we used a numerical modelling approach to evaluate the sensitivity of groundwater microbial biomass distribution and nutrient cycling to spatial heterogeneity in different scenarios accounting for various residence times.
The model results gave us an insight into domain characteristics with respect to the presence of oxic niches in predominantly anoxic zones and vice versa depending on the extent of spatial heterogeneity and the flow regime.
The obtained results show that microbial abundance, distribution, and activity are sensitive to the applied flow regime and that the mobile (i.e. observable by groundwater sampling) fraction of microbial biomass is a varying, yet only a small, fraction of the total biomass in a domain. Furthermore, spatial heterogeneity resulted in anaerobic niches in the domain and shifts in microbial biomass between active and inactive states. The lack of consideration of spatial heterogeneity, thus, can result in inaccurate estimation of microbial activity. In most cases this leads to an overestimation of nutrient removal (up to twice the actual amount) along a flow path.
We conclude that the governing factors for evaluating this are the residence time of solutes and the Damkohler number (Da) of the biogeochemical reactions in the domain. We propose a relationship to scale the impact of spatial heterogeneity on nutrient removal governed by the logioDa.
This relationship may be applied in upscaled descriptions of microbially mediated nutrient cycling dynamics in the subsurface thereby resulting in more accurate predictions of, for example, carbon and nitrogen cycling in groundwater over long periods at the catchment scale.
ABSTRACT: Structural evolution of cesium triiodide at high pressures has been revealed by synchrotron single-crystal X-ray diffraction. Cesium triiodide undergoes a first-order phase transition above 1.24(3) GPa from an orthorhombic to a trigonal system. This transition is coupled with severe reorganization of the polyiodide network from a layered to three-dimensional architecture. Quantum chemical calculations show that even though the two polymorphic phases are nearly isoenergetic under ambient conditions, the PV term is decisive in stabilizing the trigonal polymorph above the transition point. Phonon calculations using a non-local correlation functional that accounts for dispersion interactions confirm that this polymorph is dynamically unstable under ambient conditions. The high-pressure behavior of crystalline CsI3 can be correlated with other alkali metal trihalides, which undergo a similar sequence of structural changes upon load.
Alpine glacial erosion exerts a first-order control on mountain topography and sediment production, but its mechanisms are poorly understood. Observational data capable of testing glacial erosion and transport laws in glacial models are mostly lacking. New insights, however, can be gained from detrital tracer thermochronology. Detrital tracer thermochronology works on the premise that thermochronometer bedrock ages vary systematically with elevation, and that detrital downstream samples can be used to infer the source elevation sectors of sediments. We analyze six new detrital samples of different grain sizes (sand and pebbles) from glacial deposits and the modern river channel integrated with data from 18 previously analyzed bedrock samples from an elevation transect in the Leones Valley, Northern Patagonian Icefield, Chile (46.7 degrees S). We present 622 new detrital zircon (U-Th)/He (ZHe) single-grain analyses and 22 new bedrock ZHe analyses for two of the bedrock samples to determine age reproducibility. Results suggest that glacial erosion was focused at and below the Last Glacial Maximum and neoglacial equilibrium line altitudes, supporting previous modeling studies. Furthermore, grain age distributions from different grain sizes (sand, pebbles) might indicate differences in erosion mechanisms, including mass movements at steep glacial valley walls. Finally, our results highlight complications and opportunities in assessing glacigenic environments, such as dynamics of sediment production, transport, transient storage, and final deposition, that arise from settings with large glacio-fluvial catchments.
Frequency-domain electromagnetic (FDEM) data are commonly inverted to characterize subsurface geoelectrical properties using smoothness constraints in 1D inversion schemes assuming a layered medium.
Smoothness constraints are suitable for imaging gradual transitions of subsurface geoelectrical properties caused, for example, by varying sand, clay, or fluid content. However, such inversion approaches are limited in characterizing sharp interfaces. Alternative regularizations based on the minimum gradient support (MGS) stabilizers can, instead, be used to promote results with different levels of smoothness/sharpness selected by simply acting on the so-called focusing parameter.
The MGS regularization has been implemented for different kinds of geophysical data inversion strategies. However, concerning FDEM data, the MGS regularization has only been implemented for vertically constrained inversion (VCI) approaches but not for laterally constrained inversion (LCI) approaches.
We present a novel LCI approach for FDEM data using the MGS regularization for the vertical and lateral direction. Using synthetic and field data examples, we demonstrate that our approach can efficiently and automatically provide a set of model solutions characterized by different levels of sharpness and variable lateral consistencies.
In terms of data misfit, the obtained set of solutions contains equivalent models allowing us also to investigate the non-uniqueness of FDEM data inversion.
Nocardioides alcanivorans sp. nov., a novel hexadecane-degrading species isolated from plastic waste
(2022)
Strain NGK65(T), a novel hexadecane degrading, non-motile, Gram-positive, rod-to-coccus shaped, aerobic bacterium, was isolated from plastic polluted soil sampled at a landfill.
Strain NGK65(T) hydrolysed casein, gelatin, urea and was catalase-positive. It optimally grew at 28 degrees C. in 0-1% NaCl and at pH 7.5-8.0. Glycerol, D-glucose, arbutin, aesculin, salicin, potassium 5-ketogluconate. sucrose, acetate, pyruvate and hexadecane were used as sole carbon sources.
The predominant membrane fatty acids were iso-C-16:0 followed by iso-C(17:)0 and C-18:1 omega 9c. The major polar lipids were phosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol and hydroxyphosphatidylinositol.
The cell-wall peptidoglycan type was A3 gamma, with LL-diaminopimelic acid and glycine as the diagnostic amino acids. MK 8 (H-4) was the predominant menaquinone. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK65(T) belongs to the genus Nocardioides (phylum Actinobacteria). appearing most closely related to Nocardioides daejeonensis MJ31(T) (98.6%) and Nocardioides dubius KSL-104(T) (98.3%).
The genomic DNA G+C content of strain NGK65(T) was 68.2%.
Strain NGK65(T) and the type strains of species involved in the analysis had average nucleotide identity values of 78.3-71.9% as well as digital DNA-DNA hybridization values between 22.5 and 19.7%, which clearly indicated that the isolate represents a novel species within the genus Nocardioides.
Based on phenotypic and molecular characterization, strain NGK65(T) can clearly be differentiated from its phylogenetic neighbours to establish a novel species, for which the name Nocardioides alcanivorans sp. nov. is proposed.
The type strain is NGK65(T) (=DSM 113112(T)=NCCB 100846(T)).
Understanding the key factors influencing the water quality of large river systems forms an important basis for the assessment and protection of cross-regional ecosystems and the implementation of adapted water management concepts. However, identifying these factors requires in-depth comprehension of the unique environmental systems, which can only be achieved by detailed water quality monitoring.
Within the scope of the joint science and sports event "Elbschwimmstaffel" (swimming relay on the river Elbe) in June/July 2017 organized by the German Ministry of Education and Research, water quality data were acquired along a 550 km long stretch of the Elbe River in Germany. During the survey, eight physiochemical water quality parameters were recorded in high spatial and temporal resolution with the BIOFISH multisensor system. Multivariate statistical methods were applied to identify and delineate processes influencing the water quality.
The BIOFISH dataset revealed that phytoplankton activity has a major impact on the water quality of the Elbe River in the summer months. The results suggest that phytoplankton biomass constitutes a substantial proportion of the suspended particles and that photosynthetic activity of phytoplankton is closely related to significant temporal changes in pH and oxygen saturation.
An evaluation of the BIOFISH data based on the combination of statistical analysis with weather and discharge data shows that the hydrological and meteorological history of the sampled water body was the main driver of phytoplankton dynamics. This study demonstrates the capacity of longitudinal river surveys with the BIOFISH or similar systems for water quality assessment, the identification of pollution sources and their utilization for online in situ monitoring of rivers.
The authors demonstrate that a vegetation system's ability to recover from disturbances-its resilience-can be estimated from its natural variability. Global patterns of resilience loss and gains since the early 1990s reveal shifts towards widespread resilience loss since the early 2000s.
The character and health of ecosystems worldwide is tightly coupled to changes in Earth's climate. Theory suggests that ecosystem resilience-the ability of ecosystems to resist and recover from external shocks such as droughts and fires-can be inferred from their natural variability. Here, we quantify vegetation resilience globally with complementary metrics based on two independent long-term satellite records. We first empirically confirm that the recovery rates from large perturbations can be closely approximated from internal vegetation variability across vegetation types and climate zones. On the basis of this empirical relationship, we quantify vegetation resilience continuously and globally from 1992 to 2017. Long-term vegetation resilience trends are spatially heterogeneous, with overall increasing resilience in the tropics and decreasing resilience at higher latitudes. Shorter-term trends, however, reveal a marked shift towards a global decline in vegetation resilience since the early 2000s, particularly in the equatorial rainforest belt.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
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.
The Walker Circulation (WC) is an east-west trending band of atmospheric circulation cells along the equator and the predominant controller of heat and moisture transport in the tropics. Its variability is closely linked to the sea-surface temperature (SST) changes across the Pacific, the Indian and the Atlantic Oceans and can have pronounced effects on the humidity regimes of the adjacent continents. In recent years, the evolution of the WC during the Plioand Pleistocene epochs has been intensely studied in the context of the effectiveness of the tropics in modulating global climate change (e.g., the intensification of Northern Hemisphere glaciation). However, the onset of the modern WC pattern as well as its global impact during the Plioand Pleistocene is controversially assessed in the literature. For its onset, previous studies have suggested dates ranging between 2.4 and 0.8 million years ago (Myr), while its argued impact ranges from crucially influencing the increase of Northern Hemisphere ice sheet growth by channelling heat and moisture from the tropics into the high latitudes to having no effect on global ice volume changes. In order to achieve a comprehensive understanding of the spatiotemporal evolution of the WC during this time frame, we statistically analysed 30 globally distributed SST records covering the low and high latitudes between 3.5 and 1.5 Myr, encompassing the Late Pliocene to Early Pleistocene. We utilized a statistical change-point regression model to determine significant change points in the SST evolution of the (sub)-tropics and high latitudes that potentially relate to changes in the WC. We find that the WC experienced a multifaceted evolution between the Late Pliocene and the Early Pleistocene with significant transitional steps at-2.7 and-2.1 Ma. Our results suggest after the Late Pliocene, a pre-modern WC set in, which was characterized by a progressively strengthened Pacific Walker Cell alongside a weakened Indian Ocean Walker Cell. This change was potentially triggered by the constriction of the Indonesian seaway, an important transmitter between the Pacific and Indian Ocean. The ensuing mode of the WC intensified until-2.1 Myr, when SST values around the global scale signalled a progressive strengthening of the Indian Walker Cell in phase with the progressive strengthening of the Pacific and Atlantic Cells. Our findings indicate that a shift from a pre-modern to a modern-like WC potentially only occurred during the mid-Pleistocene.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.