Irrigation and warming drive the decreases in surface albedo over High Mountain Asia.

Human and climate induced land surface changes resulting from irrigation, snow cover decreases, and greening impact the surface albedo over High Mountain Asia (HMA). Here we use a partial information decomposition approach and remote sensing data to quantify the effects of the changes in leaf area index, soil moisture, and snow cover on the surface albedo in HMA, home to over a billion people, from 2003 to 2020. The study establishes strong evidence of anthropogenic agricultural water use over irrigated lands (e.g., Ganges-Brahmaputra) which causes the highest surface albedo decreases (≤ 1%/year). Greening and decreased snow cover from warming also drive changes in visible and near-infrared surface albedo in different areas of HMA. The significant role of irrigation and greening in influencing albedo suggests the potential of a positive feedback cycle where albedo decreases lead to increased evaporative demand and increased stress on water resources.

Modeling geogenic and atmospheric nitrogen through the East River Watershed, Colorado Rocky Mountains.

There is a growing understanding of the role that bedrock weathering can play as a source of nitrogen (N) to soils, groundwater and river systems. The significance is particularly apparent in mountainous environments where weathering fluxes can be large. However, our understanding of the relative contributions of rock-derived, or geogenic, N to the total N supply of mountainous watersheds remains poorly understood. In this study, we develop the High-Altitude Nitrogen Suite of Models (HAN-SoMo), a watershed-scale ensemble of process-based models to quantify the relative sources, transformations, and sinks of geogenic and atmospheric N through a mountain watershed. Our study is based in the East River Watershed (ERW) in the Upper Colorado River Basin. The East River is a near-pristine headwater watershed underlain primarily by an N-rich Mancos Shale bedrock, enabling the timing and magnitude of geogenic and atmospheric contributions to watershed scale dissolved N-exports to be quantified. Several calibration scenarios were developed to explore equifinality using >1600 N concentration measurements from streams, groundwater, and vadose zone samples collected over the course of four years across the watershed. When accounting for recycling of N through plant litter turnover, rock weathering accounts for approximately 12% of the annual dissolved N sources to the watershed in the most probable calibration scenario (0-31% in other scenarios), and 21% (0-44% in other scenarios) when considering only "new" N sources (i.e. geogenic and atmospheric). On an annual scale, instream dissolved N elimination, plant turnover (including cattle grazing) and atmospheric deposition are the most important controls on N cycling.

Benchmarking numerical codes for tracer transport with the aid of laboratory-scale experiments in 2D heterogeneous porous media.

We present a combined experimental and numerical modeling study that addresses two principal questions: (i) is any particular Eulerian-based method used to solve the classical advection-dispersion equation (ADE) clearly superior (relative to the others), in terms of yielding solutions that reproduce BTCs of the kind that are typically sampled at the outlet of a laboratory cell? and (ii) in the presence of matches of comparable quality against such BTCs, do any of these methods render different (or similar) numerical BTCs at locations within the domain? To address these questions, we obtained measurements from carefully controlled laboratory experiments, and employ them as a reference against which numerical results are benchmarked and compared. The experiments measure solute transport breakthrough curves (BTCs) through a square domain containing various configurations of coarse, medium, and fine quartz sand. The approaches to solve the ADE involve Eulerian-Lagrangian and Eulerian (finite volume, finite elements, mixed and discontinuous finite elements) numerical methods. Model calibration is not examined; permeability and porosity of each sand were determined previously through separate, standard laboratory tests, while dispersivities are assigned values proportional to mean grain size. We find that the spatial discretization of the flow field is of critical importance, due to the non-uniformity of the domain. Although simulated BTCs at the system outlet are observed to be very similar for these various numerical methods, computed local (point-wise, inside the domain) BTCs can be very different. We find that none of the numerical methods is able to fully reproduce the measured BTCs. The impact of model parameter uncertainty on the calculated BTCs is characterized through a set of numerical Monte Carlo simulations; in cases where the impact is significant, assessment of simulation matches to the experimental data can be ambiguous.