Use of incipient motion data for backward erosion piping models

Backward erosion piping involves the gradual removal of granular material under the action of water flow from the foundation of a dam or levee, whereby shallow pipes are formed that grow in the direction opposite to the flow. This pipe-forming process can ultimately lead to failure of a water-retaining structure and is considered one of the most important failure mechanisms for dikes and levees in the Netherlands and the United States. Modeling of this mechanism requires the assessment of hydraulic conditions in the pipe, which are controlled by the particle equilibrium at the pipe wall. Since the pipe's dimensions are controlled by the inflow to the pipe from the porous medium, the flow through the pipe is thought to be laminar for fine- to medium-grained sands. The literature provides data for incipient motion in laminar flow, which is reviewed here and complemented with data from backward erosion experiments. The experiments illustrate the applicability of the laminar incipient motion data to determine the erosion pipe dimensions and corresponding pipe hydraulics for fine- to medium-grained sands, for the purpose of backward erosion piping modeling.     

Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW”

The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.     

Classifying zones of suitability for manual drilling using textural and hydraulic parameters of shallow aquifers: a case study in northwestern Senegal

A method is proposed that uses analysis of borehole stratigraphic logs for the characterization of shallow aquifers and for the assessment of areas suitable for manual drilling. The model is based on available borehole-log parameters: depth to hard rock, depth to water, thickness of laterite and hydraulic transmissivity of the shallow aquifer. The model is applied to a study area in northwestern Senegal. A dataset of boreholes logs has been processed using a software package (TANGAFRIC) developed during the research. After a manual procedure to assign a standard category describing the lithological characteristics, the next step is the automated extraction of different textural parameters and the estimation of hydraulic conductivity using reference values available in the literature. The hydraulic conductivity values estimated from stratigraphic data have been partially validated, by comparing them with measured values from a series of pumping tests carried out in large-diameter wells. The results show that this method is able to produce a reliable interpretation of the shallow hydrogeological context using information generally available in the region. The research contributes to improving the identification of areas where conditions are suitable for manual drilling. This is achieved by applying the described method, based on a structured and semi-quantitative approach, to classify the zones of suitability for given manual drilling techniques using data available in most African countries. Ultimately, this work will support proposed international programs aimed at promoting low-cost water supply in Africa and enhancing access to safe drinking water for the population.     

Concentration,composition and sources of PAHs in the coastal sediments of the exclusive economic zone (EEZ) of Qatar,Arabian Gulf

Surface sediments were collected from sixteen locations in order to assess levels and sources of polycyclic aromatic hydrocarbons (PAHs) in sediments of Qatar exclusive economic zone (EEZ). Samples were analyzed for 16 parent PAHs, 18 alkyl homologs and for dibenzothiophenes. Total PAHs concentration (∑PAHs) ranged from 2.6 ng g⁻¹ to 1025 ng g⁻¹. The highest PAHs concentrations were in sediments in and adjacent to harbors. Alkylated PAHs predominated most of the sampling locations reaching up to 80% in offshore locations. Parent PAHs and parent high molecular weight PAHs dominated location adjacent to industrial activities and urban areas. The origin of PAHs sources to the sediments was elucidated using ternary plot, indices, and molecular ratios of specific compounds such as (Ant/Phe + Ant), (Flt/Flt + Pyr). PAHs inputs to most coastal sites consisted of mixture of petroleum and combustion derived sources. However, inputs to the offshore sediments were mainly of petroleum origin.     

Long-term diagnostics of precipitation estimates and the development of radar hardware monitoring within a radar product data quality management system

Abstract

Quality is key to ensuring that the potential offered by weather radar networks is realized. To achieve optimum quality, a comprehensive radar data quality management system, designed to monitor the end-to-end radar data processing chain and evaluate product quality, is being developed at the UK Met Office. Three contrasting elements of this system are described: monitoring of key radar hardware performance indicators; generation of long-term integrations of radar products; and monitoring of radar reflectivity factor using synthesized observations from numerical weather prediction model fields. Examples of each component are presented and ways in which the different types of monitoring information have been used to both identify issues with the radar product data quality and help formulate solutions are given.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Harrison, D., Georgiou, S., Gaussiat, N., and Curtis, A., 2013. Long-term diagnostics of precipitation estimates and the development of radar hardware monitoring within a radar product data quality management system. Hydrological Sciences Journal, 59 (7), 1327–1342. http://dx.doi.org/10.1080/02626667.2013.841316     

Evaluating the sensitivity of wetlands to climate change with remote sensing techniques

Wetlands are valuable ecosystems that provide many valuable services, yet many of these important ecosystems are at risk because of current trends in climate change. The Prairie Pothole Region (PPR) in the upper‐midwest of the United States and south‐central Canada, characterized by glacially sculpted landscapes and abundant wetlands, is one such vulnerable region. According to regional/global climate model predictions, drought occurrence will increase in the PPR region through the 21st century and thus will probably cause the amount of water in wetlands to decline. Water surface area (WSA) of Kidder County, ND, from 1984–2011 was measured by classifying TM/ETM+(Landsat Thematic Mapper / Enhanced Thematic Mapper Plus) images through the modified normalized difference water index. We then developed a linear model based on the WSA of these wetlands and historical climate data and used this to determine the wetland sensitivity to climate change and predict future wetlands WSA in the PPR. Our model based on Palmer drought severity index (PDSI) of the current year (PDSI_{t ? 0}) and of the previous two years (PDSI_{t ? 2}) can explain 79% of the annual wetland WSA variance, suggesting a high sensitivity of wetlands to drought/climate change. We also predicted the PPR wetlands WSA in the 21st century under A1B scenario (a mid‐carbon emission scenario) using simulated PDSI based on Intergovernmental Panel on Climate Change AR4 22‐model ensemble climate. According to our prediction, the WSA of the PPR wetlands will decrease to less than half of the baseline WSA (defined as the mean wetlands WSA of the 2000s) by the mid of the 21st century, and to less than one‐third by the 2080s, and will then slightly increase in the 2090s. This considerable future wetland loss caused only by climate change provides important implication to future wetland management and climate adaptation policy. Copyright © 2012 John Wiley & Sons, Ltd.     

Geologic and hydrologic settings for development of freshwater lenses in arid lands

Satellite observations were used to test the validity of previously identified favourable conditions for the formation of freshwater lenses, identify additional potential occurrences, and model modern potential recharge in the Raudhatain Watershed (3696) in northern Kuwait. Favourable conditions include infrequent yet intensive precipitation events, drainage depressions to collect the limited runoff, and presence of conditions (e.g. high infiltration capacity) that promote groundwater recharge and preservation (e.g. underlying saline aquifer) of infiltrating groundwater as freshwater lenses floating over saline aquifer water due to differences in density. Specifically, the following field and satellite‐based observations were noted for the Raudhatain Watershed: (1) Over ~30 precipitation events were identified from the Tropical Rainfall Measuring Mission precipitation data (1998–2009); (2) slope is gentle (2 m/km), and the surface is largely (80%) covered by alluvial deposits with high infiltration capacities (up to 9 m/day); (3) no flows and long‐term ponding were reported at the watershed outlet or detected from Landsat thematic mapper images; (4) infiltration is high based on increases in soil moisture content (from an advanced microwave scanning radiometer) and vegetation index following large precipitation events; and (5) freshwater lenses that overlie highly saline [total dissolved solids (TDS): >35 000] unconfined aquifers underlying the watershed are absent in the southern regions, where infiltrating fresh water mixes with the less saline groundwater (TDS: <10 000). Twenty potential locations (size: 1 to 75 km²) for freshwater lens development were identified in northern Kuwait, and continuous rainfall–runoff models (Soil Water and Assessment Tool) were applied to provide a first‐order estimation of the average annual recharge in the watershed (127 × 10⁶ m³) and freshwater lenses (8.17 × 10⁶ m³). Results demonstrate the settings for enhanced opportunities for groundwater recharge, outline the amounts of and preservation conditions for the groundwater feeding the freshwater lenses, and highlight potential applications and locations of freshwater lenses in similar settings elsewhere in the Arabian Peninsula and beyond. Copyright © 2013 John Wiley & Sons, Ltd.     

El Chichón volcano, April 4, 1982: volcanic cloud history and fine ash fallout

This retrospective study focuses on the fine silicate particles (<62 µm in diameter) produced in a large eruption that was otherwise well studied. Fine particles represent a potential hazard to aircraft, because as simple particles they have very low terminal velocities and could potentially stay aloft for weeks. New data were collected to describe the fine particle size distributions of distal fallout samples collected soon after eruption. Although, about half of the mass of silicate particles produced in this eruption of ~1 km³ dense rock equivalent magma were finer than 62 µm in diameter, and although these particles were in a stratospheric cloud after eruption, almost all of these fine particles fell to the ground near (<300 km) the volcano in a day or two. Particles falling out from 70 to 300 km from the volcano are mostly <62 µm in diameter. The most plausible explanation for rapid fallout is that the fine ash nucleates ice in the convective cloud and initiates a process of meteorological precipitation that efficiently removes fine silicates. These observations are similar to other eruptions and we conclude that ice formation in convective volcanic clouds is part of an effective fine ash removal process that affects all or most volcanic clouds. The existence of pyroclastic flows and surges in the El Chichón eruption increased the overall proportion of fine silicates, probably by milling larger glassy pyroclasts.     

Atmospheric circulations of terrestrial planets orbiting low-mass stars

Circulations and habitable zones of planets orbiting low-mass stars are investigated. Many of these planets are expected to rotate synchronously relative to their parent stars, thereby raising questions about their surface temperature distributions and habitability. We use a global circulation model to study idealized, synchronously rotating (tidally locked) planets of various rotation periods, with surfaces of all land or all water, but with an Earth-like atmosphere and solar insolation. The dry planets exhibit wide variations in surface temperature: >80 °C on the dayside to <−110 °C on the nightside for the 240-h rotator, for example. The water-covered aquaplanets are warmer and exhibit narrower ranges of surface temperatures, e.g., ∼40 °C to >−60 °C for the 240-h orbiter. They also have a larger habitable area, defined here as the region where average surface temperatures are between 0 °C and 50 °C. This concept has little relevance for either dry or aquaplanets, but might become relevant on a planet with both land area and oceans.The circulations on these tidally locked planets exhibit systematic changes as the rotation period is varied. However, they also reveal abrupt transitions between two different circulation regimes and multiple equilibria. For the dry planet, the transition occurs between a 4-day and a 5-day period, while for the aquaplanet, it occurs between a 3-day and a 4-day period. For both dry and aqua planets, this transition occurs when the Rossby deformation radius exceeds half the planetary radius. Further investigation on the dry planet reveals that multiple equilibria exist between 100- and 221-h periods. These multiple equilibria may be relevant for real planets within the habitable zones of late K and M stars, because these planets are expected to have rotation periods between 8 and 100 Earth days.     

Coupled convection and tidal dissipation in Europa’s ice shell using non-Newtonian grain-size-sensitive (GSS) creep rheology

We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10⁻⁵), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10⁻⁵ on Europa.