Refine
Year of publication
Document Type
- Article (96)
- Postprint (14)
- Other (2)
- Part of a Book (1)
- Conference Proceeding (1)
- Doctoral Thesis (1)
Language
- English (115)
Keywords
- Himalaya (7)
- erosion (6)
- TRMM (4)
- cosmogenic nuclides (4)
- extreme rainfall (4)
- lidar (4)
- remote sensing (4)
- InSAR (3)
- Precipitation (3)
- TanDEM-X (3)
Institute
- Institut für Geowissenschaften (115) (remove)
The mechanisms by which climate and vegetation affect erosion rates over various time scales lie at the heart of understanding landscape response to climate change. Plot-scale field experiments show that increased vegetation cover slows erosion, implying that faster erosion should occur under low to moderate vegetation cover. However, demonstrating this concept over long time scales and across landscapes has proven to be difficult, especially in settings complicated by tectonic forcing and variable slopes. We investigate this problem by measuring cosmogenic Be-10-derived catchment-mean denudation rates across a range of climate zones and hillslope gradients in the Kenya Rift, and by comparing our results with those published from the Rwenzori Mountains of Uganda. We find that denudation rates from sparsely vegetated parts of the Kenya Rift are up to 0.13 mm/yr, while those from humid and more densely vegetated parts of the Kenya Rift flanks and the Rwenzori Mountains reach a maximum of 0.08 mm/yr, despite higher median hillslope gradients. While differences in lithology and recent land-use changes likely affect the denudation rates and vegetation cover values in some of our studied catchments, hillslope gradient and vegetation cover appear to explain most of the variation in denudation rates across the study area. Our results support the idea that changing vegetation cover can contribute to complex erosional responses to climate or land-use change and that vegetation cover can play an important role in determining the steady-state slopes of mountain belts through its stabilizing effects on the land surface.
Mapping urban forest leaf area index with airborne lidar using penetration metrics and allometry
(2015)
In urban areas, leaf area index (LAI) is a key ecosystem structural attribute with implications for energy and water balance, gas exchange, and anthropogenic energy use. In this study, we estimated LAI spatially using airborne lidar in downtown Santa Barbara, California, USA. We implemented two different modeling approaches. First, we directly estimated effective LAI (LAIe) using scan angle- and clump-corrected lidar laser penetration metrics (LPM). Second, we adapted existing allometric equations to estimate crown structural metrics including tree height and crown base height using lidar. The latter approach allowed for LAI estimates at the individual tree-crown scale. The LPM method, at both high and decimated point densities, resulted in good linear agreement with estimates from ground-based hemispherical photography (r(2) = 0.82, y = 0.99x) using a model that assumed a spherical leaf angle distribution. Within individual tree crown segments, the lidar estimates of crown structure closely paralleled field measurements (e.g., r(2) = 0.87 for crown length). LAI estimates based on the lidar crown measurements corresponded well with estimates from field measurements (r(2) = 0.84, y = 0.97x + 0.10). Consistency of the LPM and allometric lidar methods was also strong at 71 validation plots (r(2) = 0.88) and at 450 additional sample locations across the entire study area (r(2) = 0.72). This level of correspondence exceeded that of the canopy hemispherical photography and allometric, ground-based estimates (r(2) = 0.53). The first-order alignment of these two disparate methods may indicate that the error bounds for mapping LAI in cities are small enough to pursue large scale, spatially explicit estimation. (C) 2015 Elsevier Inc All rights reserved.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
We modeled the two most extreme highstands of Lake Naivasha during the last 175 k.y. to estimate potential precipitation/ evaporation changes in this basin. In a first step, the bathymetry of the paleolakes at f135 and 9 k.y. BP was reconstructed from sediment cores and surface outcrops. Second, we modeled the paleohydrologic budget during the highstands using a simplified coupled energy mass-balance model. Our results show that the hydrologic and hence the climate conditions at f135 and 9 k.y. BP were similar, but significantly different from today. The main difference is a f15% higher value in precipitation compared to the present. An adaptation and migration of vegetation in the cause of climate changes would result in a f30% increase in precipitation. The most likely cause for such a wetter climate at f135 and 9 k.y. BP is a more intense intertropical convergence and increased precipitation in East Africa.
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
Extreme Rainfall of the South American Monsoon System: A Dataset Comparison Using Complex Networks
(2015)
In this study, the authors compare six different rainfall datasets for South America with a focus on their representation of extreme rainfall during the monsoon season (December February): the gauge-calibrated TRMM 3B42 V7 satellite product; the near-real-time TRMM 3B42 V7 RT, the GPCP 1 degrees daily (1DD) V1.2 satellite gauge combination product, the Interim ECMWF Re-Analysis (ERA-Interim) product; output of a high-spatial-resolution run of the ECHAM6 global circulation model; and output of the regional climate model Eta. For the latter three, this study can be understood as a model evaluation. In addition to statistical values of local rainfall distributions, the authors focus on the spatial characteristics of extreme rainfall covariability. Since traditional approaches based on principal component analysis are not applicable in the context of extreme events, they apply and further develop methods based on complex network theory. This way, the authors uncover substantial differences in extreme rainfall patterns between the different datasets: (i) The three model-derived datasets yield very different results than the satellite gauge combinations regarding the main climatological propagation pathways of extreme events as well as the main convergence zones of the monsoon system. (ii) Large discrepancies are found for the development of mesoscale convective systems in southeastern South America. (iii) Both TRMM datasets and ECHAM6 indicate a linkage of extreme rainfall events between the central Amazon basin and the eastern slopes of the central Andes, but this pattern is not reproduced by the remaining datasets. The authors' study suggests that none of the three model-derived datasets adequately captures extreme rainfall patterns in South America.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.
Birds migrating across the Himalayan region fly over the highest peaks in the world, facing immense physiological and climatic challenges. The authors show the different strategies used by birds to cope with these challenges. Many wetland avian species are seen in the high-altitude lakes of the Himalayas and the adjoining Tibetan Plateau, such as Bar-Headed Geese. Ringing programmes have generated information about origins and destinations, and this book is the first to present information on the bird's exact migratory paths. Capitalising on knowledge generated through satellite telemetry, the authors describe the migratory routes of a multitude of birds flying over or skirting the Himalayas. The myriad of threats to migratory birds and the wetland system in the Central Asian Flyway are discussed, with ways to mitigate them. This volume will inform and persuade policy-makers and conservation practitioners to take appropriate measures for the long-term survival of this unique migration
The India-Eurasia continental collision zone provides a spectacular example of active mountain building and climatic forcing. In order to quantify the critically important process of mass removal, I analyzed spatial and temporal precipitation patterns of the oscillating monsoon system and their geomorphic imprints. I processed passive microwave satellite data to derive high-resolution rainfall estimates for the last decade and identified an abnormal monsoon year in 2002. During this year, precipitation migrated far into the Sutlej Valley in the northwestern part of the Himalaya and reached regions behind orographic barriers that are normally arid. There, sediment flux, mean basin denudation rates, and channel-forming processes such as erosion by debris-flows increased significantly. Similarly, during the late Pleistocene and early Holocene, solar forcing increased the strength of the Indian summer monsoon for several millennia and presumably lead to analogous precipitation distribution as were observed during 2002. However, the persistent humid conditions in the steep, high-elevation parts of the Sutlej River resulted in deep-seated landsliding. Landslides were exceptionally large, mainly due to two processes that I infer for this time: At the onset of the intensified monsoon at 9.7 ka BP heavy rainfall and high river discharge removed material stored along the river, and lowered the baselevel. Second, enhanced discharge, sediment flux, and increased pore-water pressures along the hillslopes eventually lead to exceptionally large landslides that have not been observed in other periods. The excess sediments that were removed from the upstream parts of the Sutlej Valley were rapidly deposited in the low-gradient sectors of the lower Sutlej River. Timing of downcutting correlates with centennial-long weaker monsoon periods that were characterized by lower rainfall. I explain this relationship by taking sediment flux and rainfall dynamics into account: High sediment flux derived from the upstream parts of the Sutlej River during strong monsoon phases prevents fluvial incision due to oversaturation the fluvial sediment-transport capacity. In contrast, weaker monsoons result in a lower sediment flux that allows incision in the low-elevation parts of the Sutlej River.
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
Aluminum-26 and beryllium-10 surface exposure dating on cut-and-fill river-terrace surfaces from the lower Sutlej Valley (northwest Himalaya) documents the close link between Indian Summer Monsoon (ISM) oscillations and intervals of enhanced fluvial incision. During the early Holocene ISM optimum, precipitation was enhanced and reached far into the internal parts of the orogen. The amplified sediment flux from these usually dry but glaciated areas caused alluviation of downstream valleys up to 120 m above present grade at ca. 9.9 k.y. B.P. Terrace formation (i.e., incision) in the coarse deposits occurred during century-long weak ISM phases that resulted in reduced moisture availability and most likely in lower sediment flux. Here, we suggest that the lower sediment flux during weak ISM phases allowed rivers to incise episodically into the alluvial fill
Hydrological modelling of a Pleistocene landslide-dammed lake in the Santa Maria Basin, NW Argentina
(2001)
The tectonic and climatic boundary conditions of the broken foreland and the orogen interior of the southern Central Andes of northwestern Argentina cause strong contrasts in elevation, rainfall, and surface-process regimes. The climatic gradient in this region ranges from the wet, windward eastern flanks (similar to 2 m/yr rainfall) to progressively drier western basins and ranges (similar to 0.1 m/yr) bordering the arid Altiplano-Puna Plateau. In this study, we analyze the impact of spatiotemporal climatic gradients on surface erosion: First, we present 41 new catchment-mean erosion rates derived from cosmogenic nuclide inventories to document spatial erosion patterns. Second, we re-evaluate paleoclimatic records from the Calchaquies basin (66 W, 26 S), a large intermontane basin bordered by high (> 4.5 km) mountain ranges, to demonstrate temporal variations in erosion rates associated with changing climatic boundary conditions during the late Pleistocene and Holocene. Three key observations in this region emphasize the importance of climatic parameters on the efficiency of surface processes in space and time: (1) First-order spatial patterns of erosion rates can be explained by a simple specific stream power (SSP) approach. We explicitly account for discharge by routing high-resolution, satellite derived rainfall. This is important as the steep climatic gradient results in a highly non-linear relation between drainage area and discharge. This relation indicates that erosion rates (ER) scale with ER similar to SSP1.4 on cosmogenic-nuclide time scales. (2) We identify an intrinsic channel-slope behavior in different climatic compartments. Channel slopes in dry areas (< 0.25 m/yr rainfall) are slightly steeper than in wet areas (> 0.75 m/yr) with equal drainage areas, thus compensating lower amounts of discharge with steeper slopes. (3) Erosion rates can vary by an order of magnitude between presently dry (similar to 0.05 mm/yr) and well-defined late Pleistocene humid (similar to 0.5 mm/yr) conditions within an intemontane basin. Overall, we document a strong climatic impact on erosion rates and channel slopes. We suggest that rainfall reaching areas with steeper channel slopes in the orogen interior during wetter climate periods results in intensified sediment mass transport, which is primarily responsible for maintaining the balance between surface uplift, erosion, sediment routing and transient storage in the orogen.
The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved
Late Quaternary intensified monsoon phases control landscape evolution in the northwest Himalaya
(2005)
The intensity of the Asian summer-monsoon circulation varies over decadal to millennial time scales and is reflected in changes in surface processes, terrestrial environments, and marine sediment records. However, the mechanisms of long-lived (2-5 k.y.) intensified monsoon phases, the related changes in precipitation distribution, and their effect on landscape evolution and sedimentation rates are not yet well understood. The and high-elevation sectors of the orogen correspond to a climatically sensitive zone that currently receives rain only during abnormal (i.e., strengthened) monsoon seasons. Analogous to present-day rainfall anomalies, enhanced precipitation during an intensified monsoon phase is expected to have penetrated far into these geomorphic threshold regions where hillslopes are close to the angle of failure. We associate landslide triggering during intensified monsoon phases with enhanced precipitation, discharge, and sediment flux leading to an increase in pore-water pressure, lateral scouring of rivers, and over- steepening of hillslopes, eventually resulting in failure of slopes and exceptionally large mass movements. Here we use lacustrine deposits related to spatially and temporally clustered large landslides (>0.5 km(3)) in the Sutlej Valley region of the northwest Himalaya to calculate sedimentation rates and to infer rainfall patterns during late Pleistocene (29-24 ka) and Holocene (10-4 ka) intensified monsoon phases. Compared to present-day sediment-flux measurements, a fivefold increase in sediment-transport rates recorded by sediments in landslide-dammed lakes characterized these episodes of high climatic variability. These changes thus emphasize the pronounced imprint of millennial-scale climate change on surface processes and landscape evolution
Synergistic applications based on integrated hyperspectral and lidar data are receiving a growing interest from the remote-sensing community. A prerequisite for the optimum sensor fusion of hyperspectral and lidar data is an accurate geometric coalignment. The simple unadjusted integration of lidar elevation and hyperspectral reflectance causes a substantial loss of information and does not exploit the full potential of both sensors. This paper presents a novel approach for the geometric coalignment of hyperspectral and lidar airborne data, based on their respective adopted return intensity information. The complete approach incorporates ray tracing and subpixel procedures in order to overcome grid inherent discretization. It aims at the correction of extrinsic and intrinsic (camera resectioning) parameters of the hyperspectral sensor. In additional to a tie-point-based coregistration, we introduce a ray-tracing-based back projection of the lidar intensities for area-based cost aggregation. The approach consists of three processing steps. First is a coarse automatic tie-point-based boresight alignment. The second step coregisters the hyperspectral data to the lidar intensities. Third is a parametric coalignment refinement with an area-based cost aggregation. This hybrid approach of combining tie-point features and area-based cost aggregation methods for the parametric coregistration of hyperspectral intensity values to their corresponding lidar intensities results in a root-mean-square error of 1/3 pixel. It indicates that a highly integrated and stringent combination of different coalignment methods leads to an improvement of the multisensor coregistration.
The fusion of hyperspectral imaging (HSI) sensor and airborne lidar scanner (ALS) data provides promising potential for applications in environmental sciences. Standard fusion approaches use reflectance information from the HSI and distance measurements from the ALS to increase data dimen-sionality and geometric accuracy. However, the potential for data fusion based on the respective intensity information of the complementary active and passive sensor systems is high and not yet fully exploited. Here, an approach for the rigorous illumination correction of HSI data, based on the radiometric cross-calibrated return intensity information of ALS data, is presented. The cross calibration utilizes a ray tracing-based fusion of both sensor measurements by intersecting their particular beam shapes. The developed method is capable of compensating for the drawbacks of passive HSI systems, such as cast and cloud shadowing effects, illumination changes over time, across track illumination, and partly anisotropy effects. During processing, spatial and temporal differences in illumination patterns are detected and corrected over the entire HSI wavelength domain. The improvement in the classification accuracy of urban and vegetation surfaces demonstrates the benefit and potential of the proposed HSI illumination correction. The presented approach is the first step toward the rigorous in-flight fusion of passive and active system characteristics, enabling new capabilities for a variety of applications.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
Understanding the evolution of continental deformation zones relies on quantifying spatial and temporal changes in deformation rates of tectonic structures. Along the eastern boundary of the Pamir-Tian Shan collision zone, we constrain secular variations of rock uplift rates for a series of five Quaternary detachment- and fault-related folds from their initiation to the modern day. When combined with GPS data, decomposition of interferometric synthetic aperture radar time series constrains the spatial pattern of surface and rock uplift on the folds deforming at decadal rates of 1-5mm/yr. These data confirm the previously proposed basinward propagation of structures during the Quaternary. By fitting our geodetic rates and previously published geologic uplift rates with piecewise linear functions, we find that gradual rate changes over >100kyr can explain the interferometric synthetic aperture radar observations where changes in average uplift rates are greater than similar to 1 mm/yr among different time intervals (similar to 10(1), 10(4-5), and 10(5-6) years).
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.
Effects of topographic smoothing on the simulation of winter precipitation in High Mountain Asia
(2017)
Numerous studies have projected future changes in High Mountain Asia water resources based on temperature and precipitation from global circulation models (GCMs) under future climate scenarios. Although the potential benefit of such studies is immense, coarse grid-scale GCMs are unable to resolve High Mountain Asia's complex topography and thus have a biased representation of regional weather and climate. This study investigates biases in the simulation of physical mechanisms that generate snowfall and contribute to snowpack in High Mountain Asia in coarse topography experiments using the Weather Research and Forecasting model. Regional snowpack is event driven, thus 33 extreme winter orographic precipitation events are simulated at fine atmospheric resolution with 6.67 km resolution topography and smoothed 1.85° × 1.25° GCM topography. As with many modified topography experiments performed in other regions, the distribution of precipitation is highly dependent on first-order orographic effects, which dominate regional meteorology. However, we demonstrate that topographic smoothing enhances circulation in simulated extratropical cyclones, with significant impacts on orographic precipitation. Despite precipitation reductions of 28% over the highest ranges, due to reduced ascent on windward slopes, total precipitation over the study domain increased by an average of 9% in smoothed topography experiments on account of intensified extratropical cyclone dynamics and cross-barrier moisture flux. These findings identify an important source of bias in coarse-resolution simulated precipitation in High Mountain Asia, with important implications for the application of GCMs toward projecting future hydroclimate in the region.
Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes
(2020)
During the South-American Monsoon season, deep convective systems occur at the eastern flank of the Central Andes leading to heavy rainfall and flooding. We investigate the large- and meso-scale atmospheric dynamics associated with extreme discharge events (> 99.9th percentile) observed in two major river catchments meridionally stretching from humid to semi-arid conditions in the southern Central Andes. Based on daily gauge time series and ERA-Interim reanalysis, we made the following three key observations: (1) for the period 1940-2016 daily discharge exhibits more pronounced variability in the southern, semi-arid than in the northern, humid catchments. This is due to a smaller ratio of discharge magnitudes between intermediate (0.2 year return period) and rare events (20 year return period) in the semi-arid compared to the humid areas; (2) The climatological composites of the 40 largest discharge events showed characteristic atmospheric features of cold surges based on 5-day time-lagged sequences of geopotential height at different levels in the troposphere; (3) A subjective classification revealed that 80% of the 40 largest discharge events are mainly associated with the north-northeastward migration of frontal systems and 2/3 of these are cold fronts, i.e. cold surges. This work highlights the importance of cold surges and their related atmospheric processes for the generation of heavy rainfall events and floods in the southern Central Andes.
Recent studies have shown that the 1976-77 global climate shift strongly affected the South American climate. In our study, we observed a link between this climate shift and river-discharge variability in the subtropical Southern Central Andes. We analyzed the daily river-discharge time series between 1940 and 1999 from small to medium mountain drainage basins (10(2)-10(4) km(2) ) across a steep climatic and topographic gradient. We document that the discharge frequency distribution changed significantly, with higher percentiles exhibiting more pronounced trends. A change point between 1971 and 1977 marked an intensification of the hydrological cycle, which resulted in increased river discharge. In the upper Rio Bermejo basin of the northernmost Argentine Andes, the mean annual discharge increased by 40% over 7 years. Our findings are important for flood risk management in areas impacted by the 1976-77 climate shift; discharge frequency distribution analysis provides important insights into the variability of the hydrological cycle in the Andean realm.
The eastern flanks of the Central Andes are characterized by deep convection, exposing them to hydrometeorological extreme events, often resulting in floods and a variety of mass movements. We assessed the spatiotemporal pattern of rainfall trends and the changes in the magnitude and frequency of extreme events (ae<yen>95th percentile) along an E-W traverse across the southern Central Andes using rain-gauge and high-resolution gridded datasets (CPC-uni and TRMM 3B42 V7). We generated different climate indices and made three key observations: (1) an increase of the annual rainfall has occurred at the transition between low (< 0.5 km) and intermediate (0.5-3 km) elevations between 1950 and 2014. Also, rainfall increases during the wet season and, to a lesser degree, decreases during the dry season. Increasing trends in annual total amounts characterize the period 1979-2014 in the arid, high-elevation southern Andean Plateau, whereas trend reversals with decreasing annual total amounts were found at low elevations. (2) For all analyzed periods, we observed small or no changes in the median values of the rainfall-frequency distribution, but significant trends with intensification or attenuation in the 95th percentile. (3) In the southern Andean Plateau, extreme rainfall events exhibit trends towards increasing magnitude and, to a lesser degree, frequency during the wet season, at least since 1979. Our analysis revealed that low (< 0.5 km), intermediate (0.5-3 km), and high-elevation (> 3 km) areas respond differently to changing climate conditions, and the transition zone between low and intermediate elevations is characterized by the most significant changes.
This study analyzes the discharge variability of small to medium drainage basins (10(2)-10(4) km(2)) in the southern Central Andes of NW Argentina. The Hilbert-Huang Transform (HHT) was applied to evaluate non-stationary oscillatory modes of variability and trends, based on four time series of monthly normalized discharge anomaly between 1940 and 2015. Statistically significant trends reveal increasing discharge during the past decades and document an intensification of the hydrological cycle during this period. An Ensemble Empirical Mode Decomposition (EEMD) analysis revealed that discharge variability in this region can be best described by five quasi-periodic statistically significant oscillatory modes, with mean periods varying from 1 to 20 y. Moreover, we show that discharge variability is most likely linked to the phases of the Pacific Decadal Oscillation (PDO) at multi-decadal timescales (similar to 20 y) and, to a lesser degree, to the Tropical South Atlantic SST anomaly (TSA) variability at shorter timescales (similar to 2-5 y). Previous studies highlighted a rapid increase in discharge in the southern Central Andes during the 1970s, inferred to have been associated with the global 1976-77 climate shift. Our results suggest that the rapid discharge increase in the NW Argentine Andes coincides with the periodic enhancement of discharge, which is mainly linked to a negative to positive transition of the PDO phase and TSA variability associated with a long-term increasing trend. We therefore suggest that variations in discharge in this region are largely driven by both natural variability and the effects of global climate change. We furthermore posit that the links between atmospheric and hydrologic processes result from a combination of forcings that operate on different spatiotemporal scales. (C) 2017 Elsevier B.V. All rights reserved.
The structure and organization of river networks has been used for decades to investigate the influence of climate and tectonics on landscapes. The majority of these studies either analyze rivers in profile view by extracting channel steepness or calculate planform metrics such as drainage density. However, these techniques rely on the assumption of homogeneity: that intrinsic and external factors are spatially or temporally invariant over the measured profile. This assumption is violated for the majority of Earth's landscapes, where variations in uplift rate, rock strength, climate, and geomorphic process are almost ubiquitous. We propose a method for classifying river profiles to identify landscape regions with similar characteristics by adapting hierarchical clustering algorithms developed for time series data. We first test our clustering on two landscape evolution scenarios and find that we can successfully cluster regions with different erodibility and detect the transient response to sudden base level fall. We then test our method in two real landscapes: first in Bitterroot National Forest, Idaho, where we demonstrate that our method can detect transient incision waves and the topographic signature of fluvial and debris flow process regimes; and second, on Santa Cruz Island, California, where our technique identifies spatial patterns in lithology not detectable through normalized channel steepness analysis. By calculating channel steepness separately for each cluster, our method allows the extraction of more reliable steepness metrics than if calculated for the landscape as a whole. These examples demonstrate the method's ability to disentangle fluvial morphology in complex lithological and tectonic settings.
Deciphering the response of sediment routing systems to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Kangra Basin (northwest Sub-Himalaya, India), upper Pleistocene to Holocene alluvial fills and fluvial terraces record periodic fluctuations of sediment supply and transport capacity on timescales of 10(3) to 10(5) yr. To evaluate the potential influence of climate change on these fluctuations, we compare the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with climate archives. New surface-exposure dating of six terrace levels with in-situ cosmogenic Be-10 indicates the onset of incision phases. Two terrace surfaces from the highest level (T1) sculpted into the oldest preserved alluvial fan (AF1) date back to 53.4 +/- 3.2 ka and 43.0 +/- 2.7 ka (1 sigma). T2 surfaces sculpted into the remnants of AF1 have exposure ages of 18.6 +/- 1.2 ka and 15.3 +/- 0.9 ka, while terraces sculpted into the upper Pleistocene-Holocene fan (AF2) provide ages of 9.3 +/- 0.4 ka (T3), 7.1 +/- 0.4 ka (T4), 5.2 +/- 0.4 ka (T5) and 3.6 +/- 0.2 ka (T6). Together with previously published OSL ages yielding the timing of aggradation, we find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon. During periods of increased monsoon intensity and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of weakened monsoon intensity or lower sediment supply coincide with incision. (C) 2016 Elsevier B.V. All rights reserved.
The southernmost thrust of the Himalayan orogenic wedge that separates the foreland from the orogen, the Main Frontal Thrust, is thought to accommodate most of the ongoing crustal shortening in the Sub-Himalaya. Steepened longitudinal river profile segments, terrace offsets, and back-tilted fluvial terraces within the Kangra reentrant of the NW Sub-Himalaya suggest Holocene activity of the Jwalamukhi Thrust (JMT) and other thrust faults that may be associated with strain partitioning along the toe of the Himalayan wedge. To assess the shortening accommodated by the JMT, we combine morphometric terrain analyses with in situ Be-10-based surface-exposure dating of the deformed terraces. Incision into upper Pleistocene sediments within the Kangra Basin created two late Pleistocene terrace levels (T1 and T2). Subsequent early Holocene aggradation shortly before similar to 10ka was followed by episodic reincision, which created four cut-and-fill terrace levels, the oldest of which (T3) was formed at 10.10.9ka. A vertical offset of 445m of terrace T3 across the JMT indicates a shortening rate of 5.60.8 to 7.51.1mma(-1) over the last similar to 10ka. This result suggests that thrusting along the JMT accommodates 40-60% of the total Sub-Himalayan shortening in the Kangra reentrant over the Holocene. We speculate that this out-of-sequence shortening may have been triggered or at least enhanced by late Pleistocene and Holocene erosion of sediments from the Kangra Basin.
The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
The Greater and Lesser Caucasus mountains and their associated foreland basins contain similar rock types, experience a similar two-fold, along-strike variation in mean annual precipitation, and were affected by extreme base-level drops of the neighboring Caspian Sea. However, the two Caucasus ranges are characterized by decidedly different tectonic regimes and rates of deformation that are subject to moderate (less than an order of magnitude) gradients in climate, and thus allow for a unique opportunity to isolate the effects of climate and tectonics in the evolution of topography within active orogens. There is an apparent disconnect between modern climate, shortening rates, and topography of both the Greater Caucasus and Lesser Caucasus which exhibit remarkably similar topography along-strike despite the gradients in forcing. By combining multiple datasets, we examine plausible causes for this disconnect by presenting a detailed analysis of the topography of both ranges utilizing established relationships between catchment-mean erosion rates and topography (local relief, hillslope gradients, and channel steepness) and combining it with a synthesis of previously published low-temperature thermochronologic data. Modern climate of the Caucasus region is assessed through an analysis of remotely-sensed data (TRMM and MODIS) and historical streamflow data. Because along-strike variation in either erosional efficiency or thickness of accreted material fail to explain our observations, we suggest that the topography of both the western Lesser and Greater Caucasus are partially supported by different geodynamic forces. In the western Lesser Caucasus, high relief portions of the landscape likely reflect uplift related to ongoing mantle lithosphere delamination beneath the neighboring East Anatolian Plateau. In the Greater Caucasus, maintenance of high topography in the western portion of the range despite extremely low (<2-4 mm/y) modern convergence rates may be related to dynamic topography from detachment of the north-directed Greater Caucasus slab or to a recent slowing of convergence rates. Large-scale spatial gradients in climate are not reflected in the topography of the Caucasus and do not seem to exert any significant control on the tectonics or structure of either range. (C) 2016 Elsevier B.V. All rights reserved.
The northwest Argentine Andes constitute a premier natural laboratory to assess the complex interactions between isolated uplifts, orographic precipitation gradients, and related erosion and sedimentation patterns. Here we present new stratigraphic observations and age information from intermontane basin sediments to elucidate the Neogene to Quaternary shortening history and associated sediment dynamics of the broken Salta foreland. This part of the Andean orogen, which comprises an array of basement-cored range uplifts, is located at similar to 25 degrees S and lies to the east of the arid intraorogenic Altiplano/Puna plateau. In the Salta foreland, spatially and temporally disparate range uplift along steeply dipping inherited faults has resulted in foreland compartmentalization with steep basin-tobasin precipitation gradients. Sediment architecture and facies associations record a three-phase (similar to 10, similar to 5, and <2 Ma), east directed, yet unsystematic evolution of shortening, foreland fragmentation, and ensuing changes in precipitation and sediment transport. The provenance signatures of these deposits reflect the trapping of sediments in the intermontane basins of the Andean hinterland, as well as the evolution of a severed fluvial network. Present-day moisture supply to the hinterland is determined by range relief and basin elevation. The conspiring effects of range uplift and low rainfall help the entrapment and long-term storage of sediments, ultimately raising basin elevation in the hinterland, which may amplify aridification in the orogen interior.
Restoring degraded lands in rural environments that are heavily managed to meet subsistence needs is a challenge due to high rates of disturbance and resource extraction. This study investigates the efficacy of erosion control structures (ECSs) as restoration tools in the context of a watershed rehabilitation and wet meadow (bofedal) restoration program in the Bolivian Andes. In an effort to enhance water security and increase grazing stability, Aymara indigenous communities built over 15,000 check dams, 9,100 terraces, 5,300 infiltration ditches, and 35 pasture improvement trials. Communities built ECSs at different rates, and we compared vegetation change in the highest restoration management intensity, lowest restoration management intensity, and nonproject control communities. We used line transects to measure changes in vegetation cover and standing water in gullies with check dams and without check dams, and related these ground measurements to a time series (1986-2009) of normalized difference vegetation index derived from Landsat TM5 images. Evidence suggests that check dams increase bofedal vegetation and standing water at a local scale, and lead to increased greenness at a basin scale when combined with other ECSs. Watershed rehabilitation enhances ecosystem services significant to local communities (grazing stability, water security), which creates important synergies when conducting land restoration in rural development settings.
Geodetic and seismologic studies support a tectonic model for the central Himalaya wherein similar to 2 cm/yr of Indo-Asian convergence is accommodated along the primary decollement under the range, the Main Himalayan thrust. A steeper midcrustal ramp in the Main Himalayan thrust is commonly invoked as driving rapid rock uplift along a range-parallel band in the Greater Himalaya. This tectonic model, developed primarily from studies in central Nepal, is commonly assumed to project along strike with little lateral variation in Main Himalayan thrust geometry or associated rock uplift patterns. Here, we synthesize multiple lines of evidence for a major discontinuity in the Main Himalayan thrust in western Nepal. Analysis of topography and seismicity indicates that west of similar to 82.5 degrees E, the single band of steep topography and seismicity along the Main Himalayan thrust ramp in central Nepal bifurcates around a high-elevation, low-relief landscape, resulting in a two-step topographic front along an similar to 150 km segment of the central Himalaya. Although multiple models could explain this bifurcation, the full suite of data appears to be most consistent with a northward bend to the Main Himalayan thrust ramp and activation of a young duplex horse to the south. This poorly documented segmentation of the Main Himalayan thrust has important implications for the seismogenic potential of the western Nepal seismic gap and for models of the ongoing evolution of the orogen.
Dating growth strata and basin fill by combining 26Al/10Be burial dating and magnetostratigraphy
(2018)
Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.
With the advent of the two Sentinel-1 (S1) satellites, Synthetic Aperture Radar (SAR) data with high temporal and spatial resolution are freely available. This provides a promising framework to facilitate detailed investigations of surface instabilities and movements on large scales with high temporal resolution, but also poses substantial processing challenges because of storage and computation requirements. Methods are needed to efficiently detect short term changes in dynamic environments. Approaches considering pair-wise processing of a series of consecutive scenes to retain maximum temporal resolution in conjunction with time series analyses are required. Here we present OSARIS, the “Open Source SAR Investigation System,” as a framework to process large stacks of S1 data on high-performance computing clusters. Based on Generic Mapping Tools SAR, shell scripts, and the workload manager Slurm, OSARIS provides an open and modular framework combining parallelization of high-performance C programs, flexible processing schemes, convenient configuration, and generation of geocoded stacks of analysis-ready base data, including amplitude, phase, coherence, and unwrapped interferograms. Time series analyses can be conducted by applying automated modules to the data stacks. The capabilities of OSARIS are demonstrated in a case study from the northwestern Tien Shan, Central Asia. After merging of slices, a total of 80 scene pairs were processed from 174 total input scenes. The coherence time series exhibits pronounced seasonal variability, with relatively high coherence values prevailing during the summer months in the nival zone. As an example of a time series analysis module, we present OSARIS' “Unstable Coherence Metric” which identifies pixels affected by significant drops from high to low coherence values. Measurements of motion provided by LOSD measurements require careful evaluation because interferometric phase unwrapping is prone to errors. Here, OSARIS provides a series of modules to detect and mask unwrapping errors, correct for atmospheric disturbances, and remove large-scale trends. Wall clock processing time for the case study (area ~9,000 km2) was ~12 h 4 min on a machine with 400 cores and 2 TB RAM. In total, ~12 d 10 h 44 min (~96%) were saved through parallelization. A comparison of selected OSARIS datasets to results from two state-of-the-art SAR processing suites, ISCE and SNAP, shows that OSARIS provides products of competitive quality despite its high level of automatization. OSARIS thus facilitates efficient S1-based region-wide investigations of surface movement events over multiple years.
Glacial deposits on the high-altitude, arid southern Central Andean Plateau (CAP), the Puna in northwestern Argentina, document past changes in climate, but the associated geomorphic features have rarely been directly dated. This study provides direct age control of glacial moraine deposits from the central Puna (24 degrees S) at elevations of 3900-5000 m through surface exposure dating with cosmogenic nuclides. Our results show that the most extensive glaciations occurred before 95 ka and an additional major advance occurred between 46 and 39 ka. The latter period is synchronous with the highest lake levels in the nearby Pozuelos basin and the Minchin (Inca Huasi) wet phase on the Altiplano in the northern CAP. None of the dated moraines produced boulder ages corresponding to the Tauca wet phase (24-15 ka). Additionally, the volcanic lithologies of the deposits allow us to establish production ratios at low latitude and high elevation for five different nuclide and mineral systems: Be-10, Ne-21, and Al-26 from quartz (11 or 12 samples) and He-3 and Ne-21 from pyroxene (10 samples). We present production ratios for all combinations of the measured nuclides and cross-calibrated production rates for 21Ne in pyroxene and quartz for the high, (sub-)tropical Andes. The production rates are based on our Be-10-normalized production ratios and a weighted mean of reference 10Be production rates calibrated in the high, tropical Andes (4.02 +/- 0.12 at g(-1) yr(-1)). These are, Ne-21(qtz): 18.1 +/- 1.2 at g(-1) yr(-1) and Ne-21(px): 36.6 +/- 1.8 at g(-1) yr(-1) (En(88-94)) scaled to sea level and high latitude using the Lal/Stone scheme, with 1 sigma uncertainties. As He-3 and Al-26 have been directly calibrated in the tropical Andes, we recommend using those rates. Finally, we compare exposure ages calculated using all measured cosmogenic nuclides from each sample, including 11 feldspar samples measured for Cl-36, and a suite of previously published production rates. (C) 2018 Published by Elsevier B.V.
In this study, we provide a comprehensive analysis of trends in the extremes during the Indian summer monsoon (ISM) months (June to September) at different temporal and spatial scales. Our goal is to identify and quantify spatiotemporal patterns and trends that have emerged during the recent decades and may be associated with changing climatic conditions. Our analysis primarily relies on quantile regression that avoids making any subjective choices on spatial, temporal, or intensity pattern of extreme rainfall events. Our analysis divides the Indian monsoon region into climatic compartments that show different and partly opposing trends. These include strong trends toward intensified droughts in Northwest India, parts of Peninsular India, and Myanmar; in contrast, parts of Pakistan, Northwest Himalaya, and Central India show increased extreme daily rain intensity leading to higher flood vulnerability. Our analysis helps explain previously contradicting results of trends in average ISM rainfall.
The effect of Indian Summer Monsoon rainfall on surface water delta D values in the central Himalaya
(2018)
Stable isotope proxy records, such as speleothems, plant-wax biomarker records, and ice cores, are suitable archives for the reconstruction of regional palaeohydrologic conditions. But the interpretation of these records in the tropics, especially in the Indian Summer Monsoon (ISM) domain, is difficult due to differing moisture and water sources: precipitation from the ISM and Winter Westerlies, as well as snow- and glacial meltwater. In this study, we use interannual differences in ISM strength (2011-2012) to understand the stable isotopic composition of surface water in the Arun River catchment in eastern Nepal. We sampled main stem and tributary water (n = 204) for stable hydrogen and oxygen isotope analysis in the postmonsoon phase of two subsequent years with significantly distinct ISM intensities. In addition to the 2011/2012 sampling campaigns, we collected a 12-month time series of main stem waters (2012/2013, n = 105) in order to better quantify seasonal effects on the variability of surface water delta O-18/delta D. Furthermore, remotely sensed satellite data of rainfall, snow cover, glacial coverage, and evapotranspiration was evaluated. The comparison of datasets from both years revealed that surface waters of the main stem Arun and its tributaries were D-enriched by similar to 15 parts per thousand when ISM rainfall decreased by 20%. This strong response emphasizes the importance of the ISM for surface water run-off in the central Himalaya. However, further spatio-temporal analysis of remote sensing data in combination with stream water d-excess revealed that most high-altitude tributaries and the Tibetan part of the Arun receive high portions of glacial melt water and likely Winter Westerly Disturbances precipitation. We make the following two implications: First, palaeohydrologic archives found in high-altitude tributaries and on the southern Tibetan Plateau record a mixture of past precipitation delta D values and variable amounts of additional water sources. Second, surface water isotope ratios of lower elevated tributaries strongly reflect the isotopic composition of ISM rainfall implying a suitable region for the analysis of potential delta D value proxy records.
This work explores the control of fore-arc structure on segmentation of megathrust earthquake ruptures using coastal geomorphic markers. The Arauco-Nahuelbuta region at the south-central Chile margin constitutes an anomalous fore- arc sector in terms of topography, geology, and exhumation, located within the overlap between the Concepcion and Valdivia megathrust segments. This boundary, however, is only based on similar to 500 years of historical records. We integrate deformed marine terraces dated by cosmogenic nuclides, syntectonic sediments, published fission track data, seismic reflection profiles, and microseismicity to analyze this earthquake boundary over 10(2) -10(6) years. Rapid exhumation of Nahuelbuta's dome-like core started at 4 +/- 1.2 Ma, coeval with inversion of the adjacent Arauco basin resulting in emergence of the Arauco peninsula. Here, similarities between topography, spatiotemporal trends in fission track ages, Pliocene-Pleistocene growth strata, and folded marine terraces suggest that margin-parallel shortening has dominated since Pliocene time. This shortening likely results from translation of a fore-arc sliver or microplate, decoupled from South America by an intra-arc strike-slip fault. Microplate collision against a buttress leads to localized uplift at Arauco accrued by deep-seated reverse faults, as well as incipient oroclinal bending. The extent of the Valdivia segment, which ruptured last in 1960 with an M-w 9.5 event, equals the inferred microplate. We propose that mechanical homogeneity of the fore-arc microplate delimits the Valdivia segment and that a marked discontinuity in the continental basement at Arauco acts as an inhomogeneous barrier controlling nucleation and propagation of 1960-type ruptures. As microplate-related deformation occurs since the Pliocene, we propose that this earthquake boundary and the extent of the Valdivia segment are spatially stable seismotectonic features at million year scale.
Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014–2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet–dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust’s composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR–SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite–gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan’s moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.
Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014–2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet–dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust’s composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR–SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite–gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan’s moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.
A valley-filling ignimbrite re-exposed through subsequent river incision at the southern margin of the Andean (Puna) plateau preserves pristine geological evidence of pre-late Miocene palaeotopography in the north western Argentine Andes. Our new Ar-40/(39) Ar dating of the Las Papas Ignimbrites yields a plateau age of 9.24 +/- 0.03 Ma, indicating valley-relief and orographic-barrier conditions comparable to the present-day. A later infill of Plio-Pleistocene coarse conglomerates has been linked to wetter conditions, but resulted in no additional net incision of the Las Papas valley, considering that the base of the ignimbrite remains unexposed in the valley bottom. Our observations indicate that at least 550 m of local plateau margin relief (and likely > 2 km) existed by 9 Ma at the southern Puna margin, which likely aided the efficiency of the orographic barrier to rainfall along the eastern and south eastern flanks of the Puna and causes aridity in the plateau interior.
Resolving Earth's surface at the meter scale is essential for an improved understanding of the dynamics of mass-movement processes. In this study, we explore the applicability and potential of digital elevation models (DEMs) derived from stereophotogrammetry to detect debris-flow channels in the Quebrada del Toro in the northwestern Argentine Andes. Our analysis relies on a high-resolution (3 m) DEM created from SPOT-7 tri-stereo satellite data. We carefully validated DEM quality with ∼6,000 differential GPS points and identified optimal parameters for DEM generation in high-relief terrain. After multiple processing steps, we achieved an accuracy of 0.051 ± 1.915 m (1σ) using n = 3,139 control points with cm precision. Previous studies have used the drainage area and slope framework to identify topographic signatures of debris flows within a catchment. We built upon this and investigated individual river-channel segments using connected-component (CC) analysis on meter-scale topographic data. We define CC as segments of similar slope along the channel profile. Based on seven manually identified debris-flow catchments, we developed a debris-flow similarity index using component length and mean channel-segment slope and identified channel segments that have likely been shaped by debris flows. The presented approach has the potential to resolve intra-catchment variability of transport processes, allows to constrain the extent of debris-flow channels more precisely than slope-area analysis, and highlights the versatility of combined space- and field-based observations for natural-hazard assessments.
This study demonstrates the potential of using single-pass TanDEM-X (TDX) radar imagery to analyse inter- and intra-annual glacier changes in mountainous terrain. Based on SAR images acquired in February 2012, March 2013 and November 2013 over the Inylchek Glacier, Kyrgyzstan, we discuss in detail the processing steps required to generate three reliable digital elevation models (DEMs) with a spatial resolution of 10 m that can be used for glacial mass balance studies. We describe the interferometric processing steps and the influence of a priori elevation information that is required to model long wavelength topographic effects. We also focus on DEM alignment to allow optimal DEM comparisons and on the effects of radar signal penetration on ice and snow surface elevations. We finally compare glacier elevation changes between the three TDX DEMs and the C-band shuttle radar topography mission (SRTM) DEM from February 2000. We introduce a new approach for glacier elevation change calculations that depends on the elevation and slope of the terrain. We highlight the superior quality of the TDX DEMs compared to the SRTM DEM, describe remaining DEM uncertainties and discuss the limitations that arise due to the side-looking nature of the radar sensor. (C) 2017 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.
Streams commonly respond to base-level fall by localizing erosion within steepened, convex knickzone reaches. Localized incision causes knickzone reaches to migrate upstream. Such migrating knickzones dictate the pace of landscape response to changes in tectonics or erosional efficiency and can help quantify the timing and source of base-level fall. Identification of knickzones typically requires individual selection of steepened reaches: a process that is tedious and subjective and has no efficient means to measure knickzone size. We construct an algorithm to automate this procedure by selecting the bounds of knickzone reaches in a -space (drainage-area normalized) framework. An automated feature calibrates algorithm parameters to a subset of knickzones handpicked by the user. The algorithm uses these parameters as consistent criteria to identify knickzones objectively, and then the algorithm measures the height, length, and slope of each knickzone reach. We test the algorithm on 1, 10, and 30m resolution digital elevation models (DEMs) of six catchments (trunk-stream lengths: 2.1-5.4km) on Santa Cruz Island, southern California. On the 1m DEM, algorithm-selected knickzones confirm 93% of handpicked knickzone positions (n=178) to a spatial accuracy of 100m, 88% to an accuracy within 50m, and 46% to an accuracy within 10m. Using 10 and 30m DEMs, accuracy is similar: 88-86% to 100m and 82% to 50m (n=38 and 36, respectively). The algorithm enables efficient regional comparison of the size and location of knickzones with geologic structures, mapped landforms, and hillslope morphology, thereby facilitating approaches to characterize the dynamics of transient landscapes. Plain Language Summary The shape of rivers reflects the environments that they flow through and the environments that they link together: mountains and oceans. Anywhere along the length of a river, changes in environmental conditions are propagated upstream and downstream as the river changes its morphology to match the new environmental conditions. Commonly, rivers steepen as land uplifts faster in regions of high tectonic convergence. The steepening of river gradients is propagated upstream and can be mapped to trace zones of high tectonic activity across landscapes and estimate the source and timing of environmental change. Such insights may indicate regions where earthquakes have become more frequent in the recent past and how rivers respond to these changes. In this submission, we detail an algorithm that can use digital topographic data (similar to google earth), to automatically map and measure anomalously steep river reaches across continental scales. This technology can highlight areas that have experienced recent sustained changes in environmental conditions, evident by changes in the morphology of rivers. Such environmental conditions could be changes in tectonic uplift and earthquake activity, changes in sea level, changes in land-use, or changes in climate, all factors that can produce measurable differences in river morphology over time.