Estimation of unrecorded groundwater abstractions in Jordan through regional groundwater modelling

Jordan suffers from water scarcity and groundwater covers the majority of Jordan’s water supply. Therefore, there is an urgent need to manage this resource conscientiously. A regional numerical groundwater flow model, developed as part of a decision support system for the country of Jordan, allows for quantification of the overexploitation of groundwater resources and enables determination of the extent of unrecorded agricultural groundwater abstraction. Groundwater in Jordan is abstracted from three main aquifers partly separated by aquitards. With updated geological, structural, and hydrogeological data available in the country, a regional numerical groundwater flow model for the whole of Jordan and the southernmost part of Syria was developed using MODFLOW. It was first calibrated for a steady-state condition using data from the 1960s, when groundwater abstraction was negligible. After transient calibration using groundwater level measurements from all aquifers, model results reproduce the large groundwater-level declines experienced in the last decades, which have led to the drying out of numerous springs. They show a reversal of groundwater flow directions in some regions, due to over-abstraction, and demonstrate that documented abstractions are not sufficient to cause the observed groundwater-level decline. Only after considering irrigation water demand derived from remote sensing data, the model is able to simulate these declines. Illegal abstractions can be quantified and predictive scenarios show the potential impact of different management strategies on future groundwater resources.

     

Characterization of deep saline aquifers in the Bécancour area,St. Lawrence Lowlands,Québec,Canada: implications for CO2 geological storage

The present paper provides a case study of the assessment of the potential for CO₂ storage in the deep saline aquifers of the Bécancour region in southern Québec. This assessment was based on a hydrogeological and petrophysical characterization using existing and newly acquired core and well log data from hydrocarbon exploration wells. Analyses of data obtained from different sources provide a good understanding of the reservoir hydrogeology and petrophysics. Profiles of formation pressure, temperature, density, viscosity, porosity, permeability, and net pay were established for Lower Paleozoic sedimentary aquifers. Lateral hydraulic continuity is dominant at the regional scale, whereas vertical discontinuities are apparent for most physical and chemical properties. The Covey Hill sandstone appears as the most suitable saline aquifer for CO₂ injection/storage. This unit is found at a depth of more than 1 km and has the following properties: fluid pressures exceed 14 MPa, temperature is above 35 °C, salinity is about 108,500 mg/l, matrix permeability is in the order of 3 × 10^?16 m² (0.3 mDarcy) with expected higher values of formation-scale permeability due to the presence of natural fractures, mean porosity is 6 %, net pay reaches 282 m, available pore volume per surface area is 17 m³/m², rock compressibility is 2 × 10^?9 Pa^?1 and capillary displacement pressure of brine by CO₂ is about 0.4 MPa. While the containment for CO₂ storage in the Bécancour saline aquifers can be ensured by appropriate reservoir characteristics, the injectivity of CO₂ and the storage capacity could be limiting factors due to the overall low permeability of aquifers. This characterization offers a solid basis for the subsequent development of a numerical hydrogeological model, which will be used for CO₂ injection capacity estimation, CO₂ injection scenarios and risk assessment.     

Neoproterozoic Volhynia-Brest magmatic province in the western East European craton: Within-plate magmatism in an ancient suture zone

The reasons for the isotopic and geochemical heterogeneity of magmatism of the Neoproterozoic large Volhynia-Brest igneous province (VBP) are considered. The province was formed at 550 Ma in response to the break up of the Rodinia supercontinent and extends along the western margin of the East European craton, being discordant to the Paleoproterozoic mobile zone that separates Sarmatia and Fennoscandia and the Mesoproterozoic Volhynia-Orsha aulacogen. The basalts of VBP show prominent spatiotemporal geochemical zoning. Based on petrographic, mineralogical, geochemical, and isotopic data, the following types of basalts can be distinguished: olivine-normative subalkaline basalts consisting of low-Ti (sLT, < 1.10–2.0 wt % TiO₂; ε_Nd(550) from ?6.6 to ?2.7) and medium-Ti (sMT, 2.0–3.0 wt % TiO₂, occasionally up to 3.6 wt % TiO₂; ε_Nd(550) from ?3.55 to + 0.6) varieties; normal quartz-normative basalts (tholeiites) including low-Ti (tLT, < 1.75–2.0 wt % TiO₂) and medium-to-high-Ti (tHT1, 2.0–3.6 wt % TiO₂, ε_Nd(550) from ?1.3 to + 1.0) varieties. The hypabyssal bodies are made up of subalkaline low-Ti olivine dolerites (LT, 1.2–1.5 wt % TiO₂; ε_Nd(550) = ?5.8) and subalkaline high-Ti olivine gabbrodolerites (HT2, 3.0–4.5 wt % TiO₂; ε_Nd(550) = ?2.5). Felsic rocks of VBP are classed as volcanic rocks of normal (andesidacites, dacites, and rhyodacites) and subalkaline (trachyrhyodacites) series with TiO₂ 0.72–0.77 wt% and ε_Nd(550) of ?12. The central part of VBP is underlain by a Paleoproterozoic domain formed by continent-arc accretion and contains widespread sills of HT2 dolerites and lavas of LT basalts; the northern part of the province is underlain by the juvenile Paleoproterozoic crust dominated by MT and HT1 basalts. MT and LT basalts underwent significant AFC-style upper crustal contamination. During their long residence in the upper crustal magmatic chambers, the basaltic melts fractionated and caused notable heating of the wall rocks and, correspondingly, nonmodal melting of the upper crustal protolith containing high-Rb phase (biotite), thus producing the most felsic rocks of the province. The basalts of VBP were derived from geochemically different sources: probably, the lithosphere and a deep-seated plume (PREMA type). The HT2 dolerites were generated mainly from a lithospheric source: by 3–4% melting of the geochemically enriched garnet lherzolite mantle. LT dolerites were obtained by partial melting of the modally metasomatized mantle containing volatile-bearing phases. The concepts of VBP formation were summarized in the model of three-stage plume-lithosphere interaction.     

How Well Does Chlorophyll Explain the Seasonal Variation in Phytoplankton Activity?

The seasonal variation in phytoplankton activity is determined by analysing 1385 primary production (PP) profiles, chlorophyll a (Chl) concentration profiles and phytoplankton carbon biomass concentrations (C) from the period 1998–2012. The data was collected at six different stations in the Baltic Sea transition zone (BSTZ) which is a location with strong seasonal production patterns with light as the key parameter controlling this productivity. We show that the use of Chl as a proxy for phytoplankton activity strongly overestimates the contribution from the spring production to annual pelagic carbon flow. Spring (February and March) Chl comprised 16–30% of the total annual Chl produced, whereas spring C was much lower (8–23%) compared to the annual C. Spring PP accounted for 10–18% of the total annual PP, while the July–August production contributed 26–33%, i.e. within the time frame when zooplankton biomass and grazing pressure are highest. That is, Chl failed in this study to reflect the importance of the high summer PP. A better proxy for biomass may be C, which correlated well with the seasonal pattern of PP (Pearson correlation, p < 0.05). Thus, this study suggests to account for the strong seasonal pattern in C/Chl ratios when considering carbon flow in coastal systems. Seasonal data for PP were fitted to a simple sinusoidal wave model describing the seasonal distribution of PP in the BSTZ and were proposed to present a better parameterizaton of PP in shallow stratified temperate regions than more commonly applied proxies.     

Response of A Gravel - Bed River To Dam Closure: Insights From Sediment Transport Processes And Channel Morphodynamics

This work deals with the impacts of dams on large gravel -bed rivers in terms of altering coarse transport regimes and the relationship with river morphodynamics. Using data collected by a tracer -based monitoring programme carried out in a 4 -km -long study sector of the Parma River (Italy), located downstream from a relatively recently established dam, we applied a virtual velocity approach to estimate the coarse bed material load at four river cross -sections. Monitoring and calculation results provided new insights into the impacts of the dam on streambed material mobility and the sediment regime over the 17 -month calculation period. A longitudinal gradient of effects was observed along the study sector. Sections located closer to the dam are characterized by more evident impacts due to deficits in coarse sediment input from upstream. Sediment mobility here is strongly altered, especially in the highly armoured main channel, and the overall bed material load is extremely low. A partial recovery of sediment dynamics was observed at the sections located further from the dam, where estimates indicate higher sediment yield. The observed longitudinal trend in the coarse sediment transport regime matches the morphology, as the river shifts downstream from a sinuous configuration with alternate bars to a wandering one. The novel insights into alteration of coarse sediment dynamics and the relationship with river morphodynamics are potentially applicable to many other fluvial contexts affected by similar impoundments. © 2019 John Wiley & Sons, Ltd.     

Numerical simulation of the electric field and zonal current in the Earth’s ionosphere: The dynamo field and equatorial electrojet

The work is devoted to the numerical simulation of the dynamo electric field and its effects in the Earth’s ionosphere within the scope of the thermosphere-ionosphere-protonosphere global self-consistent model developed at WD IZMIRAN. The new electric field calculation block, which was used to obtain results of the self-consistent calculations of the electric field potential generated by the dynamo effect of the thermospheric winds (the dynamo field) and the equatorial electrojet for March 22, 1987, is briefly described in this work. A comparison of the obtained results with the experimental data showed a satisfactory agreement. Moreover, the proposed model was used to calculate the diurnal variations in the ionospheric parameters for Jicamarca equatorial station under the same conditions with the help of the new block of the electric field. The results of these calculations are also presented and discussed in this work. It has been indicated that the new model satisfactorily describes the specific features of electric field distribution at the geomagnetic equator and the well-known phenomenon of equatorial electrojet.     

Validating and assessing the sensitivity of the climate model with an ocean general circulation model developed at the Institute of Atmospheric Physics,Russian Academy of Sciences

A new version of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), climate model (CM) has been developed using an ocean general circulation model instead of the statistical-dynamical ocean model applied in the previous version. The spatial resolution of the new ocean model is 3° in latitude and 5° in longitude, with 25 unevenly spaced vertical levels. In the previous version of the oceanic model, as in the atmospheric model, the horizontal resolution was 4.5° in latitude and 6° in longitude, with four vertical levels (the upper quasi-homogeneous layer, seasonal thermocline, abyssal ocean, and bottom friction layer). There is no correction for the heat and momentum fluxes between the atmosphere and ocean in the new version of the IAP RAS CM. Numerical experiments with the IAP RAS CM have been performed under current initial and boundary conditions, as well as with an increasing concentration of atmospheric carbon dioxide. The main simulated atmospheric and oceanic fields agree quite well with observational data. The new version’s equilibrium temperature sensitivity to atmospheric CO₂ doubling was found to be 2.9 K. This value lies in the mid-range of estimates (2–4.5 K) obtained from simulations with state-of-the-art models of different complexities.     

Isotope Constraints on the Aquatic Carbon Budget: Langat Watershed,Malaysia

Langat River drains a tropical watershed in the southwest of the Malaysian Peninsula. The watershed is heavily urbanized in its downstream portion. Water samples were collected from May 2010 to December 2011, at three localities along the main stem river, 1 location at its Semenyih tributary and from an upstream groundwater source. Concentration and δ¹³C data of riverine DIC and DOC indicate the dominance of C3 plant-derived material as the primary source of carbon, with δ¹³C_DIC values enriched in ¹³C relative to that of the C3 source. This enrichment is likely due to CO₂ outgassing, as calculated concentrations of riverine CO₂ are significantly higher than ambient atmospheric values, with methanogenic activity a theoretically possible contributing factor, particularly at the upstream location. The Langat River therefore acts as a net source of CO₂, with a total sub-basin flux of 19.7 × 10³ t C year^?1. This is comparable to the sum of riverine DOC, DIC and POC loss rates from the sub-basin, calculated as 24.5 × 10³ t C year^?1, and highlights the significance of CO₂ evasion from water bodies to the atmosphere for balancing the budget of the terrestrial carbon cycle. The DIC and DOC concentration and δ¹³C data also suggests that in the more urbanized downriver areas, much of the organic carbon input may be anthropogenicaly derived due to ubiquity of sewage treatment plants and landfill sites. Such human-induced perturbations to riverine carbon cycling should be taken into account in future studies of urbanized watersheds.     

Geochemistry of the Adamello massif (northern Italy)

The Tertiary Adamello massif, outcropping over an area of more than 550 km² in the southern Alps (northern Italy) is composed mainly of granitoid rocks (granodiorite, tonalite, quartz diorite) with minor amounts of diorite and gabbro. The major and trace element composition of these rocks is comparable to calc-alkaline volcanic rocks of continental margins. The granitoid rocks display spatial and temporal variations in their composition, particularly in Na, P, Sr, La, Nb and Y contents and ⁸⁷Sr/⁸⁶Sr ratios. The variations were probably produced by concurrent contamination/wall-rock assimilation and fractional crystallization of high-alumina basaltic magma.     

U-Pb isotopic dating of zircons (SHRIMP-II) from granulites of the Ol’khon region of Western Baikal area

Doklady Earth Sciences -