Convection with heat flux prescribed on the boundaries of the system. I. The effect of temperature dependence of material properties

Abstract

We study the bifurcation to steady two-dimensional convection with the heat flux prescribed on the fluid boundaries. The fluid is weakly non-Boussinesq on account of a slight temperature dependence of its material properties. Using expansions in the spirit of shallow water theory based on the preference for large horizontal scales in fixed flux convection, we derive an evolution equation for the horizontal structure of convective cells. In the steady state, this reduces to a simple nonlinear ordinary differential equation. When the horizontal scales of the cells exceed a certain critical size, the bifurcation to steady convection is subcritical and the degree of subcriticality increases with increasing cell size.     

Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer,Cyprus

Abstract

We investigate the general methodology for an intensive development of coastal aquifers, described in a companion paper, through its application to the management of the Akrotiri aquifer, Cyprus. The Zakaki area of that aquifer, adjacent to Lemessos City, is managed such that it permits a fixed annual agricultural water demand to be met, as well as and a fraction of the water demand of Lemessos, which varies according to available surface water. Effluents of the Lemessos wastewater treatment plant are injected into the aquifer to counteract the seawater intrusion resulting from the increased pumping. The locations of pumping and injection wells are optimized based on least-cost, subject to meeting the demand. This strategy controls sea intrusion so effectively that desalting of only small volumes of slightly brackish groundwater is required over short times, while ～2.3 m³ of groundwater is produced for each 1 m³ of injected treated wastewater. The cost over the 20-year period 2000–2020 of operation is ～40 M€ and the unit production cost of potable water is under 0.2 €/m³. The comparison between the deterministic and stochastic analyses of the groundwater dynamics indicates the former as conservative, i.e. yielding higher groundwater salinity at the well. The Akrotiri case study shows that the proposed aquifer management scheme yields solutions that are preferable to the widely promoted seawater desalination, also considering the revenues from using the treated wastewater for irrigation.

Citation Koussis, A. D., Georgopoulou, E., Kotronarou, A., Mazi, K., Restrepo, P., Destouni, G., Prieto, C., Rodriguez, J. J., Rodriguez-Mirasol, J., Cordero, T., Ioannou, C., Georgiou, A., Schwartz, J. & Zacharias, I. (2010) Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer, Cyprus. Hydrol. Sci. J. 55(7), 1234–1245.     

Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

THE INFLUENCE OF SLOPE AND RAIN INTENSITY ON RUNOFF AND INFILTRATION / L'influence de l'inclinaison de terrain et de l'intensité de pluie sur l'écoulement et l'infiltration

Abstract

The design and construction of a special-purpose laboratory catchment and rainfall simulator is described. The equipment consists of a soil catchment area that can be inclined at various angles. Additional instrumentation then measures the flow of water across the surface of, and through, the soil bed. Precipitation is provided by a unit that simulates rainfall at particular rates with uniform distribution.

The equipment was used to examine infiltration, runoff and other hydrological properties of a number of soils under different rainfall intensities and with different catchment slopes. Correlations were obtained for these variables.     

Grid Refinement in Cartesian Coordinates for Groundwater Flow Models Using the Divergence Theorem and Taylor's Series

Grid refinement is introduced in a numerical groundwater model to increase the accuracy of the solution over local areas without compromising the run time of the model. Numerical methods developed for grid refinement suffered certain drawbacks, for example, deficiencies in the implemented interpolation technique; the non‐reciprocity in head calculations or flow calculations; lack of accuracy resulting from high truncation errors, and numerical problems resulting from the construction of elongated meshes. A refinement scheme based on the divergence theorem and Taylor's expansions is presented in this article. This scheme is based on the work of De Marsily (1986) but includes more terms of the Taylor's series to improve the numerical solution. In this scheme, flow reciprocity is maintained and high order of refinement was achievable. The new numerical method is applied to simulate groundwater flows in homogeneous and heterogeneous confined aquifers. It produced results with acceptable degrees of accuracy. This method shows the potential for its application to solving groundwater heads over nested meshes with irregular shapes.     

Evaluation of the response modification coefficient and collapse potential of special concentrically braced frames

Special concentrically braced frames (SCBFs) are commonly used for seismic design of buildings. Their large elastic stiffness and strength efficiently sustains the seismic demands during smaller, more frequent earthquakes. During large, infrequent earthquakes, SCBFs exhibit highly nonlinear behavior due to brace buckling and yielding and the inelastic behavior induced by secondary deformation of the framing system. These response modes reduce the system demands relative to an elastic system without supplemental damping using a response modification coefficient, commonly termed the R factor. More recently, procedures put forth in FEMAP695 have been made to quantify the R factor through a formalized procedure that accounts for collapse potential. The primary objective of the research in this paper was to evaluate the approach for SCBFs. An improved model for SCBFs that permits simulation of brace fracture was used to conduct response history analyses. A series of three‐story, nine‐story and 20‐story SCBFs were designed and evaluated. Initially, the FEMAP695 method was conducted to estimate collapse and the corresponding R factor. An alternate procedure for scaling the multiple acceleration records to the seismic design hazard was also evaluated. The results show significant variation between the two methods. Of the three variations of buildings studied, the largest vulnerability was identified for the three‐story building. To achieve a consistent margin of safety against collapse, a significantly lower R factor is required for the low‐rise SCBFs (three‐story), whereas the mid‐rise and high‐rise SCBFs (nine‐story and 20‐story) may continue to use the current value of 6, as provided in ASCE‐07. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatial interaction between collapsed pipes and landslides in hilly regions with loess‐derived soils

While most studies focus on the effect of soil pipes on hillslope stability, this present study investigates the impact of landsliding on pipe development. It is hypothesized that poorly drained active and dormant landslides change the hillslope hydrology through (i) surface flow obstruction, by changing topography, as well as (ii) subsurface flow obstruction by tilting less‐permeable clay‐rich substrates. Hence, new preferential flow paths are created at reverse slopes within the landslide zone and at the boundary of the landslide, enhancing pipe formation. This study aims at a better understanding of the interaction between collapsed pipe (CP) occurrence and landslide (LS) occurrence in the Flemish Ardennes (Belgium) by comparing their respective spatial patterns. At least 24.5% of the 139 sites with CP were related to the occurrence of an observed LS. Poorly drained LS may create favourable conditions for pipe development. Outside LS, natural and anthropogenic (e.g. broken field drains, road drainage) causes may result in concentrated subsurface flow, resulting in pipe development. No evidence was found that pipe development enhanced LS, probably because the subsurface drainage discharge generated upslope of the LS is too low. Even when pipes become blocked, it is more likely that new pipes develop and new collapses occur than they trigger or reactivate LS. A conceptual model is presented summarizing all elements that influence piping erosion in the Flemish Ardennes, including the role of LS. Copyright © 2012 John Wiley & Sons, Ltd.     

Gravel‐bed rivers: from particles to patterns

The papers in this special issue reflect several of the major themes and topics from the 7^th International Workshop on Gravel‐Bed Rivers. The papers focus primarily on aspects of bed material transport in gravel‐bed rivers and larger scale morpho‐dynamics. Research in gravel‐bed rivers is increasingly integrating processes over a wide range of temporal and spatial scales by combining field observation, lab experimentation, numerical modeling and theory testing in a range of river types, aided by new technological developments in particle tracking, computational modeling and high resolution spatial data. This is leading to greater understanding of the processes leading to distinctive morpho‐dynamics of river types and a more reliable basis for river management. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating water balance components of tropical wetland lakes in the Pantanal dry season,Brazil

Abstract

Water balance studies with stable water isotopes have rarely been conducted in remote and tropical wetland areas. As such, little is known regarding the water balance and groundwater–surface water interaction in the Pantanal, one of the largest and most pristine wetlands in the world. We applied MINA TrêS, a water balance model utilizing stable water isotopes (δ¹⁸O, δ²H) and chloride (Cl^-) to assess the dry-season hydrological processes controlling groundwater–surface water interactions and the water balance of six floodplain lakes in the northern Pantanal, Brazil. Qualitatively, all lakes exhibited similarity in hydrological controls. Quantitatively, they differed significantly due to morphological differences in controlling groundwater inflow and lake volume. Our approach is readily transferable to other remote and tropical wetland systems with minimal data input requirements, which is useful in regions with sparse hydrometric monitoring.

Editor Z.W. Kundzewicz     

Spectral method for estimating the quality of the medium from one source and two stations

Izvestiya, Physics of the Solid Earth - The method is suggested for estimating the quality (coefficient Q) of the medium from the longitudinal waves of the earthquakes in the regional zone (Q P )....