Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

The link between abrupt climate change and basin stratigraphy: a numerical approach

To use basin stratigraphy for studying past climate change, it is important to understand the influence of evolving boundary conditions (river discharge and sediment flux, initial bathymetry, sea level, subsidence) and the complex interplay of the redistribution processes (plumes, turbidity currents, debris flows). To provide understanding of this complexity, we have employed source to sink numerical models to evaluate which process dominates the observed variability in a sedimentary record of two coastal Pacific basins, Knight Inlet in British Columbia and the Eel Margin of northern California.During the last glacial period, the Eel River supplied comparatively more sediment with a less variable flux to the ocean, while today the river is dominated by episodic events. Model results show this change in the variability of sediment flux to be as important to the deposit character as is the change in the volume of sediment supply. Due to the complex interaction of flooding events and ocean storm events, the more episodic flood deposits of recent times are less well preserved than the flood deposits associated with an ice-age climate.In Knight Inlet, the evolving boundary conditions (rapidly prograding coastline, secondary transport by gravity flows from sediment failures) are a strong influence on the sedimentary record. The delta and gravity flow deposits punctuate the sedimentary record formed by hemipelagic sedimentation from river plumes. Missing time intervals due to sediment failures can take away the advantage of the otherwise amplified lithologic record of discharge events, given the enclosed nature of the fjord basin.     

Dependence on agriculture and ecosystem services for livelihood in Northeast India and Bhutan: vulnerability to climate change in the Tropical River Basins of the Upper Brahmaputra

The Upper Brahmaputra River Basin is prone to natural disasters and environmental stresses (floods, droughts and bank erosion, delayed rainfall, among others) creating an environment of uncertainty and setting the basin back in terms of socio-economic development. The climate change literature shows that agriculture and ecosystems and their services are highly climate sensitive, yet they are the main sources of livelihood that supports a large proportion of residents of the tributaries of the Brahmaputra River Basin. The continuous depletion of ecosystems and loss of agricultural outputs resulting from environmental stressors has a substantial impact on the socio-economic wellbeing of the basins residents, particularly the vulnerable rural poor. This paper uses spatially explicit data from Census, Household Surveys and Earth Observation to develop a transferable methodological approach which investigates the extent of dependence on agriculture and ecosystems as a source of livelihood in the contrasting sub-basins of the Brahmaputra River in the State of Assam, India and Bhutan, and the risk to these livelihood dependencies in these sub-basins due to potential environmental impacts of climate change. The results from this study constitute a case study in the development of a systematic and spatially explicit set of tools that inform and assist policy makers in the appropriate interventions to secure the livelihood benefits of sustainably managed agriculture in the face of environmental change.     

Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.     

Development of daily spatial heat unit mapping from monthly climatic surfaces for the Australian continent

In absence of other limitations, the growth rate of a plant is dependent upon the amount of heat it receives. Each species, whether a crop, weed or disease organism, is adapted to grow at its optimum rate within a specific temperature range. Within this range, the growing degree days (GDD) is the heat accumulation above a given base temperature for a specific time period, such as a crop's growing season or phenological stage. In this paper we detail a methodology to predict GDD for synthetically generated average growing seasons derived from long term average climate data over the Australian continent. An application of these techniques has been made using the GEODATA 9 second DEM, with temperature threshold values estimated to characterize optimum growth in citrus (Citrus sinensis (L.) Osbeck). Three major determinants of the annual growth cycle of Citrus sp. were established and predicted on a spatial basis including the starting day of the growing season, the GDD for a growing season, and the time required to accumulate an arbitrarily selected 2000 GDD from the estimated starting day. When these critical environmental factors are expressed on a spatial basis, covering the Australian continent, the combination can be used to identify locations where new crop varieties can most effectively be grown to maximize fruit quality and productivity, or to extend the harvest season. Likewise, new germplasm introduced to Australia from overseas can be horticulturally assessed at sites climatically matched to the source location.     

Evaluation of infiltration‐based stormwater management to restore hydrological processes in urban headwater streams

Urbanization threatens headwater stream ecosystems globally. Watershed restoration practices, such as infiltration‐based stormwater management, are implemented to mitigate the detrimental effects of urbanization on aquatic ecosystems. However, their effectiveness for restoring hydrologic processes and watershed storage remains poorly understood. Our study used a comparative hydrology approach to quantify the effects of urban watershed restoration on watershed hydrologic function in headwater streams within the Coastal Plain of Maryland, USA. We selected 11 headwater streams that spanned an urbanization–restoration gradient (4 forested, 4 urban‐degraded, and 3 urban‐degraded) to evaluate changes in watershed hydrologic function from both urbanization and watershed restoration. Discrete discharge and continuous, high‐frequency rainfall‐stage monitoring were conducted in each watershed. These datasets were used to develop 6 hydrologic metrics describing changes in watershed storage, flowpath connectivity, or the resultant stream flow regime. The hydrological effects of urbanization were clearly observed in all metrics, but only 1 of the 3 restored watersheds exhibited partially restored hydrologic function. At this site, a larger minimum runoff threshold was observed relative to the urban‐degraded watersheds, suggesting enhanced infiltration of stormwater runoff within the restoration structure. However, baseflow in the stream draining this watershed remained low compared to the forested reference streams, suggesting that enhanced infiltration of stormwater runoff did not recharge subsurface storage zones contributing to stream baseflow. The highly variable responses among the 3 restored watersheds were likely due to the spatial heterogeneity of urban development, including the level of impervious cover and extent of the storm sewer network. This study yielded important knowledge on how restoration strategies, such as infiltration‐based stormwater management, modulated—or failed to modulate—hydrological processes affected by urbanization, which will help improve the design of future urban watershed management strategies. More broadly, we highlighted a multimetric approach that can be used to monitor the restoration of headwater stream ecosystems in disturbed landscapes.     

Sources of sodium in the lunar exosphere: Modeling using ground-based observations of sodium emission and spacecraft data of the plasma

Observations of the equatorial lunar sodium emission are examined to quantify the effect of precipitating ions on source rates for the Moon’s exospheric volatile species. Using a model of exospheric sodium transport under lunar gravity forces, the measured emission intensity is normalized to a constant lunar phase angle to minimize the effect of different viewing geometries. Daily averages of the solar Lyman α flux and ion flux are used as the input variables for photon-stimulated desorption (PSD) and ion sputtering, respectively, while impact vaporization due to the micrometeoritic influx is assumed constant. Additionally, a proxy term proportional to both the Lyman α and to the ion flux is introduced to assess the importance of ion-enhanced diffusion and/or chemical sputtering. The combination of particle transport and constrained regression models demonstrates that, assuming sputtering yields that are typical of protons incident on lunar soils, the primary effect of ion impact on the surface of the Moon is not direct sputtering but rather an enhancement of the PSD efficiency. It is inferred that the ion-induced effects must double the PSD efficiency for flux typical of the solar wind at 1 AU. The enhancement in relative efficiency of PSD due to the bombardment of the lunar surface by the plasma sheet ions during passages through the Earth’s magnetotail is shown to be approximately two times higher than when it is due to solar wind ions. This leads to the conclusion that the priming of the surface is more efficiently carried out by the energetic plasma sheet ions.