Filtern
Volltext vorhanden
- nein (240) (entfernen)
Erscheinungsjahr
- 2019 (240) (entfernen)
Dokumenttyp
- Wissenschaftlicher Artikel (216)
- Sonstiges (9)
- Rezension (7)
- Dissertation (6)
- Habilitation (2)
Gehört zur Bibliographie
- ja (240)
Schlagworte
- climate (4)
- Holocene (3)
- paleoclimate (3)
- remote sensing (3)
- AMOC (2)
- Asia (2)
- Central Andes (2)
- Climate (2)
- Climate change (2)
- Climate change adaptation (2)
Institut
- Institut für Geowissenschaften (240) (entfernen)
Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds.
The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
Changes in the steepness of river profiles or abrupt vertical steps (i.e. waterfalls) are thought to be indicative of changes in erosion rates, lithology or other factors that affect landscape evolution. These changes are referred to as knickpoints or knickzones and are pervasive in bedrock river systems. Such features are thought to reveal information about landscape evolution and patterns of erosion, and therefore their locations are often reported in the geomorphic literature. It is imperative that studies reporting knickpoints and knickzones use a reproducible method of quantifying their locations, as their number and spatial distribution play an important role in interpreting tectonically active landscapes. In this contribution we introduce a reproducible knickpoint and knickzone extraction algorithm that uses river profiles transformed by integrating drainage area along channel length (the so-called integral or chi method). The profile is then statistically segmented and the differing slopes and step changes in the elevations of these segments are used to identify knickpoints, knickzones and their relative magnitudes. The output locations of identified knickpoints and knickzones compare favourably with human mapping: we test the method on Santa Cruz Island, CA, using previously reported knickzones and also test the method against a new dataset from the Quadrilatero Ferrifero in Brazil. The algorithm allows for the extraction of varying knickpoint morphologies, including stepped, positive slope-break (concave upward) and negative slope-break knickpoints. We identify parameters that most affect the resulting knickpoint and knickzone locations and provide guidance for both usage and outputs of the method to produce reproducible knickpoint datasets.
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
We present a new algorithm for solving the common problem of flow trapped in closed depressions within digital elevation models, as encountered in many applications relying on flow routing. Unlike other approaches (e.g., the Priority-Flood depression filling algorithm), this solution is based on the explicit computation of the flow paths both within and across the depressions through the construction of a graph connecting together all adjacent drainage basins. Although this represents many operations, a linear time complexity can be reached for the whole computation, making it very efficient. Compared to the most optimized solutions proposed so far, we show that this algorithm of flow path enforcement yields the best performance when used in landscape evolution models. In addition to its efficiency, our proposed method also has the advantage of letting the user choose among different strategies of flow path enforcement within the depressions (i.e., filling vs. carving). Furthermore, the computed graph of basins is a generic structure that has the potential to be reused for solving other problems as well, such as the simulation of erosion. This sequential algorithm may be helpful for those who need to, e.g., process digital elevation models of moderate size on single computers or run batches of simulations as part of an inference study.
Vermetid reefs and rocky shores are hot spots of biodiversity, often referred to as the subtropical equivalent of coral reefs. The development of the ecosystem depends on the activity of several reef builders, including red crustose coralline algae (CCA) such as Neogoniolithon brassica-florida. Despite its importance, little is known about Neogoniolithon sp. acclimation to rapid changes in light intensity and corresponding photosynthetic activity. To overcome the large spatial variability in the light field (due to location and the porous nature of the rocks) we grew Neogoniolithon sp. on glass slides and characterized its photosynthetic performance in response to various light intensities by following O-2 exchange and fluorescence parameters. This was also performed on rock-inhabiting thalli collected from the east Mediterranean basin. Generally, maximal photosynthetic rate was reached when Neogoniolithon sp. thalli grown under low illumination (such as in protected niches where the light intensity can be as low as 1% of surface illumination) were examined. When exposed to light intensities higher than those experienced during growth, Neogoniolithon sp. activates adaptive/protective mechanisms such as state transition and nonphotochemical fluorescence quenching and increases the dark respiration thereafter. We find that the Fv/Fm parameter (variable/maximal fluorescence) is not suitable to assess photosynthetic performance in Neogoniolithon sp. and propose using instead an alternative parameter recently developed. Our findings help to clarify why Neogoniolithon sp. is usually observed in shaded niches along the reef surfaces.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
In Magnetotellurics (MT) natural electromagnetic field variations are recorded to study the electrical conductivity structure of the subsurface. Thereby long time-series of electromagnetic data are subdivided into smaller segments, which are Fourier transformed and typically averaged in a statistically robust manner to obtain MT transfer functions. Unfortunately, nowadays the presence of man-made electromagnetic noise sources often deteriorates a significant fraction of the recorded time-series by overprinting the desired natural field variations. Available approaches to obtain undisturbed and high quality MT results include, for example robust statistics, remote reference or multi-station analyses which aim at the removal of outliers or uncorrelated noise. However, we have observed that intermittent noise often affects a certain time span resulting in a second cluster of transfer functions in addition to the expected true MT distribution. In this paper, we present a novel criterion for the detection and pre-selection of EM noise in form of outliers or additional clusters based on a distance measure of each data segment with regard to the centre of the data distribution. For this purpose, we utilize the Mahalanobis distance (MD) which computes the distance between two multivariate points considering the covariance matrix of the data that quantifies the shape and the size of multivariate data distributions. As the MD considers the covariance matrix, it corrects not only for different variances but also for any correlation between the data. The computation of both, the mean value and covariance matrix, is susceptible to ouliers (e.g. noise) and requires a statistically robust estimation. We tested several robust estimators, for example median absolute deviation or minimum covariance determinant algorithm and finally implemented an automatic criterion using a deterministic minimum covariance determinant algorithm. We will present results using MT data from various field experiments all over the world, which illustrate successfull data improvement. This approach is able to remove scattered data points as well as to reject complete data cluster originating from noise sources. However, like all purely statistical algorithms the criterion is limited to cases where the majority of the recorded data is well-behaved, that is noise content is below 50 per cent. If the majority of data points originates from noise sources, the new criterion will fail if used in an automatic way. In these cases, additional input by the user either manually or in an automated fashion can be utilized. We therefore suggest to use an add-on criterion to back the MD selection and subsequent robust stacking in form of a physically motivated constraint based on the magnetic incidence direction. This property indicates whether the magnetic field originates from various sources in the far field or from a strong and well defined source in the near field.
The intangible impacts of floods on welfare are not well investigated, even though they are important aspects of welfare. Moreover, flooding has gender based impacts on welfare. These differing impacts create a gender based flood risk resilience gap. We study the intangible impacts of flood risk on the subjective well-being of residents in central Vietnam. The measurement of intangible impacts through subjective well-being is a growing field within flood risk research. We find an initial drop in welfare through subjective well-being across genders when a flood is experienced. Male respondents tended to recover their welfare losses by around 80% within 5 years while female respondents were associated with a welfare recovery of around 70%. A monetization of the impacts floods have on an individual’s subjective well-being shows that for the average female respondent, between 41% to 86% of annual income would be required to compensate subjective well-being losses after 5 years of experiencing a flood. The corresponding value for males is 30% to 57% of annual income. This shows that the intangible impacts of flood risk are important (across genders) and need to be integrated into flood (or climate) risk assessments to develop more socially appropriate risk management strategies.
An Overview of Using Weather Radar for Climatological Studies: Successes, Challenges, and Potential
(2019)
Weather radars have been widely used to detect and quantify precipitation and nowcast severe weather for more than 50 years. Operational weather radars generate huge three-dimensional datasets that can accumulate to terabytes per day. So it is essential to review what can be done with existing vast amounts of data, and how we should manage the present datasets for the future climatologists. All weather radars provide the reflectivity factor, and this is the main parameter to be archived. Saving reflectivity as volumetric data in the original spherical coordinates allows for studies of the three-dimensional structure of precipitation, which can be applied to understand a number of processes, for example, analyzing hail or thunderstorm modes. Doppler velocity and polarimetric moments also have numerous applications for climate studies, for example, quality improvement of reflectivity and rain rate retrievals, and for interrogating microphysical and dynamical processes. However, observational data alone are not useful if they are not accompanied by sufficient metadata. Since the lifetime of a radar ranges between 10 and 20 years, instruments are typically replaced or upgraded during climatologically relevant time periods. As a result, present metadata often do not apply to past data. This paper outlines the work of the Radar Task Team set by the Atmospheric Observation Panel for Climate (AOPC) and summarizes results from a recent survey on the existence and availability of long time series. We also provide recommendations for archiving current and future data and examples of climatological studies in which radar data have already been used.
Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high-resolution seismicity catalogs framing the occurrence of 20 M-L > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events.
We investigate the inclusions hosted in peritectic garnet from metapelitic migmatites of the Kinzigite Formation (Ivrea Zone, NW Italy) to evaluate the starting composition of the anatectic melt and fluid regime during anatexis throughout the upper amphibolite facies, transition, and granulite facies zones. Inclusions have negative crystal shapes, sizes from 2 to 10 mu m and are regularly distributed in the core of the garnet. Microstructural and micro-Raman investigations indicate the presence of two types of inclusions: crystallized silicate melt inclusions (i.e., nanogranitoids, NI), and fluid inclusions (FI). Microstructural evidence suggests that FI and NI coexist in the same cluster and are primary (i.e., were trapped simultaneously during garnet growth). FI have similar compositions in the three zones and comprise variable proportions of CO2, CH4, and N-2, commonly with siderite, pyrophyllite, and kaolinite, suggesting a COHN composition of the trapped fluid. The mineral assemblage in the NI contains K-feldspar, plagioclase, quartz, biotite, muscovite, chlorite, graphite and, rarely, calcite. Polymorphs such as kumdykolite, cristobalite, tridymite, and less commonly kokchetavite, were also found. Rehomogenized NI from the different zones show that all the melts are leucogranitic but have slightly different compositions. In samples from the upper amphibolite facies, melts are less mafic (FeO + MgO = 2.0-3.4 wt%), contain 860-1700 ppm CO2 and reach the highest H2O contents (6.5-10 wt%). In the transition zone melts have intermediate H2O (4.8-8.5 wt%), CO2 (457-1534 ppm) and maficity (FeO + MgO = 2.3-3.9 wt%). In contrast, melts at granulite facies reach highest CaO, FeO + MgO (3.2-4.7 wt%), and CO2 (up to 2,400 ppm), with H2O contents comparable (5.4-8.3 wt%) to the other two zones. Our results represent the first clear evidence for carbonic fluid-present melting in the Ivrea Zone. Anatexis of metapelites occurred through muscovite and biotite breakdown melting in the presence of a COH fluid, in a situation of fluid-melt immiscibility. The fluid is assumed to have been internally derived, produced initially by devolatilization of hydrous silicates in the graphitic protolith, then as a result of oxidation of carbon by consumption of Fe3+-bearing biotite during melting. Variations in the compositions of the melts are interpreted to result from higher T of melting. The H2O contents of the melts throughout the three zones are higher than usually assumed for initial H2O contents of anatectic melts. The CO2 contents are highest at granulite facies, and show that carbon-contents of crustal magmas are not negligible at high T. The activity of H2O of the fluid dissolved in granitic melts decreases with increasing metamorphic grade. Carbonic fluid-present melting of the deep continental crust represents, together with hydrate-breakdown melting reactions, an important process in the origin of crustal anatectic granitoids.
Arboreal epiphytes (plants residing in forest canopies) are present across all major climate zones and play important roles in forest biogeochemistry. The substantial water storage capacity per unit area of the epiphyte "bucket" is a key attribute underlying their capability to influence forest hydrological processes and their related mass and energy flows. It is commonly assumed that the epiphyte bucket remains saturated, or near-saturated, most of the time; thus, epiphytes (particularly vascular epiphytes) can store little precipitation, limiting their impact on the forest canopy water budget. We present evidence that contradicts this common assumption from (i) an examination of past research; (ii) new datasets on vascular epiphyte and epi-soil water relations at a tropical montane cloud forest (Monteverde, Costa Rica); and (iii) a global evaluation of non-vascular epiphyte saturation state using a process-based vegetation model, LiBry. All analyses found that the external and internal water storage capacity of epiphyte communities is highly dynamic and frequently available to intercept precipitation. Globally, non-vascular epiphytes spend <20% of their time near saturation and regionally, including the humid tropics, model results found that non-vascular epiphytes spend similar to 1/3 of their time in the dry state (0-10% of water storage capacity). Even data from Costa Rican cloud forest sites found the epiphyte community was saturated only 1/3 of the time and that internal leaf water storage was temporally dynamic enough to aid in precipitation interception. Analysis of the epi-soils associated with epiphytes further revealed the extent to which the epiphyte bucket emptied-as even the canopy soils were often <50% saturated (29-53% of all days observed). Results clearly show that the epiphyte bucket is more dynamic than currently assumed, meriting further research on epiphyte roles in precipitation interception, redistribution to the surface and chemical composition of "net" precipitation waters reaching the surface.
The occurrence of mounds dominated by siliceous sponges and microbialites is often related to distal, deep settings of middle ramps and shelves. This paper presents evidence for Bajocian (Garanliana garantiana Zone) microbial-siliceous sponge mounds formed in open marine but relatively shallow settings of a ramp from the Iberian Basin of eastern Spain. Marked differences in mound spacing, morphology, and composition of the related intermound facies are observed from distal to more proximal settings. The distal (below storm wave base) settings are characterized by alternating tabular-bedded marls and limestones rich in pelagic fossils (ammonites, belemnites), open-marine thin-shelled bivalves (Bositra-like), as well as peloids, which include widely or randomly spaced isolated, small (up to 0.4 m high) and larger (up to 2.5 m high) mounds with upward accretion. The intermediate (near to above storm wave base) settings show tabular, thickened beds of peloidal and/or intraclastic limestones with closely spaced mounds (similar to 1 m high), which often coalesce laterally, forming extensive lenticular structures (up to 10 m wide). The proximal (above storm wave base) depositional settings consist of tabular to irregular beds of intraclastic limestones with widely spaced small (up to 0.4 m high) mounds with mainly tabular geometries. The mound framework contains variable proportions of microbialites (dense to clotted peloidal thrombolitic fabrics) and siliceous sponges (hexactinellids and lithistids in similar proportion) ranging from planar to conic shapes. These morphological and compositional changes allow characterizing three shallowing-upward sequences (sequences 1-3) developed in the overall regressive trend of a basin-wide, upper Bajocian T-R cycle. Episodic wave reworking of the early-cemented mounds resulted in the formation of peloids, small rounded intraclasts, and large, rounded or subangular intraclasts. These nonskeletal micritic grains show internal fabrics related to those of the mound and/or microbialites. A progressive textural gradation towards greater size and lesser roundness of the nonskeletal grains in the areas in the vicinity of the main mound factory is documented (i.e., from large, subangular intraclasts in the areas close to the main mound factory to peloids in the areas that are far from it). We discuss the alternative model of internal waves (instead of storm-induced waves) as the hydrodynamic agent providing the high-energy events needed to explain the origin of the peloidal-intraclastic intermound facies and, likely, also the nutrients needed by the microbialites and siliceous sponges to grow.
In Germany, the irrigation sector accounts for only 1% of water use. In recent years, however, this sector has attracted more attention due to the occurrence of severe drought periods. Irrigation scheduling systems could support adaptation strategies but little is known about current providers, performance and users. In this study we aimed to depict the current situation of the existence and functioning of irrigation scheduling systems available in Germany. Six methods were identified and assessed based on direct interviews with end-users and a comparative analysis. The results showed a positive feedback from the users. However, the recommendations were rarely implemented, while only the seasonal irrigation requirement was considered to support actual water abstraction. These results were corroborated by the comparative analysis. Five of the six irrigation scheduling systems estimated the seasonal irrigation amount consistently, while wider differences were found by looking at the irrigation season and at the number of irrigations. Overall, it is found that irrigation support systems are valuable tools for supporting adaptation strategies to fast changes in agro-environmental conditions. However, specific assessments based on real measurements should be considered in order to improve the performance of the systems and provide more consistent support to end-users. (c) 2019 John Wiley & Sons, Ltd.
High spectral resolution (hyperspectral) remote sensing has already demonstrated its capabilities for soil constituent mapping based on absorption feature parameters. This paper tests different parametrizations of the 1.75 μm gypsum feature for the determination of gypsum abundances, from the laboratory to remote sensing applications of recent as well as upcoming hyperspectral sensors. In particular, this study focuses on remote sensing imagery over the large body of the Omongwa pan located in the Namibian Kalahari. Four common absorption feature parameters are compared: band ratio through the introduction of the Normalized Differenced Gypsum Index (NDGI), the shape-based parameters Slope, and Half-Area, and the Continuum Removed Absorption Depth (CRAD). On laboratory soil samples from the pan, CRAD and NDGI approaches perform best to determine gypsum content tested in cross validated regression models with XRD mineralogical data (R² = 0.84 for NDGI and R² = 0.86 for CRAD). Subsequently the laboratory prediction functions are transferred to remote sensing imagery of spaceborne Hyperion, airborne HySpex and simulated spaceborne EnMAP sensor. Variable results were obtained depending on sensor characteristics, data quality, preprocessing and spectral parameters. Overall, the CRAD parameter in this wavelength region proved not to be robust for remote sensing applications, and the simple band ratio based parameter, the NDGI, proved robust and is recommended for future use for the determination of gypsum content in bare soils based on remote sensing hyperspectral imagery.
Water fluxes in highly impounded regions are heavily dependent on reservoir properties. However, for large and remote areas, this information is often unavailable. In this study, the geometry and volume of small surface reservoirs in the semi-arid region of Brazil were estimated using terrain and shape attributes extracted by remote sensing. Regression models and data classification were used to predict the volumes, at different water stages, of 312 reservoirs for which topographic information is available. The power function used to describe the reservoir shapes tends to overestimate the volumes; therefore, a modified shape equation was proposed. Among the methods tested, four were recommended based on performance and simplicity, for which the mean absolute percentage errors varied from 24 to 39%, in contrast to the 94% error achieved with the traditional method. Despite the challenge of precisely deriving the flooded areas of reservoirs, water management in highly reservoir-dense environments should benefit from volume prediction based on remote sensing.
As a consequence of the rapid growing worldwide seismic data set, a huge variety of automatized data-processing methods have been developed. To perform automatized waveform-based seismological studies aiming for magnitudes or source process inversion, it is crucial to identify network stations with erroneous transfer functions, gain factors, or component orientations. We developed a new tool dedicated to automated station quality control of dense seismic networks and arrays. The python-based AutoStatsQ toolbox uses the pyrocko seismic data-processing environment. The toolbox automatically downloads data and metadata for selected teleseismic events and performs different tests. As a result, relative gain factors, sensor orientation corrections, and reliable frequency bands are computed for all stations in a chosen time period. Relative gain factors are calculated for all stations and events in a time domain based on maximum P-phase amplitudes. A Rayleigh-wave polarization analysis is used to identify deviating sensor orientations. The power spectra of all stations in a given frequency range are compared with synthetic ones, accessing Global Centroid Moment Tensor (CMT) solutions. Frequency ranges of coinciding synthetic and recorded power spectral densities (PSDs) may serve as guidelines for choosing band-pass filters for moment tensor (MT) inversion and help confirm the corner frequency of the instrument. The toolbox was applied to the permanent and temporary AlpArray networks as well as to the denser SWATH-D network, a total of over 750 stations. Stations with significantly deviating gain factors were identified, as well as stations with inverse polarity and misorientations of the horizontal components. The tool can be used to quickly access network quality and to omit or correct stations before MT inversion. Electronic Supplement: List of teleseismic events and tables of median, mean, and standard deviation of relative gain factors, and figures of relative gain factors of all event-station pairs, waveform example showing inverse polarity of horizontal components on ZS.D125, histograms of median, mean, and standard deviation of the correction angles, examples of synthetic and recorded frequency spectra of ZS.D046 and NI.VINO.
Quantifying erosion rates, and how they compare to rock uplift rates, is fundamental for understanding landscape response to tectonics and associated sediment fluxes from upland areas. The erosional response to uplift is well-represented by river incision and the associated landslide activity. However, characterising the relationship between these processes remains a major challenge in tectonically active areas, in some cases because landslides can preclude obtaining reliable erosion rates from cosmogenic radionuclide (CRN) concentrations. Here, we quantify the control of tectonics and its coupled geomorphic response on the erosion rates of catchments in southern Italy that are experiencing a transient response to normal faulting. We analyse in-situ Be-10 concentrations for detrital sediment samples, collected along the strike of faults with excellent tectonic constraints and landslide inventories. We demonstrate that Be-10-derived erosion rates are controlled by fault throw rates and the extent of transient incision and associated landsliding in the catchments. We show that the low-relief sub-catchments above knickpoints erode at uniform background rates of similar to 0.10 mm/yr, while downstream of knickpoints, erosion removes similar to 50% of the rock uplifted by the faults, at rates of 0.10-0.64 mm/yr. Despite widespread landsliding, CRN samples provide relatively consistent and accurate erosion rates, most likely because landslides are frequent, small, and shallow, and represent the integrated record of landsliding over several seismic cycles. Consequently, we combine these validated Be-10 erosion rates and data from a geomorphological landslide inventory in a published numerical model, to gain further insight into the long-term landslide rates and sediment mixing, highlighting the potential of CRN data to study landslide dynamics. (C) 2018 Elsevier B.V. All rights reserved.
Breakup Without Borders
(2019)
Relative plate motions during continental rifting result from the interplay of local with far-field forces. Here we study the dynamics of rifting and breakup using large-scale numerical simulations of mantle convection with self-consistent evolution of plate boundaries. We show that continental separation follows a characteristic evolution with four distinctive phases: (1) an initial slow rifting phase with low divergence velocities and maximum tensional stresses, (2) a synrift speed-up phase featuring an abrupt increase of extension rate with a simultaneous drop of tensional stress, (3) the breakup phase with inception of fast sea-floor spreading, and (4) a deceleration phase occurring in most but not all models where extensional velocities decrease. We find that the speed-up during rifting is compensated by subduction acceleration or subduction initiation even in distant localities. Our study illustrates new links between local rift dynamics, plate motions, and subduction kinematics during times of continental separation.
Convergence between the Indian and Asian plates has reshaped large parts of Asia, changing regional climate and biodiversity, yet geodynamic models fundamentally diverge on how convergence was accommodated since the India-Asia collision. Here we report palaeomagnetic data from the Burma Terrane, which is at the eastern edge of the collision zone and is famous for its Cretaceous amber biota, to better determine the evolution of the India-Asia collision. The Burma Terrane was part of a Trans-Tethyan island arc and stood at a near-equatorial southern latitude at similar to 95 Ma, suggesting island endemism for the Burmese amber biota. The Burma Terrane underwent significant clockwise rotation between similar to 80 and 50 Ma, causing its subduction margin to become hyper-oblique. Subsequently, it was translated northward on the Indian Plate by an exceptional distance of at least 2,000 km along a dextral strike-slip fault system in the east. Our reconstructions are only compatible with geodynamic models involving an initial collision of India with a near-equatorial Trans-Tethyan subduction system at similar to 60 Ma, followed by a later collision with the Asian margin.
A proper assessment of seismic reference site conditions has important applications as they represent the basis on which ground motions and amplifications are generally computed. Besides accounting for the average S-wave velocity over the uppermost 30 m (V-S30), the parameterization of high-frequency ground motions beyond source-corner frequency received significant attention. kappa, an empirical parameter introduced by Anderson and Hough (1984), is often used to represent the spectral decay of the acceleration spectrum at high frequencies. The lack of hard-rock records and the poor understanding of the physics of kappa introduced significant epistemic uncertainty in the final seismic hazard of recent projects. Thus, determining precise and accurate regional hard-rock kappa(0) values is critical. We propose an alternative procedure for capturing the reference kappa(0) on regional scales by linking thewell-known high-frequency attenuation parameter kappa and the properties of multiple-scattered coda waves. Using geological and geophysical data around more than 1300 stations for separating reference and soft soil sites and based on more than 10,000 crustal earthquake recordings, we observe that kappa(0) from multiple-scattered coda waves seems to be independent of the soil type but correlated with the hard-rock kappa(0), showing significant regional variations across Europe. The values range between 0.004 s for northern Europe and 0.020 s for the southern and southeastern parts. On the other hand, measuring kappa (and correspondingly kappa(0)) on the S-wave window (as classically proposed), the results are strongly affected by transmitted (reflected, refracted, and scattered) waves included in the analyzed window biasing the proper assessment of kappa(0). This effect is more pronounced for soft soil sites. In this way, kappa(coda)(0) can serve as a proxy for the regional hard-rock kappa(0) at the reference sites.
A wide variety of processes controls the time of occurrence, duration, extent, and severity of river floods. Classifying flood events by their causative processes may assist in enhancing the accuracy of local and regional flood frequency estimates and support the detection and interpretation of any changes in flood occurrence and magnitudes. This paper provides a critical review of existing causative classifications of instrumental and preinstrumental series of flood events, discusses their validity and applications, and identifies opportunities for moving toward more comprehensive approaches. So far no unified definition of causative mechanisms of flood events exists. Existing frameworks for classification of instrumental and preinstrumental series of flood events adopt different perspectives: hydroclimatic (large-scale circulation patterns and atmospheric state at the time of the event), hydrological (catchment scale precipitation patterns and antecedent catchment state), and hydrograph-based (indirectly considering generating mechanisms through their effects on hydrograph characteristics). All of these approaches intend to capture the flood generating mechanisms and are useful for characterizing the flood processes at various spatial and temporal scales. However, uncertainty analyses with respect to indicators, classification methods, and data to assess the robustness of the classification are rarely performed which limits the transferability across different geographic regions. It is argued that more rigorous testing is needed. There are opportunities for extending classification methods to include indicators of space-time dynamics of rainfall, antecedent wetness, and routing effects, which will make the classification schemes even more useful for understanding and estimating floods. This article is categorized under: Science of Water > Water Extremes Science of Water > Hydrological Processes Science of Water > Methods
The establishment and evolution of the Asian monsoons and arid interior have been linked to uplift of the Tibetan Plateau, retreat of the inland proto-Paratethys Sea and global cooling during the Cenozoic. However, the respective role of these driving mechanisms remains poorly constrained. This is partly due to a lack of continental records covering the key Eocene epoch marked by the onset of Tibetan Plateau uplift, proto-Paratethys Sea incursions and long-term global cooling. In this study, we reconstruct paleoenvironments in the Xining Basin, NE Tibet, to show a long-term drying of the Asian continental interior from the early Eocene to the Oligocene. Superimposed on this trend are three alternations between arid mudflat and wetter saline lake intervals, which are interpreted to reflect atmospheric moisture fluctuations in the basin. We date these fluctuations using magnetostratigraphy and the radiometric age of an intercalated tuff layer. The first saline lake interval is tentatively constrained to the late Paleocene-early Eocene. The other two are firmly dated between similar to 46 Ma (top magnetochron C21n) and similar to 41 Ma (base C18r) and between similar to 40 Ma (base C18n) and similar to 37 Ma (top C17n). Remarkably, these phases correlate in time with highstands of the proto-Paratethys Sea. This strongly suggests that these sea incursions enhanced westerly moisture supply as far inland as the Xining Basin. We conclude that the proto-Paratethys Sea constituted a key driver of Asian climate and should be considered in model and proxy interpretations. (C) 2019 Elsevier B.V. All rights reserved.
Permafrost thaw subjects previously frozen soil organic carbon (SOC) to microbial degradation to the greenhouse gases carbon dioxide (CO2) and methane (CH4). Emission of these gases constitutes a positive feedback to climate warming. Among numerous uncertainties in estimating the strength of this permafrost carbon feedback (PCF), two are: (i) how mineralization of permafrost SOC thawed in saturated anaerobic conditions responds to changes in temperature and (ii) how microbial communities and temperature sensitivities change over time since thaw. To address these uncertainties, we utilized a thermokarst-lake sediment core as a natural chronosequence where SOC thawed and incubated in situ under saturated anaerobic conditions for up to 400 years following permafrost thaw. Initial microbial communities were characterized, and sediments were anaerobically incubated in the lab at four temperatures (0 °C, 3 °C, 10 °C, and 25 °C) bracketing those observed in the lake's talik. Net CH4 production in freshly-thawed sediments near the downward-expanding thaw boundary at the base of the talik were most sensitive to warming at the lower incubation temperatures (0 °C to 3 °C), while the overlying sediments which had been thawed for centuries had initial low abundant methanogenic communities (< 0.02%) and did not experience statistically significant increases in net CH4 production potentials until higher incubation temperatures (10 °C to 25 °C). We propose these observed differences in temperature sensitivities are due to differences in SOM quality and functional microbial community composition that evolve over time; however further research is necessary to better constrain the roles of these factors in determining temperature controls on anaerobic C mineralization.
Britholite group minerals (REE,Ca)(5)[(Si,P)O-4](3)(OH,F) are widespread rare-earth minerals in alkaline rocks and their associated metasomatic zones, where they usually are minor accessory phases. An exception is the REE deposit Rodeo de los Molles, Central Argentina, where fluorbritholite-(Ce) (FBri) is the main carrier of REE and is closely intergrown with fluorapatite (FAp). These minerals reach an abundance of locally up to 75 modal% (FBri) and 20 modal% (FAp) in the vein mineralizations. The Rodeo de los Molles deposit is hosted by a fenitized monzogranite of the Middle Devonian Las Chacras-Potrerillos batholith. The REE mineralization consists of fluorbritholite-(Ce), britholite-(Ce), fluorapatite, allanite-(Ce), and REE fluorcarbonates, and is associated with hydrothermal fluorite, quartz, albite, zircon, and titanite. The REE assemblage takes two forms: irregular patchy shaped REE-rich composites and discrete cross-cutting veins. The irregular composites are more common, but here fluorbritholite-(Ce) is mostly replaced by REE carbonates. The vein mineralization has more abundant and better-preserved britholite phases. The majority of britholite grains at Rodeo de los Molles are hydrothermally altered, and alteration is strongly enhanced by metamictization, which is indicated by darkening of the mineral, loss of birefringence, porosity, and volume changes leading to polygonal cracks in and around altered grains. A detailed electron microprobe study of apatite-britholite minerals from Rodeo de los Molles revealed compositional variations in fluorapatite and fluorbritholite-(Ce) consistent with the coupled substitution of REE3+ + Si4+ = Ca2+ + P5+ and a compositional gap of similar to 4 apfu between the two phases, which we interpret as a miscibility gap. Micrometer-scale intergrowths of fluorapatite in fluorbritholite-(Ce) minerals and vice versa are chemically characterized here for the first time and interpreted as exsolution textures that formed during cooling below the proposed solvus.
The Last Interglacial (Eemian, MIS 5e) can be considered a test-bed for climate dynamics under a warmer-than-present climate. In this study we present a chironomid record from the high latitude Sokli site (N Finland), where a long continuous sediment sequence from the last interglacial has been preserved from glacial erosion. The chironomid-analysis shows a diverse fauna, with dominance of warm-water indicators and shifts in assemblage composition that can be attributed to temperature, lake depth, productivity and habitat availability. Quantitative mean July paleotemperature estimates based on the chironomid data indicate overall mean July air temperatures up to 1 degrees C warmer than present. Two cooling events can be discerned, the Tunturi event, dated to about 127.5kaBP, in the lower part of the sequence, and the Varrio event, dated to about 119kaBP, associated with the beginning of a cooling trend in the upper part of the record. Warm conditions already at the onset of the interglacial contrast with a recent chironomid-based last interglacial temperature reconstruction from Denmark, which suggests a late onset of Eemian warming. The relatively small increase in inferred temperatures compared to present day temperatures at Sokli differs from other high latitude Eemian sites, and likely reflects the influence of the Atlantic Meridional Overturning Circulation in maintaining already elevated temperatures in Fennoscandia during interglacials.
The Olorgesailie Drilling Project and the related Hominin Sites and Paleolakes Drilling Project in East Africa were initiated to test hypotheses and models linking environmental change to hominin evolution by drilling lake basin sediments adjacent to important archeological and paleoanthropological sites. Drill core OL012-1A recovered 139 m of sedimentary and volcaniclastic strata from the Koora paleolake basin, southern Kenya Rift, providing the opportunity to compare paleoenvironmental influences over the past million years with the parallel record exposed at the nearby Olorgesailie archeological site. To refine our ability to link core-to-outcrop paleoenvironmental records, we institute here a methodological framework for deriving a robust age model for the complex lithostratigraphy of OL012-1A. Firstly, chronostratigraphic control points for the core were established based on 4 Ar/39Ar ages from intercalated tephra deposits and a basal trachyte flow, as well as the stratigraphic position of the Brunhes-Matuyama geomagnetic reversal. This dataset was combined with the position and duration of paleosols, and analyzed using a new Bayesian algorithm for high-resolution age-depth modeling of hiatus-bearing stratigraphic sections. This model addresses three important aspects relevant to highly dynamic, nonlinear depositional environments: 1) correcting for variable rates of deposition, 2) accommodating hiatuses, and 3) quantifying realistic age uncertainty with centimetric resolution. Our method is applicable to typical depositional systems in extensional rifts as well as to drill cores from other dynamic terrestrial or aquatic environments. We use the core age model and lithostratigraphy to examine the inter connectivity of the Koora Basin to adjacent areas and sources of volcanism. (C) 2019 Elsevier Ltd. All rights reserved.
The Baringo-Tugen-Barsemoi 2013 drillcore (BTB13), acquired as part of the Hominin Sites and Paleolakes Drilling Project, recovered 228 m of fluviolacustrine sedimentary rocks and tuffs spanning a similar to 3.29-2.56 Ma interval of the highly fossiliferous and hominin-bearing Chemeron Formation, Tugen Hills, Kenya. Here we present a Bayesian stratigraphic age model for the core employing chronostratigraphic control points derived from Ar-40/Ar-39 dating of tuffs from core and outcrop, Ar-40/Ar-39 age calibration of related outcrop diatomaceous units, and core magnetostratigraphy. The age model reveals three main intervals with distinct sediment accumulation rates: an early rapid phase from 3.2 to 2.9 Ma; a relatively slow phase from 2.9 to 2.7 Ma; and the highest rate of accumulation from 2.7 to 2.6 Ma. The intervals of rapid accumulation correspond to periods of high Earth orbital eccentricity, whereas the slow accumulation interval corresponds to low eccentricity at 2.9-2.7 Ma, suggesting that astronomically mediated climate processes may be responsible for the observed changes in sediment accumulation rate. Lacustrine transgression-regression events, as delineated using sequence stratigraphy, dominantly operate on precession scale, particularly within the high eccentricity periods. A set of erosively based fluvial conglomerates correspond to the 2.9-2.7 Ma interval, which could be related to either the depositional response to low eccentricity or to the development of unconformities due to local tectonic activity. Age calibration of core magnetic susceptibility and gamma density logs indicates a close temporal correspondence between a shift from high- to low-frequency signal variability at similar to 3 Ma, approximately coincident the end of the mid-Piacenzian Warm Period, and the beginning of the cooling of world climate leading to the initiation of Northern Hemispheric glaciation c. 2.7 Ma. BTB13 and the Baringo Basin records may thus provide evidence of a connection between high-latitude glaciation and equatorial terrestrial climate toward the end of the Pliocene.
The Chew Bahir Drilling Project (CBDP) aims to test possible linkages between climate and evolution in Africa through the analysis of sediment cores that have recorded environmental changes in the Chew Bahir basin. In this statistical project we consider the Chew Bahir palaeolake to be a dynamical system consisting of interactions between its different components, such as the waterbody, the sediment beneath lake, and the organisms living within and around the lake. Recurrence is a common feature of such dynamical systems, with recurring patterns in the state of the system reflecting typical influences. Identifying and defining these influences contributes significantly to our understanding of the dynamics of the system. Different recurring changes in precipitation, evaporation, and wind speed in the Chew Bahir basin could result in similar (but not identical) conditions in the lake (e.g., depth and area of the lake, alkalinity and salinity of the lake water, species assemblages in the water body, and diagenesis in the sediments). Recurrence plots (RPs) are graphic displays of such recurring states within a system. Measures of complexity were subsequently introduced to complement the visual inspection of recurrence plots, and provide quantitative descriptions for use in recurrence quantification analysis (RQA). We present and discuss herein results from an RQA on the environmental record from six short (< 17 m) sediment cores collected during the CBDP, spanning the last 45 kyrs. The different types of variability and transitions in these records were classified to improve our understanding of the response of the biosphere to climate change, and especially the response of humans in the area.
The processes that control long term landscape evolution in continental interiors and, in particular, along passive margins such as in southern Africa, are still the subject of much debate (e.g. Braun, 2018). Although today the Namibian margin is characterized by an arid climate, it has experienced climatic fluctuations during the Cenozoic and, yet, to date no study has documented the potential role of climate on its erosion history. In western Namibia, the Brandberg Massif, an erosional remnant or inselberg, provides a good opportunity to document the Cenozoic denudation history of the margin using the relationship between rock cooling or exhumation ages and their elevation. Here we provide new apatite (UThSm)/He dates on the Brandberg Inselberg that range from 151 +/- 12 to 30 +/- 2 Ma. Combined with existing apatite fission track data, they yield new constraints on the denudation history of the margin. These data document two main cooling phases since continental break-up 130 Myr ago, a rapid one (similar to 10 degrees C/Myr) following break-up and a slower one (similar to 12 degrees C/Myr) between 65 and 35 Ma. We interpret them respectively to be related to escarpment erosion following rifting and continental break-up and as a phase of enhanced denudation during the Early Eocene Climatic Optimum. We propose that during the Early Eocene Climatic Optimum chemical weathering was important and contributed significantly to the denudation of the Namibian margin and the formation of a pediplain around the Brandberg and enhanced valley incision within the massif. Additionally, aridification of the region since 35 Ma has resulted in negligible denudation rates since that time. (C) 2019 Elsevier B.V. All rights reserved.
The structure and organization of river networks has been used for decades to investigate the influence of climate and tectonics on landscapes. The majority of these studies either analyze rivers in profile view by extracting channel steepness or calculate planform metrics such as drainage density. However, these techniques rely on the assumption of homogeneity: that intrinsic and external factors are spatially or temporally invariant over the measured profile. This assumption is violated for the majority of Earth's landscapes, where variations in uplift rate, rock strength, climate, and geomorphic process are almost ubiquitous. We propose a method for classifying river profiles to identify landscape regions with similar characteristics by adapting hierarchical clustering algorithms developed for time series data. We first test our clustering on two landscape evolution scenarios and find that we can successfully cluster regions with different erodibility and detect the transient response to sudden base level fall. We then test our method in two real landscapes: first in Bitterroot National Forest, Idaho, where we demonstrate that our method can detect transient incision waves and the topographic signature of fluvial and debris flow process regimes; and second, on Santa Cruz Island, California, where our technique identifies spatial patterns in lithology not detectable through normalized channel steepness analysis. By calculating channel steepness separately for each cluster, our method allows the extraction of more reliable steepness metrics than if calculated for the landscape as a whole. These examples demonstrate the method's ability to disentangle fluvial morphology in complex lithological and tectonic settings.
Debate persists concerning the timing and geodynamics of intercontinental collision, style of syncollisional deformation, and development of topography and fold-and-thrust belts along the >1,700-km-long Izmir-Ankara-Erzincan suture zone (IAESZ) in Turkey. Resolving this debate is a necessary precursor to evaluating the integrity of convergent margin models and kinematic, topographic, and biogeographic reconstructions of the Mediterranean domain. Geodynamic models argue either for a synchronous or diachronous collision during either the Late Cretaceous and/or Eocene, followed by Eocene slab breakoff and postcollisional magmatism. We investigate the collision chronology in western Anatolia as recorded in the sedimentary archives of the 90-km-long Saricakaya Basin perched at shallow structural levels along the IAESZ. Based on new zircon U-Pb geochronology and depositional environment and sedimentary provenance results, we demonstrate that the Saricakaya Basin is an Eocene sedimentary basin with sediment sourced from both the IAESZ and Sogut Thrust fault to the south and north, respectively, and formed primarily by flexural loading from north-south shortening along the syncollisional Sogut Thrust. Our results refine the timing of collision between the Anatolides and Pontide terranes in western Anatolia to Maastrichtian-Middle Paleocene and Early Eocene crustal shortening and basin formation. Furthermore, we demonstrate contemporaneous collision, deformation, and magmatism across the IAESZ, supporting synchronous collision models. We show that regional postcollisional magmatism can be explained by renewed underthrusting instead of slab breakoff. This new IAESZ chronology provides additional constraints for kinematic, geodynamic, and biogeographic reconstructions of the Mediterranean domain.
Encroachment of shrubs into the unique pastoral grassland ecosystems of the Tibetan Plateau has significant impact on ecosystem services, especially forage production. We developed a process-based ecohydrological model to identify the relative importance of the main drivers of shrub encroachment for the alpine meadows within the Qinghai province. Specifically, we explored the effects of summer livestock grazing (intensity and type of livestock) together with the effects of climate warming, including interactions between herbaceous and woody vegetation and feedback loops between soil, water and vegetation. Under current climatic conditions and a traditional herd composition, an increasing grazing intensity above a threshold value of 0.32 +/- 0.10 large stock units (LSU) ha(-1) day(-1) changes the vegetation composition from herbaceous towards a woody and bare soil dominated system. Very high grazing intensity (above 0.8 LSU ha(-1) day(-1)) leads to a complete loss of any vegetation. Under warmer conditions, the vegetation showed a higher resilience against livestock farming. This resilience is enhanced when the herd has a higher browser : grazer ratio. A cooler climate has a shrub encroaching effect, whereas warmer conditions increase the cover of the herbaceous vegetation. This effect was primarily due to season length and an accompanied competitive loss of slower growing shrubs, rather than evaporative water loss leading to less soil water in deeper soil layers for deeper rooting shrubs. If climate warming is driving current shrub encroachment, we conclude it is only indirectly so. It would be manifest by an advancing shrubline and could be regarded as a climatic escape of specific shrub species such as Potentilla fruticosa. Under the recent high intensity of grazing, only herding by more browsing animals can potentially prevent both shrub encroachment and the complete loss of herbaceous vegetation.
The Arctic-Boreal regions experience strong changes of air temperature and precipitation regimes, which affect the thermal state of the permafrost. This results in widespread permafrost-thaw disturbances, some unfolding slowly and over long periods, others occurring rapidly and abruptly. Despite optical remote sensing offering a variety of techniques to assess and monitor landscape changes, a persistent cloud cover decreases the amount of usable images considerably. However, combining data from multiple platforms promises to increase the number of images drastically. We therefore assess the comparability of Landsat-8 and Sentinel-2 imagery and the possibility to use both Landsat and Sentinel-2 images together in time series analyses, achieving a temporally-dense data coverage in Arctic-Boreal regions. We determined overlapping same-day acquisitions of Landsat-8 and Sentinel-2 images for three representative study sites in Eastern Siberia. We then compared the Landsat-8 and Sentinel-2 pixel-pairs, downscaled to 60 m, of corresponding bands and derived the ordinary least squares regression for every band combination. The acquired coefficients were used for spectral bandpass adjustment between the two sensors. The spectral band comparisons showed an overall good fit between Landsat-8 and Sentinel-2 images already. The ordinary least squares regression analyses underline the generally good spectral fit with intercept values between 0.0031 and 0.056 and slope values between 0.531 and 0.877. A spectral comparison after spectral bandpass adjustment of Sentinel-2 values to Landsat-8 shows a nearly perfect alignment between the same-day images. The spectral band adjustment succeeds in adjusting Sentinel-2 spectral values to Landsat-8 very well in Eastern Siberian Arctic-Boreal landscapes. After spectral adjustment, Landsat and Sentinel-2 data can be used to create temporally-dense time series and be applied to assess permafrost landscape changes in Eastern Siberia. Remaining differences between the sensors can be attributed to several factors including heterogeneous terrain, poor cloud and cloud shadow masking, and mixed pixels.
Lineament mapping, which is an important part of any structural geological investigation, is made more efficient and easier by the availability of optical as well as radar remote sensing data, such as Landsat and Sentinel with medium and high spatial resolutions. However, the results from these multi-resolution data vary due to their difference in spatial resolution and sensitivity to soil occupation. The accuracy and quality of extracted lineaments depend strongly on the spatial resolution of the imagery. Therefore, the aim of this study was to compare the optical Landsat-8, Sentinel-2A, and radar Sentinel-1A satellite data for automatic lineament extraction. The framework of automatic approach includes defining the optimal parameters for automatic lineament extraction with a combination of edge detection and line-linking algorithms and determining suitable bands from optical data suited for lineament mapping in the study area. For the result validation, the extracted lineaments are compared against the manually obtained lineaments through the application of directional filtering and edge enhancement as well as to the lineaments digitized from the existing geological maps of the study area. In addition, a digital elevation model (DEM) has been utilized for an accuracy assessment followed by the field verification. The obtained results show that the best correlation between automatically extracted lineaments, manual interpretation, and the preexisting lineament map is achieved from the radar Sentinel-1A images. The tests indicate that the radar data used in this study, with 5872 and 5865 lineaments extracted from VH and VV polarizations respectively, is more efficient for structural lineament mapping than the Landsat-8 and Sentinel-2A optical imagery, from which 2338 and 4745 lineaments were extracted respectively.
Obtaining representative hydrometric values is essential for characterizing extreme events, hydrological dynamics and detecting possible changes on the long-term hydrology. Reliability of streamflow data requires a temporal continuity and a maintenance of the gauging stations, which data are affected by epistemic and random sources of error. An assessment of discharge meterings' and stage-discharge rating curves' uncertainties were carried out by comparing the accuracy of the measuring instruments of two different hydrometric networks (i.e., one analogical and one digital) established in the same river location at the Mediterranean island of Mallorca. Furthermore, the effects of such uncertainties were assessed on the hydrological dynamics, considering the significant global change impacts beset this island. Evaluation was developed at four representative gauging stations of the hydrographic network with analogic (≈40 years) and digital (≈10 years) data series. The study revealed that the largest source of uncertainty in the analogical (28 to 274%) and in the digital (17–37%) networks were the stage-discharge rating curves. Their impact on the water resources was also evaluated at the event and annual scales, resulting in an average difference of water yields of 183% and 142% respectively. Such improvement on the comprehension of hydrometric networks uncertainties will dramatically benefit the interpretation of the long-term streamflow by providing better insights into the hydrologic and flood hazard planning, management and modelling.
Climate change is affecting the rate of carbon cycling, particularly in the Arctic. Permafrost degradation through deeper thaw and physical disturbances results in the release of carbon dioxide and methane to the atmosphere and to an increase in lateral dissolved organic matter (DOM) fluxes. Whereas riverine DOM fluxes of the large Arctic rivers are well assessed, knowledge is limited with regard to small catchments that cover more than 40% of the Arctic drainage basin. Here, we use absorption measurements to characterize changes in DOM quantity and quality in a low Arctic (Herschel Island, Yukon, Canada) and a high Arctic (Cape Bounty, Melville Island, Nunavut, Canada) setting with regard to geographical differences, impacts of permafrost degradation, and rainfall events. We find that DOM quantity and quality is controlled by differences in vegetation cover and soil organic carbon content (SOCC). The low Arctic site has higher SOCC and greater abundance of plant material resulting in higher chromophoric dissolved organic matter (cDOM) and dissolved organic carbon (DOC) than in the high Arctic. DOC concentration and cDOM in surface waters at both sites show strong linear relationships similar to the one for the great Arctic rivers. We used the optical characteristics of DOM such as cDOM absorption, specific ultraviolet absorbance (SUVA), ultraviolet (UV) spectral slopes (S275-295), and slope ratio (SR) for assessing quality changes downstream, at base flow and storm flow conditions, and in relation to permafrost disturbance. DOM in streams at both sites demonstrated optical signatures indicative of photodegradation downstream processes, even over short distances of 2000 m. Flow pathways and the connected hydrological residence time control DOM quality. Deeper flow pathways allow the export of permafrost-derived DOM (i.e. from deeper in the active layer), whereas shallow pathways with shorter residence times lead to the export of fresh surface- and near-surface-derived DOM. Compared to the large Arctic rivers, DOM quality exported from the small catchments studied here is much fresher and therefore prone to degradation. Assessing optical properties of DOM and linking them to catchment properties will be a useful tool for understanding changing DOM fluxes and quality at a pan-Arctic scale.
Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.
Water stable isotope signatures can provide valuable insights into the catchment internal runoff processes. However, the ability of the water isotope data to constrain the internal apportionments of runoff components in hydrological models for glacierized basins is not well understood. This study developed an approach to simultaneously model the water stable isotopic compositions and runoff processes in a glacierized basin in Central Asia. The fractionation and mixing processes of water stable isotopes in and from the various water sources were integrated into a glacio-hydrological model. The model parameters were calibrated on discharge, snow cover and glacier mass balance data, and additionally isotopic composition of streamflow. We investigated the value of water isotopic compositions for the calibration of model parameters, in comparison to calibration methods without using such measurements. Results indicate that: (1) The proposed isotope-hydrological integrated modeling approach was able to reproduce the isotopic composition of streamflow, and improved the model performance in the evaluation period; (2) Involving water isotopic composition for model calibration reduced the model parameter uncertainty, and helped to reduce the uncertainty in the quantification of runoff components; (3) The isotope-hydrological integrated modeling approach quantified the contributions of runoff components comparably to a three-component tracer-based end-member mixing analysis method for summer peak flows, and required less water tracer data. Our findings demonstrate the value of water isotopic compositions to improve the quantification of runoff components using hydrological models in glacierized basins.
The Schwarzenberg mining district in the western Erzgebirge hosts numerous skarn-hosted tin-polymetallic deposits, such as Breitenbrunn. The St. Christoph mine is located in the Breitenbrunn deposit and is the locus typicus of christophite, an iron-rich sphalerite variety, which can be associated with indium enrichment. This study presents a revision of the paragenetic scheme, a contribution to the indium behavior and potential, and discussion on the origin of the sulfur. This was achieved through reflected light microscopy, SEM-based MLA, EPMA, and bulk mineral sulfur isotope analysis on 37 sulfide-rich skarn samples from a mineral collection. The paragenetic scheme includes: a pre-mineralization stage of anhydrous calc-silicates and hydrous minerals; an oxide stage, dominated by magnetite; a sulfide stage of predominantly sphalerite, minor pyrite, chalcopyrite, arsenopyrite, and galena. Some sphalerite samples present elevated indium contents of up to 0.44 wt%. Elevated iron contents (4-10 wt%) in sphalerite can be tentatively linked to increased indium incorporation, but further analyses are required. Analyzed sulfides exhibit homogeneous delta S-34 values (-1 to +2 parts per thousand VCDT), assumed to be post-magmatic. They correlate with other Fe-Sn-Zn-Cu-In skarn deposits in the western Erzgebirge, and Permian vein-hosted associations throughout the Erzgebirge region.
We show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects.
Aquifer thermal energy storage (ATES) as a complement to fluctuating renewable energy systems is a reliable technology to guarantee continuous energy supply for heating and air conditioning. We investigated a high-temperature (HT) mono-well system (c. 100 degrees C), where the well screens are separated vertically within the aquifer, as an alternative to conventional doublet ATES systems for an underground storage in northern Oman. We analysed the impact of thermal inference between injection and extraction well screens on the heat recovery factor (HRF) in order to define the optimal screento-screen distance for best possible systems efficiency. Two controlling interference parameters were considered: the vertical screen-to-screen distance and aquifer heterogeneities. The sensitivity study shows that with decreasing screen-to-screen distances, thermal interference increases storage performance. A turning point is reached if the screen distance is too close, causing either water breakthrough or negative thermal interference between the screens. Our simulations show that a combined heat plume with spherical geometry results in the highest heat recovery factors due to the lowest surface area to volume ratios. Thick aquifers for mono-well HT-ATES are thus not mandatory Our study shows that a HT-ATES mono-well system is a feasible storage design with high heat recovery factors for continuous cooling or heating purposes.
Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights.
In response to collision and convergence between India and Asia during the Cenozoic, convergence took place between the Pamir and South Tian Shan. Here we present new detrital zircon U-Pb ages coupled with conglomerate clast counting and sedimentary data from the late Cenozoic Wuheshalu section in the convergence zone, to shed light on the convergence process of the Pamir and South Tian Shan. Large Triassic zircon U-Pb age populations in all seven samples suggest that Triassic igneous rocks from the North Pamir were the major source area for the late Cenozoic Wuheshalu section. In the Miocene, large populations of the North Pamir component supports rapid exhumation in the North Pamir and suggest that topography already existed there since the early Miocene. Exhumation of the South Tian Shan was relatively less important in the Miocene and its detritus could only reach a limited area in the foreland area. Gradually increasing sediment loading and convergence of the Pamir and South Tian Shan caused rapid subsidence in the convergence area. Since ca. 6-5.3 Ma, the combination of a major North Pamir component and a minor South Tian Shan component at the Wuheshalu section is consistent with active deformation of the South Tian Shan and the North Pamir. During deposition of the upper Atushi Formation, a larger proportion of North Pamir-derived sediments was deposited in the Wuheshalu section, maybe because faulting and northward propagation of the North Pamir caused northward displacement of the depocenter to north of the Wuheshalu section.
One commonly proposed method to limit flood risk is land-use or zoning policies which regulates construction in high-risk areas, in order to reduce economic exposure and its vulnerability to flood events. Although such zoning regulations can be effective in limiting trends in flood risk, they also have adverse impacts on society, for instance by limiting local development of areas near the water. In order to judge whether proposed land-use or zoning policies are a net benefit to society, they should be accepted or rejected based on a societal cost-benefit analysis (CBA). However, conducting a CBA of zoning regulation is complex and comprehensive guidelines of how to do such an analysis are lacking. We offer guidelines for good practice. In order to assess the costs and benefits of zoning as a climate change adaption strategy, they should be assessed at a societal level in order to account for public good features of flood risk reduction strategies, and because costs in one area can be benefits in another region. We propose a multistep process: first, determine the spatial extent of the zoning policy and how interconnected the zoned area is to other locations; second, conduct a CBA using monetary costs and benefits estimated from an integrated hydro-economic model to investigate if total benefits exceed total costs; third, conduct a sensitivity analysis regarding the main assumptions; fourth, conduct a multicriteria analysis (MCA) of the normative outcomes of a zoning policy. A desirable policy is preferred in both the CBA and MCA. This article is categorized under: Engineering Water > Planning Water Human Water > Value of Water Science of Water > Water Extremes Human Water > Methods
Garnet of eclogite (formerly termed garnet clinopyroxenite) hosted in lenses of orogenic garnet peridotite from the Granulitgebirge, NW Bohemian Massif, contains unique inclusions of granitic melt, now either glassy or crystallized. Analysed glasses and re‐homogenized inclusions are hydrous, peraluminous, and enriched in highly incompatible elements characteristic of the continental crust such as Cs, Li, B, Pb, Rb, Th, and U. The original melt thus represents a pristine, chemically evolved metasomatic agent, which infiltrated the mantle via deep continental subduction during the Variscan orogeny. The bulk chemical composition of the studied eclogites is similar to that of Fe‐rich basalt and the enrichment in LILE and U suggest a subduction‐related component. All these geochemical features confirm metasomatism. In comparison with many other garnet+clinopyroxene‐bearing lenses in peridotites of the Bohemian Massif, the studied samples from Rubinberg and Klatschmühle are more akin to eclogite than pyroxenites, as reflected in high jadeite content in clinopyroxene, relatively low Mg, Cr, and Ni but relatively high Ti. However, trace elements of both bulk rock and individual mineral phases show also important differences making these samples rather unique. Metasomatism involving a melt requiring a trace element pattern very similar to the composition reported here has been suggested for the source region of rocks of the so‐called durbachite suite, that is, ultrapotassic melanosyenites, which are found throughout the high‐grade Variscan basement. Moreover, the Th, U, Pb, Nb, Ta, and Ti patterns of these newly studied melt inclusions (MI) strongly resemble those observed for peridotite and its enclosed pyroxenite from the T‐7 borehole (Staré, České Středhoři Mountains) in N Bohemia. This suggests that a similar kind of crustal‐derived melt also occurred here. This study of granitic MI in eclogites from peridotites has provided the first direct characterization of a preserved metasomatic melt, possibly responsible for the metasomatism of several parts of the mantle in the Variscides.
Oceanic Anoxic Event-2 (OAE-2; Cenomanian-Turonian) is characterized by extensive deposition of organic carbon-rich deposits (black shales) in ocean basins worldwide as result of a major perturbation of the global carbon cycle. While the sedimentological, geochemical, and paleontological aspects of deep water expressions of OAE-2 have been intensively studied in the last few decades, much less attention has been given to the coeval shallow water deposits. In this study, we present the results of a detailed facies and petrographic (optical microscope and scanning electron microscopy) and geochemical studies (delta C-13(carb), delta C-13 (org), delta N-15(bulk), TOC, and Rock-Eval pyrolysis) on two key shallow marine sections from the Apennine Carbonate Platform (ACP; Italy). Here a continuous record of shallow water carbonates through the OAE-2 interval is preserved, offering the unique opportunity to document the archive of paleoenvironmental changes in a neritic setting, at a tropical latitude and far from the influence of a large continental block. Two conspicuous intervals are characterized by abundant and closely spaced dark microbial laminites found at correlative stratigraphic horizons in the two stratigraphic sections. These laminites contain elevated concentrations of TOC (up to 1%) relative to microbial capping cycles laminites stratigraphically above and below. The organic matter preserved in these fine-grained laminites is dominated by cyanobacteria remains, which accounted for most of the organic matter produced on the ACP in these intervals. Our study suggests that Tethyan carbonate platforms experienced significant biological changes during OAE-2, alternating, in few kiloyears, between eutrophic phases dominated by microbial communities and mesotrophic/oligotrophic conditions favoring normal carbonate production/sedimentation. The synchronous occurrence of microbialite facies at different locations across the ACP and on other platforms worldwide (e.g., Mexico and Croatia) suggests a causal link between Large Igneous Province volcanism and the environmental conditions necessary to trigger cyanobacterial proliferation on shallow carbonate platforms.