TY - JOUR A1 - Krstulović, Marija A1 - Rosa, Angelika D. A1 - Ferreira Sanchez, Dario A1 - Libon, Lélia A1 - Albers, Christian A1 - Merkulova, Margarita A1 - Grolimund, Daniel A1 - Irifune, Tetsuo A1 - Wilke, Max T1 - Effect of temperature on the densification of silicate melts to lower earth's mantle conditions JF - Physics of the earth and planetary interiors N2 - Physical properties of silicate melts play a key role for global planetary dynamics, controlling for example volcanic eruption styles, mantle convection and elemental cycling in the deep Earth. They are significantly modified by structural changes at the atomic scale due to external parameters such as pressure and temperature or due to chemistry. Structural rearrangements such as 4- to 6-fold coordination change of Si with increasing depth may profoundly influence melt properties, but have so far mostly been studied at ambient temperature due to experimental difficulties. In order to investigate the structural properties of silicate melts and their densification mechanisms at conditions relevant to the deep Earth's interior, we studied haplo basaltic glasses and melts (albite-diopside composition) at high pressure and temperature conditions in resistively and laser-heated diamond anvil cells using X-ray absorption near edge structure spectroscopy. Samples were doped with 10 wt% of Ge, which is accessible with this experimental technique and which commonly serves as a structural analogue for the network forming cation Si. We acquired spectra on the Ge K edge up to 48 GPa and 5000 K and derived the average Ge-O coordination number NGe-O, and bond distance RGe-O as functions of pressure. Our results demonstrate a continuous transformation from tetrahedral to octahedral coordination between ca. 5 and 30 GPa at ambient temperature. Above 1600 K the data reveal a reduction of the pressure needed to complete conversion to octahedral coordination by ca. 30 %. The results allow us to determine the influence of temperature on the Si coordination number changes in natural melts in the Earth's interior. We propose that the complete transition to octahedral coordination in basaltic melts is reached at about 40 GPa, corresponding to a depth of ca. 1200 km in the uppermost lower mantle. At the core-mantle boundary (2900 km, 130 GPa, 3000 K) the existence of non-buoyant melts has been proposed to explain observed low seismic wave velocity features. Our results highlight that the melt composition can affect the melt density at such extreme conditions and may strongly influence the structural response. KW - Silicate melts KW - Densification KW - High pressure and high temperature; KW - XANES KW - Coordination number KW - Ultra-low velocity zones Y1 - 2022 U6 - https://doi.org/10.1016/j.pepi.2021.106823 SN - 0031-9201 SN - 1872-7395 VL - 323 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Nooshiri, Nima A1 - Bean, Christopher J. A1 - Dahm, Torsten A1 - Grigoli, Francesco A1 - Kristjansdottir, Sigriour A1 - Obermann, Anne A1 - Wiemer, Stefan T1 - A multibranch, multitarget neural network for rapid point-source inversion in a microseismic environment BT - examples from the Hengill Geothermal Field, Iceland JF - Geophysical journal international N2 - Despite advanced seismological techniques, automatic source characterization for microseismic earthquakes remains difficult and challenging since current inversion and modelling of high-frequency signals are complex and time consuming. For real-time applications such as induced seismicity monitoring, the application of standard methods is often not fast enough for true complete real-time information on seismic sources. In this paper, we present an alternative approach based on recent advances in deep learning for rapid source-parameter estimation of microseismic earthquakes. The seismic inversion is represented in compact form by two convolutional neural networks, with individual feature extraction, and a fully connected neural network, for feature aggregation, to simultaneously obtain full moment tensor and spatial location of microseismic sources. Specifically, a multibranch neural network algorithm is trained to encapsulate the information about the relationship between seismic waveforms and underlying point-source mechanisms and locations. The learning-based model allows rapid inversion (within a fraction of second) once input data are available. A key advantage of the algorithm is that it can be trained using synthetic seismic data only, so it is directly applicable to scenarios where there are insufficient real data for training. Moreover, we find that the method is robust with respect to perturbations such as observational noise and data incompleteness (missing stations). We apply the new approach on synthesized and example recorded small magnitude (M <= 1.6) earthquakes at the Hellisheioi geothermal field in the Hengill area, Iceland. For the examined events, the model achieves excellent performance and shows very good agreement with the inverted solutions determined through standard methodology. In this study, we seek to demonstrate that this approach is viable for microseismicity real-time estimation of source parameters and can be integrated into advanced decision-support tools for controlling induced seismicity. KW - Neural networks KW - fuzzy logic KW - Computational seismology KW - Induced seismicity KW - Earthquake source observations Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab511 SN - 0956-540X SN - 1365-246X VL - 229 IS - 2 SP - 999 EP - 1016 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Ezpeleta, Miguel A1 - Parra, Mauricio A1 - Collo, Gilda A1 - Wunderlin, Cecilia A1 - Borrego, Angeles G. A1 - Sobel, Edward A1 - Glodny, Johannes T1 - Thermochronometry unveils ancient thermal regimes in the NW Pampean Ranges, Argentina BT - from Mesozoic rifting to Miocene flat-slab subduction JF - Basin research N2 - Reconstructing thermal histories in thrust belts is commonly used to infer the age and rates of thrusting and hence the driving mechanisms of orogenesis. In areas where ancient basins have been incorporated into the orogenic wedge, a quantitative reconstruction of the thermal history helps distinguish among potential mechanisms responsible for heating events. We present such a reconstruction for the Ischigualasto-Villa Union basin in the western Pampean Ranges of Argentina, where Triassic rifting and late Cretaceous-Cenozoic retroarc foreland basin development has been widely documented, including Miocene flat-slab subduction. We report results of organic and inorganic thermal indicators acquired along three stratigraphic sections, including vitrinite reflectance and X-ray diffractometry in claystones and new thermochronological [(apatite fission-track and apatite and zircon [U-Th]/He)] analyses. Despite up to 5 km-thick Cenozoic overburden and unlike previously thought, the thermal peak in the basin is not due to Cenozoic burial but occurred in the Triassic, associated with a high heat flow of up to 90 mWm(-2) and <2 km of burial, which heated the base of the Triassic strata to similar to 160 degrees C. Following exhumation, attested by the development of an unconformity between the Triassic and Late-Cretaceous-Cenozoic sequences, Cenozoic re-burial increased the temperature to similar to 110 degrees C at the base of the Triassic section and only similar to 50 degrees C 7 km upsection, suggesting a dramatic decrease in the thermal gradient. The onset of Cenozoic cooling occurred at similar to 10(-8) Ma, concomitant with sediment accumulation and thus preceding the latest Miocene onset of thrusting that has been independently documented by stratigraphic-cross-cutting relationships. We argue that the onset of cooling is associated with lithospheric refrigeration following establishment of flat-slab subduction, leading to the eastward displacement of the asthenospheric wedge beneath the South American plate. Our study places time and temperature constraints on flat-slab cooling that calls for a careful interpretation of exhumation signals in thrustbelts inferred from thermochronology only. KW - %Ro KW - Cenozoic flat-slab KW - Ischigualasto-Villa Union Basin KW - thermochronological modelling (AFT, AHe and ZHe) KW - Triassic rifting KW - XRD in the clay fraction KW - heat flow KW - burial Y1 - 2022 U6 - https://doi.org/10.1111/bre.12693 SN - 0950-091X SN - 1365-2117 VL - 34 IS - 6 SP - 1983 EP - 2012 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Fernandez-Palomino, Carlos Antonio A1 - Hattermann, Fred A1 - Krysanova, Valentina A1 - Lobanova, Anastasia A1 - Vega-Jacome, Fiorella A1 - Lavado, Waldo A1 - Santini, William A1 - Aybar, Cesar A1 - Bronstert, Axel T1 - A novel high-resolution gridded precipitation dataset for peruvian and ecuadorian watersheds BT - development and hydrological evaluation JF - Journal of hydrometeorology N2 - A novel approach for estimating precipitation patterns is developed here and applied to generate a new hydrologically corrected daily precipitation dataset, called RAIN4PE (Rain for Peru and Ecuador), at 0.1 degrees spatial resolution for the period 1981-2015 covering Peru and Ecuador. It is based on the application of 1) the random forest method to merge multisource precipitation estimates (gauge, satellite, and reanalysis) with terrain elevation, and 2) observed and modeled streamflow data to first detect biases and second further adjust gridded precipitation by inversely applying the simulated results of the ecohydrological model SWAT (Soil and Water Assessment Tool). Hydrological results using RAIN4PE as input for the Peruvian and Ecuadorian catchments were compared against the ones when feeding other uncorrected (CHIRP and ERA5) and gauge-corrected (CHIRPS, MSWEP, and PISCO) precipitation datasets into the model. For that, SWAT was calibrated and validated at 72 river sections for each dataset using a range of performance metrics, including hydrograph goodness of fit and flow duration curve signatures. Results showed that gauge-corrected precipitation datasets outperformed uncorrected ones for streamflow simulation. However, CHIRPS, MSWEP, and PISCO showed limitations for streamflow simulation in several catchments draining into the Pacific Ocean and the Amazon River. RAIN4PE provided the best overall performance for streamflow simulation, including flow variability (low, high, and peak flows) and water budget closure. The overall good performance of RAIN4PE as input for hydrological modeling provides a valuable criterion of its applicability for robust countrywide hydrometeorological applications, including hydroclimatic extremes such as droughts and floods. Significance StatementWe developed a novel precipitation dataset RAIN4PE for Peru and Ecuador by merging multisource precipitation data (satellite, reanalysis, and ground-based precipitation) with terrain elevation using the random forest method. Furthermore, RAIN4PE was hydrologically corrected using streamflow data in watersheds with precipitation underestimation through reverse hydrology. The results of a comprehensive hydrological evaluation showed that RAIN4PE outperformed state-of-the-art precipitation datasets such as CHIRP, ERA5, CHIRPS, MSWEP, and PISCO in terms of daily and monthly streamflow simulations, including extremely low and high flows in almost all Peruvian and Ecuadorian catchments. This underlines the suitability of RAIN4PE for hydrometeorological applications in this region. Furthermore, our approach for the generation of RAIN4PE can be used in other data-scarce regions. KW - Amazon region KW - Complex terrain KW - South America KW - Streamflow KW - Precipitation KW - Hydrology KW - Water budget / balance KW - Inverse methods KW - Mountain meteorology KW - Machine learning Y1 - 2022 U6 - https://doi.org/10.1175/JHM-D-20-0285.1 SN - 1525-755X SN - 1525-7541 VL - 23 IS - 3 SP - 309 EP - 336 PB - American Meteorological Soc. CY - Boston ER - TY - JOUR A1 - Schüler, Lennart A1 - Calabrese, Justin M. A1 - Attinger, Sabine T1 - Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany JF - PLoS one N2 - The SARS-CoV-2 virus has spread around the world with over 100 million infections to date, and currently many countries are fighting the second wave of infections. With neither sufficient vaccination capacity nor effective medication, non-pharmaceutical interventions (NPIs) remain the measure of choice. However, NPIs place a great burden on society, the mental health of individuals, and economics. Therefore the cost/benefit ratio must be carefully balanced and a target-oriented small-scale implementation of these NPIs could help achieve this balance. To this end, we introduce a modified SEIRD-class compartment model and parametrize it locally for all 412 districts of Germany. The NPIs are modeled at district level by time varying contact rates. This high spatial resolution makes it possible to apply geostatistical methods to analyse the spatial patterns of the pandemic in Germany and to compare the results of different spatial resolutions. We find that the modified SEIRD model can successfully be fitted to the COVID-19 cases in German districts, states, and also nationwide. We propose the correlation length as a further measure, besides the weekly incidence rates, to describe the current situation of the epidemic. Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0254660 SN - 1932-6203 VL - 16 IS - 8 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Huang, Sichao A1 - Stoof-Leichsenring, Kathleen R. A1 - Liu, Sisi A1 - Courtin, Jeremy A1 - Andreev, Andrej A. A1 - Pestryakova, Luidmila. A. A1 - Herzschuh, Ulrike T1 - Plant sedimentary ancient DNA from Far East Russia covering the last 28,000 years reveals different assembly rules in cold and warm climates JF - Frontiers in Ecology and Evolution N2 - Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules. Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years. Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest. We reveal highest richness during 28-23 ka and a second richness peak during 13-9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed. This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community. KW - sedimentary ancient DNA (sedaDNA) KW - metabarcoding KW - phylogenetic and taxonomic plant diversity KW - Arctic Russia KW - Siberia KW - holocene KW - glacial KW - treeline Y1 - 2021 U6 - https://doi.org/10.3389/fevo.2021.763747 SN - 2296-701X VL - 9 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - West, Charles A1 - Rosolem, Rafael A1 - MacDonald, Alan M. A1 - Cuthbert, Mark O. A1 - Wagener, Thorsten T1 - Understanding process controls on groundwater recharge variability across Africa through recharge landscapes JF - Journal of hydrology N2 - Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information. KW - Groundwater recharge KW - Africa KW - Recharge controls KW - Ground-based estimates KW - Landscapes KW - Comparative hydrology Y1 - 2022 U6 - https://doi.org/10.1016/j.jhydrol.2022.127967 SN - 0022-1694 SN - 1879-2707 VL - 612 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Ghaffar, Salman A1 - Jomaa, Seifeddine A1 - Meon, Günter A1 - Rode, Michael T1 - Spatial validation of a semi-distributed hydrological nutrient transport model JF - Journal of hydrology N2 - Semi-distributed hydrological and water quality models are increasingly used as innovative and scientific-based management tools. However, their application is usually restricted to the gauging stations where they are originally calibrated, limiting their spatial capability. In this study, the semi-distributed hydrological water quality model HYPE (HYdrological Predictions for the Environment) was tested spatially to represent nitrate-N (NO3- N) and total phosphorus (TP) concentrations and loads of the nested and heterogeneous Selke catchment (463 km(2)) in central Germany. First, an automatic calibration procedure and uncertainty analysis were conducted using the DiffeRential Evolution Adaptive Metropolis (DREAM) tool to simulate discharge, NO3--N and TP concentrations. A multi-site and multi-objective calibration approach was applied using three main gauging stations, covering the most important hydro-meteorological and physiographical characteristics of the whole catchment. Second, the model's capability was tested to represent further internal stations, which were not initially considered for calibration. Results showed that discharge was well represented by the model at all three main stations during both calibration (1994-1998) and validation (1999-2014) periods with lowest Nash-Sutcliffe Efficiency (NSE) of 0.71 and maximum Percentage BIAS (PBIAS) of 18.0%. The model was able to reproduce the seasonal dynamics of NO3--N and TP concentrations with low predictive uncertainty at the three main stations, reflected by PBIAS values in the ranges from 16.1% to 6.4% and from 20.0% to 11.5% for NO3--N and TP load simulations, respectively. At internal stations, the model could represent reasonably well the seasonal variation of nutrient concentrations with PBIAS values in the ranges from 9.0% to 14.2% for NO3--N and from 25.3% to 34.3% for TP concentration simulations. Overall, results suggested that the spatial validation of a nutrient transport model can be better ensured when a multi-site and multi-objective calibration approach using archetypical gauging stations is implemented. Further, results revealed that the delineation of sub-catchments should put more focus on hydro-meteorological conditions than on land-use features. KW - HYPE model KW - Nitrate-N KW - Phosphorus KW - internal validation KW - uncertainty KW - analysis KW - archetypical gauging station Y1 - 2021 U6 - https://doi.org/10.1016/j.jhydrol.2020.125818 SN - 0022-1694 SN - 1879-2707 VL - 593 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Janocha, Julian A1 - Smyrak-Sikora, Aleksandra A1 - Senger, Kim A1 - Birchall, Thomas T1 - Seeing beyond the outcrop BT - integration of ground-penetrating radar with digital outcrop models of a paleokarst system JF - Marine and petroleum geology N2 - Paleokarst breccias are a common feature of sedimentary rift basins. The Billefjorden Trough in the High Arctic archipelago of Svalbard is an example of such a rift. Here the Carboniferous stratigraphy exhibits intervals of paleokarst breccias formed by gypsum dissolution. In this study we integrate digital outcrop models (DOMs) with a 2D ground penetrating radar (GPR) survey to extrapolate external irregular paleokarst geometries beyond the 2D outcrops. DOMs are obtained through combining a series of overlapping photographs with structure-frommotion photogrammetry, to create mmto dm-resolution georeferenced DOMs. GPR is typically used for surveying the shallow subsurface and relies on detecting the contrasts in electro-magnetic permittivity. We defined three geophysical facies based on their appearance in GPR. By integrating subsurface geophysical data with DOMs we were able to correlate reflection patterns in GPR with outcrop features. The chaotic nature of paleokarst breccias is seen both in outcrop and GPR. Key horizons in outcrop and the GPR profiles allow tying together observations between these methods. Furthermore, we show that this technique expands the twodimensional outcrop surface into a three-dimensional domain, thus complementing, strengthening and extending outcrop interpretations. KW - Digital geology KW - Svalbard KW - Billefjorden trough KW - Upper Paleozoic KW - Fortet member KW - Dissolution collapse breccia KW - Syn-rift breccia Y1 - 2021 U6 - https://doi.org/10.1016/j.marpetgeo.2020.104833 SN - 0264-8172 SN - 1873-4073 VL - 125 PB - Elsevier Science CY - Amsterdam [u.a.] ER - TY - JOUR A1 - Gomez Zapata, Juan Camilo A1 - Zafrir, Raquel A1 - Pittore, Massimiliano A1 - Merino, Yvonne T1 - Towards a sensitivity analysis in seismic risk with probabilistic building exposure models BT - an application in Valparaiso, Chile using ancillary open-source data and parametric ground motions JF - ISPRS International Journal of Geo-Information N2 - Efforts have been made in the past to enhance building exposure models on a regional scale with increasing spatial resolutions by integrating different data sources. This work follows a similar path and focuses on the downscaling of the existing SARA exposure model that was proposed for the residential building stock of the communes of Valparaiso and Vina del Mar (Chile). Although this model allowed great progress in harmonising building classes and characterising their differential physical vulnerabilities, it is now outdated, and in any case, it is spatially aggregated over large administrative units. Hence, to more accurately consider the impact of future earthquakes on these cities, it is necessary to employ more reliable exposure models. For such a purpose, we propose updating this existing model through a Bayesian approach by integrating ancillary data that has been made increasingly available from Volunteering Geo-Information (VGI) activities. Its spatial representation is also optimised in higher resolution aggregation units that avoid the inconvenience of having incomplete building-by-building footprints. A worst-case earthquake scenario is presented to calculate direct economic losses and highlight the degree of uncertainty imposed by exposure models in comparison with other parameters used to generate the seismic ground motions within a sensitivity analysis. This example study shows the great potential of using increasingly available VGI to update worldwide building exposure models as well as its importance in scenario-based seismic risk assessment. KW - exposure KW - buildings KW - Bayesian model KW - downscaling KW - OpenStreetMap KW - ground motion fields KW - sensitivity KW - earthquake KW - vulnerability KW - risk Y1 - 2022 U6 - https://doi.org/10.3390/ijgi11020113 SN - 2220-9964 VL - 11 IS - 2 PB - MDPI CY - Basel ER - TY - JOUR A1 - Hannemann, Katrin A1 - Eulenfeld, Tom A1 - Krüger, Frank A1 - Dahm, Torsten T1 - Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic JF - Geophysical journal international N2 - The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere. KW - Body waves KW - Seismic attenuation KW - Wave scattering and diffraction Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab493 SN - 0956-540X SN - 1365-246X VL - 229 IS - 2 SP - 948 EP - 961 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Ben Nsir, Siwar A1 - Jomaa, Seifeddine A1 - Yildirim, Umit A1 - Zhou, Xiangqian A1 - D'Oria, Marco A1 - Rode, Michael A1 - Khlifi, Slaheddine T1 - Assessment of climate change impact on discharge of the lakhmass catchment (Northwest Tunisia) JF - Water N2 - The Mediterranean region is increasingly recognized as a climate change hotspot but is highly underrepresented in hydrological climate change studies. This study aims to investigate the climate change effects on the hydrology of Lakhmass catchment in Tunisia. Lakhmass catchment is a part of the Medium Valley of Medjerda in northwestern Tunisia that drains an area of 126 km(2). First, the Hydrologiska Byrans Vattenbalansavdelning light (HBV-light) model was calibrated and validated successfully at a daily time step to simulate discharge during the 1981-1986 period. The Nash Sutcliffe Efficiency and Percent bias (NSE, PBIAS) were (0.80, +2.0%) and (0.53, -9.5%) for calibration (September 1982-August 1984) and validation (September 1984-August 1986) periods, respectively. Second, HBV-light model was considered as a predictive tool to simulate discharge in a baseline period (1981-2009) and future projections using data (precipitation and temperature) from thirteen combinations of General Circulation Models (GCMs) and Regional Climatic Models (RCMs). We used two trajectories of Representative Concentration Pathways, RCP4.5 and RCP8.5, suggested by the Intergovernmental Panel on Climate Change (IPCC). Each RCP is divided into three projection periods: near-term (2010-2039), mid-term (2040-2069) and long-term (2070-2099). For both scenarios, a decrease in precipitation and discharge will be expected with an increase in air temperature and a reduction in precipitation with almost 5% for every +1 degrees C of global warming. By long-term (2070-2099) projection period, results suggested an increase in temperature with about 2.7 degrees C and 4 degrees C, and a decrease in precipitation of approximately 7.5% and 15% under RCP4.5 and RCP8.5, respectively. This will likely result in a reduction of discharge of 12.5% and 36.6% under RCP4.5 and RCP8.5, respectively. This situation calls for early climate change adaptation measures under a participatory approach, including multiple stakeholders and water users. KW - hydrological modeling KW - HBV-light model KW - Mediterranean KW - discharge KW - climate change KW - RCP4,5 and 8,5 Y1 - 2022 U6 - https://doi.org/10.3390/w14142242 SN - 2073-4441 VL - 14 IS - 14 PB - MDPI CY - Basel ER - TY - JOUR A1 - McCool, Weston C. A1 - Codding, Brian F. A1 - Vernon, Kenneth B. A1 - Wilson, Kurt M. A1 - Yaworsky, Peter M. A1 - Marwan, Norbert A1 - Kennett, Douglas J. T1 - Climate change-induced population pressure drives high rates of lethal violence in the Prehispanic central Andes JF - Proceedings of the National Academy of Sciences of the United States of America : PNAS N2 - Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of-and the interaction between-climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity-whether driven by reduced resource abundance or increased competition-can lead to violence in subsistence societies when the outcome is lower per capita resource availability. KW - climate change KW - population pressure KW - warfare KW - lethal violence KW - Andes Y1 - 2022 U6 - https://doi.org/10.1073/pnas.2117556119 SN - 0027-8424 SN - 1091-6490 VL - 119 IS - 17 PB - National Acad. of Sciences CY - Washington ER - TY - JOUR A1 - Yuan, Xiaoping P. A1 - Jiao, Ruohong A1 - Dupont-Nivet, Guillaume A1 - Shen, Xiaoming T1 - Southeastern Tibetan Plateau growth revealed by inverse analysis of landscape evolution model JF - Geophysical research letters N2 - The Cenozoic history of the Tibetan Plateau topography is critical for understanding the evolution of the Indian-Eurasian collision, climate, and biodiversity. However, the long-term growth and landscape evolution of the Tibetan Plateau remain ambiguous, it remains unclear if plateau uplift occurred soon after the India-Asia collision in the Paleogene (similar to 50-25 Ma) or later in the Neogene (similar to 20-5 Ma). Here, we reproduce the uplift history of the southeastern Tibetan Plateau using a 2D landscape evolution model, which simultaneously solves fluvial erosion and sediment transport processes in the drainage basins of the Three Rivers region (Yangtze, Mekong, and Salween Rivers). Our model was optimized through a formal inverse analysis with 20,000 forward simulations, which aims to reconcile the transient states of the present-day river profiles. The results, compared to existing paleoelevation and thermochronologic data, suggest initially low elevations (similar to 300-500 m) during the Paleogene, followed by a gradual southeastward propagation of topographic uplift of the plateau margin. KW - Tibetan Plateau KW - landscape evolution KW - fluvial erosion KW - inverse analysis KW - mountain growth KW - propagating uplift Y1 - 2022 U6 - https://doi.org/10.1029/2021GL097623 SN - 0094-8276 SN - 1944-8007 VL - 49 IS - 10 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Drozdov, Alexander A1 - Allison, Hayley J. A1 - Shprits, Yuri A1 - Usanova, Maria E. A1 - Saikin, Anthony A1 - Wang, Dedong T1 - Depletions of Multi-MeV Electrons and their association to Minima in Phase Space Density JF - Geophysical research letters N2 - Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations. KW - radiation belts KW - EMIC KW - VERB KW - PSD Y1 - 2022 U6 - https://doi.org/10.1029/2021GL097620 SN - 0094-8276 SN - 1944-8007 VL - 49 IS - 8 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Dahm, Torsten A1 - Heimann, Sebastian A1 - Metz, Malte A1 - Isken, Marius Paul T1 - A self-similar dynamic rupture model based on the simplified wave-rupture analogy JF - Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society N2 - The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions. KW - Earthquake dynamics KW - Earthquake ground motions KW - Earthquake hazards KW - Earthquake source observations Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab045 SN - 0956-540X SN - 1365-246X VL - 225 IS - 3 SP - 1586 EP - 1604 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Krstulović, Marija A1 - Rosa, Angelika D. A1 - Biedermann, Nicole A1 - Irifune, Tetsuo A1 - Wilke, Max T1 - Structural changes in aluminosilicate glasses up to 164 GPa and the role of alkali, alkaline earth cations and alumina in the densification mechanism JF - Chemical geology : official journal of the European Association for Geochemistry N2 - Pressure induced structural changes in silicate melts have a great impact on their physico-chemical properties and hence on their behaviour in the deep Earth's interior. In order to gain a deeper understanding we have studied the densification mechanism in multicomponent aluminosilicate glasses (albitic and albit-diopside composition) by means of extended X-ray absorption fine structure spectroscopy coupled to a diamond anvil cell up to 164 GPa. We have monitored the structural modifications from the network-former Ge as well as the network-modifier Sr. Notably, we tracked the evolution of Ge-O and Sr-O bond lengths (RGe-O, RSr-O) and their coordination number with pressure. We show that RGe-O increases strongly up to about 32 GPa, whereas RSr-O increases only slightly up to similar to 26 GPa. We assign these extensions to the increase of the coordination number from 4 to 6 (Ge) and from similar to 6 to at least 9 (Sr). Upon further compression RGe-O and RSr-O exhibit a continuous decrease to the highest probed pressure. These bond contractions, notably of RGe-O, that are continuous and exceed the one observed in pure SiO2 and GeO2, reflect a higher structural flexibility of multi-component glasses compared to those simple systems. Particularly, the high fraction of non-bridging oxygen atoms due to the presence of Na, Sr, Ca, Mg in the studied glasses, favours the simple compression of the highly-coordinated polyhedra of Si and Ge at pressure greater than 30 GPa. This is in strong contrast to pure oxides where cation polyhedral distortions govern the densification mechanism of the glass. The results of this study demonstrate that low field-strength alkali and alkaline earth cations, ubiquitous in deep Earth's melts, have a profound influence on the densification mechanism of glasses. Our results provide important constrains for interpreting the observed low velocity anomalies at the Earth's core-mantle boundary that have been, beyond others, referred to the presence of high-density melts. The hypothesis that non-buoyant melts at the Earth's core-mantle boundary can be formed by peculiar structural transformations in melts leading to higher coordination numbers compared to their crystalline equivalents is not supported from the present observations. The present results rather suggest that if velocity anomalies are to be explained by melts, these likely have considerable differences in chemical composition to the surrounding crystalline phase assemblage. Y1 - 2020 U6 - https://doi.org/10.1016/j.chemgeo.2020.119980 SN - 0009-2541 SN - 1872-6836 VL - 560 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Zuhr, Alexandra M. A1 - Dolman, Andrew M. A1 - Ho, Sze Ling A1 - Groeneveld, Jeroen A1 - Loewemark, Ludvig A1 - Grotheer, Hendrik A1 - Su, Chih-Chieh A1 - Laepple, Thomas T1 - Age-heterogeneity in marine sediments revealed by three-dimensional high-resolution radiocarbon measurements JF - Frontiers in Earth Science N2 - Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from. However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions. To shed light on 1) the age-depth relationship and its full uncertainty, 2) the magnitude of mixing processes affecting the downcore proxy variations, and 3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle. Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth. Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met. Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples. We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records. KW - paleoceanography KW - radiocarbon KW - age-heterogeneity KW - marine sediments KW - planktonic foraminifera KW - bioturbation KW - agemodeling KW - South China Sea Y1 - 2022 U6 - https://doi.org/10.3389/feart.2022.871902 SN - 2296-6463 VL - 10 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Zhu, Chuanbin A1 - Cotton, Fabrice A1 - Kwak, Dong-Youp A1 - Ji, Kun A1 - Kawase, Hiroshi A1 - Pilz, Marco T1 - Within-site variability in earthquake site response JF - Geophysical journal international N2 - The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses. In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings. We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models. Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites. Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s). However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models. KW - earthquake ground motions KW - earthquake hazards KW - site effects Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab481 SN - 0956-540X SN - 1365-246X VL - 229 IS - 2 SP - 1268 EP - 1281 PB - Oxford Univ. Press CY - Oxford ER - TY - CHAP A1 - Simon, François-Xavier A1 - Papadopoulos, Nikos A1 - Guillemoteau, Julien A1 - Oikonomou, Dimitris A1 - Simirdanis, Kleanthis T1 - Multi-frequency loop electromagnetic system measurement on shallow offshore archaeological site of Oulos T2 - ArcheoSciences : revue d'archéométrie / Groupe des Méthodes Pluridisciplinaires Contribuant à l'Archéologie (GMPCA) KW - hallow offshore KW - multi-frequency KW - electromagnetic KW - modelling KW - case study Y1 - 2021 U6 - https://doi.org/10.4000/archeosciences.9690 SN - 1960-1360 SN - 2104-3728 VL - 45 IS - 1 SP - 215 EP - 218 PB - Presses Universitaires de Rennes CY - Rennes ER -