Soil is a critical natural resource - the foundation of conventional agriculture and terrestrially foraged food, battery for below-ground carbon and nutrient storage, filter for water, and support system for most other ecosystem services. Fundamental questions remain regarding soil sustainability under rapidly-changing climatic and environmental conditions, including massive human impacts and climate change. Accelerated rates of soil erosion threaten the stability of ecosystems, nutrient and carbon cycles, clean water reserves and global food supplies if the processes that produce soil cannot keep pace. Weathering of underlying bedrock substrate helps dislodge rock fragments and prepares particles for movement and incorporation of minerals into the mobile regolith layer, but soil is produced from below and above. The flux of weathered minerals into the near-surface zone acts from below, but the incorporation of dust, transported soil and organic matter near the surface of the soil amends soil from above. Thus, the total amount of soil represents the balance between these soil formation mechanisms and organic matter degradation from physical erosion and chemical transformations in the soil. Sustainability of soil over long timescales requires a positive balance of inputs and outputs. A mismatch of erosion rate and soil production rate can persist for a while before badlands are created, but the buffering timescale depends on the thickness of the existing soil stocks. Fortunately, over millennial timescales, the rate of soil production is thought to keep pace with the rate of surface erosion through negative feedbacks between soil thickness and the rate at which soil is produced from the underlying mineral substrate. This paradigm in the Earth Sciences holds that an underlying mechanism lowers the rate of soil production when soil is thick and increases the rate of soil production when soils are thin. This dynamic balance lends support to two observations: first, soil covers >90% of Earth‘s ice-free surface despite natural soil erosion rates that vary by three orders of magnitude, and second, the thickness of soils on Earth exists within a relatively narrow range even in old and deeply weathered landscapes. However, the actual coupling mechanism between soil thickness and soil production is still unknown, and the environmental controls governing the maximum soil depth or maximum production rate in a given landscape are also yet to be discovered. Unfortunately, a new view is emerging that previously reported rates of soil production were subject to an artifact of the method and that increases in erosion may not give rise to a concomitant, matched rise in soil production. Instead, rates of soil production seem to be set by soil ‘age‘ and the aridity index, i.e. the amount of water left from precipitation after evapotranspiration is accounted for. These controls on soil production are either set by landscape position or features of the Earth system that may be swinging out of balance due to climate change or climate change-induced erosion processes. If correct, this result would have drastic consequences for soil sustainability in the context of anthropogenically accelerated soil erosion. If we fast-forward with current rates of erosion and new soil production rates, how long until soils themselves are endangered?

Rocky coast cliff retreat presents a hazard to coastal communities and infrastructure that is potentially amplified under rising sea level conditions, among other factors. Unfortunately, constraints on retreat rates are typically limited to those derived from imagery and maps spanning the last ∼100 years. Here, we use a newly developed coupled model of shore platform profile evolution and cosmogenic radionuclide production that considers the influence of relative sea level (RSL) rise, weathering, material resistance, and wave height decay to model cliff retreat over millennial timescales and its potential drivers in Del Mar, California. We demonstrate the ability to use topographic and bathymetric measurements from a narrow shore platform along with a limited data set of nine cosmogenic 10 Be concentrations extending ∼125 m from the cliff base to obtain modeled cliff retreat rates that steadily range from 5.0 to 12.5 cm yr −1 over the last two millennia until 100 years before present. These rates are consistent with modern retreat rates of about 2–19 cm yr −1 here. RSL rise in Southern California remained relatively constant during the late Holocene, potentially explaining the relatively stable modeled cliff retreat rate over this time. We also explore the relative influence of weathering, material resistance, and wave erosion efficacy and find that both weathering and wave‐driven erosion are necessary to replicate the measured data at this location, with the latter exerting a stronger control on model acceptance, suggesting that waves may provide a possible mechanism by which RSL rise may influence coastal cliff erosion in southern California.

Datasets, regression data, and descriptions of minor corrections to the datasets, and a detailed description of the statistical methods.

The lateral migration of a river meander is driven by erosion on the outer bank and deposition on the inner bank, both of which are affected by shear stress (and therefore channel slope) through complex morphodynamic feedbacks. To test the sensitivity of lateral migration to channel slope, we quantify slope change induced by glacial isostatic adjustment along the Red River (North Dakota, USA and Manitoba, Canada) and two of its tributaries over the past 8.5 ka. We demonstrate a statistically significant, positive relationship between normalized cutoff count, which we interpret as a proxy for channel lateral migration rate, and slope change. We interpret this relationship as the signature of slope change modulating the magnitude of shear stress on riverbanks, suggesting that slope changes that occur over thousands of years are recorded in river floodplain morphology.

Topography is commonly viewed as a passive backdrop on which vegetation grows. Yet, in certain circumstances, a bidirectional feedback may develop between the control of topography and the spatial distribution of vegetation and landform development, because vegetation modulates the erosion of the land surface. Therefore, if reinforcing feedbacks are established between erosion and land cover distribution over timescales relevant to landform development, then the interactions between vegetation and topography may create distinctive landforms, shaped by vegetation. We expose here a strong correlation between the spatial distribution of vegetation, erosion rates, and topography at a characteristic length scale of 102-103m (mesoscale topography) in the Luquillo Experimental forest (LEF) of Puerto Rico. We use high-resolution LiDAR topography to characterize landforms, satellite images to classify the vegetation into forest types, and in-situ produced cosmogenic 10Be in the quartz extracted from soils and stream sediments to document spatial variations in soil erosion. The data document a strong correlation between forest type and topographic position (hilltop vs. valleys), and a correlation between topographic position and 10Be-derived erosion rates over 103−104 years. Erosion is faster in valleys, which are mostly covered by monocot Palm Forest, and slower on surrounding hills mostly covered by the dicot Palo Colorado Forest. Transition from one forest type to the next occurs across a break-in-slope that separates shallowly convex hilltops from deeply concave valleys (coves). The break-in-slope is the consequence of a longer-lasting erosional imbalance whereby coves erode faster than hills over landscape-shaping timescales. Such a deepening of the coves is usually spurred by external drivers, but such drivers are here absent. This implies that cove erosion is driven by a process originating within the coves themselves. We propose that vegetation is the primary driver of this imbalance, soil erosion being faster under Palm forest than under Palo Colorado forest. Concentration of the Palm forest in the deepening coves is reinforced by the better adaptation of Palm trees to the erosive processes that take place in the coves, once these develop steep slopes. At the current rate of landscape development, we find that the imbalance started within the past 0.1–1.5 My. The initiation of the process could correspond to time of settlement of these mountain slopes by the Palm and Palo Colorado forests.

The exceptional concentration of vertebrate diversity in continental freshwaters has been termed the "freshwater fish paradox," with > 15,000 fish species representing more than 20% of all vertebrate species compressed into tiny fractions of the Earth's land surface area (<0.5%) or total aquatic habitat volume (<0.001%). This study asks if the fish species richness of the world's river basins is explainable in terms of river captures using topographic metrics as proxies. The River Capture Hypothesis posits that drainage-network rearrangements have accelerated biotic diversification through their combined effects on dispersal, speciation, and extinction. Yet rates of river capture are poorly constrained at the basin scale worldwide. Here we assess correlations between fish species density (data for 14,953 obligate freshwater fish species) and basin-wide metrics of landscape evolution (data for 3,119 river basins), including: topography (elevation, average relief, slope, drainage area) and climate (average rainfall and air temperature). We assess the results in the context of both static landscapes (e.g., species-area and habitat heterogeneity relationships) and transient landscapes (e.g., river capture, tectonic activity, landscape disequilibrium). We also relax assumptions of functional neutrality of basins (tropical vs. extratropical, tectonically stable vs. active terrains). We found a disproportionate number of freshwater species in large, lowland river basins of tropical South America, Africa, and Southeast Asia, under predictable conditions of large geographic area, tropical climate, low topographic relief, and high habitat volume (i.e., high rainfall rates). However, our results show that these conditions are only necessary, but not fully sufficient, to explain the basins with the highest diversity. Basins with highest diversity are all located on tectonically stable regions, places where river capture is predicted to be most conducive to the formation of high fish species richness over evolutionary timescales. Our results are consistent with predictions of several landscape evolution models, including the River Capture Hypothesis, Mega Capture Hypothesis, and Intermediate Capture Rate Hypothesis, and support conclusions of numerical modeling studies indicating landscape transience as a mechanistic driver of net diversification in riverine and riparian organisms with widespread continental distributions.

Measuring rates of valley head migration and determining the timing of canyon-opening are insightful for the evolution of planetary surfaces. Spatial gradients of in situ-produced cosmogenic nuclide concentrations along horizontal transects provide a framework for assessing the migration of valley networks and similar topographic features. We developed a new derivation for valley head retreat rates from the concentration of in situ produced cosmogenic radionuclides in valley walls. The retreat rate is inversely proportional to the magnitude of the spatial concentration gradient and proportional to local nuclide production rates. By solving for a spatial gradient in concentration along a valley parallel transect, we created an expression for the explicit determination of valley head retreat, which we refer to herein as unzipping. We applied this expression to a seepage-derived drainage network developing along the Apalachicola River, Florida, USA. Sample concentrations along a valley margin transect varied systematically from 2.9x105 atoms/g to 3.5x105 atoms/g resulting in a gradient of 160 atoms/g/m, and from this value a valley head retreat rate of 0.025 m/y was found. The discrepancy between overall network age and current rates of valley head migration suggests intermittent network growth which is consistent with glacial-interglacial precipitation variations during the Pleistocene. This method can be applied to a wide range of Earth-surface environments. For the 10Be system, this method should be sensitive to unzipping rates bounded between 10-6 m/y and 100 m/y.

Over the last decade, an increasing number of studies have used soundscapes to address diverse ecological questions. Sound represents one of the few sources of information capable of providing in situ insights into processes occurring within opaque soil matrices. To date, the use of soundscapes for soil macrofauna monitoring has been experimentally tested only in controlled laboratory environments. Here we assess the validity of laboratory predictions and explore the use of soil soundscape proxies for monitoring soil macrofauna (i.e., earthworm) activities in an outdoor context. In a common garden experiment in northern Sweden, we constructed outdoor mesocosm plots (N = 36) containing two different Arctic vegetation types (meadow and heath) and introduced earthworms to half of these plots. Earthworms substantially altered the ambient soil soundscape under both vegetation types, as measured by both traditional soundscape indices and frequency band power levels, although their acoustic impacts were expressed differently in heath versus meadow soils. While these findings support the as-of-yet untapped promise of using belowground soundscape analyses to monitor soil ecosystem health, direct acoustic emissions from earthworm activities appear to be an unlikely proxy for tracking worm activities at daily timescales. Instead, earthworms indirectly altered the soil soundscape by 're-engineering' the soil matrix: an effect that was dependent on vegetation type. Our findings suggest that long-term (i.e., seasonal) earthworm activities in natural soil settings can likely be monitored indirectly via their impacts on soundscape measures and acoustic indices. Analyzing soil soundscapes may enable larger-scale monitoring of high-latitude soils and is directly applicable to the specific case of earthworm invasions within Arctic soils, which has recently been identified as a potential threat to the resilience of high-latitude ecosystems. Soil soundscapes could also offer a novel means to monitor soils and soil-plant-faunal interactions in situ across diverse pedogenic, agronomic, and ecological systems.

Active, Cordilleran-type margins are distributed from Patagonia to Alaska and the western limits of South, Central, and North America and host some of the tallest and longest mountain belts on Earth. Cordilleran-type margins form a highly dynamic subduction zone system in which plate velocities and slab dynamics (i.e., rollback and break offs) affect coupling and decoupling of the interacting plates and this, in turn, controls deformation and mountain building in the overriding plate and quasi-cyclic mountain belt evolution. Cyclicity in the orogenic cycle in subduction zones and processes related to climate and erosion create dynamic feedbacks that remain outstanding questions in tectonic geomorphology.

Over the last decade, an increasing number of studies have used soundscapes to address diverse ecological questions. Sound represents one of the few sources of information capable of providing in situ insights into processes occurring within opaque soil matrices. To date, the use of soundscapes for soil macrofauna monitoring has been experimentally tested only in controlled laboratory environments. Here we assess the validity of laboratory predictions and explore the use of soil soundscape proxies for monitoring soil macrofauna ( i.e. , earthworm) activities in an outdoor context. In a common garden experiment in northern Sweden, we constructed outdoor mesocosm plots (N= 36) containing two different Arctic vegetation types (meadow and heath) and introduced earthworms to half of these plots. Earthworms substantially altered the ambient soil soundscape under both vegetation types, although their acoustic impacts were expressed differently in heath versus meadow soils. While these findings support the as-of-yet untapped promise of using belowground soundscape analyses to monitor soil ecosystem health, direct acoustic emissions from earthworm activities appear to be an unlikely proxy for tracking worm activities at daily timescales. Instead, earthworms indirectly altered the soil soundscape by ‘re-engineering’ the soil matrix: an effect that was dependent on vegetation type. Our findings suggest that long-term ( i.e. , seasonal) earthworm activities in natural soil settings can likely be monitored indirectly via their impacts on soundscape measures and acoustic indices. Analyzing soil soundscapes may enable larger-scale monitoring of high-latitude soils and is directly applicable to the specific case of earthworm invasions within Arctic soils, which has recently been identified as a potential threat to the resilience of high-latitude ecosystems. Soil soundscapes could also offer a novel means to monitor soils and soil-plant-faunal interactions in situ across diverse pedogenic, agronomic, and ecological systems.

Cosmogenic nuclide techniques have advanced the geosciences by providing tools for exposure age dating, burial dating, quantification of denudation rates and more. Advances in geochemistry, accelerator mass spectrometry and atom trap trace analyses are ushering in a new cosmogenic nuclide era, by improving the sensitivity of measurements to ultra-trace levels that now allow new applications of these techniques to numerous Earth surface processes. The advances in cosmogenic nuclide techniques have equipped the next generation of geoscientists with invaluable tools for understanding the planet, but addressing pressing needs requires rising to an even greater challenge: imbuing within the cosmogenic community, and the geosciences as a whole, a commitment to justice, equity, diversity and inclusion that matches our dedication to scientific research. In this Primer, we review the state of the art and recent exciting breakthroughs in the use of cosmogenic nuclide techniques, focusing on erosion factories over space and time, and new perspectives on ice sheet stability. We also highlight promising ways forward in enhancing inclusion in the field, as well as obstacles that remain to be overcome.

The rise of a mountain range affects moisture circulation in the atmosphere and water runoff across the land surface, modifying the distribution of precipitation and drainage patterns in its vicinity. Water routing in turn affects erosion on hillslopes and incision in river channels on surrounding mountain ranges. In central Guatemala, two parallel, closely spaced mountain ranges formed during two consecutive pulses of uplift, the first between 12 and 7 Ma (Sierra de Chuacús–Sierra de las Minas), and the second after 7 Ma (Altos de Cuchumatanes). We explore the climatic and tectonic processes through which the rise of the most recent range drove the slowing of river incision and hillslope erosion over the previously uplifted range. The 40Ar/39Ar dating of perched volcanic deposits documents the sequential rise and incision of these mountain ranges. Terrestrial cosmogenic 10Be in river sediments indicates that currently hillslopes in the older range erode more slowly than in the younger range (20–150 vs. 300 m Myr−1). These differences mimic the current distribution of precipitation, with the younger range intercepting the atmospheric moisture before it reaches the older range. River channel steepness and deformation of paleovalleys in the new range further indicate that the younger range has been rising faster than the older range up to today. We review how atmospheric moisture interception and river long-profile adjustment to the rise of the new range have contributed to the decline of erosion rates over the old range. We also explore the consequences of this decline and of aridification on the topographic evolution of the older range. The older range undergoes a slow topographic decay, dominated by backwearing, by the stacking of slowly migrating erosion waves along the mountain flanks, and by the formation of pediments around its base. The morphology of the old range is therefore transitioning from that of a front range to that of a dry interior range.

Accurately inferring erosion rates from cosmogenic isotope concentrations in river sand assumes temporally steady concentrations; few studies test this assumption. Following Hurricane Maria in Puerto Rico, we quantified temporal variability in meteoric and in situ Be-10 (Be-10(m), Be-10(i)) on sand-sized grains of riverine transported material in landslide-prone basins. We analyzed 20 samples collected over 18 months from the channels of two nested watersheds: the Icacos (3.14 km(2), 0.09% active landslide area) and Guaba basins (0.11 km(2), 1.23% active landslide area). Be-10(i) concentrations in Icacos basin sediment remained steady over time whereas concentrations in Guaba basin sediment were initially half those in the Icacos basin and increased linearly over 18 months, constraining recovery time to <2 yrs for this basin. Be-10(m) concentrations in both drainages did not change consistently over time and were not related to precipitation events; Pb-210(ex) and Cs-137 were below detection limits in all samples. Our data demonstrate that Be-10(i) concentrations in river sand can be lowered for months to years after major landscape disturbing events, such as large or extensive mass movements. Sampling soon after a landslide will result in over-estimates of long-term erosion rates. Such bias can be reduced by repeated sampling over time and by sampling numerous similar watersheds of different sizes and different concentrations of landslides in a study area. (C) 2021 Elsevier B.V. All rights reserved.

Meteoric 10Be and 7Be produced in the atmosphere from high-energy spallation reactions are deposited onto the Earth's surface through wet and dry deposition and are sorbed onto the surfaces of particles. On land, the sorbed concentrations scale with the residence time of sediments in a landscape-offset by slow (10Be) and fast (7Be) radioactive decay. Additionally, the amount of native 9Be, leached from minerals, correlates with the chemical weathering of soils. However, previous work has shown that chemical and physical properties of soils and river sediments affects sorption of beryllium. Therefore, the magnitude of sorbed beryllium concentrations may be more representative of the sorption capacity of the system rather than its erosional or weathering history. Although previous work has examined the physical and chemical properties of soil that influence beryllium sorption, these studies either lack consensus or exclude potentially important variables. In this work, we provide a thorough examination of variables previously reported to have influence on beryllium chemistry as well as new variables such as nitrogen, phosphorus and sulfur concentrations in order to determine which factors best predict beryllium sorption. We selected two soil endmembers with differing compositions, separated them into different size fractions, and characterized the surface area, cation exchange capacity (CEC), mineralogy, sulfur, carbon, nitrogen and phosphorus concentrations. We determined that the inverse percent abundance of quartz and the CEC best predict beryllium sorption potential in these soils. By deriving a model that relates these two variables to the percent sorbed beryllium, we were able to predict the sorption capacity of our system and reduced the error in sorbed beryllium amounts due to differences in soil properties by about 42%. From these results, we provide insight as to why there is inconsistency in the literature with regards to the physio-chemical controls on the environmental behavior of beryllium.

Abstract. We report previously unpublished evidence for a Marine Isotope Stage 3 (MIS 3; 60–26 ka) glacial outburst flood in the Torngat Mountains (northern Quebec/Labrador, Canada). We present 10Be cosmogenic exposure ages from legacy fieldwork for a glacial lake shoreline with evidence for outburst flooding in the Torngat Mountains, with a minimum age of 36 ± 3 ka (we consider the most likely age, corrected for burial, to be ~56 ± 3 ka). This shoreline position and age can potentially constrain the Laurentide Ice Sheet margin in the Torngat Mountains. This region, considered a site of glacial inception, has no published dated geologic constraints for high-elevation MIS 3 ice margins. We estimate the freshwater flux associated with the inferred glacial outburst flood using high-resolution digital elevation maps corrected for glacial isostatic adjustment. Using assumptions about the ice-dammed locations we find that a freshwater flood volume of 1.14 × 1012 m3 could have entered the Hudson Strait. This glacial outburst flood volume could have contributed to surface ocean freshening to cause a measurable meltwater signal in δ18O records, but would not necessarily have been associated with substantial ice rafted debris. Future work is required to refine estimates of the size and timing of such a glacial outburst flood. Nevertheless, we outline testable hypotheses about the Laurentide Ice Sheet and glacial outburst floods, including possible implications for Heinrich events and glacial inception in North America, that can be assessed with additional fieldwork and cosmogenic measurements.

We present cospatial, contemporary archival records of biotic and abiotic terrestrial processes operating over the past similar to 25 ky within the Rio Fajardo watershed, in northeastern Puerto Rico. The proxy records were derived from a 5-m-thick stratigraphic section exposed by cut bank incision. We interpreted ecosystem dynamics from changes in the stable carbon isotopic ratio of sedimentary organic material compared to delta C-13 ratios of contemporary carbon sources. Sedimentary organic material had delta C-13 values ranging from -29.715 to -15.291. We derived a record of paleo-erosion rates in the catchment from the concentration of meteoric Be-10 in layers of the floodplain sediments. Paleo-erosion rates ranged from 13 to 356 mm ky(-1). The chronology of the sediments was constrained with the radiocarbon ages of organic deposits, the oldest age was calibrated to similar to 22.4 ky BP (thousand years before present) and retrieved at 440 cm depth. We collected grain size data, clay mineralogies, and analyzed geochemical indices including the chemical weathering index, salinization, and base cation loss down profile. This stratigraphic sequence captures major shifts in the Caribbean climate, the intensification of the El Nino-Southern Oscillation, and the arrival of humans on the island. During the last glacial and early Holocene epochs both biotic (delta C-13) and abiotic proxies (Be-10(met) and geochemical data) indicated dynamic equilibrium with climate. The past five thousand years (ky) of record are characterized instead by pulsed responses to disturbances in both systems. Colonial-era land use drove changes that significantly exceeded natural variability in any proxy over the period of record. (C) 2020 Elsevier Ltd. All rights reserved.

To prevent exposure of millions of people living near thousands of contaminated sites to asbestos, the sites are typically capped with soil and amendments, assuming asbestos mobility underground is negligible in all conditions. Here, we disprove this long-standing assumption and showed that the presence of certain types of dissolved organic matter (DOM) can facilitate the transport of the asbestos fibers through packed sand or soil columns. The column scale simulates asbestos transport without and without DOM, although asbestos transport at a longer scale through aquifer may vary because of soil heterogeneity. Our data shows that the extent of fiber transport in soil is affected by DOM quality: fulvic acid > humic acid > natural organic matter isolated from Suwanee River, a model terrestrial organic carbon. We attribute the results to the ability of DOM to affect aggregation of asbestos fibers and limit their attachment on soil. Thus, types of organic amendments applied on asbestos-contaminated sites may have an unintended consequence: transport of asbestos fibers to shallow groundwater to receiving streams or rivers, from where they could be resuspended in the air via irrigation or drying of the riverbed.

Continental sedimentation was widespread across the Central Anatolian Plateau in Miocene–Pliocene time, during the early stages of plateau uplift. Today, however, most sediment produced on the plateau is dispersed by a well-integrated drainage and released into surrounding marine depocenters. Residual long-term (106–107 yr) sediment storage on the plateau is now restricted to a few closed catchments. Lacustrine sedimentation was widespread in the Miocene–Pliocene depocenters. Today, it is also restricted to the residual closed catchments. The present-day association of closed catchments, long-term sediment storage, and lacustrine sedimentation suggests that the Miocene–Pliocene sedimentation also occurred in closed catchments. The termination of sedimentation across the plateau would therefore mark the opening of these closed catchments, their integration, and the formation of the present-day drainage. By combining newly dated volcanic markers with previously dated sedimentary sequences, we show that this drainage integration occurred remarkably rapidly, within 1.5 m.y., at the turn of the Pliocene. The evolution of stream incision documented by these markers and newly obtained 10Be erosion rates allow us to discriminate the respective contributions of three potential processes to drainage integration, namely, the capture of closed catchments by rivers draining the outer slopes of the plateau, the overflow of closed lakes, and the avulsion of closed catchments. Along the southern plateau margin, rivers draining the southern slope of the Central Anatolian Plateau expanded into the plateau interior; however, only a small amount of drainage integration was achieved by this process. Instead, avulsion and/or overflow between closed catchments achieved most of the integration, and these top-down processes left a distinctive sedimentary signal in the form of terminal lacustrine limestone sequences. In the absence of substantial regional climate wetting during the early Pliocene, we propose that two major tectonic events triggered drainage integration, separately or in tandem: the uplift of the Central Anatolian Plateau and the tectonic completion of the Anatolian microplate. Higher surface uplift of the eastern Central Anatolian Plateau relative to the western Central Anatolian Plateau promoted more positive water balances in the eastern catchments, higher water discharge, and larger sediment fluxes. Overflow/avulsion in some of the eastern catchments triggered a chain of avulsions and/or overflows, sparking sweeping integration across the plateau. Around 5 Ma, the inception of the full escape of the Anatolian microplate led to the disruption of the plateau surface by normal and strike-slip faults. Fault scarps partitioned large catchments fed by widely averaged sediment and water influxes into smaller catchments with more contrasted water balances and sediment fluxes. The evolution of the Central Anatolian Plateau shows that top-down processes of integration can outcompete erosion of outer plateau slopes to reintegrate plateau interior drainages, and this is overlooked in current models, in which drainage evolution is dominated by bottom-up integration. Top-down integration has the advantage that it can be driven by more subtle changes in climatic and tectonic boundary conditions than bottom-up integration.