Refine
Year of publication
Document Type
- Article (2810) (remove)
Keywords
- Holocene (35)
- climate change (30)
- Earthquake source observations (21)
- Climate change (18)
- Pollen (18)
- erosion (18)
- permafrost (18)
- Himalaya (17)
- Tibetan Plateau (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2810) (remove)
Crustal earthquake swarms are an expression of intensive cracking and rock damaging over periods of days, weeks or month in a small source region in the crust. They are caused by longer lasting stress changes in the source region. Often, the localized stressing of the crust is associated with fluid or gas migration, possibly in combination with pre-existing zones of weaknesses. However, verifying and quantifying localized fluid movement at depth remains difficult since the area affected is small and geophysical prospecting methods often cannot reach the required resolution.
We apply a simple and robust method to estimate the velocity ratio between compressional (P) and shear (S) waves (upsilon(P)/upsilon(S)-ratio) in the source region of an earthquake swarm. The upsilon(P)/upsilon(S)-ratio may be unusual small if the swarm is related to gas in a porous or fractured rock. The method uses arrival time difference between P and S waves observed at surface seismic stations, and the associated double differences between pairs of earthquakes. An advantage is that earthquake locations are not required and the method seems lesser dependent on unknown velocity variations in the crust outside the source region. It is, thus, suited for monitoring purposes.
Applications comprise three natural, mid-crustal (8-10 km) earthquake swarms between 1997 and 2008 from the NW-Bohemia swarm region. We resolve a strong temporal decrease of upsilon(P)/upsilon(S) before and during the main activity of the swarm, and a recovery of upsilon(P)/upsilon(S) to background levels at the end of the swarms. The anomalies are interpreted in terms of the Biot-Gassman equations, assuming the presence of oversaturated fluids degassing during the beginning phase of the swarm activity.
QuestionHow important is the effect of micro-relief and vegetation type on the characteristics of vascular plants and bryophytes in low-centred polygons? LocationSiberian Arctic, Russia. MethodsEight low-centred polygons in northern Siberia were surveyed for vegetation along transects running from the rim to the pond via the rim-pond transition of each polygon and across a vegetation type gradient from open forest to tundra. ResultsThe cover of vascular plants and bryophytes displays no significant differences between the rim and rim-pond transition but is significantly lower in the pond section of the polygons. Alpha-diversity of vascular plants decreases strongly from rim to pond, whereas bryophyte diversity in pond plots is significantly distinct from the rim and the rim-pond transition. There is no clear trend in cover for either plant group along the vegetation type transect and only a weak trend in -diversity. However, both gradients are reflected in the compositional turnover. The applied indicator species analysis identified taxa characteristic of certain environmental conditions. Among others, we found vascular plants primarily characteristic of the rim and bryophyte taxa characteristic of each micro-relief level and vegetation type. ConclusionsThe observed gradual pattern in -diversity and composition of polygonal vegetation suggests that micro-relief is the main driver of changes in the vegetation composition, while vegetation type and the related forest cover change are of subordinate importance for polygonal vegetation patterns along the Siberian tree line.
African savannas are primarily used as pastures and are subject to changes in climate and management strategies. For sustainable management of these landscapes ecological knowledge on seasonal and long-term variability in plant community composition and the availability of green biomass is essential. In this study, we assessed the effects of dry and wet season on species richness and beta diversity for three sites along a gradient of increasing vegetation cover and precipitation in northwest Namibia. A hexagonal systematic sampling design was used to record floristic data. The Simple Matching, Soerensen, and multi-plot similarity coefficient and distance decay analyses were applied for examining beta diversity. Analyses were repeated while separating the plots according to the presence of woody vegetation. Species richness nearly doubled from dry to wet season; compositional similarity increased from dry to wet season and with increasing aridity of the study sites: distance decay was more pronounced in the dry season without any link to the precipitation gradient. Woody elements in the landscape, which occur along drainage lines or as tree islands, govern spatial and seasonal plant diversity fluctuations. Monitoring them is important for conservation strategies and for establishing grazing rules that ensure a sustainable use of savanna ecosystems.
We assessed tropical montane cloud forest (TMCF) sensitivity to natural disturbance by drought, fire, and dieback with a 7300-year-long paleorecord. We analyzed pollen assemblages, charcoal accumulation rates, and higher plant biomarker compounds (average chain length [ACL] of n-alkanes) in sediments from Wai 'anapanapa, a small lake near the upper forest limit and the mean trade wind inversion ('IWI) in Hawai`i. The paleorecord of ACL suggests increased drought frequency and a lower awl elevation from 2555-1323 cal yr B.P. and 606-334 cal yr B.P. Charcoal began to accumulate and a novel fire regime was initiated ca. 880 cal yr B.P., followed by a decreased fire return interval at ca. 550 cal yr B.P. Diebacks occurred at 2931, 2161, 1162, and 306 cal yr B.P., and two of these were independent of drought or fire. Pollen assemblages indicate that on average species composition changed only 2.8% per decade. These dynamics, though slight, were significantly associated with disturbance. The direction of species composition change varied with disturbance type. Drought was associated with significantly more vines and lianas; fire was associated with an increase in the tree fern Sadleria and indicators of open, disturbed landscapes at the expense of epiphytic ferns; whereas stand-scale dieback was associated with an increase in the tree fern Cibotium. Though this cloud forest was dynamic in response to past disturbance, it has recovered, suggesting a resilient TMCF with no evidence of state change in vegetation type (e.g., grassland or shrubland).
By reducing current velocity, tidal marsh vegetation can diminish storm surges and storm waves. Conversely, currents often exert high mechanical stresses onto the plants and hence affect vegetation structure and plant characteristics. In our study, we aim at analysing this interaction from both angles. On the one hand, we quantify the reduction of current velocity by Bolboschoenus maritimus, and on the other hand, we identify functional traits of B. maritimus’ ramets along environmental gradients. Our results show that tidal marsh vegetation is able to buffer a large proportion of the flow velocity at currents under normal conditions. Cross-shore current velocity decreased with distance from the marsh edge and was reduced by more than 50% after 15 m of vegetation. We were furthermore able to show that plants growing at the marsh edge had a significantly larger diameter than plants from inside the vegetation. We found a positive correlation between plant thickness and cross-shore current which could provide an adaptive value in habitats with high mechanical stress. With the adapted morphology of plants growing at the highly exposed marsh edge, the entire vegetation belt is able to better resist the mechanical stress of high current velocities. This self-adaptive effect thus increases the ability of B. maritimus to grow and persist in the pioneer zone and may hence better contribute to ecosystem-based coastal protection by reducing current velocity.
Siberian arctic vegetation and lake water communities, known for their temperature dependence, are expected to be particularly impacted by recent climate change and high warming rates. However, decadal information on the nature and strength of recent vegetation change and its time lag to climate signals are rare. In this study, we present a Pb-210/Cs-137 dated pollen and Pediastrum species record from a unnamed lake in the south of the Taymyr peninsula covering the period from AD 1706 to 2011. Thirty-nine palynomorphs and 10 morphotypes of Pediastrum species were studied to assess changes in vegetation and lake conditions as probable responses to climate change. We compared the pollen record with Pediastrum species, which we consider to be important proxies of climate changes. Three pollen assemblage zones characterised by Betula nana, Alnus viridis and Larix gmelinii (1706-1808); herbs such as Cyperaceae, Artemisia or Senecio (1808-1879), and higher abundance of Larix pollen (1955-2011) are visible. Also, three Pediastrum assemblage zones show changes of aquatic conditions: higher abundances of Pediastrum boryanum var. brevicorne (1706-1802); medium abundances of P. kawraiskyi and P. integrum (1802-1840 and 1920-1980), indicating cooler conditions while less eutrophic conditions are indicated by P. boryanum, and a mainly balanced composition with only small changes of cold- and warm-adapted Pediastrum species (1965-2011). In general, compositional Pediastrum species turnover is slightly higher than that indicated by pollen data (0.54 vs 0.34 SD), but both are only minor for this treeline location. In conclusion, the relevance of differentiation of Pediastrum species is promising and can give further insights into the relationship between lakes and their surrounding vegetation transferred onto climatic conditions.
The sedimentary record of the Dead Sea provides an exceptional high-resolution archive of past climate changes in the drought-sensitive eastern Mediterranean-Levant, a key region for the development of humankind at the boundary of global climate belts. Moreover, it is the only deep hypersaline lake known to have deposited long sequences of finely laminated, annually deposited sediments (i.e. varves) of varied compositions, including aragonite, gypsum, halite and clastic sediments. Vast efforts have been made over the years to decipher the environmental information stored in these evaporitic-clastic sequences spanning from the Pleistocene Lake Amora to the Holocene Dead Sea. A general characterisation of sediment facies has been derived from exposed sediment sections, as well as from shallow- and deep-water sediment cores. During high lake stands and episodes of positive water budget, mostly during glacial times, alternating aragonite and detritus laminae (‘aad’ facies) were accumulated, whereas during low lake stands and droughts, prevailing during interglacials, laminated detritus (‘ld’ facies) and laminated halite (‘lh’ facies) dominate the sequence. In this paper, we (i) review the three types of laminated sediments of the Dead Sea sedimentary record (‘aad’, ‘ld’ and ‘lh’ facies), (ii) discuss their modes of formation, deposition and accumulation, and their interpretation as varves, and (iii) illustrate how Dead Sea varves are utilized for palaeoclimate reconstructions and for establishing floating chronologies.
Fluvial planation surfaces, such as straths, commonly serve as recorders of climatic and tectonic changes and are formed by the lateral erosion of rivers, a process that remains poorly understood. Here we present a study of kilometer-wide, fluvially eroded, low-relief surfaces on rapidly uplifting folds in the foreland of the southwestern Tian Shan. A combination of field work, digital elevation model analysis, and dating of fluvial deposits reveals that despite an arid climate and rapid average rock-uplift rates of 1-3mm/yr, rivers cut extensive (>1-2km wide) surfaces with typical height variations of <6m over periods of >2-6kyr. The extent of this beveling varies in space and time, such that different beveling episodes affect individual structures. Between times of planation, beveled surfaces are abandoned, incised, and deformed across the folds. In a challenge to models that link strath cutting and abandonment primarily to changes in river incision rates, we demonstrate that lateral erosion rates of antecedent streams crossing the folds have to vary by more than 1 order of magnitude to explain the creation of beveled platforms in the past and their incision at the present day. These variations do not appear to covary with climate variability and might be caused by relatively small (much less than an order of magnitude) changes in sediment or water fluxes. It remains uncertain in which settings variations in lateral bedrock erosion rates predominate over changes in vertical erosion rates. Therefore, when studying fluvial planation and strath terraces, variability of both lateral and vertical erosion rates should be considered.
Rivers transfer particulate organic carbon (POC) from eroding mountains into geological sinks. Organic carbon source composition and selective mobilization have been shown to affect the type and quantity of POC export, but their combined effects across complex mountain ranges remain underexplored. Here, we examine the variation in organic carbon sourcing and transport in the trans-Himalayan Kali Gandaki River catchment, along strong gradients in precipitation, rock type and vegetation. Combining bulk stable nitrogen, and stable and radioactive organic carbon isotopic composition of bedrock, litter, soil and river sediment samples with a Bayesian end-member mixing approach, we differentiate POC sources along the river and quantify their export. Our analysis shows that POC export from the Tibetan segment of the catchment, where carbon bearing shales are partially covered by aged and modern soils, is dominated by petrogenic POC. Based on our data we re-assess the presence of aged biospheric OC in this part of the catchment, and its contribution to the river load. In the High Himalayan segment, we observed low inputs of petrogenic and biospheric POC, likely due to very low organic carbon concentrations in the metamorphic bedrock, combined with erosion dominated by deep-seated landslides. Our findings show that along the Kali Gandaki River, the sourcing of sediment and organic carbon are decoupled, due to differences in rock organic carbon content, soil and above ground carbon stocks, and geomorphic process activity. While the fast eroding High Himalayas are the principal source of river sediment, the Tibetan headwaters, where erosion rates are lower, are the principal source of organic carbon. To robustly estimate organic carbon export from the Himalayas, the mountain range should be divided into tectono-physiographic zones with distinct organic carbon yields due to differences in substrate and erosion processes and rates.
Small-scale variations in mineral chemistry, textures, and platinum group element (PGE) mineralization were investigated in the Lower and Middle Group chromitite layers LG6, LG6a, MG1, MG2, and MG2 II from vertical drill core profiles at the Thaba mine in the northwestern limb of the Bushveld Complex. We present detailed geochemical profiles of chromite composition and chromite crystal size distribution curves to shed light on the processes of chromite accumulation and textural modification as well as mineralization. Multiple samples within each layer were assayed for PGE concentrations, and the respective platinum group mineral association was determined by mineral liberation analysis (MLA). There is strong evidence for postcumulus changes in the chromitites. The crystal size distribution curves suggest that the primary chromite texture was coarsened by a combination of adcumulus growth and textural equilibration, while compaction of the crystal mush played only a minor role. Mineral compositions were also modified by postcumulus processes, but because of the very high modal amount of chromite and its local preservation in orthopyroxene oikocrysts, that phase retained much primary information. Vertical variations of chromite composition within chromitite layers and from one layer to another do not support the idea of chromite accumulation from crystal-rich slurries or crystal settling from a large magma chamber. Instead, we favor a successive buildup of chromitite layers by repeated injections of relatively thin layers of chromite-saturated magmas, with in situ crystallization occurring at the crystal mush-magma interface. The adcumulus growth of chromite grains to form massive chromitite required addition of Cr to the layers, which we attribute to downward percolation from the overlying magma. The PGE concentrations are elevated in all chromitite layers compared to adjacent silicate rocks and show a systematic increase upward from LG6 (avg 807 ppb Ir + Ru + Rh + Pt + Pd + Au) to MG2 II (avg 2,062 ppb). There are also significant internal variations in all layers, with enrichments at hanging and/or footwalls. The enriched nature of chromitites in PGEs compared to host pyroxenites is a general feature, independent of the layer thickness. The MLA results distinguish two principal groups of PGE mineral associations: the LG6, LG6, and MG1 are dominated by the malanite series, laurite, and PGE sulfarsenides, while the MG2 and MG2 II layers are characterized by laurite and PGE sulfides as well as Pt-Fe-Sn and PGE-Sb-Bi-Pb alloys. Differences in the PGE associations are attributed to postcumulus alteration of the MG2 and MG2 II layer, while the chromitites below, particularly LG6 and LG6a, contain a more pristine association.
With the growing size and use of night light time series from the Visible Infrared Imaging Radiometer Suite Day/Night Band (DNB), it is important to understand the stability of the dataset. All satellites observe differences in pixel values during repeat observations. In the case of night light data, these changes can be due to both environmental effects and changes in light emission. Here we examine the stability of individual locations of particular large scale light sources (e.g., airports and prisons) in the monthly composites of DNB data from April 2012 to September 2017. The radiances for individual pixels of most large light emitters are approximately normally distributed, with a standard deviation of typically 15-20% of the mean. Greenhouses and flares, however, are not stable sources. We observe geospatial autocorrelation in the monthly variations for nearby sites, while the correlation for sites separated by large distances is small. This suggests that local factors contribute most to the variation in the pixel radiances and furthermore that averaging radiances over large areas will reduce the total variation. A better understanding of the causes of temporal variation would improve the sensitivity of DNB to lighting changes.
Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal
(2020)
The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.
The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
Variability of the Cold Season Climate in Central Asia. Part II: Hydroclimatic Predictability
(2019)
Central Asia (CA) is subjected to a large variability of precipitation. This study presents a statistical model, relating precipitation anomalies in three subregions of CA in the cold season (November-March) with various predictors in the preceding October. Promising forecast skill is achieved for two subregions covering 1) Uzbekistan, Turkmenistan, Kyrgyzstan, Tajikistan, and southern Kazakhstan and 2) Iran, Afghanistan, and Pakistan. ENSO in October is identified as the major predictor. Eurasian snow cover and the quasi-biennial oscillation further improve the forecast performance. To understand the physical mechanisms, an analysis of teleconnections between these predictors and the wintertime circulation over CA is conducted. The correlation analysis of predictors and large-scale circulation indices suggests a seasonal persistence of tropical circulation modes and a dynamical forcing of the westerly circulation by snow cover variations over Eurasia. An EOF analysis of pressure and humidity patterns allows separating the circulation variability over CA into westerly and tropical modes and confirms that the identified predictors affect the respective circulation characteristics. Based on the previously established weather type classification for CA, the predictors are investigated with regard to their effect on the regional circulation. The results suggest a modification of the Hadley cell due to ENSO variations, with enhanced moisture supply from the Arabian Gulf during El Nino. They further indicate an influence of Eurasian snow cover on the wintertime Arctic Oscillation (AO) and Northern Hemispheric Rossby wave tracks. Positive anomalies favor weather types associated with dry conditions, while negative anomalies promote the formation of a quasi-stationary trough over CA, which typically occurs during positive AO conditions.
A reconstruction of Milankovitch to millennial-scale variability of sea surface temperature (SST) and sea surface productivity in the Pleistocene midlatitude North Atlantic Ocean (marine isotope stage (MIS) 16-9) and its relationship to ice sheet instability was carried out on sediments from Integrated Ocean Drilling Program (IODP) Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, U-37(K)' index, total organic carbon (TOC) and carbonate contents, X-ray diffraction data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggests increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19 degrees C during MIS 11.3 and absolute minima of <10 degrees C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2 are related to Heinrich-like event meltwater pulses, as supported by the coincidence of SST minima and maxima in C-37:4 alkenones and dolomite. These sudden meltwater pulses, especially during terminations IV and V, probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the northern high latitudes.
To better understand the reaction of Arctic coasts to increasing environmental pressure, coastal changes along a 210-km length of the Yukon Territory coast in north-west Canada were investigated. Shoreline positions were acquired from aerial and satellite images between 1951 and 2011. Shoreline change rates were calculated for multiple time periods along the entire coast and at six key sites. Additionally, Differential Global Positioning System (DGPS) measurements of shoreline positions from seven field sites were used to analyze coastal dynamics from 1991 to 2015 at higher spatial resolution. The whole coast has a consistent, spatially averaged mean rate of shoreline change of 0.7 +/- 0.2 m/a with a general trend of decreasing erosion from west to east. Additional data from six key sites shows that the mean shoreline change rate decreased from -1.3 +/- 0.8 (1950s-1970s) to -0.5 +/- 0.6 m/a (1970s-1990s). This was followed by a significant increase in shoreline change to -1.3 +/- 0.3 m/a in the 1990s to 2011. This increase is confirmed by DGPS measurements that indicate increased erosion rates at local rates up to -8.9 m/a since 2006. Ground surveys and observations with remote sensing data indicate that the current rate of shoreline retreat along some parts of the Yukon coast is higher than at any time before in the 64-year-long observation record. Enhanced availability of material in turn might favor the buildup of gravel features, which have been growing in extent throughout the last six decades. Plain Language Summary The Arctic is warming, but the impacts on its coasts are not well documented. To better understand the reaction of Arctic coasts to increasing environmental pressure, shoreline position changes along a 210-km length of the Yukon Territory coast in northwest Canada were investigated for the time period from 1951 to 2015. Shoreline positions were extracted from historical aerial images from the 1950s, 1970s, and 1990s and from satellite images from 2011. Additionally, measurements of shoreline positions from field sites were used to analyze coastal dynamics from 1991 to 2015. The mean shoreline change rate was -1.3 m/a between the 1950s and 1970s and followed by a decrease to -0.5 m/a between the 1970s to 1990s. This was followed by a significant increase in mean shoreline change rates again to -1.3 m/a in the 1990s to 2011 time period. This acceleration in erosion is confirmed by field measurements that indicate increased erosion rates at high local rates up to -8.9 m/a since 2006. Enhanced coastal erosion might, in turn, favor the buildup of gravel features, which have been growing in extent throughout the last six decades.
In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12–30 m), and ALOS World 3D (5–30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12–30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m=n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2–8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e.g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.
In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12-30 m), and ALOS World 3D (5-30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12-30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m = n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2-8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e. g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.
The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.
The apparent isotope enrichment factor epsilon(macrophyte) of submerged plants (epsilon(macrophyte-DIC) = delta C-13(macrophyte) - delta C-13(DIC)) is indicative of dissolved inorganic carbon (DIC) supply in neutral to alkaline waters and is related to variations in aquatic productivity (Papadimitriou et al. in Limnol Oceanogr 50:1084-1095, 2005). This paper aims to evaluate the usage of epsilon(macrophyte) inferred from isotopic analyses of submerged plant fossils in addition to analyses of lake carbonate as a palaeolimnological proxy for former HCO3 (-) concentrations. Stable carbon isotopic analysis of modern Potamogeton pectinatus leaves and its host water DIC from the Tibetan Plateau and Central Yakutia (Russia) yielded values between -23.3 and +0.4aEuro degrees and between +14.0 and +6.5aEuro degrees, respectively. Values of epsilon (Potamogeton-DIC) (range -15.4 to +1.1aEuro degrees) from these lakes are significantly correlated with host water HCO3 (-) concentration (range 78-2,200 mg/l) (r = -0.86; P < 0.001), thus allowing for the development of a transfer function. Palaeo-epsilon (Potamogeton-ostracods) values from Luanhaizi Lake on the NE Tibetan Plateau, as inferred from the stable carbon isotope measurement of fossil Potamogeton pectinatus seeds (range -24 to +2.8aEuro degrees) and ostracods (range -7.8 to +7.5%) range between -14.8 and 1.6aEuro degrees. Phases of assumed disequilibrium between delta C-13(DIC) and delta C-13(ostracods) known to occur in charophyte swards (as indicated by the deposition of charophyte fossils) were excluded from the analysis of palaeo-epsilon. The application of the epsilon (Potamogeton-DIC)-HCO3 (-) transfer function yielded a median palaeo-HCO3 (-) -concentration of 290 mg/l. Variations in the dissolved organic carbon supply compare well with aquatic plant productivity changes and lake level variability as inferred from a multiproxy study of the same record including analyses of plant macrofossils, ostracods, carbonate and organic content.
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations, and biological pCO(2) drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using Th-230-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections, and problems with the reference data sets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.
Due to increasing demands and competition for high quality groundwater resources in many parts of the world, there is an urgent need for efficient methods that shed light on the interplay between complex natural settings and anthropogenic impacts. Thus a new approach is introduced, that aims to identify and quantify the predominant processes or factors of influence that drive groundwater and lake water dynamics on a catchment scale. The approach involves a non-linear dimension reduction method called Isometric feature mapping (Isomap). This method is applied to time series of groundwater head and lake water level data from a complex geological setting in Northeastern Germany. Two factors explaining more than 95% of the observed spatial variations are identified: (1) the anthropogenic impact of a waterworks in the study area and (2) natural groundwater recharge with different degrees of dampening at the respective sites of observation. The approach enables a presumption-free assessment to be made of the existing geological conception in the catchment, leading to an extension of the conception. Previously unknown hydraulic connections between two aquifers are identified, and connections revealed between surface water bodies and groundwater. (C) 2014 Elsevier B.V. All rights reserved.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
The aims of this study were to identify areas of wind erosion and dust deposition and to quantify the effects of different grazing intensities on soil redistribution rates in grasslands based on the Cs-137 technique. Because the method uses a reference inventory as threshold for erosion or deposition, the classification of any other site as source or sink for dust depends on the accurate selection of this reference site.
Measurements of Cs-137 inventories and depth distributions were carried out at pasture sites with predominant species of Stipa grandis and Leymus chinensis which are grazed with different intensities. Additional measurements were made at arable land, plant-covered sand dunes and alluvial plains. Wind-induced soil erosion and dust deposition rates were calculated from Cs-137 inventories by means of the "Profile-Distribution" and the "Mass Balance II" models.
The selection of the reference site was based on fluid dynamical and process-determining parameters. The chosen site should meet the following four conditions: (i) located at a summit position with obviously low deposition rates, (ii) sufficient vegetation cover to prevent wind erosion, (iii) plane to exclude water erosion and (iv) in the wind/dust shadow of a higher elevation. The measured reference inventory of Cs-137 was 1967(+/- 102) Bqm(-2) located at a summit position of moderately grazed Leymus chinensis steppe. The Cs-137 inventories at other sites ranged from 1330 Bqm(-2) at heavily grazed sites to 5119 Bqm(-2) at river deposits, representing annual average soil losses of up to 130 tkm(-2) and deposits of up to 540 tkm(-2), respectively. The calculated annual averages of dust depositions at ungrazed Leymus chinensis sites were related to the dust storm frequencies of the last 50 years resulting in a description of the temporal variability of annual dust depositions from about 154 tkm(-2) in the 1960s to 26 tkm(-2) at recent times. Based on this quantification already 80% of the total dust depositions can be related to the 20 years between the 1960s and the end of the 1970s and only 20% to the time between 1980 and 2001.
Cs-137 technique is a promising method to assess the effect of grazing intensity and land use types on the spatial variability of wind-induced soil and dust redistribution processes in semi-arid grasslands. However, considerable efforts are needed to identify a reliable reference site, because erosion and deposition induced by wind may occur at the same places. The combination of the dust deposition rates derived from Cs-137 profile data with the dust storm frequencies is helpful for a better reconstruction of the temporal variability of dust deposition and wind erosion in this region. The calculated recent deposition rates of about 20 tkm(-2) are in good agreement with data of other authors.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
Cities will play a key role in the grand challenge of nourishing a growing global population, because, due to their population density, they set the demand. To ensure that food systems are sustainable, as well as nourishing, one solution often suggested is to shorten their supply chains toward a regional rather than a global basis. While such regional systems may have a range of costs and benefits, we investigate the mitigation potential of regionalized urban food systems by examining the greenhouse gas emissions associated with food transport. Using data on food consumption for 7108 urban administrative units (UAUs), we simulate total transport emissions for both regionalized and globalized supply chains. In regionalized systems, the UAUs' demands are fulfilled by peripheral food production, whereas to simulate global supply chains, food demand is met from an international pool (where the origin can be any location globally). We estimate that regionalized systems could reduce current emissions from food transport. However, because longer supply chains benefit from maximizing comparative advantage, this emission reduction would require closing yield gaps, reducing food waste, shifting toward diversified farming, and consuming seasonal produce. Regionalization of food systems will be an essential component to limit global warming to well below 2 degrees C in the future.
The chapter illustrates how cities can be understood from a system–theory perspective as complex social systems. It argues that the classical and linear intervention methods are often no (longer) suitable for the complex structure, temporal dynamics and multifaceted processuality of urban development. It offers a systemic and systems theory-inspired method as an alternative approach, which allows for extended possibilities that are more appropriate for dealing with urban development processes. The method was developed on the basis of practical experience and theoretical insights. The approach should demonstrate for local decision-makers potential areas of activity for organising urban changes through co-design.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Passive continental margins offer the unique opportunity to study the processes involved in continental extension and break-up. Within the LISPWAL (Llthospheric Structure of the Namibian continental Passive margin at the intersection with the Walvis Ridge from amphibious seismic investigations) project, combined on- and offshore seismic experiments were designed to characterize the Southern African passive margin at the Walvis Ridge in northern Namibia. In addition to extensive analysis of the crustal structures, we carried out seismic investigations targeting the velocity structure of the upper mantle in the landfall region of the Walvis Ridge with the Namibian coast. Upper mantle P-n travel time tomography from controlled source, amphibious seismic data was used to investigate the sub-Moho upper mantle seismic velocity. We succeeded in imaging upper mantle structures potentially associated with continental break-up and/or the Tristan da Cunha hotspot track. We found mostly coast-parallel sub-Moho velocity anomalies, interpreted as structures which were created during Gondwana break-up. (C)2016 Published by Elsevier B.V.
We present a high-resolution study of the upper mantle structure of Central Europe, including the western part of the East European Platform, based on S-receiver functions of 345 stations. A distinct contrast is found between Phanerozoic Europe and the East European Craton across the Trans-European Suture Zone. To the west, a pronounced velocity reduction with depth interpreted as lithosphere-asthenosphere boundary (LAB) is found at an average depth of 90 km. Beneath the craton, no strong and continuous LAB conversion is observed. Instead we find a distinct velocity reduction within the lithosphere, at 80-120 km depth. This mid-lithospheric discontinuity (MLD) is attributed to a compositional boundary between depleted and more fertile lithosphere created by late Proterozoic metasomatism. A potential LAB phase beneath the craton is very weak and varies in depth between 180 and 250 km, consistent with a reduced velocity contrast between the lower lithosphere and the asthenosphere. Within the Trans-European Suture Zone, lithospheric structure is characterized by strong heterogeneity. A dipping or step-wise increase to LAB depth of 150 km is imaged from Phanerozoic Europe to 20-22 degrees E, whereas no direct connection to the cratonic LAB or MLD to the east is apparent. At larger depths, a positive conversion associated with the lower boundary of the asthenosphere is imaged at 210-250 km depth beneath Phanerozoic Europe, continuing down to 300 km depth beneath the craton. Conversions from both 410 km and 660 km discontinuities are found at their nominal depth beneath Phanerozoic Europe, and the discontinuity at 410 km depth can also be traced into the craton. A potential negative conversion on top of the 410 km discontinuity found in migrated images is analyzed by modeling and attributed to interference with other converted phases.
We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.
This study presents an unsupervised feature selection and learning approach for the discovery and intuitive imaging of significant temporal patterns in seismic single-station or network recordings. For this purpose, the data are parametrized by real-valued feature vectors for short time windows using standard analysis tools for seismic data, such as frequency-wavenumber, polarization, and spectral analysis. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure, which is in particular suitable for high-dimensional data sets. Our feature selection method is based on significance testing using the Wald-Wolfowitz runs test for-individual features and on correlation hunting with SOMs in feature subsets. Using synthetics composed of Rayleigh and Love waves and real-world data, we show the robustness and the improved discriminative power of that approach compared to feature subsets manually selected from individual wavefield parametrization methods. Furthermore, the capability of the clustering and visualization techniques to investigate the discrimination of wave phases is shown by means of synthetic waveforms and regional earthquake recordings.
Salt marshes filter pollutants, protect coastlines against storm surges, and sequester carbon, yet are under threat from sea level rise and anthropogenic modification. The sustained existence of the salt marsh ecosystem depends on the topographic evolution of marsh platforms. Quantifying marsh platform topography is vital for improving the management of these valuable landscapes. The determination of platform boundaries currently relies on supervised classification methods requiring near-infrared data to detect vegetation, or demands labour-intensive field surveys and digitisation. We propose a novel, unsupervised method to reproducibly isolate salt marsh scarps and platforms from a digital elevation model (DEM), referred to as Topographic Identification of Platforms (TIP). Field observations and numerical models show that salt marshes mature into subhorizontal platforms delineated by subvertical scarps. Based on this premise, we identify scarps as lines of local maxima on a slope raster, then fill landmasses from the scarps upward, thus isolating mature marsh platforms. We test the TIP method using lidar-derived DEMs from six salt marshes in England with varying tidal ranges and geometries, for which topographic platforms were manually isolated from tidal flats. Agreement between manual and unsupervised classification exceeds 94% for DEM resolutions of 1 m, with all but one site maintaining an accuracy superior to 90% for resolutions up to 3 m. For resolutions of 1 m, platforms detected with the TIP method are comparable in surface area to digitised platforms and have similar elevation distributions. We also find that our method allows for the accurate detection of local block failures as small as 3 times the DEM resolution. Detailed inspection reveals that although tidal creeks were digitised as part of the marsh platform, unsupervised classification categorises them as part of the tidal flat, causing an increase in false negatives and overall platform perimeter. This suggests our method may benefit from combination with existing creek detection algorithms. Fallen blocks and high tidal flat portions, associated with potential pioneer zones, can also lead to differences between our method and supervised mapping. Although pioneer zones prove difficult to classify using a topographic method, we suggest that these transition areas should be considered when analysing erosion and accretion processes, particularly in the case of incipient marsh platforms. Ultimately, we have shown that unsupervised classification of marsh platforms from high-resolution topography is possible and sufficient to monitor and analyse topographic evolution.
The Subandean fold and thrust belt of Bolivia constitutes the easternmost part of the Andean orogen that reflects thin-skinned shortening and eastward propagation of the Andean deformation front. The exact interplay of tectonics, climate, and erosion in the deposition of up to 7.5 km of late Cenozoic strata exposed in the Subandes remains unclear. To better constrain these relationships, we use four W-E industry seismic reflection profiles, eight new zircon U-Pb ages from Mio-Pliocene sedimentary strata, and cross-section balancing to evaluate the rates of thrust propagation, shortening, and deposition pinch-out migration. Eastward thrusting arrived in the Subandean belt at similar to 12.4 +/- 0.5 Ma and propagated rapidly toward the foreland unit approximately 6 Ma. This was followed by out-of- sequence deformation from ca. 4 to 2.1 Ma and by renewed eastward propagation thereafter. Our results show that the thrust-front propagation- and deposition pinch-out migration rates mimic the sediment accumulation rate. The rates of deposition pinchout migration and thrust propagation increased three- and two fold, respectively (8 mm/a; 3.3 mm/a) at 86 Ma. The three-fold increase in deposition pinch-out migration rate at this time is an indication of enhanced erosional efficiency in the hinterland, probably coupled with flexural rebound of the basin. Following the pulse of pinch-out migration, the Subandean belt witnessed rapid similar to 80 km eastward propagation of thrusting to the La Vertiente structure at 6 Ma. As there is no evidence for this event of thrust front migration being linked to an increase in shortening rate, the enhanced frontal accretion suggests a shift to supercritical wedge taper conditions. We propose that the supercritical state was due to a drop in basal strength, caused by sediment loading and pore fluid overpressure. This scenario implies that climate-controlled variation in erosional efficiency was the driver of late Miocene mass redistribution, which induced flexural rebound of the Subandean thrust belt, spreading of a large clastic wedge across the basin, and subsequent thrust-front propagation.
The Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
The architecture of coastal sequences in tectonically-active regions results mostly from a combination of sea-level and land-level changes. The objective of this study is to unravel these signals by combining sequence stratigraphy and sedimentology of near-shore sedimentary sequences in wave-built terraces. We focus on Santa Maria Island at the south-central Chile margin, which hosts excellent exposures of coastal sediments from Marine Isotope Stage 3. A novel method based on statistical analysis of grain-size distributions coupled with fades descriptions provided a detailed account of transgressive-regressive cycles. Radiocarbon ages from paleosols constrain the chronology between >53 and similar to 31 cal ka BP. Because the influence of glaciations can be neglected, we calculated relative sea-level curves by tying the onset of deposition on a bedrock abrasion platform to a global sea-level curve. The observed depositional cycles match those predicted for uplift rates between 1.2 and 1.8 m/ka. The studied sedimentary units represent depositional cycles that resulted in reoccupation events of an existing marine terrace. Our study demonstrates wave-built marine terrace deposits along clastic shorelines in temperate regions can be used to distinguish between tectonic uplift and climate-induced sea-level changes.
Reductions in emissions have successfully led to a regional decline in atmospheric nitrogen depositions over the past 20 years. By analyzing long-term data from 110 mountainous streams draining into German drinking water reservoirs, nitrate concentrations indeed declined in the majority of catchments. Furthermore, our meta-analysis indicates that the declining nitrate levels are linked to the release of dissolved iron to streams likely due to a reductive dissolution of iron(III) minerals in riparian wetland soils. This dissolution process mobilized adsorbed compounds, such as phosphate, dissolved organic carbon and arsenic, resulting in concentration increases in the streams and higher inputs to receiving drinking water reservoirs. Reductive mobilization was most significant in catchments with stream nitrate concentrations < 6 mg L-1. Here, nitrate, as a competing electron acceptor, was too low in concentration to inhibit microbial iron(III) reduction. Consequently, observed trends were strongest in forested catchments, where nitrate concentrations were unaffected by agricultural and urban sources and which were therefore sensitive to reductions of atmospheric nitrogen depositions. We conclude that there is strong evidence that the decline in nitrogen deposition toward pre-industrial conditions lowers the redox buffer in riparian soils, destabilizing formerly fixed problematic compounds, and results in serious implications for water quality.
Due to the challenges in upscaling daily climatic forcing to geological time, physically realistic models describing how rainfall drives fluvial erosion are lacking. To bridge this gap between short-term hydrology and long-term geomorphology, we derive a theoretical framework for long-term fluvial erosion rates driven by realistic climate by integrating an established stochastic-mechanistic model of hydrology into a threshold-stochastic formulation of stream power. The hydrological theory provides equations for the daily streamflow probability distribution as a function of climatic boundary conditions. The new parameters introduced are rooted firmly in established climatic and hydrological theory. This allows us to account for how fluvial erosion rates respond to changes in rainfall intensity, frequency, evapotranspiration rates, and soil moisture dynamics in a way that is consistent with existing theories. We use this framework to demonstrate how hydroclimatic conditions and erosion threshold magnitude control the degree of nonlinearity between steepness index and erosion rate. We find that hydrological processes can have a significant influence on how erosive a particular climatic forcing will be. By accounting for the influence of hydrology on fluvial erosion, we conclude that climate is an important control on erosion rates and long-term landscape evolution.
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
Understanding the geologic evolution of Northern Tibetan Plateau with multiple thermochronometers
(2018)
The early onset of deformation following the India-Asia collision, Neogene expanse of uplift, and complex systems that comprise strike-slip faults, thrust faults, and intermontane basins characterize the Cenozoic tectonism of Northern Tibetan Plateau and raise two prominent questions in orogenic geodynamics: 1) What mechanism(s) control(s) the transfer of stress related to the India-Asia collision across the distance of >2000 km; and 2) Why the development of high topography was delayed in the Northern Tibetan Plateau and what does it reveal about how the internal forces and external boundary conditions evolved. To address these two questions, we reconstruct a holistic spatial-temporal deformation history of the Northern Tibetan Plateau by using a range of thermochronometers, with closure temperature spanning from 350 degrees C to-60-70 degrees C. This multi-thermochronometer study reveals three stages of faulting related cooling, in the early Cretaceous, in Paleocene-Eocene and in middle-late Miocene. We observe that Paleocene-Eocene deformation was spatially restricted and mostly occurred on reactivated Cretaceous structures, indicating a control of pre-existing weakness on early Cenozoic deformation. Extensive Neogene deformation contrasts with restricted Paleocene-Eocene deformation and relatively quiescent shortening during the Oligocene-early Miocene, which implies a change in the regional tectonics regime. Global plate reconstructions show that this tectonic reorganization is coeval with an increase in Pacific-Asia plate convergence rates. We argue that this change in regional tectonics is a result of increasing constrictive environment of the eastern plate boundary, which changed the behavior of the Altyn Tagh fault the boundary fault of Northern Tibetan Plateau, causing it to change from feeding slip into structures out of the plateau to feeding slip into structures at plateau margins.
Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30–50%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ S2S (20–30%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset’s advantages is the large number of recordings per station (9–90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock.
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.
Statistical distributions of flood peak discharge often show heavy tail behavior, that is, extreme floods are more likely to occur than would be predicted by commonly used distributions that have exponential asymptotic behavior.
This heavy tail behavior may surprise flood managers and citizens, as human intuition tends to expect light tail behavior, and the heaviness of the tails is very difficult to predict, which may lead to unnecessarily high flood damage.
Despite its high importance, the literature on the heavy tail behavior of flood distributions is rather fragmented.
In this review, we provide a coherent overview of the processes causing heavy flood tails and the implications for science and practice.
Specifically, we propose nine hypotheses on the mechanisms causing heavy tails in flood peak distributions related to processes in the atmosphere, the catchment, and the river system.
We then discuss to which extent the current knowledge supports or contradicts these hypotheses.
We also discuss the statistical conditions for the emergence of heavy tail behavior based on derived distribution theory and relate them to the hypotheses and flood generation mechanisms.
We review the degree to which the heaviness of the tails can be predicted from process knowledge and data. Finally, we recommend further research toward testing the hypotheses and improving the prediction of heavy tails.
Understanding the rates and pattern of erosion is a key aspect of deciphering the impacts of climate and tectonics on landscape evolution. Denudation rates derived from terrestrial cosmogenic nuclides (TCNs) are commonly used to quantify erosion and bridge tectonic (Myr) and climatic (up to several kiloyears) time scales. However, how the processes of erosion in active orogens are ultimately reflected in Be-10 TCN samples remains a topic of discussion. We investigate this problem in the Arun Valley of eastern Nepal with 34 new Be-10-derived catchment-mean denudation rates. The Arun Valley is characterized by steep north-south gradients in topography and climate. Locally, denudation rates increase northward, from <0.2mmyr(-1) to similar to 1.5mmyr(-1) in tributary samples, while main stem samples appear to increase downstream from similar to 0.2mmyr(-1) at the border with Tibet to 0.91mmyr(-1) in the foreland. Denudation rates most strongly correlate with normalized channel steepness (R-2=0.67), which has been commonly interpreted to indicate tectonic activity. Significant downstream decrease of Be-10 concentration in the main stem Arun suggests that upstream sediment grains are fining to the point that they are operationally excluded from the processed sample. This results in Be-10 concentrations and denudation rates that do not uniformly represent the upstream catchment area. We observe strong impacts on Be-10 concentrations from local, nonfluvial geomorphic processes, such as glaciation and landsliding coinciding with areas of peak rainfall rates, pointing toward climatic modulation of predominantly tectonically driven denudation rates.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.
Saturated hydraulic conductivity (K-s) is an important soil characteristic affecting soil water storage, runoff generation and erosion processes. In some areas where high-intensity rainfall coincides with low K-s values at shallow soil depths, frequent overland flow entails dense drainage networks. Consequently, linear structures such as flowlines alternate with inter-flowline areas. So far, investigations of the spatial variability of K-s mainly relied on isotropic covariance models which are unsuitable to reveal patterns resulting from linear structures. In the present study, we applied two sampling approaches so as to adequately characterize K-s spatial variability in a tropical forest catchment that features a high density of flowlines: A classical nested sampling survey and a purposive sampling strategy adapted to the presence of flowlines. The nested sampling approach revealed the dominance of small-scale variability, which is in line with previous findings. Our purposive sampling, however, detected a strong spatial gradient: surface K-s increased substantially as a function of distance to flowline; 10 m off flowlines, values were similar to the spatial mean of K-s. This deterministic trend can be included as a fixed effect in a linear mixed modeling framework to obtain realistic spatial fields of K-s. In a next step we used probability maps based on those fields and prevailing rainfall intensities to assess the hydrological relevance of the detected pattern. This approach suggests a particularly good agreement between the probability statements of K-s exceedance and observed overland flow occurrence during wet stages of the rainy season.
The Late Permian Zechstein Group in northeastern Germany is characterized by shelf and slope carbonates that rimmed a basin extending from eastern England through the Netherlands and Germany to Poland. Conventional reservoirs are found in grainstones rimming islands created by pre-existing paleohighs and platform-rimming shoals that compose steep margins in the north and ramp deposits in the southern part. The slope and basin deposits are characterized by debris flows and organic-rich mudstones. Lagoonal and basinal evaporites formed the seal for these carbonate and underlying sandstone reservoirs. The objective of this investigation is to evaluate potential unconventional reservoirs in organic-rich, fine-grained and/or tight mudrocks in slope and basin as well as platform carbonates occurring in this stratigraphic interval. Therefore, a comprehensive study was conducted that included sedimentology, sequence stratigraphy, petrography, and geochemistry. Sequence stratigraphic correlations from shelf to basin are crucial in establishing a framework that allows correlation of potential productive facies in fine-grained, organic-rich basinal siliceous and calcareous mudstones or interfingering tight carbonates and siltstones, ranging from the lagoon, to slope to basin, which might be candidates for forming an unconventional reservoir. Most organic-rich shales worldwide are associated with eustatic transgressions. The basal Zechstein cycles, Z1 and Z2, contain organic-rich siliceous and calcareous mudstones and carbonates that form major transgressive deposits in the basin. Maturities range from over-mature (gas) in the basin to oil-generation on the slope with variable TOC contents. This sequence stratigraphic and sedimentologic evaluation of the transgressive facies in the Z1 and Z2 assesses the potential for shale-gas/oil and hybrid unconventional plays. Potential unconventional reservoirs might be explored in laminated organic-rich mudstones within the oil window along the northern and southern slopes of the basin. Although the Zechstein Z1 and Z2 cycles might have limited shale-gas potential because of low thickness and deep burial depth to be economic at this point, unconventional reservoir opportunities that include hybrid and shale-oil potential are possible in the study area.
Shrinking glaciers in the Hindu Kush-Karakoram-Himalaya-Nyainqentanglha (HKKHN) region have formed several thousand moraine-dammed glacial lakes(1-3), some of these having grown rapidly in past decades(3,4). This growth may promote more frequent and potentially destructive glacial lake outburst floods (GLOFs)(5-7). Testing this hypothesis, however, is confounded by incomplete databases of the few reliable, though selective, case studies. Here we present a consistent Himalayan GLOF inventory derived automatically from all available Landsat imagery since the late 1980s. We more than double the known GLOF count and identify the southern Himalayas as a hotspot region, compared to the more rarely affected Hindu Kush-Karakoram ranges. Nevertheless, the average annual frequency of 1.3 GLOFs has no credible posterior trend despite reported increases in glacial lake areas in most of the HKKHN3,8, so that GLOF activity per unit lake area has decreased since the late 1980s. We conclude that learning more about the frequency and magnitude of outburst triggers, rather than focusing solely on rapidly growing glacial lakes, might improve the appraisal of GLOF hazards.
Cosmic-ray neutron sensing (CRNS) is a promising proximal soil sensing technique to estimate soil moisture at intermediate scale and high temporal resolution. However, the signal shows complex and non-unique response to all hydrogen pools near the land surface, providing some challenges for soil moisture estimation in practical applications. Aims of the study were 1) to assess the uncertainty of CRNS as a stand-alone approach to estimate volumetric soil moisture in cropped field 2) to identify the causes of this uncertainty 3) and possible improvements. Two experimental sites in Germany were equipped with a CRNS probe and point-scale soil moisture network. Additional monitoring activities were conducted during the crop growing season to characterize the soil-plant systems. This data is used to identify and quantify the different sources of uncertainty (factors). An uncertainty analysis, based on Monte Carlo approach, is applied to propagate these uncertainties to CRNS soil moisture estimations. In addition, a sensitivity analysis based on the Sobol’ method is performed to identify the most important factors explaining this uncertainty. Results show that CRNS soil moisture compares well to the soil moisture network when these point-scale values are weighted to account for the spatial sensitivity of the signal and other sources of hydrogen (lattice water and organic carbon) are added to the water content. However, the performance decreases when CRNS is considered as a stand-alone method to retrieve the actual (non-weighted) volumetric soil moisture. The support volume (penetration depth and radius) shows also a considerable uncertainty, especially in relatively dry soil moisture conditions. Four of the seven factors analyzed (the vertical soil moisture profile, bulk density, incoming neutron correction and the calibrated parameter N0) were found to play an important role. Among the possible improvements identified, a simple correction factor based on vertical point-scale soil moisture profiles shows to be a promising approach to account for the sensitivity of the CRNS signal to the upper soil layers.
The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Empirical species distribution models (SDMs) constitute often the tool of choice for the assessment of rapid climate change effects on species vulnerability. Conclusions regarding extinction risks might be misleading, however, because SDMs do not explicitly incorporate dispersal or other demographic processes. Here, we supplement SDMs with a dynamic population model 1) to predict climate-induced range dynamics for black grouse in Switzerland, 2) to compare direct and indirect measures of extinction risks, and 3) to quantify uncertainty in predictions as well as the sources of that uncertainty. To this end, we linked models of habitat suitability to a spatially explicit, individual-based model. In an extensive sensitivity analysis, we quantified uncertainty in various model outputs introduced by different SDM algorithms, by different climate scenarios and by demographic model parameters. Potentially suitable habitats were predicted to shift uphill and eastwards. By the end of the 21st century, abrupt habitat losses were predicted in the western Prealps for some climate scenarios. In contrast, population size and occupied area were primarily controlled by currently negative population growth and gradually declined from the beginning of the century across all climate scenarios and SDM algorithms. However, predictions of population dynamic features were highly variable across simulations. Results indicate that inferring extinction probabilities simply from the quantity of suitable habitat may underestimate extinction risks because this may ignore important interactions between life history traits and available habitat. Also, in dynamic range predictions uncertainty in SDM algorithms and climate scenarios can become secondary to uncertainty in dynamic model components. Our study emphasises the need for principal evaluation tools like sensitivity analysis in order to assess uncertainty and robustness in dynamic range predictions. A more direct benefit of such robustness analysis is an improved mechanistic understanding of dynamic species responses to climate change.
Tillage erosion on arable land is a very important process leading to a net downslope movement of soil and soil constitutes. Tillage erosion rates are commonly in the same order of magnitude as water erosion rates and can be even higher, especially under highly mechanized agricultural soil management. Despite its prevalence and magnitude, tillage erosion is still understudied compared to water erosion. The goal of this study was to bring together experts using different techniques to determine tillage erosion and use the different results to discuss and quantify uncertainties associated with tillage erosion measurements. The study was performed in northeastern Germany on a 10 m by 50 m plot with a mean slope of 8%. Tillage erosion was determined after two sequences of seven tillage operations. Two different micro-tracers (magnetic iron oxide mixed with soil and fluorescent sand) and one macro-tracer (passive radio-frequency identification transponders (RFIDs), size: 4 x 22 mm) were used to directly determine soil fluxes. Moreover, tillage induced changes in topography were measured for the entire plot with two different terrestrial laser scanners and an unmanned aerial system for structure from motion topography analysis. Based on these elevation differences, corresponding soil fluxes were calculated. The mean translocation distance of all techniques was 0.57 m per tillage pass, with a relatively wide range of mean soil translocation distances ranging from 039 to 0.72 m per pass. A benchmark technique could not be identified as all used techniques have individual error sources, which could not be quantified. However, the translocation distances of the macro-tracers used were consistently smaller than the translocation distances of the micro-tracers (mean difference = -26 +/- 12%), which questions the widely used assumption of non-selective soil transport via tillage operations. This study points out that tillage erosion measurements, carried out under almost optimal conditions, are subject to major uncertainties that are far from negligible. (C) 2018 Elsevier B.V. All rights reserved.
The advantages of remote sensing using Unmanned Aerial Vehicles (UAVs) are a high spatial resolution of images, temporal flexibility and narrow-band spectral data from different wavelengths domains. This enables the detection of spatio-temporal dynamics of environmental variables, like plant-related carbon dynamics in agricultural landscapes. In this paper, we quantify spatial patterns of fresh phytomass and related carbon (C) export using imagery captured by a 12-band multispectral camera mounted on the fixed wing UAV Carolo P360. The study was performed in 2014 at the experimental area CarboZALF-D in NE Germany. From radiometrically corrected and calibrated images of lucerne (Medicago sativa), the performance of four commonly used vegetation indices (VIs) was tested using band combinations of six near-infrared bands. The highest correlation between ground-based measurements of fresh phytomass of lucerne and VIs was obtained for the Enhanced Vegetation Index (EVI) using near-infrared band b(899). The resulting map was transformed into dry phytomass and finally upscaled to total C export by harvest. The observed spatial variability at field- and plot-scale could be attributed to small-scale soil heterogeneity in part.
The age of Proterozoic granulite facies metamorphism and deformation in the Strangways Metamorphic Complex (SMC) of central Australia is determined on zircon grown in syn-metamorphic and syn-deformational orthopyroxene-bearing, enderbitic, veins. SHRIMP zircon studies suggest that M1-M2 and the correlated periods of intense deformation (D1-D2) are part of a single tectonothermal event between 1,717-2 and 1,732-7 Ma. It is considered unlikely that the two metamorphic phases (M1, M2) suggested by earlier work represent separate events occurring within 10-25 Ma of each other. Previous higher estimates for the age of M1 granulite metamorphism in the SMC (Early Strangways event at ca. 1,770 Ma) based on U-Pb zircon dating of granitic, intrusive rocks, are not believed to relate to the metamorphism, but to represent pre-metamorphic intrusion ages. Conventional multi-grain U-Pb monazite analyses on high-grade metasediments from three widely spaced localities in the western SMC yield 207Pb/235U ages between 1,728-11 and 1,712-2 Ma. The age range of the monazites corresponds to the SHRIMP zircon ages in the granulitic veins and is interpreted to record monazite growth (prograde in the metasedimentary rocks). The data imply a maximum time-span of 30 Ma for high-grade metamorphism and deformation in the SMC. There is, thus, no evidence for an extremely long period of continuous high- temperature conditions from 1,770 to ca. 1,720 Ma as previously proposed. The results firmly establish that the SMC has a very different high-grade metamorphic history than the neighbouring Harts Range, where upper amphibolite facies metamorphism in the Palaeozoic caused widespread growth or recrystallization of monazite.
The type-locality granulites from the Granulitgebirge of Saxony, Germany, are rocks of broadly granitic composition containing minor garnet and kyanite within a commonly mylonitised matrix of feldspars and quartz. Petrographic evidence indicates a primary assemblage of ternary feldspar + quartz + garnet + kyanite + rutile, most likely resulting from partial melting of a granitic protolith, for which equilibrium temperature and pressure conditions of > 1000 degrees C and > 1.5 GPa have been deduced. These extreme (for crustal rocks) conditions, and the inferred peak assemblage, are supported by the newly-developed Zr-in-rutile geothermometer and experimental studies on the same bulk composition, respectively. As these conditions lie above those required for plagioclase stability in quartz tholeiites, they are thus in the eclogite facies. Widespread modification of the peak assemblage, for example mesoperthite formation after ternary feldspar, deformation-induced recrystallisation of perthites to two-feldspar + quartz aggregates, biotite replacing garnet, Ca-loss at garnet rims, sillimanite replacing kyanite or secondary garnet growth, makes reliable interpretation of equilibrium assemblages and compositions very difficult and explains the spread of published pressure- temperature values and consequent confusion about formation depths and the validity of tectonometamorphic models. Such extreme metamorphic conditions in rock compositions typical for the upper continental crust, reflecting a hot subduction environment, has important consequences for understanding some collisional orogens
TURBO2 - a MATLAB simulation to study the effects of bioturbation on paleoceanographic time series
(2013)
Bioturbation (or benthic mixing) causes significant distortions in marine stable isotope signals and other palaeoceanographic records. Although the influence of bioturbation on these records is well known it has rarely been dealt systematically. The MATLAB program called TURBO2 can be used to simulate the effect of bioturbation on individual sediment particles. It can therefore be used to model the distortion of all physical, chemical, and biological signals in deep-sea sediments, such as Mg/Ca ratios and UK37-based sea-surface temperature (SST) variations. In particular, it can be used to study the distortions in paleoceanographic records that are based on individual sediment particles, such as SST records based on foraminifera assemblages. Furthermore. TURBO2 provides a tool to study the effect of benthic mixing of isotope signals such as C-14, delta O-18, and delta C-13, measured in a stratigraphic carrier such as foraminifera shells.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Process life cycle assessment (PLCA) is widely used to quantify environmental flows associated with the manufacturing of products and other processes. As PLCA always depends on defining a system boundary, its application involves truncation errors. Different methods of estimating truncation errors are proposed in the literature; most of these are based on artificially constructed system complete counterfactuals. In this article, we review the literature on truncation errors and their estimates and systematically explore factors that influence truncation error estimates. We classify estimation approaches, together with underlying factors influencing estimation results according to where in the estimation procedure they occur. By contrasting different PLCA truncation/error modeling frameworks using the same underlying input-output (I-O) data set and varying cut-off criteria, we show that modeling choices can significantly influence estimates for PLCA truncation errors. In addition, we find that differences in I-O and process inventory databases, such as missing service sector activities, can significantly affect estimates of PLCA truncation errors. Our results expose the challenges related to explicit statements on the magnitude of PLCA truncation errors. They also indicate that increasing the strictness of cut-off criteria in PLCA has only limited influence on the resulting truncation errors. We conclude that applying an additional I-O life cycle assessment or a path exchange hybrid life cycle assessment to identify where significant contributions are located in upstream layers could significantly reduce PLCA truncation errors.
The remaining carbon stocks in wet tropical forests are currently at risk because of anthropogenic deforestation, but also because of the possibility of release driven by climate change. To identify the relative roles of CO2 increase, changing temperature and rainfall, and deforestation in the future, and the magnitude of their impact on atmospheric CO2 concentrations, we have applied a dynamic global vegetation model, using multiple scenarios of tropical deforestation (extrapolated from two estimates of current rates) and multiple scenarios of changing climate (derived from four independent offline general circulation model simulations). Results show that deforestation will probably produce large losses of carbon, despite the uncertainty about the deforestation rates. Some climate models produce additional large fluxes due to increased drought stress caused by rising temperature and decreasing rainfall. One climate model, however, produces an additional carbon sink. Taken together, our estimates of additional carbon emissions during the twenty-first century, for all climate and deforestation scenarios, range from 101 to 367 Gt C, resulting in CO2 concentration increases above background values between 29 and 129 p.p.m. An evaluation of the method indicates that better estimates of tropical carbon sources and sinks require improved assessments of current and future deforestation, and more consistent precipitation scenarios from climate models. Notwithstanding the uncertainties, continued tropical deforestation will most certainly play a very large role in the build-up of future greenhouse gas concentrations
Triassic Latemar cycle tops - Subaerial exposure of platform carbonates under tropical arid climate
(2012)
The Triassic Latemar platform in the Dolomites, Italy, is the site of several ongoing controversies. Perhaps the most interesting debate focuses on apparent cyclic deposition within the Latemar platform, whose nature and duration are still open to debate. Further disagreement concerns the lack of meteoric diagenesis-related isotope shifts at cycle tops that bear circumstantial petrographic evidence for subaerial emergence. Here, an evaluation of the nature of Latemar cycle tops is presented combining evidence from previous work and new field, petrographic and geochemical data. Cycle tops are ranked according to increasing exposure duration and spatial extent: type I surfaces lacking unequivocal evidence of prolonged supratidal conditions; type II dolomite caps formed in warm, evaporitic, intertidal lagoonal waters followed by exposure of perhaps intermediate duration; type III clastic-rich, red calcareous horizons with some showing platform-wide extent, representing prolonged supratidal conditions, and type IV discontinuities in tepee belts, genetically related to type II and III surfaces, but likely representing shorter-lived exposure stages. Petrographic and geochemical criteria indicate that most diagenesis occurred in the shallow marine and burial domain whilst an extensive meteoric overprint of cycle tops is lacking. This is underlined by the scarcity of meteoric diagenetic fabrics such as gravitational cements that, where present, are here interpreted as marine-vadose in origin. The scarcity of carbon and oxygen isotope signatures commonly assigned to subaerial exposure stages is best explained in the context of mid-Triassic climate. The low latitude, tropical but arid setting of the Latemar, situated in the western extension of the Tethys ocean, its isolation from nearby continental areas and overall short-term emergence episodes are in agreement with a limited degree of meteoric alteration of most cycle tops. High amounts of aeolian clastic material beneath some cycle tops, along with high Fe and Mn elemental abundances argue for intermittent subaerial conditions. This study proposes an enhancement of the classical Allan and Matthews (1982) isotope model for subaerial exposure under strongly arid climates. As the subaerial exposure nature of Latemar cycle tops, and therefore eustasy as the cause for cyclicity, have been previously challenged due to the lack of meteoric-induced isotopic signatures, the outcome of this study is of significance for the ongoing Latemar stratigraphic controversy.
We analyzed published records of terrigenous dust flux from marine sediments off subtropical West Africa, the eastern Mediterranean Sea, and the Arabian Sea, and lake records from East Africa using statistical methods to detect trends, rhythms and events in Plio-Pleistocene African climate. The critical reassessment of the environmental significance of dust flux and lake records removes the apparent inconsistencies between marine vs. terrestrial records of African climate variability. Based on these results, major steps in mammalian and hominin evolution occurred during episodes of a wetter, but highly variable climate largely controlled by orbitally induced insolation changes in the low latitudes.
Thousands of glacier lakes have been forming behind natural dams in high mountains following glacier retreat since the early 20th century. Some of these lakes abruptly released pulses of water and sediment with disastrous downstream consequences. Yet it remains unclear whether the reported rise of these glacier lake outburst floods (GLOFs) has been fueled by a warming atmosphere and enhanced meltwater production, or simply a growing research effort. Here we estimate trends and biases in GLOF reporting based on the largest global catalog of 1,997 dated glacier-related floods in six major mountain ranges from 1901 to 2017. We find that the positive trend in the number of reported GLOFs has decayed distinctly after a break in the 1970s, coinciding with independently detected trend changes in annual air temperatures and in the annual number of field-based glacier surveys (a proxy of scientific reporting). We observe that GLOF reports and glacier surveys decelerated, while temperature rise accelerated in the past five decades. Enhanced warming alone can thus hardly explain the annual number of reported GLOFs, suggesting that temperature-driven glacier lake formation, growth, and failure are weakly coupled, or that outbursts have been overlooked. Indeed, our analysis emphasizes a distinct geographic and temporal bias in GLOF reporting, and we project that between two to four out of five GLOFs on average might have gone unnoticed in the early to mid-20th century. We recommend that such biases should be considered, or better corrected for, when attributing the frequency of reported GLOFs to atmospheric warming.
We analyze trends in compound flooding resulting from high coastal water levels (HCWLs) and peak river discharge over northwestern Europe during 1901-2014. Compound peak discharge associated with 37 stream gauges with at least 70 years of record availability near the North and Baltic Sea coasts is used. Compound flooding is assessed using a newly developed index, compound hazard ratio, that compares the severity of river flooding associated with HCWL with the at-site, T-year (a flood with 1/T chance of being exceeded in any given year) fluvial peak discharge. Our findings suggest a spatially coherent pattern in the dependence between HCWL and river peaks and in compound flood magnitudes and frequency. For higher return levels, we find upward trends in compound hazard ratio frequency at midlatitudes (gauges from 47 degrees N to 60 degrees N) and downward trends along the high latitude (>60 degrees N) regions of northwestern Europe. Plain Language Summary Compound floods in delta areas, that is, the co-occurrence of high coastal water levels (HCWLs) and high river discharge, are a particular challenge for disaster management. Such events are caused by two distinct mechanisms: (1) HCWLs may affect river flows and water levels by backwater effects or by reversing the seaward flow of rivers, particularly in regions with elevation less than 10 m in northwestern Europe. (2) The correlation between HCWL and river flow peaks may also stem from a common meteorological driver. Severe storm periods may be associated with high winds leading to storm surges, and at the same time with high precipitation followed by inland flooding. Understanding the historical trends in compound flooding, owing to changes in relative sea levels, in river flooding and in the dependence between these two drivers, is essential for projecting future changes and disaster management. The risk assessment frameworks are often limited to assessing flood risk from a single driver only. We present a new approach to assess compound flood severity resulting from extreme coastal water level and peak river discharge. We find upward trends in compound flooding for midlatitude regions and downward trends for high latitudes in northwestern Europe.
Reliable flood risk analyses, including the estimation of damage, are an important prerequisite for efficient risk management. However, not much is known about flood damage processes affecting companies. Thus, we conduct a flood damage assessment of companies in Germany with regard to two aspects. First, we identify relevant damage-influencing variables. Second, we assess the prediction performance of the developed damage models with respect to the gain by using an increasing amount of training data and a sector-specific evaluation of the data. Random forests are trained with data from two postevent surveys after flood events occurring in the years 2002 and 2013. For a sector-specific consideration, the data set is split into four subsets corresponding to the manufacturing, commercial, financial, and service sectors. Further, separate models are derived for three different company assets: buildings, equipment, and goods and stock. Calculated variable importance values reveal different variable sets relevant for the damage estimation, indicating significant differences in the damage process for various company sectors and assets. With an increasing number of data used to build the models, prediction errors decrease. Yet the effect is rather small and seems to saturate for a data set size of several hundred observations. In contrast, the prediction improvement achieved by a sector-specific consideration is more distinct, especially for damage to equipment and goods and stock. Consequently, sector-specific data acquisition and a consideration of sector-specific company characteristics in future flood damage assessments is expected to improve the model performance more than a mere increase in data.
Vertical flow filters and vertical flow constructed wetlands are established wastewater treatment systems and have also been proposed for the treatment of contaminated groundwater. This study investigates the removal processes of volatile organic compounds in a pilot-scale vertical flow filter. The filter is intermittently irrigated with contaminated groundwater containing benzene, MTBE and ammonium as the main contaminants. The system is characterized by unsaturated conditions and high contaminant removal efficiency. The aim of the present study is to evaluate the contribution of biodegradation and volatilization to the overall removal of benzene and MTBE. Tracer tests and flow rate measurements showed a highly transient flow and heterogeneous transport regime. Radon-222, naturally occurring in the treated groundwater, was used as a gas tracer and indicated a high volatilization potential. Radon-222 behavior was reproduced by numerical simulations and extrapolated for benzene and MTBE, and indicated these compounds also have a high volatilization potential. In contrast, passive sampler measurements on top of the filter detected only low benzene and MTBE concentrations. Biodegradation potential was evaluated by the analysis of catabolic genes involved in organic compound degradation and a quantitative estimation of biodegradation was derived from stable isotope fractionation analysis. Results suggest that despite the high volatilization potential, biodegradation is the predominant mass removal process in the filter system, which indicates that the volatilized fraction of the contaminants is still subject to subsequent biodegradation. In particular, the upper filter layer located between the injection tubes and the surface of the system might also contribute to biodegradation, and might play a crucial role in avoiding the emission of volatilized contaminants into the atmosphere.
Traveltime residuals for worldwide seismic stations are calculated. We use P and S waves from earthquakes in SE-Asia at teleseismic and regional distances. The obtained station residuals help to enhance earthquake localisation. Furthermore we calculated regional source dependent station residuals. They show a systematic dependence of the locality of the source. These source dependent residuals reflect heterogenities along the path and can be used for a refinement of earthquake localisation.
On 7 January 2020, an M-w 6.4 earthquake occurred in the northeastern Caribbean, a few kilometers offshore of the island of Puerto Rico. It was the mainshock of a complex seismic sequence, characterized by a large number of energetic earthquakes illuminating an east-west elongated area along the southwestern coast of Puerto Rico. Deformation fields constrained by Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data indicate that the coseismic movements affected only the western part of the island. To assess the mainshock's source fault parameters, we combined the geodetically derived coseismic deformation with teleseismic waveforms using Bayesian inference. The results indicate a roughly east-west oriented fault, dipping northward and accommodating similar to 1.4 m of transtensional motion. Besides, the determined location and orientation parameters suggest an offshore continuation of the recently mapped North Boqueron Bay-Punta Montalva fault in southwest Puerto Rico. This highlights the existence of unmapped faults with moderate-to-large earthquake potential within the Puerto Rico region.
Different tectonic episodes from Late Triassic to recent times in the eastern Binalud Mountains have resulted from convergence and transpression between the Turan and Central Iran plates. Heterogeneous deformation and variable portions of pure and simple shear, demonstrated by finite strain and vorticity analysis in the Mashhad metamorphic rocks, indicate strain partitioning during the first tectonic episode. Modern strain partitioning is characterized by reverse and strike-slip faulting along the Neyshabur fault system and Shandiz fault zone in the southern and northern flanks of the eastern Binalud, respectively. Time-transgressive regional deformation migrated from the hinterland of the belt into the foreland basin, from northeast to the southwest of the mountains. Different generations of deformation resulted in obliteration of the subduction-related accretionary wedge, and growth of an orogenic wedge resulted from collision between the Central Iran and Turan plates.
Injection of nanoscale zero-valent iron (nZVI) has recently gained great interest as emerging technology for in-situ remediation of chlorinated organic compounds from groundwater systems. Zero-valent iron (ZVI) is able to reduce organic compounds and to render it to less harmful substances. The use of nanoscale particles instead of granular or microscale particles can increase dechlorination rates by-orders of magnitude due to its high surface area. However, classical nZVI appears to be hampered in its environmental application by its limited mobility. One approach is colloid supported transport of nZVI, where the nZVI gets transported by a Mobile colloid. In this study transport properties of activated carbon colloid supported nZVI (c-nZVI; d(50) = 2.4 mu m) are investigated in column tests using columns of 40 cm length, which were filled with porous media. A suspension was pumped through the column under different physicochemical conditions (addition of a polyanionic stabilizer and changes in pH and ionic strength). Highest observed breakthrough was 62% of the injected concentration in glass beads with addition of stabilizer. Addition of mono- and bivalent salt, e.g. more than 0.5 mM/L CaCl2, can decrease mobility and changes in pH to values below six can inhibit mobility at all. Measurements of colloid sizes and zeta potentials show changes in the mean particle size by a factor of ten and an increase of zeta potential from -62 mV to -80 mV during the transport experiment. However, results suggest potential applicability of c-nZVI under field conditions. (C) 2014 Elsevier B.V. All rights reserved.
Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.
A geological transect across the suture separating northwestern South America from the Panama Arc helps document the provenance and thermal history of both crustal domains and the suture zone. During middle Miocene, strata were being accumulated over the suture zone between the Panama Arc and the continental margin. Integrated provenance analyses of those middle Miocene strata show the presence of mixed sources that includes material derived from the two major crustal domains: the old northwestern South American orogens and the younger Panama Arc. Coeval moderately rapid exhumation of Upper Cretaceous to Paleogene sediments forming the reference continental margin is suggested from our inverse thermal modeling. Strata within the suture zone are intruded by similar to 12 Ma magmatic arc-related plutons, marking the transition from a collisional orogen to a subduction-related one. Renewed late Miocene to Pliocene acceleration of the exhumation rates is the consequence of a second tectonic pulse, which is likely to be triggered by the onset of a flat-slab subduction of the Nazca plate underneath the northernmost Andes of Colombia, suggesting that late Miocene to Pliocene orogeny in the Northern Andes is controlled by at least two different tectonic mechanisms.
In this paper we report on a series of replicated tracer experiments with deuterium conducted under controlled, steady stormflow conditions at the Gardsjon G1 catchment in south-western Sweden. In five different years, these experiments were carried out in a subcatchment of G1. Deuterium was applied as a narrow pulse so that distributions of water transit times could be directly inferred from the observed tracer breakthrough curves. Significantly different transit times of water were observed under similar experimental conditions. Coefficients of variation for estimated mean transit times were greater than 60%, which can be understood as a measure of the interannual variability for this type of experiments. Implications for water transit times under more natural flow conditions as wells as for future experimentation are discussed. Copyright (c) 2015 John Wiley & Sons, Ltd.
Earthquake-triggered landslide dams are potentially dangerous disrupters of water and sediment flux in mountain rivers, and capable of releasing catastrophic outburst flows to downstream areas. We analyze an inventory of 828 landslide dams in the Longmen Shan mountains, China, triggered by the M-w 7.9 2008 Wenchuan earthquake. This database is unique in that it is the largest of its kind attributable to a single regional-scale triggering event: 501 of the spatially clustered landslides fully blocked rivers, while the remainder only partially obstructed or diverted channels in steep watersheds of the hanging wall of the Yingxiu-Beichuan Fault Zone. The size distributions of the earthquake-triggered landslides, landslide dams, and associated lakes (a) can be modeled by an inverse gamma distribution; (b) show that moderate-size slope failures caused the majority of blockages; and (c) allow a detailed assessment of seismically induced river-blockage effects on regional water and sediment storage. Monte Carlo simulations based on volumetric scaling relationships for soil and bedrock failures respectively indicate that 14% (18%) of the estimated total coseismic landslide volume of 6.4 (14.6) x 10(9) m(3) was contained in landslide dams, representing only 1.4% of the >60,000 slope failures attributed to the earthquake. These dams have created storage capacity of similar to 0.6x 10(9) m(3) for incoming water and sediment. About 25% of the dams containing 2% of the total river-blocking debris volume failed one week after the earthquake; these figures had risen to 60% (similar to 20%), and >90% (>90%) within one month, and one:year, respectively, thus also emptying similar to 92% of the total potential water and sediment storage behind these, dams within one year following the earthquake. Currently only similar to 0.08 x 10(9) m(3) remain available as natural reservoirs for storing water and sediment, while similar to 0.19 x 10(9) m(3), i.e. about a third of the total river-blocking debris volume, has been eroded by rivers. Dam volume and upstream catchment area control to first order the longevity of the barriers, and bivariate domain plots are consistent with the observation that most earthquake-triggered landslide dams were ephemeral. We conclude that the river-blocking portion of coseismic slope failures disproportionately modulates the post-seismic sediment flux in the Longmen Shan on annual to decadal timescales.
Thermokarst lakes are prevalent in Arctic coastal lowland regions and sublake permafrost degradation and talik development contributes to greenhouse gas emissions by tapping the large permafrost carbon pool. Whereas lateral thermokarst lake expansion is readily apparent through remote sensing and shoreline measurements, sublake thawed sediment conditions and talik growth are difficult to measure. Here we combine transient electromagnetic surveys with thermal modeling, backed up by measured permafrost properties and radiocarbon ages, to reveal closed-talik geometry associated with a thermokarst lake in continuous permafrost. To improve access to talik geometry data, we conducted surveys along three transient electromagnetic transects perpendicular to lakeshores with different decadal-scale expansion rates of 0.16, 0.38, and 0.58m/year. We modeled thermal development of the talik using boundary conditions based on field data from the lake, surrounding permafrost and a borehole, independent of the transient electromagnetics. A talik depth of 91m was determined from analysis of the transient electromagnetic surveys. Using a lake initiation age of 1400years before present and available subsurface properties the results from thermal modeling of the lake center arrived at a best estimate talk depth of 80m, which is on the same order of magnitude as the results from the transient electromagnetic survey. Our approach has provided a noninvasive estimate of talik geometry suitable for comparable settings throughout circum-Arctic coastal lowland regions.
Earthquakes impart an impressive force on epicentral landscapes, with immediate catastrophic hillslope response. However, their legacy on geomorphic process rates remains poorly constrained. We have determined the evolution of landslide rates in the epicentral areas of four intermediate to large earthquakes (M-w, 6.6-7.6). In each area, landsliding correlates with the cumulative precipitation during a given interval. Normalizing for this meteorological forcing, landslide rates have been found to peak after an earthquake and decay to background values in 1-4 yr, with the decay time scale probably proportional to the earthquake magnitude. The transient pulse of landsliding is not related to external forcing such as rainfall or aftershocks, and we tentatively attribute it to the reduction and subsequent recovery of ground strength. Observed geomorphic trends are not linked with groundwater level changes or root system damage, both of which could affect substrate strength. We propose that they are caused by reversible damage of rock mass and/or loosening of regolith. Qualitative accounts of ground cracking due to strong ground motion abound, and our observations are circumstantial evidence of its potential importance in setting landscape sensitivity to meteorological forcing after large earthquakes.
Tomato (Solanum lycopersicum L.) being an important vegetable is cultivated and used throughout the world. It not only contributes in fulfilling the basic nutritional requirements of the human body but also has many health benefits due to its rich biochemical composition. However, its production at large scale is hampered by many limiting factors such as biotic and abiotic stresses. Among the different abiotic stresses, drought poses drastic impact on tomato yield. Drought stress is genetically regulated by many transcription factors that not only regulate the stress responsive mechanism but also facilitate the growth and development of tomato plants. NAC is an important stress related transcription factor genes family, and the ATAF1 gene, a member of this family, is involved in ABA signaling and stress response. In this study, tomato variety Rio Drande was transformed with drought resistant ATAF1 gene via Agrobacterium mediated gene transformation method. The ATAF1 gene was first cloned in the pK7WFG2 vector having kanamycin selectable marker and then it was introduced in the Agrobacterium tumefaciens strain GV3101 through heat shock method. The tomato cotyledon and hypocotyl ex-plants of variety "Rio Ggrande" were cultured on callus induction medium (MS + 2.5 mg/L IAA + 2 mg/L BAP). The calli were then infected with Agrobacterium tumefaciens strain GV3101 containing ATAF1 gene and selection was carried out on the kanamycin selectable medium (MS + 100 mg/L Kan), and were regenerated on MS medium with 1 mg/L IAA + 1 mg/L BAP. Out of 216 putative transformed calli, 13 calli were able to regenerate on the selection medium. Of the 13 calli, three transgenic tomato plantlets were recovered, and these were confirmed through PCR analysis for the presence of 432 bp fragment of ATAF1 gene. The transformation protocol reported here can be used to generate drought resistant tomato plants in future.
On 26 December 2004, a moment magnitude M-w = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region(1). The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave- height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN)(2,3), located approximately 9,000 km from the event. Here we present a modification of a standard array- seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.
Geothermobarometric, radiogenic isotopic and thermochronologic data are used to track the influence of an ancient continental margin (Western Province) on development of an adjacent Carboniferous-Cretaceous magmatic arc (Outboard Median Batholith) in Fiordland, New Zealand. The data show a record of complicated Mesozoic Gondwana margin growth. Paragneiss within the Outboard Median Batholith is of Carboniferous to Jurassic age and records burial to middle crustal depths in Late Jurassic-Early Cretaceous during subduction-related plutonism and arc thickening. In contrast, Western Province metasedimentary rocks in the area of study immediately west of the Outboard Median Batholith are Late Cambrian-Early Ordovician in age, recrystallized at the amphibolite facies in the Late Devonian-Early Carboniferous and exhibit no evidence for Mesozoic textural or isotopic reequilibration. A phase of deformation, between 128 and 116 Ma deformed, exhumed, and cooled the Outboard Median Batholith to greenschist facies temperatures, while large parts of the Western Province underwent >= 9 kbar metamorphic conditions. Zircon grains from Mesozoic inboard plutons are isotopically more evolved (epsilon Hf(t) = +2.3 to +4.0) than those in the Outboard Median Batholith (epsilon Hf(t) = +9.4 to +11.1). The contrasting zircon Hf isotope ratios, absence of S-type plutons or Proterozoic-Early Paleozoic inherited zircon, and the apparent absence of Early Paleozoic metasedimentary rocks indicates that the Outboard Median Batholith is unlikely to be underlain by the Western Province continental lithosphere. The new data are consistent with the Outboard Median Batholith representing an allochthonous (although not necessarily exotic) arc that was juxtaposed onto the Gondwana continental margin along the intervening Grebe Mylonite Zone.