@article{WulfBookhagenScherler2016, author = {Wulf, Hendrik and Bookhagen, Bodo and Scherler, Dirk}, title = {Differentiating between rain, snow, and glacier contributions to river discharge in the western Himalaya using remote-sensing data and distributed hydrological modeling}, series = {Advances in water resources}, volume = {88}, journal = {Advances in water resources}, publisher = {Elsevier}, address = {Oxford}, issn = {0309-1708}, doi = {10.1016/j.advwatres.2015.12.004}, pages = {152 -- 169}, year = {2016}, abstract = {Rivers draining the southern Himalaya provide most of the water supply for the densely populated Indo-Gangetic plains. Despite the importance of water resources in light of climate change, the relative contributions of rainfall, snow and glacier melt to discharge are not well understood, due to the scarcity of ground-based data in this complex terrain. Here, we quantify discharge sources in the Sutlej Valley, western Himalaya, from 2000 to 2012 with a distributed hydrological model that is based on daily, ground-calibrated remote-sensing observation. Based on the consistently good model performance, we analyzed the spatiotemporal distribution of hydrologic components and quantified their contribution to river discharge. Our results indicate that the Sutlej River's annual discharge at the mountain front is sourced to 55\% by effective rainfall (rainfall reduced by evapotranspiration), 35\% by snow melt and 10\% by glacier melt. In the high-elevation orogenic interior glacial runoff contributes \&\#8764;30\% to annual river discharge. These glacier melt contributions are especially important during years with substantially reduced rainfall and snowmelt runoff, as during 2004, to compensate for low river discharge and ensure sustained water supply and hydropower generation. In 2004, discharge of the Sutlej River totaled only half the maximum annual discharge; with 17.3\% being sourced by glacier melt. Our findings underscore the importance of calibrating remote-sensing data with ground-based data to constrain hydrological models with reasonable accuracy. For instance, we found that TRMM (Tropical Rainfall Measuring Mission) product 3B42 V7 systematically overestimates rainfall in arid regions of our study area by a factor of up to 5. By quantifying the spatiotemporal distribution of water resources we provide an important assessment of the potential impact of global warming on river discharge in the western Himalaya. Given the near-global coverage of the utilized remote-sensing datasets this hydrological modeling approach can be readily transferred to other data-sparse regions.}, language = {en} } @article{StolbovaSurovyatkinaBookhagenetal.2016, author = {Stolbova, Veronika and Surovyatkina, Elena and Bookhagen, Bodo and Kurths, J{\"u}rgen}, title = {Tipping elements of the Indian monsoon: Prediction of onset and withdrawal}, series = {Geophysical research letters}, volume = {43}, journal = {Geophysical research letters}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL068392}, pages = {3982 -- 3990}, year = {2016}, abstract = {Forecasting the onset and withdrawal of the Indian summer monsoon is crucial for the life and prosperity of more than one billion inhabitants of the Indian subcontinent. However, accurate prediction of monsoon timing remains a challenge, despite numerous efforts. Here we present a method for prediction of monsoon timing based on a critical transition precursor. We identify geographic regions-tipping elements of the monsoon-and use them as observation locations for predicting onset and withdrawal dates. Unlike most predictability methods, our approach does not rely on precipitation analysis but on air temperature and relative humidity, which are well represented both in models and observations. The proposed method allows to predict onset 2 weeks earlier and withdrawal dates 1.5 months earlier than existing methods. In addition, it enables to correctly forecast monsoon duration for some anomalous years, often associated with El Nino-Southern Oscillation.}, language = {en} } @article{SmithBookhagen2016, author = {Smith, Taylor and Bookhagen, Bodo}, title = {Assessing uncertainty and sensor biases in passive microwave data across High Mountain Asia}, series = {Remote sensing of environment : an interdisciplinary journal}, volume = {181}, journal = {Remote sensing of environment : an interdisciplinary journal}, publisher = {Elsevier}, address = {New York}, issn = {0034-4257}, doi = {10.1016/j.rse.2016.03.037}, pages = {174 -- 185}, year = {2016}, abstract = {Snowfall comprises a significant percentage of the annual water budget in High Mountain Asia (HMA), but snow water equivalent (SWE) is poorly constrained due to lack of in-situ measurements and complex terrain that limits the efficacy of modeling and observations. Over the past few decades, SWE has been estimated with passive microwave (PM) sensors with generally good results in wide, flat, terrain, and lower reliability in densely forested, complex, or high-elevation areas. In this study, we use raw swath data from five satellite - sensors the Special Sensor Microwave/Imager (SSMI) and Special Sensor Microwave Imager/Sounder (SSMIS) (1987-2015, F08, F11, F13, F17), Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E, 2002-2011), AMSR2 (2012-2015), and the Global Precipitation Measurement (GPM, 2014-2015) - in order to understand the spatial and temporal structure of native sensor, topographic, and land cover biases in SWE estimates in HMA. We develop a thorough understanding of the uncertainties in our SWE estimates by examining the impacts of topographic parameters (aspect, relief, hillslope angle, and elevation), land cover, native sensor biases, and climate parameters (precipitation, temperature, and wind speed). HMA, with its high seasonality, large topographic gradients and low relief at high elevations provides an excellent context to examine a wide range of climatic, land-cover, and topographic settings to better constrain SWE uncertainties and potential sensor bias. Using a multi-parameter regression, we compare long-term SWE variability to forest fraction, maximal multiyear snow depth, topographic parameters, and long-term average wind speed across both individual sensor time series and a merged multi-sensor dataset. In regions where forest cover is extensive, it is the strongest control on SWE variability. In those regions where forest density is low (<5\%), maximal snow depth dominates the uncertainty signal. In our regression across HMA, we find that forest fraction is the strongest control on SWE variability (75.8\%), followed by maximal multi-year snow depth (7.82\%), 90th percentile 10-m wind speed of a 10-year December-January-February (DJF) time series (5.64\%), 25th percentile DJF 10-m wind speed (5.44\%), and hillslope angle (5.24\%). Elevation, relief, and terrain aspect show very low influence on SWE variability (<1\%). We find that the GPM sensor provides the most robust regression results, and can be reliably used to estimate SWE in our study region. While forest cover and elevation have been integrated into many SWE algorithms, wind speed and long-term maximal snow depth have not. Our results show that wind redistribution of snow can have impacts on SWE, especially over large, flat, areas. Using our regression results, we have developed an understanding of sensor specific SWE uncertainties and their spatial patterns. The uncertainty maps developed in this study provide a first-order approximation of SWE-estimate reliability for much of HMA, and imply that high-fidelity SWE estimates can be produced for many high-elevation areas. (C) 2016 Elsevier Inc. All rights reserved.}, language = {en} } @article{SchildgenRobinsonSavietal.2016, author = {Schildgen, Taylor F. and Robinson, Ruth A. J. and Savi, Sara and Phillips, William M. and Spencer, Joel Q. G. and Bookhagen, Bodo and Scherler, Dirk and Tofelde, Stefanie and Alonso, Ricardo N. and Kubik, Peter W. and Binnie, Steven A. and Strecker, Manfred}, title = {Landscape response to late Pleistocene climate change in NW Argentina: Sediment flux modulated by basin geometry and connectivity}, series = {Journal of geophysical research : Earth surface}, volume = {121}, journal = {Journal of geophysical research : Earth surface}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1002/2015JF003607}, pages = {392 -- 414}, year = {2016}, abstract = {Fluvial fill terraces preserve sedimentary archives of landscape responses to climate change, typically over millennial timescales. In the Humahuaca Basin of NW Argentina (Eastern Cordillera, southern Central Andes), our 29 new optically stimulated luminescence ages of late Pleistocene fill terrace sediments demonstrate that the timing of past river aggradation occurred over different intervals on the western and eastern sides of the valley, despite their similar bedrock lithology, mean slopes, and precipitation. In the west, aggradation coincided with periods of increasing precipitation, while in the east, aggradation coincided with decreasing precipitation or more variable conditions. Erosion rates and grain size dependencies in our cosmogenic Be-10 analyses of modern and fill terrace sediments reveal an increased importance of landsliding compared to today on the west side during aggradation, but of similar importance during aggradation on the east side. Differences in the timing of aggradation and the Be-10 data likely result from differences in valley geometry, which causes sediment to be temporarily stored in perched basins on the east side. It appears as if periods of increasing precipitation triggered landslides throughout the region, which induced aggradation in the west, but blockage of the narrow bedrock gorges downstream from the perched basins in the east. As such, basin geometry and fluvial connectivity appear to strongly influence the timing of sediment movement through the system. For larger basins that integrate subbasins with differing geometries or degrees of connectivity (like Humahuaca), sedimentary responses to climate forcing are likely attenuated.}, language = {en} } @article{OlenBookhagenStrecker2016, author = {Olen, Stephanie M. and Bookhagen, Bodo and Strecker, Manfred}, title = {Role of climate and vegetation density in modulating denudation rates in the Himalaya}, series = {Earth \& planetary science letters}, volume = {445}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.03.047}, pages = {57 -- 67}, year = {2016}, abstract = {Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new Be-10 denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially-dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation. (C) 2016 Elsevier B.V. All rights reserved.}, language = {en} } @article{MalikBookhagenMucha2016, author = {Malik, Nishant and Bookhagen, Bodo and Mucha, Peter J.}, title = {Spatiotemporal patterns and trends of Indian monsoonal rainfall extremes}, series = {Geophysical research letters}, volume = {43}, journal = {Geophysical research letters}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL067841}, pages = {1710 -- 1717}, year = {2016}, abstract = {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.}, language = {en} } @article{HoffmannFeakinsBookhagenetal.2016, author = {Hoffmann, Bernd and Feakins, Sarah J. and Bookhagen, Bodo and Olen, Stephanie M. and Adhikari, Danda P. and Mainali, Janardan and Sachse, Dirk}, title = {Climatic and geomorphic drivers of plant organic matter transport in the Arun River, E Nepal}, series = {Earth \& planetary science letters}, volume = {452}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.07.008}, pages = {104 -- 114}, year = {2016}, language = {en} } @article{HartmanBookhagenChadwick2016, author = {Hartman, Brett D. and Bookhagen, Bodo and Chadwick, Oliver A.}, title = {The effects of check dams and other erosion control structures on the restoration of Andean bofedal ecosystems}, series = {Restoration Ecology}, volume = {24}, journal = {Restoration Ecology}, publisher = {Wiley-Blackwell}, address = {Hoboken}, issn = {1061-2971}, doi = {10.1111/rec.12402}, pages = {761 -- 772}, year = {2016}, abstract = {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.}, language = {en} } @article{ForteWhippleBookhagenetal.2016, author = {Forte, Adam M. and Whipple, Kelin X. and Bookhagen, Bodo and Rossi, Matthew W.}, title = {Decoupling of modern shortening rates, climate, and topography in the Caucasus}, series = {Earth \& planetary science letters}, volume = {449}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.06.013}, pages = {282 -- 294}, year = {2016}, abstract = {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.}, language = {en} } @article{DeyThiedeSchildgenetal.2016, author = {Dey, Saptarshi and Thiede, Rasmus Christoph and Schildgen, Taylor F. and Wittmann, Hella and Bookhagen, Bodo and Scherler, Dirk and Strecker, Manfred}, title = {Holocene internal shortening within the northwest Sub-Himalaya: Out-of-sequence faulting of the Jwalamukhi Thrust, India}, series = {Tectonics}, volume = {35}, journal = {Tectonics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1002/2015TC004002}, pages = {2677 -- 2697}, year = {2016}, abstract = {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.}, language = {en} } @article{DeyThiedeSchildgenetal.2016, author = {Dey, Saptarshi and Thiede, Rasmus Christoph and Schildgen, Taylor F. and Wittmann, Hella and Bookhagen, Bodo and Scherler, Dirk and Jain, Vikrant and Strecker, Manfred}, title = {Climate-driven sediment aggradation and incision since the late Pleistocene in the NW Himalaya, India}, series = {Earth \& planetary science letters}, volume = {449}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.05.050}, pages = {321 -- 331}, year = {2016}, abstract = {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.}, language = {en} } @article{CastinoBookhagenStrecker2016, author = {Castino, Fabiana and Bookhagen, Bodo and Strecker, Manfred}, title = {River-discharge dynamics in the Southern Central Andes and the 1976-77 global climate shift}, series = {Geophysical research letters}, volume = {43}, journal = {Geophysical research letters}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL070868}, pages = {11679 -- 11687}, year = {2016}, abstract = {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.}, language = {en} } @article{BrellRogassSegletal.2016, author = {Brell, Maximilian and Rogass, Christian and Segl, Karl and Bookhagen, Bodo and Guanter, Luis}, title = {Improving Sensor Fusion: A Parametric Method for the Geometric Coalignment of Airborne Hyperspectral and Lidar Data}, series = {IEEE transactions on geoscience and remote sensing}, volume = {54}, journal = {IEEE transactions on geoscience and remote sensing}, publisher = {Inst. of Electr. and Electronics Engineers}, address = {Piscataway}, issn = {0196-2892}, doi = {10.1109/TGRS.2016.2518930}, pages = {3460 -- 3474}, year = {2016}, abstract = {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.}, language = {en} } @article{BoersBookhagenMarwanetal.2016, author = {Boers, Niklas and Bookhagen, Bodo and Marwan, Norbert and Kurths, J{\"u}rgen}, title = {Spatiotemporal characteristics and synchronization of extreme rainfall in South America with focus on the Andes Mountain range}, series = {Climate dynamics : observational, theoretical and computational research on the climate system}, volume = {46}, journal = {Climate dynamics : observational, theoretical and computational research on the climate system}, publisher = {Springer}, address = {New York}, issn = {0930-7575}, doi = {10.1007/s00382-015-2601-6}, pages = {601 -- 617}, year = {2016}, abstract = {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.}, language = {en} }