Refine
Year of publication
- 2014 (226) (remove)
Document Type
- Article (196)
- Doctoral Thesis (17)
- Preprint (8)
- Review (4)
- Postprint (1)
Is part of the Bibliography
- yes (226) (remove)
Keywords
- Earthquake source observations (7)
- Holocene (6)
- Tibetan Plateau (4)
- Erosion (3)
- Himalaya (3)
- Seismicity and tectonics (3)
- Aleatory variability (2)
- Andes (2)
- Baseline shift (2)
- Body waves (2)
Institute
- Institut für Geowissenschaften (226) (remove)
Surface displacement at volcanic edifices is related to subsurface processes associated with magma movements, fluid transfers within the volcano edifice and gravity-driven deformation processes. Understanding of associated ground displacements is of importance for assessment of volcanic hazards. For example, volcanic unrest is often preceded by surface uplift, caused by magma intrusion and followed by subsidence, after the withdrawal of magma. Continuous monitoring of the surface displacement at volcanoes therefore might allow the forecasting of upcoming eruptions to some extent. In geophysics, the measured surface displacements allow the parameters of possible deformation sources to be estimated through analytical or numerical modeling. This is one way to improve the understanding of subsurface processes acting at volcanoes. Although the monitoring of volcanoes has significantly improved in the last decades (in terms of technical advancements and number of monitored volcanoes), the forecasting of volcanic eruptions remains puzzling. In this work I contribute towards the understanding of the subsurface processes at volcanoes and thus to the improvement of volcano eruption forecasting. I have investigated the displacement field of Llaima volcano in Chile and of Tendürek volcano in East Turkey by using synthetic aperture radar interferometry (InSAR). Through modeling of the deformation sources with the extracted displacement data, it was possible to gain insights into potential subsurface processes occurring at these two volcanoes that had been barely studied before. The two volcanoes, although of very different origin, composition and geometry, both show a complexity of interacting deformation sources. At Llaima volcano, the InSAR technique was difficult to apply, due to the large decorrelation of the radar signal between the acquisition of images. I developed a model-based unwrapping scheme, which allows the production of reliable displacement maps at the volcano that I used for deformation source modeling. The modeling results show significant differences in pre- and post-eruptive magmatic deformation source parameters. Therefore, I conjecture that two magma chambers exist below Llaima volcano: a post-eruptive deep one and a shallow one possibly due to the pre-eruptive ascent of magma. Similar reservoir depths at Llaima have been confirmed by independent petrologic studies. These reservoirs are interpreted to be temporally coupled. At Tendürek volcano I have found long-term subsidence of the volcanic edifice, which can be described by a large, magmatic, sill-like source that is subject to cooling contraction. The displacement data in conjunction with high-resolution optical images, however, reveal arcuate fractures at the eastern and western flank of the volcano. These are most likely the surface expressions of concentric ring-faults around the volcanic edifice that show low magnitudes of slip over a long time. This might be an alternative mechanism for the development of large caldera structures, which are so far assumed to be generated during large catastrophic collapse events. To investigate the potential subsurface geometry and relation of the two proposed interacting sources at Tendürek, a sill-like magmatic source and ring-faults, I have performed a more sophisticated numerical modeling approach. The optimum source geometries show, that the size of the sill-like source was overestimated in the simple models and that it is difficult to determine the dip angle of the ring-faults with surface displacement data only. However, considering physical and geological criteria a combination of outward-dipping reverse faults in the west and inward-dipping normal faults in the east seem to be the most likely. Consequently, the underground structure at the Tendürek volcano consists of a small, sill-like, contracting, magmatic source below the western summit crater that causes a trapdoor-like faulting along the ring-faults around the volcanic edifice. Therefore, the magmatic source and the ring-faults are also interpreted to be temporally coupled. In addition, a method for data reduction has been improved. The modeling of subsurface deformation sources requires only a relatively small number of well distributed InSAR observations at the earth’s surface. Satellite radar images, however, consist of several millions of these observations. Therefore, the large amount of data needs to be reduced by several orders of magnitude for source modeling, to save computation time and increase model flexibility. I have introduced a model-based subsampling approach in particular for heterogeneously-distributed observations. It allows a fast calculation of the data error variance-covariance matrix, also supports the modeling of time dependent displacement data and is, therefore, an alternative to existing methods.
Changes in the Atlantic overturning circulation have a strong influence on European temperatures, North American sea level and other climate phenomena worldwide. A meaningful assessment of associated societal impacts needs to be based on the full range of its possible future evolution. This requires capturing both the uncertainty in future warming pathways and the inherently long-term response of the ocean circulation. While probabilistic projections of the global mean and regional temperatures exist, process-based probabilistic assessments of large-scale dynamical systems such as the Atlantic overturning are still missing. Here we present such an assessment and find that a reduction of more than 50 % in Atlantic overturning strength by the end of the 21 (s t) century is within the likely range under an unmitigated climate change scenario (RCP8.5). By combining linear response functions derived from comprehensive climate simulations with the full range of possible future warming pathways, we provide probability estimates of overturning changes by the year 2100. A weakening of more than 25 % is found to be very unlikely under a climate protection scenario (RCP2.6), but likely for unmitigated climate change. The method is able to reproduce the modelled recovery caused by climatic equilibration under climate protection scenarios which provides confidence in the approach. Within this century, a reduction of the Atlantic overturning is a robust climatic phenomena that intensifies with global warming and needs to be accounted for in global adaptation strategies.
Coastal flood damage and adaptation costs under 21st century sea-level rise are assessed on a global scale taking into account a wide range of uncertainties in continental topography data, population data, protection strategies, socioeconomic development and sea-level rise. Uncertainty in global mean and regional sea level was derived from four different climate models from the Coupled Model Intercomparison Project Phase 5, each combined with three land-ice scenarios based on the published range of contributions from ice sheets and glaciers. Without adaptation, 0.2-4.6% of global population is expected to be flooded annually in 2100 under 25-123 cm of global mean sea-level rise, with expected annual losses of 0.3-9.3% of global gross domestic product. Damages of this magnitude are very unlikely to be tolerated by society and adaptation will be widespread. The global costs of protecting the coast with dikes are significant with annual investment and maintenance costs of US$ 12-71 billion in 2100, but much smaller than the global cost of avoided damages even without accounting for indirect costs of damage to regional production supply. Flood damages by the end of this century are much more sensitive to the applied protection strategy than to variations in climate and socioeconomic scenarios as well as in physical data sources (topography and climate model). Our results emphasize the central role of long-term coastal adaptation strategies. These should also take into account that protecting large parts of the developed coast increases the risk of catastrophic consequences in the case of defense failure.
Making confident statements about the evolution of an ice-sheet shelf system with a numerical model requires the capability to reproduce the migration of the grounding line. Here we show that the shallow-ice approximation/shallow-shelf approximation hybrid-type Parallel Ice Sheet Model (PISM), with its recent improvements, is capable of modeling the grounding line motion in a perturbed ice-sheet shelf system. The model is set up according to the three-dimensional Marine Ice-Sheet Model Intercomparison Project (MISMIP3d), and simulations are carried out across a broad range of spatial resolutions. Using (1) a linear interpolation of the grounding line with locally interpolated basal friction and (2) an improved driving-stress computation across the grounding line, the reversibility of the grounding line (i.e. its retreat after an advance forced by a local perturbation of basal resistance) is captured by the model even at medium and low resolutions (Delta x > 10 km). The transient model response is qualitatively similar to that of higher-order models but reveals a higher initial sensitivity to perturbations on very short timescales. Our findings support the application of PISM to the Antarctic ice sheet from regional up to continental scales and on relatively low spatial resolutions.
The Isabena catchment (445 km(2)), Spain, features highly diverse spatial heterogeneity in land use, lithology and rainfall. Consequently, the relative contribution in terms of water and sediment yield varies immensely between its subcatchments, and also temporally. This study presents the synthesis of similar to 2.5 years of monitoring rainfall, discharge and suspended sediment concentration (SSC) in the five main subcatchments of the Isabena and its outlet.
Continuous discharge at the subcatchment outlets, nine tipping bucket rainfall and automatic SSC samplers (complemented by manual samples), were collected from June 2011 until November 2013. The water stage records were converted to discharge using a rating curve derived with Bayesian regression. For reconstructing sediment yields, the data from the intermittent SSC sampling needed to be interpolated. We employed non-parametric multivariate regression (Quantile Regression Forests, QRF) using the discharge and rainfall data plus different aggregation levels of these as ancillary predictors. The subsequent Monte Carlo simulations allowed the determination of monthly sediment yields and their uncertainty.
The Isabena catchment shows high erosion dynamics with great variability in space and time, with stark contrasts even between adjacent subcatchments. The natural conditions make water and sediment monitoring and instrumentation very challenging; the measurement of discharge is particularly prone to considerable uncertainties. The QRF method employed for reconstructing sedigraphs and monthly yields proved well suited for the task.
The southern foreland basin of the Alborz Mountains of northern Iran is characterized by an approximately 7.3-km-thick sequence of Miocene sedimentary rocks, constituting three basin-wde coarsening-upward units spanning a period of 10(6)years. We assess available magnetostratigraphy, paleoclimatic reconstructions, stratal architecture, records of depositional environments, and sediment-provenance data to characterize the relationships between tectonically-generated accommodation space (A) and sediment supply (S). Our analysis allows an inversion of the stratigraphy for particular forcing mechanisms, documenting causal relationships, and providing a basis to decipher the relative contributions of tectonics and climate (inferred changes in precipitation) in controlling sediment supply to the foreland basin. Specifically, A/S>1, typical of each basal unit (17.5-16.0, 13.8-13.1 and 10.3-9.6Ma), is associated with sharp facies retrogradation and reflects substantial tectonic subsidence. Within these time intervals, arid climatic conditions, changes in sediment provenance, and accelerated exhumation in the orogen suggest that sediment supply was most likely driven by high uplift rates. Conversely, A/S<1 (13.8 and 13.8-11Ma, units 1, and 2) reflects facies progradation during a sharp decline in tectonic subsidence caused by localized intra-basinal uplift. During these time intervals, climate continued to be arid and exhumation active, suggesting that sediment supply was again controlled by tectonics. A/S<1, at 11-10.3Ma and 9-6-7.6Ma (and possibly 6.2; top of units 2 and 3), is also associated with two episodes of extensive progradation, but during wetter phases. The first episode appears to have been linked to a pulse in sediment supply driven by an increase in precipitation. The second episode reflects a balance between a climatically-induced increase in sediment supply and a reduction of subsidence through the incorporation of the proximal foreland into the orogenic wedge. This in turn caused an expansion of the catchment and a consequent further increase in sediment supply.
Some of the largest and most erosive floods on Earth result from the failure of glacial dams. While potentially cataclysmic ice dams are recognized to have repeatedly formed along ice-sheet margins, much less is known about the frequency and longevity of ice dams caused by mountain glaciers, and their impact on landscape evolution. Here we present field observations and results from cosmogenic nuclide dating that allow reconstructing a > 100-k.y.-long history of glacial damming in the Shyok Valley, eastern Karakoram (South Asia). Our field observations provide evidence that Asia's second-longest glacier, the Siachen, once extended for over 180 km and blocked the Shyok River during the penultimate glacial period, leading to upstream deposition of a more than 400-m-thick fluvio-lacustrine valley fill. Be-10-depth profile modeling indicates that glacial damming ended with the onset of the Eemian interglacial and that the Shyok River subsequently incised the valley fill at an average rate of similar to 4-7 m k.y.(-1). Comparison with contemporary ice-dammed lakes in the Karakoram and elsewhere suggests recurring outburst floods during the aggradation period, while over 25 cycles of fining-upward lake deposits within the valley fill indicate impounding of floods from farther upstream. Despite prolonged damming, the net effect of this and probably earlier damming episodes by the Siachen Glacier is dominated by glacial erosion in excess of fluvial incision, as evidenced by a pronounced overdeepening that follows the glaciated valley reach. Strikingly similar overdeepened valleys at all major confluences of the Shyok and Indus Rivers with Karakoram tributaries indicate that glacial dams and subsequent outburst floods have been widespread and frequent in this region during the Quaternary. Our study suggests that the interaction of Karakoram glaciers with the Shyok and Indus Rivers promoted valley incision and headward erosion into the western margin of the Tibetan Plateau.
Sr isotope records from marginal marine basins track the mixing between seawater and local continental runoff, potentially recording the effects of sea level, tectonic, and climate forcing in marine fossils and sediments. Our 110 new Sr-87/Sr-86 analyses on oyster and foraminifera samples from six late Miocene stratigraphic sections in southern Turkey, Crete, and Sicily show that Sr-87/Sr-86 fell below global seawater values in the basins several million years before the Messinian Salinity Crisis, coinciding with tectonic uplift and basin shallowing. 87Sr/86Sr from more centrally located basins (away from the Mediterranean coast) drop below global seawater values only during the Messinian Salinity Crisis. In addition to this general trend, 55 new Sr-87/Sr-86 analyses from the astronomically tuned Lower Evaporites in the central Apennines (Italy) allow us to explore the effect of glacio-eustatic sea level and precipitation changes on Sr-87/Sr-86. Most variation in our data can be explained by changes in sea level, with greatest negative excursions from global seawater values occurring during relative sea level lowstands, which generally coincided with arid conditions in the Mediterranean realm. We suggest that this greater sensitivity to lowered sea level compared with higher runoff could relate to the inverse relationship between Sr concentration and river discharge. Variations in the residence time of groundwater within the karst terrain of the circum-Mediterranean region during arid and wet phases may help to explain the single (robust) occurrence of a negative excursion during a sea level highstand, but this explanation remains speculative without more detailed paleoclimatic data for the region.
Erosion in the Himalaya is responsible for one of the greatest mass redistributions on Earth and has fueled models of feedback loops between climate and tectonics. Although the general trends of erosion across the Himalaya are reasonably well known, the relative importance of factors controlling erosion is less well constrained. Here we present 25 Be-10-derived catchment-averaged erosion rates from the Yamuna catchment in the Garhwal Himalaya, northern India. Tributary erosion rates range between similar to 0.1 and 0.5mmyr(-1) in the Lesser Himalaya and similar to 1 and 2mmyr(-1) in the High Himalaya, despite uniform hillslope angles. The erosion-rate data correlate with catchment-averaged values of 5 km radius relief, channel steepness indices, and specific stream power but to varying degrees of nonlinearity. Similar nonlinear relationships and coefficients of determination suggest that topographic steepness is the major control on the spatial variability of erosion and that twofold to threefold differences in annual runoff are of minor importance in this area. Instead, the spatial distribution of erosion in the study area is consistent with a tectonic model in which the rock uplift pattern is largely controlled by the shortening rate and the geometry of the Main Himalayan Thrust fault (MHT). Our data support a shallow dip of the MHT underneath the Lesser Himalaya, followed by a midcrustal ramp underneath the High Himalaya, as indicated by geophysical data. Finally, analysis of sample results from larger main stem rivers indicates significant variability of Be-10-derived erosion rates, possibly related to nonproportional sediment supply from different tributaries and incomplete mixing in main stem channels.
Injection of nanoscale zero-valent iron (nZVI) has recently gained great interest as emerging technology for in-situ remediation of chlorinated organic compounds from groundwater systems. Zero-valent iron (ZVI) is able to reduce organic compounds and to render it to less harmful substances. The use of nanoscale particles instead of granular or microscale particles can increase dechlorination rates by-orders of magnitude due to its high surface area. However, classical nZVI appears to be hampered in its environmental application by its limited mobility. One approach is colloid supported transport of nZVI, where the nZVI gets transported by a Mobile colloid. In this study transport properties of activated carbon colloid supported nZVI (c-nZVI; d(50) = 2.4 mu m) are investigated in column tests using columns of 40 cm length, which were filled with porous media. A suspension was pumped through the column under different physicochemical conditions (addition of a polyanionic stabilizer and changes in pH and ionic strength). Highest observed breakthrough was 62% of the injected concentration in glass beads with addition of stabilizer. Addition of mono- and bivalent salt, e.g. more than 0.5 mM/L CaCl2, can decrease mobility and changes in pH to values below six can inhibit mobility at all. Measurements of colloid sizes and zeta potentials show changes in the mean particle size by a factor of ten and an increase of zeta potential from -62 mV to -80 mV during the transport experiment. However, results suggest potential applicability of c-nZVI under field conditions. (C) 2014 Elsevier B.V. All rights reserved.
We present the results of biogeochemical and mineralogical analyses on a sediment core that covers the Holocene sedimentation history of the climatically sensitive, closed, saline, and alkaline Lonar Lake in the core monsoon zone in central India. We compare our results of C/N ratios, stable carbon and nitrogen isotopes, grain-size, as well as amino acid derived degradation proxies with climatically sensitive proxies of other records from South Asia and the North Atlantic region. The comparison reveals some more or less contemporaneous climate shifts. At Lonar Lake, a general long term climate transition from wet conditions during the early Holocene to drier conditions during the late Holocene, delineating the insolation curve, can be reconstructed. In addition to the previously identified periods of prolonged drought during 4.6-3.9 and 2.0-0.6 cal ka that have been attributed to temperature changes in the Indo Pacific Warm Pool, several additional phases of shorter term climate alteration superimposed upon the general climate trend can be identified. These correlate with cold phases in the North Atlantic region. The most pronounced climate deteriorations indicated by our data occurred during 62-5.2,4.6-3.9, and 2.0-0.6 cal ka BP. The strong dry phase between 4.6 and 3.9 cal ka BP at Lonar Lake corroborates the hypothesis that severe climate deterioration contributed to the decline of the Indus Civilisation about 3.9 ka BP. (C) 2014 Elsevier B.V. All rights reserved.
Background and Aims Dynamic processes occurring at the soil-root interface crucially influence soil physical, chemical and biological properties at a local scale around the roots, and are technically challenging to capture in situ. This study presents a novel multi-imaging approach combining fluorescence and neutron radiography that is able to simultaneously monitor root growth, water content distribution, root respiration and root exudation.
Methods Germinated seeds of white lupins (Lupinus albus) were planted in boron-free glass rhizotrons. After 11 d, the rhizotrons were wetted from the bottom and time series of fluorescence and neutron images were taken during the subsequent day and night cycles for 13 d. The following day (i.e. 25 d after planting) the rhizotrons were again wetted from the bottom and the measurements were repeated. Fluorescence sensor foils were attached to the inner sides of the glass and measurements of oxygen and pH were made on the basis of fluorescence intensity. The experimental set-up allowed for simultaneous fluorescence imaging and neutron radiography.
Key Results The interrelated patterns of root growth and distribution in the soil, root respiration, exudation and water uptake could all be studied non-destructively and at high temporal and spatial resolution. The older parts of the root system with greater root-length density were associated with fast decreases of water content and rapid changes in oxygen concentration. pH values around the roots located in areas with low soil water content were significantly lower than the rest of the root system.
Conclusions The results suggest that the combined imaging set-up developed here, incorporating fluorescence intensity measurements, is able to map important biogeochemical parameters in the soil around living plants with a spatial resolution that is sufficiently high enough to relate the patterns observed to the root system.
We used inverse modelling techniques and soil moisture measured by the cosmic-ray neutron sensing (CRS) to estimate root-zone soil hydraulic properties at the field scale. A HYDRUS-1D model was developed for inverse modelling and calibrated with parameter estimation software (PEST) using a global optimizer. Integral CRS measurements recorded from a sunflower farm in Germany comprised the model input. Data were transformed to soil water storage to enable direct model calibration with a HYDRUS soil-water balance. Effective properties at the CRS scale were compared against local measurements and other inversely estimated soil properties from independent soil moisture profiles. Moreover, CRS-scale soil properties were tested on the basis of how field soil moisture (vertical distribution) and soil water storage were reproduced. This framework provided good estimates of effective soil properties at the CRS scale. Simulated soil moisture at different depths at the CRS scale agreed with field observations. Moreover, simulated soil water storage at the CRS scale compared well with calculations from point-scale profiles, despite their different support volumes. The CRS-scale soil properties estimated with the inverse model were within the range of variation of properties identified from all inverse simulations at the local scale. This study demonstrates the potential of CRS for inverse estimation of soil hydraulic properties.
Nanoscale zero-valent iron (nZVI) has recently gained great interest in the scientific community as in situ reagent for installation of permeable reactive barriers in aquifer systems, since nZVI is highly reactive with chlorinated compounds and may render them to harmless substances. However, nZVI has a high tendency to agglomerate and sediment; therefore it shows very limited transport ranges. One new approach to overcome the limited transport of nZVI in porous media is using a suited carrier colloid. In this study we tested mobility of a carbon colloid supported nZVI particle "Carbo-Iron Colloids" (CIC) with a mean size of 0.63 mu m in a column experiment of 40 cm length and an experiment in a two-dimensional (2D) aquifer test system with dimensions of 110x40x5 cm. Results show a breakthrough maximum of 82 % of the input concentration in the column experiment and 58 % in the 2D-aquifer test system. Detected residuals in porous media suggest a strong particle deposition in the first centimeters and few depositions in the porous media in the further travel path. Overall, this suggests a high mobility in porous media which might be a significant enhancement compared to bare or polyanionic stabilized nZVI.
Sedimentary proxy records constitute a significant portion of the recorded evidence that allows us to investigate paleoclimatic conditions and variability. However, uncertainties in the dating of proxy archives limit our ability to fix the timing of past events and interpret proxy record intercomparisons. While there are various age-modeling approaches to improve the estimation of the age-depth relations of archives, relatively little focus has been placed on the propagation of the age (and radiocarbon calibration) uncertainties into the final proxy record.
We present a generic Bayesian framework to estimate proxy records along with their associated uncertainty, starting with the radiometric age-depth and proxy-depth measurements, and a radiometric calibration curve if required. We provide analytical expressions for the posterior proxy probability distributions at any given calendar age, from which the expected proxy values and their uncertainty can be estimated. We illustrate our method using two synthetic data sets and then use it to construct the proxy records for groundwater inflow and surface erosion from Lonar lake in central India.
Our analysis reveals interrelations between the uncertainty of the proxy record over time and the variance of proxies along the depth of the archive. For the Lonar lake proxies, we show that, rather than the age uncertainties, it is the proxy variance combined with calibration uncertainty that accounts for most of the final uncertainty. We represent the proxy records as probability distributions on a precise, error-free timescale that makes further time series analyses and intercomparisons of proxies relatively simple and clear. Our approach provides a coherent understanding of age uncertainties within sedimentary proxy records that involve radiometric dating. It can be potentially used within existing age modeling structures to bring forth a reliable and consistent framework for proxy record estimation.
This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
Deforestation is a prominent anthropogenic cause of erosive overland flow and slope instability, boosting rates of soil erosion and concomitant sediment flux. Conventional methods of gauging or estimating post-logging sediment flux often focus on annual timescales but overlook potentially important process response on shorter intervals immediately following timber harvest. We resolve such dynamics with non-parametric quantile regression forests (QRF) based on high-frequency (3 min) discharge measurements and sediment concentration data sampled every 30-60 min in similar-sized (similar to 0.1 km(2)) forested Chilean catchments that were logged during either the rainy or the dry season. The method of QRF builds on the random forest algorithm, and combines quantile regression with repeated random sub-sampling of both cases and predictors. The algorithm belongs to the family of decision-tree classifiers, which allow quantifying relevant predictors in high-dimensional parameter space. We find that, where no logging occurred, similar to 80% of the total sediment load was transported during extremely variable runoff events during only 5% of the monitoring period. In particular, dry-season logging dampened the relative role of these rare, extreme sediment-transport events by increasing load efficiency during more efficient moderate events. We show that QRFs outperform traditional sediment rating curves (SRCs) in terms of accurately simulating short-term dynamics of sediment flux, and conclude that QRF may reliably support forest management recommendations by providing robust simulations of post-logging response of water and sediment fluxes at high temporal resolution.