Global climate is thought to be modulated by the supply of minerals to Earth’s surface. Whereas silicate weathering removes carbon dioxide (CO₂) from the atmosphere, weathering of accessory carbonate and sulfide minerals is a geologically relevant source of CO₂. Although these weathering pathways commonly operate side by side, we lack quantitative constraints on their co-variation across erosion rate gradients. Here we use stream-water chemistry across an erosion rate gradient of three orders of magnitude in shales and sandstones of southern Taiwan, and find that sulfide and carbonate weathering rates rise with increasing erosion, while silicate weathering rates remain steady. As a result, on timescales shorter than marine sulfide compensation (approximately 10⁶–10⁷ years), weathering in rapidly eroding terrain leads to net CO₂ emission rates that are about twice as fast as CO₂ sequestration rates in slow-eroding terrain. We propose that these weathering reactions are linked and that sulfuric acid generated from sulfide oxidation boosts carbonate solubility, whereas silicate weathering kinetics remain unaffected, possibly due to efficient buffering of the pH. We expect that these patterns are broadly applicable to many Cenozoic mountain ranges that expose marine metasediments.

Experimental delta lobe size is controlled by bed adjustment to transient floods within the backwater zone.

River deltas grow in large part through repeated cycles of lobe construction and channel avulsion. Understanding avulsion cycles is important for coastal restoration and ecology, land management, and flood hazard mitigation. Emerging theories suggest that river avulsions on lowland deltas are controlled by backwater hydrodynamics; however, our knowledge of backwater‐controlled avulsion cycles is limited. Here we present results from an experimental delta that evolved under persistent backwater hydrodynamics achieved through variable flood discharges, shallow bed slopes, and subcritical flows. The experimental avulsion cycles consisted of an initial phase of avulsion setup, an avulsion trigger, selection of a new flow path, and abandonment of the parent channel. Avulsions were triggered during the largest floods (78% of avulsions) after the channel was filled by a fraction (0.3 ± 0.13) of its characteristic flow depth at the avulsion site, which occurred in the upstream part of the backwater zone. The new flow path following avulsion was consistently one of the shortest paths to the shoreline, and channel abandonment occurred through temporal decline in water flow and sediment delivery to the parent channel. Experimental synthetic stratigraphy indicates that bed thicknesses were maximum at the avulsion sites, consistent with our morphologic measurements of avulsion setup and the idea that there is a record of avulsion locations and thresholds in sedimentary rocks. Finally, we discuss the implications of our findings within the context of sustainable management of deltas, their stratigraphic record, and predicting avulsions on deltas.

Among the most urgent challenges in future climate change scenarios is accurately predicting the magnitude to which precipitation extremes will intensify. Analogous changes have been reported for an episode of millennial-scale 5 °C warming, termed the Palaeocene-Eocene Thermal Maximum (PETM; 56 Ma), providing independent constraints on hydrological response to global warming. However, quantifying hydrologic extremes during geologic global warming analogs has proven difficult. Here we show that water discharge increased by at least 1.35 and potentially up to 14 times during the early phase of the PETM in northern Spain. We base these estimates on analyses of channel dimensions, sediment grain size, and palaeochannel gradients across the early PETM, which is regionally marked by an abrupt transition from overbank palaeosol deposits to conglomeratic fluvial sequences. We infer that extreme floods and channel mobility quickly denuded surrounding soil-mantled landscapes, plausibly enhanced by regional vegetation decline, and exported enormous quantities of terrigenous material towards the ocean. These results support hypotheses that extreme rainfall events and associated risks of flooding increase with global warming at similar, but potentially at much higher, magnitudes than currently predicted.

Hydrologic connectivity controls the lateral exchange of water, solids, and solutes between rivers and floodplains, and is critical to ecosystem function, water treatment, flood attenuation, and geomorphic processes. This connectivity has been well‐studied, typically through the lens of fluvial flooding. In regions prone to heavy rainfall, the timing and magnitude of lateral exchange may be altered by pluvial flooding on the floodplain. We collected measurements of flow depth and velocity in the Trinity River floodplain in coastal Texas (USA) during Tropical Storm Imelda (2019), which produced up to 75 cm of rainfall locally. We developed a two‐dimensional hydrodynamic model at high resolution for a section of the Trinity River floodplain inspired by the compound flooding of Imelda. We then employed Lagrangian particle routing to quantify how residence times and particle velocities changed as flooding shifted from rainfall‐driven to river‐driven. Results show that heavy rainfall initiated lateral exchange before river discharge reached flood levels. The presence of rainwater also reduced floodplain storage, causing river water to be confined to a narrow corridor on the floodplain, while rainwater residence times were increased from the effect of high river flow. Finally, we analyzed the role of floodplain channels in facilitating surface‐water connectivity by varying model resolution in the floodplain. While the resolution of floodplain channels was important locally, it did not affect as much the overall floodplain behavior. This study demonstrates the complexity of floodplain hydrodynamics under conditions of heavy rainfall, with implications for sediment deposition and nutrient removal during floods.

The Himalayan mountain range produces one of the steepest and largest rainfall gradients on Earth, with >3 m/yr rainfall difference over a ∼100 km distance. The Indian Summer Monsoon (ISM) contributes more than 80% to the annual precipitation budget of the central Himalayas. The remaining 20% falls mainly during pre‐ISM season. Understanding the seasonal cycle and the transfer pathways of moisture from precipitation to the rivers is crucial for constraining water availability in a warming climate. However, the partitioning of moisture into the different storage systems such as snow, glacier, and groundwater and their relative contribution to river discharge throughout the year remains under‐constrained. Here, we present novel field data from the Kali Gandaki, a trans‐Himalayan river, and use 4‐year time series of river and rain water stable isotope composition (δ¹⁸O and δ²H values) as well as river discharge, satellite Global Precipitation Measurement amounts, and moisture source trajectories to constrain hydrological variability. We find that rainfall before the onset of the ISM is isotopically distinct and that ISM rain and groundwater have similar isotopic values. Our study lays the groundwork for using isotopic measurements to track changes in precipitation sources during the pre‐ISM to ISM transition in this key region of orographic precipitation. Specifically, we highlight the role of pre‐ISM precipitation, derived from the Gangetic plain, to define the seasonal river isotopic variability across the central Himalayas. Lastly, isotopic values across the catchment document the importance of a large well‐mixed groundwater reservoir supplying river discharge, especially during the non‐ISM season.