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.     

Oxyanion Concentrations in Eastern Sierra Nevada Rivers – 2. Arsenic and Phosphate

Water samples were collected from the Truckee River-Pyramid Lake system, the Walker River-Walker Lake system, and the Carson River, all located in eastern California and western Nevada, U.S.A., at three different times (i.e., summer 1991, spring 1992, and autumn 1992) over a two year period. The concentrations of As, Na, Cl, PO4, and pH were measured in these river samples and the associated terminal lakes. Arsenic values ranged from below 13 nmol/kg near Truckee, California to 160 nmol/kg at Nixon, Nevada in the Truckee River, from 40 nmol/kg in the headwaters of both West and East Walker Rivers to 270 nmol/kg below Weber Reservoir on the main branch of the Walker River, and from <27 nmol/kg to 234 nmol/kg for the lower Carson River system. Arsenic concentrations in Steamboat Creek (0.91 mol/kg–1.80 mol/kg) in the Truckee River catchment are above the U.S. EPA drinking water maximum contaminant level of 0.67 mol/kg, as are the As concentrations in both Pyramid Lake (1.33 mol/kg–1.57 mol/kg ) and Walker Lake (13.7 mol/kg–18.7 mol/kg). Sources of As for all three rivers include weathering of As-rich rocks and/or regolith and input of high-As geothermal spring waters, both processes primarily, although not exclusively, adding As to the headwater regions of these rivers. Steamboat Hot Springs (29 mol/kg As 54.5 mol/kg), for example, is identified as a source of As to the Truckee River via Steamboat Creek. The high As concentrations in Pyramid and Walker Lakes are likely due to (1) desorption of arsenate from aquatic particulate matter in these high pH waters (9.0 pH 9.5), (2) limited biologic uptake of arsenate, and (3) evaporative concentration of the lake waters. Evaluation of molar PO4}/As ratios of river waters and geothermal spring waters (e.g., Steamboat Hot Springs), indicates that phosphate is substantially enriched in Steamboat Creek as well as the mid to lower reaches of the Walker and Carson Rivers. These regions of each river are dominated by agricultural interests and, additionally, in the case of Steamboat Creek, residential areas and golf courses. Our data strongly imply that phosphate-rich agricultural return flow has likely added P to these streams and, consequently, increased their respective P:As ratios.     

Probing the Quiet Solar Atmosphere from the Photosphere to the Corona

We investigate the morphology and temporal variability of a quiet-Sun network region in different solar layers. The emission in several extreme ultraviolet (EUV) spectral lines through both raster and slot time-series, recorded by the EUV Imaging Spectrometer (EIS) on board the Hinode spacecraft is studied along with \(\mbox{H}\upalpha\) observations and high-resolution spectropolarimetric observations of the photospheric magnetic field. The photospheric magnetic field is extrapolated up to the corona, showing a multitude of large- and small-scale structures. We show for the first time that the smallest magnetic structures at both the network and internetwork contribute significantly to the emission in EUV lines, with temperatures ranging from \(8\times 10^{4}~\mbox{K}\) to \(6\times 10^{5}~\mbox{K}\). Two components of transition region emission are present, one associated with small-scale loops that do not reach coronal temperatures, and another component that acts as an interface between coronal and chromospheric plasma. Both components are associated with persistent chromospheric structures. The temporal variability of the EUV intensity at the network region is also associated with chromospheric motions, pointing to a connection between transition region and chromospheric features. Intensity enhancements in the EUV transition region lines are preferentially produced by \(\mbox{H}\upalpha\) upflows. Examination of two individual chromospheric jets shows that their evolution is associated with intensity variations in transition region and coronal temperatures.     

Controls on regional variability in marine pore-water diagenesis below the seafloor in Upper Jurassic–Lower Cretaceous prodeltaic sandstone and shales,Scotian Basin,Eastern Canada

Diagenesis in the uppermost Jurassic to Lower Cretaceous deltaic sandstones and shales of the Scotian Basin is an important control on reservoir quality. Ferruginous zone (sub-oxic) marine pore-water diagenesis controls the initial formation of Fe²⁺-silicates that are the precursors of grain-rimming chlorite that preserves porosity. This study assesses the regional controls on the type of marine pore-water diagenesis by studying the sedimentology, mineralogy, and geochemistry of the retrogradational units and underlying progradational units in parasequences from conventional cores in two wells in different parts of the basin. Coated grains preserve a record of whether marine pore-water diagenesis below the seafloor was dominantly in the ferruginous or sulphidic geochemical zone. Four types of coated grain were distinguished, each with a different mineral paragenesis. Mineralogical and chemical evidence of ferruginous zone diagenesis includes the presence of diagenetic chlorite and siderite, and the correlation of P with Fe or Ti. Pyrite and Fe-calcite are found where the sulphidic zone is more significant than the ferruginous zone. Ferruginous zone diagenesis was common in low-sedimentation rate retrogradational sediments with low organic carbon, and in delta-front turbidites and river-mouth sandstones. Estuarine, tidal flat and prodeltaic facies that are directly supplied by riverine sediments have a lower Fe:Ti ratio than do fully marine shoreface and open shelf facies as a result of input of detrital ilmenite and its alteration products. The relative contribution of colloidal iron (hydr)-oxides appears greater in distal low-sedimentation rate environments. Where large changes in sedimentation rate occurred at ravinement surfaces, the underlying progradational rocks have evidence of ferruginous zone diagenesis, whatever their facies. Rapid upward migration of the pore-water profile resulting from the change in sedimentation rate reduced the time available for mineral products to form in the deeper pore-water zones. This study has shown that the availability of Fe and organic carbon varying in a complex manner in marine deltaic sediments, but that the resulting diagenesis by marine pore-water can be predicted from facies and paleogeographic setting.     

Geochemistry, tectonic setting and metamorphism of mid-triassic volcanic rocks of Greece

Mid-Triassic volcanic rocks are preserved in about eight localities in the nappes of the external Hellenides, and are also widespread in the Pelagonian zone of central Greece. Most have suffered low-grade metamorphism. Immobile trace elements from 35 samples, 80 major element analyses, and pyroxene analyses show that these rocks are the product of an eastward dipping subduction zone. The subducted ocean was perhaps analogous to the Pennine zone of the Alps. The Othris and Pindos ophiolites are probably contemporaneous products of related back-arc extension.Metamorphic mineral assemblages in volcanic rocks of the Phyllite Series of the external Hellenides indicate high-pressure, low-temperature conditions, probably related to a second subduction that occurred prior to the main phase of Hellenide nappe emplacement in the Miocene.     

The relationship between length and width of plutons within the crustal-scale Cobequid Shear Zone, northern Appalachians, Canada

The lengths and widths have been measured for 69 component bodies of composite plutons along the Cobequid Shear Zone. Plutons on major fault strands, those with mylonite zones >0.1 km wide, exhibit evidence of multiple intrusion of magma batches. Small plutons along short faults in stepover zones appear related to rapid emplacement of magma in bodies 1.5–4 km long by 0.1–2 km wide. Such small plutons show low enrichment in incompatible elements in older component bodies, but increasing amounts in younger bodies as a result of progressive magma expulsion from crystal mush during crystallization and shear-enhanced compaction in fault zones. Wider plutons generally occur along longer fault strands accommodating more strain and penetrating deeper into the crust and show enrichment in incompatible elements. The width of the mylonitic fault zone is about 15% of the width of these plutons. The length-to-width ratio of component bodies and composite plutons varies between 2 and 11. The best-fit line describing these data has a slope of 1.056, which implies scaling behavior between plutonism and tectonic processes. Scalar properties of plutonic bodies are similar to those of faults, but scalar relationships observed in component bodies do not apply to composite plutons.     

Distinguishing active and legacy source contributions to stream water quality: Comparative quantification for chloride and metals

Hydrochemical constituents in streams may originate from currently active sources at the surface and/or legacy sources from earlier surface inputs, waste deposits and land contamination. Distinction and quantification of these source contributions are needed for improved interpretation of tracer data and effective reduction of waterborne environmental pollutants. This article develops a methodology that recognizes and quantifies some general mechanistic differences in stream concentration and load behavior versus discharge between such source contributions. The methodology is applied to comparative analysis of stream concentration data for chloride (Cl⁻), copper (Cu), lead (Pb), and zinc (Zn), and corresponding data for water discharge, measured over the period 1990–2018 in multiple hydrological catchments (19 for Cl⁻, 11 for Cu and Zn, 10 for Pb) around the major Lake Mälaren in Sweden. For Cl⁻, the average load fraction of active sources is quantified to be 19%, and the average active and legacy concentration contributions as 2.9 and 11 mg/L, respectively. For the metals, the average active load fractions at outlets are 1%–3% over all catchments and 9%–14% in the relatively few catchments with mixed metal sources. Average active and legacy concentration contributions are 0.14 and 3.2 μg/L for Cu, 0.05 and 1.5 μg/L for Pb, and 1.4 and 12 μg/L for Zn, respectively. This multi-catchment analysis thus indicates a widespread prevalence of legacy sources, with greater legacy than active concentration contributions for both Cl⁻ and the metals, and active contributions playing a greater role for chloride than for the metals. The relatively simple first-order methodology developed and applied in the study can be used to screen commonly available stream monitoring data for possible distinction of active and legacy contributions of any hydrochemical constituent in and across various hydrological catchment settings.     

Nutrient source attribution: Quantitative typology distinction of active and legacy source contributions to waterborne loads

Distinction between active and legacy sources of nutrients is needed for effective reduction of waterborne nutrient loads and associated eutrophication. This study quantifies main typological differences in nutrient load behaviour versus water discharge for active and legacy sources. This quantitative typology is used for source attribution based on monitoring data for water discharge and concentrations of total nitrogen (TN) and total phosphorous (TP) from 37 catchments draining into the Baltic Sea along the coastline of Sweden over the period 2003–2013. Results indicate dominant legacy source contributions to the monitored loads of TN and TP in most (33 of the total 37) study catchments. Dominant active sources are indicated in 1 catchment for TN, and mixed sources are indicated in 3 catchments for TN, and 4 catchments for TP. The TN and TP concentration contributions are quantified to be overall higher from the legacy than the active sources. Legacy concentrations also correlate well with key indicators of human activity in the catchments, agricultural land share for TN (R² = 0.65) and population density for TP (R² = 0.56). Legacy-dominated nutrient concentrations also change more slowly than in catchments with dominant active or mixed sources. Various data-based results and indications converge in indicating legacy source contributions as largely dominant, mainly anthropogenic, and with near-zero average change trends in the present study of catchments draining into the Baltic Sea along the coastline of Sweden, as in other parts of the world. These convergent indications emphasize needs to identify and map the different types of sources in each catchment, and differentiate strategies and measures to target each source type for possible achievement of shorter- and longer-term goals of water quality improvement.     

Observational Investigation of the Statistical Properties of Surface-Layer Turbulence in a Suburban Area of São Paulo,Brazil: Objective Analysis of Scaling-Parameter Accuracy and Uncertainties

Boundary-Layer Meteorology - Statistical properties of turbulence, specifically variances of velocity components, temperature, water vapor, and carbon dioxide densities, are observationally...     

Permafrost degradation and subsurface-flow changes caused by surface warming trends

Change dynamics of permafrost thaw, and associated changes in subsurface flow and seepage into surface water, are analysed for different warming trends in soil temperature at the ground surface with a three-phase two-component flow system coupled to heat transport. Changes in annual, seasonal and extreme flows are analysed for three warming-temperature trends, representing simplified climate-change scenarios. The results support previous studies of reduced temporal variability of groundwater flow across all investigated trends. Decreased intra-annual flow variability may thus serve as an early indicator of permafrost degradation before longer-term changes in mean flows are notable. This is advantageous since hydrological data are considerably easier to obtain, may be available in longer time series, and generally reflect larger-scale conditions than direct permafrost observations. The results further show that permafrost degradation first leads to increasing water discharge, which then decreases as the permafrost degradation progresses further to total thaw. The most pronounced changes occur for minimum annual flows. The configuration considered represents subsurface discharge from a generic heterogeneous soil-type domain.