Using CRD method for quantification of groundwater recharge in the Gaza Strip, Palestine

Rainfall is the main source of groundwater recharge in the Gaza Strip area in Palestine. The area is located in the semi-arid zone and there is no source of recharge other than rainfall. Estimation of groundwater recharge from rainfall is not an easy task since it depends on many uncertain parameters. The cumulative rainfall departure (CRD) method, which depends on the water balance principle, was used in this study to estimate the net groundwater recharge from rainfall. This method does not require much data as is the case with other classical recharge estimation methods. The CRD method was carried out using optimisation approach to minimise the root mean square error (RMSE) between the measured and the simulated groundwater head. The results of this method were compared with the results of other recharge estimation methods from literature. It was found that the results of the CRD method are very close to the results of the other methods, but with less data requirements and greater ease of application. Based on the CRD method, the annual amount of groundwater recharge from rainfall in the Gaza Strip is about 43 million m³. An erratum to this article can be found at     

Localization of mining-induced horizontal fractures along rock layer interfaces in overburden: field measurements and prediction

The locations of mining-induced horizontal fractures along rock interfaces in the overburden of Donetsk Coal Basin were identified using an original experimental device. The device traps methane from horizontal fracture zone (100–fold coal seam thickness) over an active longwall mining excavation. Presence or absence of horizontal fractures along rock layer interfaces is correlated with physical characteristics of the overburden, such as thickness, uniaxial compressive strength of overburden rock layers, location of rock layer interfaces and thickness of extracted coal seams. As a result, a combined criterion based on these physical characteristics is proposed to predict the presence of overburden horizontal fracturing in coal mine operations.     

Multi‐year water balance assessment of a newly constructed wetland,Fort McMurray,Alberta

Oil sands mining in Alberta transforms the boreal landscape of forests and wetlands into open pits, tailings ponds and overburden piles. Whereas reclamation efforts have primarily focused on upland forests, rebuilding wetland systems has recently become a motivation for research. Wetland creation and sustainability in this region is complicated by the sub‐humid climate and salinity of underlying mining material. In 2012, Syncrude Canada Ltd. completed the construction of the Sandhill Fen Watershed (SFW), a 52‐ha upland‐wetland system to evaluate wetland reclamation strategies on soft tailings. SFW includes an active pumping system, upland hummocks, a fen wetland and underdrains. To evaluate the influence of management practices on the hydrology of the system, this study reports the water balance from January 2013 to December 2014, the first 2 years after commissioning. A semi‐distributed approach was taken to examine the fluxes and stores of water in uplands and lowlands. Natural and artificial inputs and outputs were measured using a series of precipitation gauges and pumps, and evapotranspiration was quantified using three eddy covariance towers. A series of near surface wells recorded water table position. Both 2013 and 2014 were normal rainfall years, with 2013 having more and 2014 less snow than normal. In 2013, inflow/outflow from pumping was the predominant hydrological fluxes, resulting in considerable variability in water table position and storage changes throughout the summer. In 2014, the artificial addition of water was negligible, yet the water table remained near the surface in lowland locations, suggesting that wetland conditions could be maintained under current conditions. Evapotranspiration rates between uplands and lowlands were similar between years and sites, ranging from 2.2 ± 1.8 to 2.5 ± 1.2 mm/day and were largely controlled by climate. These rates were less than nearby older upland systems, suggesting that water balance partitioning will change as vegetation develops. Comparison between years and with natural systems provides insight on how management practices influence hydrologic dynamics and the overall water balance of the SFW. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Comparison of runoff characteristics of two adjacent basins in a tropical rainforest using a modified hydrologic cycle model with outflow

We propose a new runoff model including an outflow process that was applied to two adjacent basins (CL, TL) located in Lambir Hills National Park in north‐central Sarawak, Malaysia. Rainfall, runoff, topography, and soil layer thickness were observed. About 19% of annual runoff was observed in the CL basin (21.97 ha), whereas about 46% was observed in the TL basin (23.25 ha). It was inferred that the CL basin has an outflow because of low base flow, small runoff peak, and excessive water loss. By incorporating the outflow process into the HYdrological CYcle MODEL, good agreement between the data generated by the model and that observed was shown, with the exception of the data from the rainless period. Then, the fitting parameters for each basin were exchanged, except for the outflow parameter, and the characteristics of each basin were compared by calculating virtual runoff. As a result, the low base flow of the CL basin was estimated by the movement of the rainwater that escaped from the basin as deep percolation or lateral flow (11% of rainfall). The potential of the CL basin for mitigating flood and drought appeared to be higher than that of the TL basin. This is consistent with the topographic characteristics of the CL basin, which has a gentler slope than the TL basin. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Nine thousand years of upper montane soil/vegetation dynamics from the summit of Caratuva Peak,Southern Brazil

Biodiversity loss, climate change, and increased freshwater consumption are some of the main environmental problems on Earth. Mountain ecosystems can reduce these threats by providing several positive influences, such as the maintenance of biodiversity, water regulation, and carbon storage, amongst others. The knowledge of the history of these environments and their response to climate change is very important for management, conservation, and environmental monitoring programs. The genesis of the soil organic matter of the current upper montane vegetation remains unclear and seems to be quite variable depending on location. Some upper montane sites in the very extensive coastal Sea Mountain Range present considerable organic matter from the late Pleistocene and other from only the Holocene. Our study was carried out on three soil profiles (two cores in grassland and one in forest) on the Caratuva Peak of the Serra do Ibitiraquire (a sub-range of Sea Mountain Range – Serra do Mar) in Southern Brazil. The δ¹³C isotopic analyses of organic matter in soil horizons were conducted to detect whether C₃ or C₄ plants dominated the past communities. Complementarily, we performed a pollen analysis and ¹⁴C dating of the humin fraction to obtain the age of the studied horizons. Except for a short and probably drier period (between 6000 and 4500 cal yr BP), C₃ plants, including ombrophilous grasses and trees, have dominated the highlands of the Caratuva Peak (Pico Caratuva), as well as the other uppermost summits of the Serra do Ibitiraquire, since around 9000 cal yr BP. The Caratuva region represents a landscape of high altitude grasslands (campos de altitude altomontanos or campos altomontanos) and upper montane rain/cloud forests with soils that most likely contain some organic matter from the late Pleistocene, as has been reported in Southern and Southeastern Brazil for other sites. However, our results indicate that the studied deposits (near the summit) are from the early to late Holocene, when somewhat wetter and warmer conditions (since around 9000 cal yr BP) enabled a stronger colonization of the ridge of Pico Caratuva by mainly C₃ plants, especially grassland species. However, at the same time, even near the summit, the soil core from the forest site already presented the current physiognomy (or a shrubby/elfin or successional forest), indicating that the colonization of the neighboring uppermost saddles and valleys were probably populated mainly by upper montane forest species.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.