Water balance prediction in stormwater infiltration basins using 2-D modeling: An application to evaluate the clogging process

Sustainable urban drainage systems are built along roads and in urban areas to collect urban runoff and avoid flooding, and to filter water pollutants. Sediment collected by runoff is deposited in the stormwater basin and progressively reduces water infiltration efficiency, leading to the clogging of the basin. To help stormwater basin managers and stakeholders better understand and predict clogging rates in order to elaborate maintenance plans and schedules, water transport prediction models are necessary. However,because of the heterogeneous sediment hydrodynamic properties inside the stormwater basin, a twodimensional(2-D) water flow model is required to predict water levels and possible overflow as accurately as possible. Saturated hydraulic conductivity(Ks) and sediment water retention curves were measured in the overall sediment layer of the stormwater basin, in addition to sediment layer thickness and organic matter content(11 sampling points). Sediment depth was used to predict organic matter(OM) content, and the OM was used to predict Ks. Water height in the basin was modeled with the HYDRUS-2 D model by taking into account the sediment hydrodynamic properties distribution. The HYDRUS-2 D model gave a satisfactory representation of the measured data. Scenarios of the hydraulic properties of stormwater basin sediment were tested over time, and hydraulic resistance, R, was calculated to assess the stormwater basin performance. Presently, after 20 years of functioning, the stormwater basin still ensures efficient water infiltration, but the first outflow(Hydraulic resistance,R 24 h)) is expected to appear in the next 5 years, and clogging(R 47 h) in the next 13 years. This 2-D water balance model makes it possible to integrate the hydrodynamic heterogeneity of a stormwater basin. It gives interesting perspectives to better predict 2-D/3-D contaminant transport.     

Simulation of Turbulent Flow in An Urban Forested Park Damaged by a Windstorm

The damage caused by windstorms to forest ecosystems is often very heterogeneous. In order to improve the stability of forested landscapes, it is of great importance to identify the factors responsible for this spatial variability. The structure of the landscape itself may play a role, through possible influences of canopy heterogeneities on the development of turbulence. For the purpose of investigating the role of landscape fragmentation on turbulence development, we used a numerical flow model with a k–ε turbulence scheme model, previously validated in simple cases with well-defined surface changes (roughness change and forest edge flow). A series of two- and three-dimensional simulations were performed over a heterogeneous urban forested park in Europe, which was severely damaged in various places by the Lothar windstorm in December 1999. The model shows the development of a region of strong turbulence, resulting from the generation of large wind shear at the top of the canopy. A sensitivity study shows how the location, extension and intensity of the region depend on canopy characteristics such as the leaf density, the nature of the edge or the presence of gaps and clearings. Simulations performed in conditions representative of the windstorm show that the location of the damaged areas corresponds very closely to the regions where the turbulent kinetic energy was above a certain threshold.     

Turbulent Structures in a Pine Forest with a Deep and Sparse Trunk Space: Stand and Edge Regions

Forested landscapes often exhibit large spatial variability in vertical and horizontal foliage distributions. This variability may affect canopy-atmosphere exchanges through its action on the development of turbulent structures. Here we investigate in neutral stratification the turbulent structures encountered in a maritime pine forest characterized by a high, dense foliated layer associated with a deep and sparse trunk space. Both stand and edge regions are considered. In situ measurements and the results of large-eddy simulations are used and analyzed together. In stand conditions, far from the edge, canopy-top structures appear strongly damped by the dense crown layer. Turbulent wind fluctuations within the trunk space, where the momentum flux vanishes, are closely related to these canopy-top structures through pressure diffusion. Consequently, autocorrelation and spectral analyses are not quite appropriate to characterize the vertical scale of coherent structures in this type of canopy, as pressure diffusion enhances the actual scale of structures. At frequencies higher than those associated with canopy-top structures, wind fluctuations related to wake structures developing behind tree stems are observed within the trunk space. They manifest themselves in wind velocity spectra as secondary peaks in the inertial subrange region, confirming the hypothesis of spectral short-cuts in vegetation canopies. In the edge region specific turbulent structures develop just below the crown layer, in addition to canopy-top structures. They are generated by the wind shear induced by the sub-canopy wind jet that forms at the edge. These structures provide a momentum exchange mechanism similar to that observed at the canopy top but in the opposite direction and with a lower magnitude. They may develop as in plane mixing-layer flows, with some perturbations induced by canopy-top structures. Wake structures are also observed within the trunk space in the edge region.     

Interseismic Strain Localization in the San Jacinto Fault Zone

We investigate interseismic deformation across the San Jacinto fault at Anza, California where previous geodetic observations have indicated an anomalously high shear strain rate. We present an updated set of secular velocities from GPS and InSAR observations that reveal a 2–3 km wide shear zone deforming at a rate that exceeds the background strain rate by more than a factor of two. GPS occupations of an alignment array installed in 1990 across the fault trace at Anza allow us to rule out shallow creep as a possible contributor to the observed strain rate. Using a dislocation model in a heterogeneous elastic half space, we show that a reduction in shear modulus within the fault zone by a factor of 1.2–1.6 as imaged tomographically by Allam and Ben-Zion (Geophys J Int 190:1181–1196, 2012) can explain about 50 % of the observed anomalous strain rate. However, the best-fitting locking depth in this case (10.4 ± 1.3 km) is significantly less than the local depth extent of seismicity (14–18 km). We show that a deep fault zone with a shear modulus reduction of at least a factor of 2.4 would be required to explain fully the geodetic strain rate, assuming the locking depth is 15 km. Two alternative possibilities include fault creep at a substantial fraction of the long-term slip rate within the region of deep microseismicity, or a reduced yield strength within the upper fault zone leading to distributed plastic failure during the interseismic period.     

The Buday’ah Formation, Sultanate of Oman: A Middle Permian to Early Triassic oceanic record of the Neotethys and the late Induan microsphere bloom

The Middle Permian to Lower Triassic Buday’ah section, exposed in the Oman Mountains, is the first deep-sea section to be described in the Neotethys. The oceanic sediments were deposited along the southern Tethys margin in the newly formed Hawasina Basin. It is one of the few places where true Tethyan Permian radiolarites are exposed that allow the documentation of CCD evolution through time. The succession begins as oceanic crust pillow basalt with red ammonoid-rich pelagic limestone occurring both above and within inter-pillow cavities; the new occurrence of Clarkina postbitteri hongshuiensis indicates a late Capitanian age for the carbonate. The sharp change to overlying late Capitanian to Changhsingian radiolarite reflects rapid subsidence about 10 Myrs after initial continental breakup that resulted in the formation of the Neotethys Ocean. New conodonts indicate that the Permian-Triassic boundary succession occurs in the first platy lime mudstone beds above a Changhsingian siliceous to calcareous shale unit. The platy lime mudstone beds include an Upper Griesbachian bloom of calcite filled spheres (radiolarians?) that marks a potential world-wide event. New conodonts indicate an early Olenekian age for overlying grey papery limestone that are devoid of both macrofossils and trace fossils indicating that recovery from the Late Permian extinction has not yet progressed within this deep-water environment.δ¹³C_org, isotope values have not been disturbed and they show a negative shift just below the Permian-Triassic transition and a second one at the parvus zone level above. The Buday’ah succession may represent the most distal and probably deepest Permian and Lower Triassic depositional sequence within the basin.     

Exploring the spatiotemporal variability of the snow water equivalent in a small boreal forest catchment through observation and modelling

In snow-fed catchments, it is crucial to monitor and model the snow water equivalent (SWE), particularly when simulating the melt water runoff. SWE distribution can, however, be highly heterogeneous, particularly in forested environments. Within these locations, scant studies have explored the spatiotemporal variability in SWE in relation with vegetation characteristics, with only few successful attempts. The aim of this paper is to fill this knowledge gap, through a detailed monitoring at nine locations within a 3.49 km² forested catchment in southern Québec, Canada (47°N, 71°W). The catchment receives an annual average of 633 mm of solid precipitation and is predominantly covered with balsam fir stands. Extracted from intensive field campaign and high-resolution LiDAR data, this study explores the effect of fine scale forest features (tree height, tree diameter, canopy density, leaf area index [LAI], tree density and gap fraction) on the spatiotemporal variability in the SWE distribution. A nested stratified random sampling design was adopted to quantify small-scale variability across the catchment and 1810 manual snow samples were collected throughout the consecutive winters of 2016–17 and 2017–18. This study explored the variability of SWE using coefficients of variation (CV) and relating to the LAI. We also present existing spatiotemporal differences in maximum snow depth across different stands and its relationship with average tree diameter. Furthermore, exploiting key vegetation characteristics, this paper explores different approaches to model SWE, such as multiple linear regression, binary regression tree and neural networks (NN). We were unable to establish any relationship between the CV of SWE and the LAI. However, we observed an increase in maximum snow depth with decreasing tree diameter, suggesting an association between these variables. NN modelling (Nash-Sutcliffe efficiency [NSE] = 0.71) revealed that, snow depth, snowpack age and forest characteristics (tree diameter and tree density) are key controlling variables on SWE. Using only variables that are deemed to be more readily available (snow depth, tree height, snowpack age and elevation), NN performance falls by only 7% (NSE = 0.66).     

Considering groundwater use to improve the assessment of groundwater pumping for irrigation in North Africa

Groundwater resources in semi-arid areas and especially in the Mediterranean face a growing demand for irrigated agriculture and, to a lesser extent, for domestic uses. Consequently, groundwater reserves are affected and water-table drops are widely observed. This leads to strong constraints on groundwater access for farmers, while managers worry about the future evolution of the water resources. A common problem for building proper groundwater management plans is the difficulty in assessing individual groundwater withdrawals at regional scale. Predicting future trends of these groundwater withdrawals is even more challenging. The basic question is how to assess the water budget variables and their evolution when they are deeply linked to human activities, themselves driven by countless factors (access to natural resources, public policies, market, etc.). This study provides some possible answers by focusing on the assessment of groundwater withdrawals for irrigated agriculture at three sites in North Africa (Morocco, Tunisia and Algeria). Efforts were made to understand the different features that influence irrigation practices, and an adaptive user-oriented methodology was used to monitor groundwater withdrawals. For each site, different key factors affecting the regional groundwater abstraction and its past evolution were identified by involving farmers’ knowledge. Factors such as farmer access to land and groundwater or development of public infrastructures (electrical distribution network) are crucial to decode the results of well inventories and assess the regional groundwater abstraction and its future trend. This leads one to look with caution at the number of wells cited in the literature, which could be oversimplified.     

A combined modelling and geochemical study of the fate of terrigenous inputs from mixed natural and mining sources in a coral reef lagoon (New Caledonia)

Open-cast mining for Ni, Cr and Co was conducted in the south-west part of New Caledonia during the 20th century. Abandoned mining and prospecting sites were severely affected by erosion, resulting in an increase in the load of terrigenous particles transported to the coral reef lagoon. This article assesses the impact of a typical small catchment area (La Coulée River, 85 km2 watershed) on two bays (Boulari and Sainte Marie) located near Noumea, New Caledonia's main city. This multi-disciplinary study combines geochemical, sedimentological, and hydrodynamic approaches. Ni and Cr concentrations contained in the geochemical matrix phase of the pelitic fraction were determined. The study of the geochemical signatures together with sedimentological data and 3D numerical simulations of the transport of non-settling particles throughout the lagoon demonstrated that terrigenous inputs from the Coulée River were mainly transported and deposited along the shoreline, reaching areas as distant as Sainte Marie Bay. Although quantitatively low (about 3% of the pelite mass of the bay sediments), the terrigenous inputs in Sainte Marie Bay originating from the Coulée River could be traced. The metal content in suspended matter was over 7000 mg kg(-1) for Ni and 4200 mg kg(-1) for Cr.