Particle manifold method (PMM): A new continuum‐discontinuum numerical model for geomechanics

Particle manifold method (PMM) is a new extension of the numerical manifold method (NMM). PMM uses a mathematical cover system to describe the motion and deformation of a particle‐based physical domain. By introducing the concept of particle into NMM, PMM takes the advantages of easy topological and contact operations with particles. In this article, the methodology, formulations and implementation of the method are presented, together with modelling examples for validation. It is found that good solutions for both continuous and discontinuous problems are obtained by the new developed PMM. Due to the underlying coupled continuum‐discontinuum property of PMM, it has great potential for modelling of geomechanical problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessing the potential of reservoir outflow management to reduce sedimentation using continuous turbidity monitoring and reservoir modelling

In‐stream sensors are increasingly deployed as part of ambient water quality‐monitoring networks. Temporally dense data from these networks can be used to better understand the transport of constituents through streams, lakes or reservoirs. Data from existing, continuously recording in‐stream flow and water quality monitoring stations were coupled with the two‐dimensional hydrodynamic CE‐QUAL‐W2 model to assess the potential of altered reservoir outflow management to reduce sediment trapping in John Redmond Reservoir, located in east‐central Kansas. Monitoring stations upstream and downstream from the reservoir were used to estimate 5.6 million metric tons of sediment transported to John Redmond Reservoir from 2007 through 2010, 88% of which was trapped within the reservoir. The two‐dimensional model was used to estimate the residence time of 55 equal‐volume releases from the reservoir; sediment trapping for these releases varied from 48% to 97%. Smaller trapping efficiencies were observed when the reservoir was maintained near the normal operating capacity (relative to higher flood pool levels) and when average residence times were relatively short. An idealized, alternative outflow management scenario was constructed, which minimized reservoir elevations and the length of time water was in the reservoir, while continuing to meet downstream flood control end points identified in the reservoir water control manual. The alternative scenario is projected to reduce sediment trapping in the reservoir by approximately 3%, preventing approximately 45 000 metric tons of sediment from being deposited within the reservoir annually. This article presents an approach to quantify the potential of reservoir management using existing in‐stream data; actual management decisions need to consider the effects on other reservoir benefits, such as downstream flood control and aquatic life. Copyright © 2012 John Wiley & Sons, Ltd.     

Geodynamic modelling of the Century deposit,Mt Isa Province,Queensland

This paper is concerned with an example of quantitative modelling of orebody formation as a guide to reducing the risk for future mineral exploration. Specifically, the paper presents a detailed 3–D numerical model for the formation of the Century zinc deposit in northern Queensland. The model couples fluid flow with deformation, thermal transport and chemical reactions. The emphasis of the study is a systems approach where the holistic mineralising system is considered rather than concentrating solely on the mineral deposit. In so doing the complete plumbing system for mineralisation is considered with a view to specifying the critical conditions responsible for the ore deposit occurring where it does and having the size and metal grades that are observed. The numerical model is based on detailed geological, tectonic, isotopic and mineralogical data collected over the past 20 years. The conclusions are that the Century zinc deposit is located where it is because of the following factors: (i) a thermal anomaly is associated with the Termite Range Fault due to advection of heat from depth by fluid flow up the Termite Range Fault; (ii) bedding‐plane fissility in the shale rocks hosting the Century zinc deposit has controlled the wavelength and nature of D₁ folding in the vicinity of the deposit and has also controlled increases in permeability due to hydrofracture of the shales; such hydrofracture is also associated with the production of hydrocarbons as these shales passed through the ‘oil‐window’; (iii) Pb–Zn leached from crustal rocks in the stratigraphic column migrated up along faults normal to the Termite Range Fault driven by topographic relief associated with inversion at the end of the Isan Orogeny; these fluids mixed with H₂S derived at depth moving up the Termite Range Fault to mix with the crustal fluids to precipitate Pb–Zn in a plume downstream from the point of mixing. Critical factors to be used as exploration guides are high temperatures, carbonaceous fissile shales now folded into relatively tight D₁ folds, fault‐controlled plumbing systems that enable fluid mixing, depletion of metals upstream of the deposit and,in particular,a very wide Fe‐depletion halo upstream of the deposit.     

Snow cover and snowmelt of an extensive High Arctic wetland: spatial and temporal seasonal patterns

Abstract

This study examined the end-of-winter snow storage, its distribution and the spatial and temporal melt patterns of a large, low gradient wetland at Polar Bear Pass, Bathurst Island, Nunavut, Canada. The project utilized a combination of field observations and a physically-based snowmelt model. Topography and wind were the major controls on snow distribution in the region, and snow was routinely scoured from the hilltop regions and deposited into hillslopes and valleys. Timing and duration of snowmelt at Polar Bear Pass were similar in 2008 and 2009. The snowmelt was initiated by an increase in air temperature and net radiation receipt. Inter-annual variability in spatial snowmelt patterns was evident at Polar Bear Pass and was attributed to a non-uniform snow cover distribution and local microclimate conditions. In situ field studies and modelling remain important in High Arctic regions for assessing wetland water budgets and runoff, in addition to model parameterization and validation of satellite imagery.

Editor Z.W. Kundzewicz

Citation Assini, J. and Young, K.L., 2012. Snow cover and snowmelt of an extensive High Arctic wetland: spatial and temporal seasonal patterns. Hydrological Sciences Journal, 57 (4), 738–755.     

A dual-porosity model for nitrogen leaching from a watershed/Un modèle à double porosité pour le lessivage de l’azote d’un bassin versant

Abstract

Abstract The MASONW (MACRO + SOILN + Watershed) model describing nitrogen leaching in watersheds was developed and tested. The model is based on the MACRO and SOILN models. The dual-porosity model MACRO simulates water flow on the field scale. The SOILN model describes turnover and leaching of nitrogen. Two main features of a watershed have been added into these two models: (a) the existence of a river system, and (b) variable thickness of the aeration zone within a watershed. Good agreement between the output of the MASONW model and observed data for water discharge and nitrate concentrations were achieved in the Odense watershed (496 km²) in Denmark.     

Long‐term final void salinity prediction for a post‐mining landscape in the Hunter Valley,New South Wales,Australia

Opencast mining alters surface and subsurface hydrology of a landscape both during and post‐mining. At mine closure, following opencast mining in mines with low overburden to coal ratios, a void is left in the final landform. This final void is the location of the active mine pit at closure. Voids are generally not infilled within the mines' lifetime, because of the prohibitive cost of earthwork operations, and they become post‐mining water bodies or pit lakes. Water quality is a significant issue for pit lakes. Groundwater within coal seams and associated rocks can be saline, depending on the nature of the strata and groundwater circulation patterns. This groundwater may be preferentially drawn to and collected in the final void. Surface runoff to the void will not only collect salts from rainfall and atmospheric fallout, but also from the ground surface and the weathering of fresh rock. As the void water level rises, its evaporative surface area increases, concentrating salts that are held in solution. This paper presents a study of the long term, water quality trends in a post‐mining final void in the Hunter Valley, New South Wales, Australia. This process is complex and occurs long term, and modelling offers the only method of evaluating water quality. Using available geochemical, climate and hydrogeological data as inputs into a mass‐balance model, water quality in the final void was found to increase rapidly in salinity through time (2452 to 8909 mg l⁻¹ over 500 years) as evaporation concentrates the salt in the void and regional groundwater containing high loads of salt continues to flow into the void. Copyright © 2004 John Wiley & Sons, Ltd.     

Snow interception modelling: Isolated observations have led to many land surface models lacking appropriate temperature sensitivities

When formulating a hydrologic model, scientists rely on parameterizations of multiple processes based on field data, but literature review suggests that more frequently people select parameterizations that were included in pre-existing models rather than re-evaluating the underlying field experiments. Problems arise when limited field data exist, when “trusted” approaches do not get reevaluated, and when sensitivities fundamentally change in different environments. The physics and dynamics of snow interception by conifers is just such a case, and it is critical to simulation of the water budget and surface albedo. The most commonly used interception parameterization is based on data from four trees from one site, but results from this field study are not directly transferable to locations with relatively warmer winters, where the dominant processes differ dramatically. Here, we combine a literature review with model experiments to demonstrate needed improvements. Our results show that the choice of model form and parameters can vary the fraction of snow lost through interception by as much as 30%. In most simulations, the warming of mean winter temperatures from −7 to 0°C reduces the modelled fraction of snow under the canopy compared to the open, but the magnitude of simulated decrease varies from about 10% to 40%. The range of results is even larger when considering models that neglect the melting of in-canopy snow in higher-humidity environments where canopy sublimation plays less of a role. Thus, we recommend that all models represent canopy snowmelt and include representation of increased loading due to increased adhesion and cohesion when temperatures rise from −3 to 0°C. In addition to model improvements, field experiments across climates and forest types are needed to investigate how to best model the combination of dynamically changing forest cover and snow cover to better understand and predict changes to albedo and water supplies.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Stream nitrate uptake and transient storage over a gradient of geomorphic complexity,north‐central Colorado,USA

The understanding of nutrient uptake in streams is impeded by a limited understanding of how geomorphic setting and flow regime interact with biogeochemical processing. This study investigated these interactions as they relate to transient storage and nitrate uptake in small agricultural and urban streams. Sites were selected across a gradient of channel conditions and management modifications and included three 180‐m long geomorphically distinct reaches on each of two streams in north‐central Colorado. The agricultural stream has been subject to historically variable cattle‐grazing practices, and the urban stream exhibits various levels of stabilisation and planform alteration. Reach‐scale geomorphic complexity was characterised using highly detailed surveys of channel morphology, substrate, hydraulics and habitat units. Breakthrough‐curve modelling of conservative bromide (Br^?) and nonconservative nitrate (NO₃^?) tracer injections characterised transient storage and nitrate uptake along each reach. Longitudinal roughness and flow depth were positively associated with transient storage, which was related to nitrate uptake, thus underscoring the importance of geomorphic influences on stream biogeochemical processes. In addition, changes in geomorphic characteristics due to temporal discharge variation led to complex responses in nitrate uptake. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological response of a Brazilian semi‐arid catchment to different land use and climate change scenarios: a modelling study

Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R² of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.