Numerical model of convection in sills

The numerical model of convection in magma sills is developed. The model is based on a full system of equations of fluid dynamics and includes heat transfer, buoyancy effects and diffusion of some minor component (marker). Solidification is treated as a phase transition. The results indicate that there are some qualitative differences between very thin sills with Rayleigh number Ra = 10⁵ and thin sills with Ra = 10⁶. For a basaltic magma the first case corresponds to the thickness of the sills of approximately 30 cm and the second case corresponds to the thickness of 60 cm. In the first case mixing is inefficient and conduction is the dominant form of heat transfer. In the second case mixing is efficient and convection is the dominant form of heat transfer. Some of the results can be scaled for the more viscous magmas in thicker sills.     

Salinisation des nappes côtières : cas de la nappe nord du Sahel de Sfax,Tunisie

The intensive agricultural and economic activities induce the increase of the risk of groundwater quality degradation through high groundwater pumping rates. The salinization and contamination are the main sources of this pollution, especially in coastal aquifers. The explanation of the origin of salinity for the shallow aquifer of Northern Sahel of Sfax was analysed by a chemical study of the groundwater main compounds. The partitioning of groundwaters into homogenous groups is undertaken by graphical techniques, including a Stiff pattern diagram, an expanded Durov diagram and several binary diagrams. The study indicates the presence of various salinization processes. In the recharge area, salinization is the result of dissolution/precipitation of the aquifer formation material (group I). The irrigation water return and the intensive pumping have been identified as major sources of salinization in the south by direct cation exchange and mixing reactions (groups II and III). The anomaly of high groundwater salinity observed near the Hazeg zone was explained by the presence of a seawater intrusion in this area. This hypothesis is related to the high chloride concentration, to the presence of inverse cation exchange reactions (group IV), and to the piezometric level inferior to sea level. To cite this article: R. Trabelsi et al., C. R. Geoscience 337 (2005).     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of flow paths and confluences for saltwater intrusion in a deltaic river network

Saltwater intrusion is a serious issue in estuarine deltas all over the world due to rapid urban sprawl and water shortage. Therefore, detecting the major flow paths or locations at risk of saltwater intrusion in estuarine ecosystems is important for mitigating saltwater intrusion. In this paper, we introduce a centrality index, the betweenness centrality (BC), to address this problem. Using the BC as the weighted attribute of the river network, we identify the critical confluences for saltwater intrusion and detect the preferential flow paths for saltwater intrusion through the least‐cost‐path algorithm from a graph theory approach. Moreover, we analyse the responses of the BC values of confluences calculated in the river network to salinity. Our results show that the major flow paths and critical confluences for saltwater intrusion in a deltaic river network can be represented by the least cost paths and the BC values of confluences, respectively. In addition, a significant positive correlation between the BC values of confluences and salinity is determined in the Pearl River Delta. Changes in the salinity can produce significant variation in the BC values of confluences. Therefore, freshwater can be diverted into these major flow paths and critical confluences to improve river network management under saltwater intrusion. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of daily stream water temperatures with a Bayesian regression approach

Stream water temperature plays a significant role in aquatic ecosystems where it controls many important biological and physical processes. Reliable estimates of water temperature at the daily time step are critical in managing water resources. We developed a parsimonious piecewise Bayesian model for estimating daily stream water temperatures that account for temporal autocorrelation and both linear and nonlinear relationships with air temperature and discharge. The model was tested at 8 climatically different basins of the USA and at 34 sites within the mountainous Boise River Basin (Idaho, USA). The results show that the proposed model is robust with an average root mean square error of 1.25 °C and Nash–Sutcliffe coefficient of 0.92 over a 2‐year period. Our approach can be used to predict historic daily stream water temperatures in any location using observed daily stream temperature and regional air temperature data.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

Prediction of river water temperature: a comparison between a new family of hybrid models and statistical approaches

River water temperature is a key physical variable controlling several chemical, biological and ecological processes. Its reliable prediction is a main issue in many environmental applications, which however is hampered by data scarcity, when using data‐demanding deterministic models, and modelling limitations, when using simpler statistical models. In this work we test a suite of models belonging to air2stream family, which are characterized by a hybrid formulation that combines a physical derivation of the key equation with a stochastic calibration of parameters. The air2stream models rely solely on air temperature and streamflow, and are of similar complexity as standard statistical models. The performances of the different versions of air2stream in predicting river water temperature are compared with those of the most common statistical models typically used in the literature. To this aim, a dataset of 38 Swiss rivers is used, which includes rivers classified into four different categories according to their hydrological characteristics: low‐land natural rivers, lake outlets, snow‐fed rivers and regulated rivers. The results of the analysis provide practical indications regarding the type of model that is most suitable to simulate river water temperature across different time scales (from daily to seasonal) and for different hydrological regimes. A model intercomparison exercise suggests that the family of air2stream hybrid models generally outperforms statistical models, while cross‐validation conducted over a 30‐year period indicates that they can be suitably adopted for long‐term analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.