The effect of demand, tank parameters, and pumping stations on energy use in municipal drinking water distribution systems

Two of the indicators of the UN Millennium Development Goals ensuring environmental sustainability are energy use and per capita carbon dioxide emissions. The increasing urbanization and increasing world population may require increased energy use in order to transport enough safe drinking water to communities. In addition, the increase in water use would result in increased energy consumption, thereby resulting in increased green-house gas emissions that promote global climate change. The study of multiple Municipal Drinking Water Distribution Systems (MDWDSs) that relates various MDWDS aspects--system components and properties--to energy use is strongly desirable. The understanding of the relationship between system aspects and energy use aids in energy-efficient design. In this study, components of a MDWDS, and/or the characteristics associated with the component are termed as MDWDS aspects (hereafter--system aspects). There are many aspects of MDWDSs that affect the energy usage. Three system aspects (1) system-wide water demand, (2) storage tank parameters, and (3) pumping stations were analyzed in this study. The study involved seven MDWDSs to understand the relationship between the above-mentioned system aspects in relation with energy use. A MDWDSs model, EPANET 2.0, was utilized to analyze the seven systems. Six of the systems were real and one was a hypothetical system. The study presented here is unique in its statistical approach using seven municipal water distribution systems. The first system aspect studied was system-wide water demand. The analysis involved analyzing seven systems for the variation of water demand and its impact on energy use. To quantify the effects of water use reduction on energy use in a municipal water distribution system, the seven systems were modeled and the energy usage quantified for various amounts of water conservation. It was found that the effect of water conservation on energy use was linear for all seven systems and that all the average values of all the systems' energy use plotted on the same line with a high R 2 value. From this relationship, it can be ascertained that a 20% reduction in water demand results in approximately a 13% savings in energy use for all seven systems analyzed. This figure might hold true for many similar systems that are dominated by pumping and not gravity driven. The second system aspect analyzed was storage tank(s) parameters. Various tank parameters: (1) tank maximum water levels, (2) tank elevation, and (3) tank diameter were considered in this part of the study. MDWDSs use a significant amount of electrical energy for the pumping of water from low elevations (usually a source) to higher ones (usually storage tanks). The use of electrical energy has an effect on pollution emissions and, therefore, potential global climate change as well. Various values of these tank parameters were modeled on seven MDWDSs of various sizes using a network solver and the energy usage recorded. It was found that when averaged over all seven analyzed systems (1) the reduction of maximum tank water level by 50% results in a 2% energy reduction, (2) energy use for a change in tank elevation is system specific, and (2) a reduction of tank diameter of 50% results in approximately a 7% energy savings. The third system aspect analyzed in this study was pumping station parameters. A pumping station consists of one or more pumps. The seven systems were analyzed to understand the effect of the variation of pump horsepower and the number of booster stations on energy use. It was found that adding booster stations could save energy depending upon the system characteristics. For systems with flat topography, a single main pumping station was found to use less energy. In systems with a higher-elevation neighborhood, however, one or more booster pumps with a reduced main pumping station capacity used less energy. The energy savings for the seven systems was dependent on the number of boosters and ranged from 5% to 66% for the analyzed five systems with higher elevation neighborhoods (S3, S4, S5, S6, and S7). No energy savings was realized for the remaining two flat topography systems, S1, and S2. The present study analyzed and established the relationship between various system aspects and energy use in seven MDWDSs. This aids in estimating the amount of energy savings in MDWDSs. This energy savings would ultimately help reduce Greenhouse gases (GHGs) emissions including per capita CO 2 emissions thereby potentially lowering the global climate change effect. This will in turn contribute to meeting the MDG of ensuring environmental sustainability.

Gas-phase Dissociation of homo-DNA Oligonucleotides

Synthetic modified oligonucleotides are of interest for diagnostic and therapeutic applications, as their biological stability, pairing selectivity, and binding strength can be considerably increased by the incorporation of unnatural structural elements. Homo-DNA is an oligonucleotide homologue based on dideoxy-hexopyranosyl sugar moieties, which follows the Watson-Crick A-T and G-C base pairing system, but does not hybridize with complementary natural DNA and RNA. Homo-DNA has found application as a bioorthogonal element in templated chemistry applications. The gas-phase dissociation of homo-DNA has been investigated by ESI-MS/MS and MALDI-MS/MS, and mechanistic aspects of its gas-phase dissociation are discussed. Experiments revealed a charge state dependent preference for the loss of nucleobases, which are released either as neutrals or as anions. In contrast to DNA, nucleobase loss from homo-DNA was found to be decoupled from backbone cleavage, thus resulting in stable products. This renders an additional stage of ion activation necessary in order to generate sequence-defining fragment ions. Upon MS(3) of the primary base-loss ion, homo-DNA was found to exhibit unspecific backbone dissociation resulting in a balanced distribution of all fragment ion series.

In Situ Monitoring of Free-Phase Gas Accumulation and Release in Peatlands Using Ground Penetrating Radar (GPR)

We tested a set of surface common mid-point (CMP) ground penetrating radar (GPR) surveys combined with elevation rods ( to monitor surface deformation) and gas flux measurements to investigate in-situ biogenic gas dynamics and ebullition events in a northern peatland ( raised bog). The main findings are: ( 1) changes in the two-way travel time from the surface to prominent reflectors allow estimation of average gas contents and evolution of free-phase gas (FPG); ( 2) peat surface deformation and gas flux measurements are strongly consistent with GPR estimated changes in FPG content over time; ( 3) rapid decreases in atmospheric pressure are associated with increased gas flux; and ( 4) single ebullition events can induce releases of methane much larger ( up to 192 g/m(2)) than fluxes reported by others. These results indicate that GPR is a useful tool for assessing the spatial distribution, temporal variation, and volume of biogenic gas deposits in peatlands.

Spatial Variability in Biogenic Gas Accumulations in Peat Soils Is Revealed By Ground Penetrating Radar (GPR)

We performed surface and borehole ground penetrating radar (GPR) tests, together with moisture probe measurements and direct gas sampling to detect areas of biogenic gas accumulation in a northern peatland. The main findings are: (1) shadow zones (signal scattering) observed in surface GPR correlate with areas of elevated CH4 and CO2 concentration; (2) high velocities in zero offset profiles and lower water content inferred from moisture probes correlate with surface GPR shadow zones; (3) zero offset profiles depict depth variable gas accumulation from 0-10% by volume; (4) strong reflectors may represent confining layers restricting upward gas migration. Our results have implications for defining the spatial distribution, volume and movement of biogenic gas in peatlands at multiple scales.

Determination of respiratory gas flow by electrical impedance tomography in an animal model of mechanical ventilation

Background A recent method determines regional gas flow of the lung by electrical impedance tomography (EIT). The aim of this study is to show the applicability of this method in a porcine model of mechanical ventilation in healthy and diseased lungs. Our primary hypothesis is that global gas flow measured by EIT can be correlated with spirometry. Our secondary hypothesis is that regional analysis of respiratory gas flow delivers physiologically meaningful results. Methods In two sets of experiments n = 7 healthy pigs and n = 6 pigs before and after induction of lavage lung injury were investigated. EIT of the lung and spirometry were registered synchronously during ongoing mechanical ventilation. In-vivo aeration of the lung was analysed in four regions-of-interest (ROI) by EIT: 1) global, 2) ventral (non-dependent), 3) middle and 4) dorsal (dependent) ROI. Respiratory gas flow was calculated by the first derivative of the regional aeration curve. Four phases of the respiratory cycle were discriminated. They delivered peak and late inspiratory and expiratory gas flow (PIF, LIF, PEF, LEF) characterizing early or late inspiration or expiration. Results Linear regression analysis of EIT and spirometry in healthy pigs revealed a very good correlation measuring peak flow and a good correlation detecting late flow. PIFEIT = 0.702 · PIFspiro + 117.4, r2 = 0.809; PEFEIT = 0.690 · PEFspiro-124.2, r2 = 0.760; LIFEIT = 0.909 · LIFspiro + 27.32, r2 = 0.572 and LEFEIT = 0.858 · LEFspiro-10.94, r2 = 0.647. EIT derived absolute gas flow was generally smaller than data from spirometry. Regional gas flow was distributed heterogeneously during different phases of the respiratory cycle. But, the regional distribution of gas flow stayed stable during different ventilator settings. Moderate lung injury changed the regional pattern of gas flow. Conclusions We conclude that the presented method is able to determine global respiratory gas flow of the lung in different phases of the respiratory cycle. Additionally, it delivers meaningful insight into regional pulmonary characteristics, i.e. the regional ability of the lung to take up and to release air.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Identification, quantification, spatiotemporal distribution and genetic variation of major latex secondary metabolites in the common dandelion (Taraxacum officinale agg.)

The secondary metabolites in the roots, leaves and flowers of the common dandelion (Taraxacum officinale agg.) have been studied in detail. However, little is known about the specific constituents of the plant’s highly specialized laticifer cells. Using a combination of liquid and gas chromatography, mass spectrometry and nuclear magnetic resonance spectrometry, we identified and quantified the major secondary metabolites in the latex of different organs across different growth stages in three genotypes, and tested the activity of the metabolites against the generalist root herbivore Diabrotica balteata. We found that common dandelion latex is dominated by three classes of secondary metabolites: phenolic inositol esters (PIEs), triterpene acetates (TritAc) and the sesquiterpene lactone taraxinic acid β-d-glucopyranosyl ester (TA-G). Purification and absolute quantification revealed concentrations in the upper mg g−1 range for all compound classes with up to 6% PIEs, 5% TritAc and 7% TA-G per gram latex fresh weight. Contrary to typical secondary metabolite patterns, concentrations of all three classes increased with plant age. The highest concentrations were measured in the main root. PIE profiles differed both quantitatively and qualitatively between plant genotypes, whereas TritAc and TA-G differed only quantitatively. Metabolite concentrations were positively correlated within and between the different compound classes, indicating tight biosynthetic co-regulation. Latex metabolite extracts strongly repelled D. balteata larvae, suggesting that the latex constituents are biologically active.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

Distribution of cutin and suberin biomarkers under forest trees with different root systems

Background and aims Differences in chemical composition of root compounds and root systems among tree species may affect organic matter (OM) distribution, source and composition in forest soils. The objective of this study was to elucidate the contribution of species specific cutin and suberin biomarkers as proxies for shoot- and root-derived organic carbon (OC) to soil OM at different depths with increasing distance to the stems of four different tree species. Methods The contribution of cutin- and suberin-derived lipids to OM in a Cutanic Alisol was analyzed with increasing soil depth and distance to the stems of Fagus sylvatica L., Picea abies (L.) Karst., Quercus robur L. and Pseudotsuga menziesii (Mirb.) Franco. Cutin and suberin monomers of plants and soils were analyzed by alkaline hydrolysis and subsequent gas chromatography–mass spectrometry. Results The amount and distribution of suberin-derived lipids in soil clearly reflected the specific root system of the different tree species. The amount of cutin-derived lipids decreased strongly with soil depth, indicating that the input of leaf/needle material is restricted to the topsoil. In contrast to the suberin-derived lipids, the spatial pattern of cutin monomer contribution to soil OM did not depend on tree species. Conclusions Our results document the importance of tree species as a main factor controlling the composition and distribution of OM in forest soils. They reveal the impact of tree species on root-derived OM distribution and the necessity to distinguish among different zones when studying soil OM storage in forests.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.

Narrow dust jets in a diffuse gas coma: a natural product of small active regions on comets

Comets often display narrow dust jets but more diffuse gas comae when their eccentric orbits bring them into the inner solar system and sunlight sublimates the ice on the nucleus. Comets are also understood to have one or more active areas covering only a fraction of the total surface active with sublimating volatile ices. Calculations of the gas and dust distribution from a small active area on a comet’s nucleus show that as the gas moves out radially into the vacuum of space it expands tangentially, filling much of the hemisphere centered on the active region. The dust dragged by the gas remains more concentrated over the active area. This explains some puzzling appearances of comets having collimated dust jets but more diffuse gaseous atmospheres. Our test case is 67P/Churyumov–Gerasimenko, the Rosetta mission target comet, whose activity is dominated by a single area covering only 4% of its surface.

Effect of the Tiger Stripes on the water vapor distribution in Enceladus' exosphere

The jet activity emanating from Enceladus' exosphere south pole region observed by Cassini is a subject of intensive study. The in situ and remote sensing observations performed since 2005 triggered an active modeling campaign. Such modeling is essential for better understanding of the measurements performed by individual instruments as well as to link them for a more complete picture of the volatile and ice grain distribution in Enceladus' exosphere. This paper is focused on the investigation of the effect that diffuse gas sources along the Tiger Stripes have on distribution of the water vapor in Enceladus' exosphere using the updated version of our multiplume model. We have found that accounting for the gas production by Tiger Stripes is critical for interpretation of the Cassini data. According to our calculations, sources along the Tiger Stripes (apart from those originally identified by Spitale and Porco (2007)) must contribute about 23–32% to the total plume source rate, which varies in the range of (6.4–29) ×1027 s−1. The effect of the previously unidentified source suggested in the paper is found to be critical for explaining the Ultraviolet Imaging Spectrograph 2007 and 2010 observations in the whole range of the elapsed times.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Distribution of organic carbon in surface sediments covered during SONNE cruise SO184

Sediments were sampled and oxygen profiles of the water column were determined in the Indian Ocean off west and south Indonesia in order to obtain information on the production, transformation, and accumulation of organic matter (OM). The stable carbon isotope composition (d13Corg) in combination with C/N ratios depicts the almost exclusively marine origin of sedimentary organic matter in the entire study area. Maximum concentrations of organic carbon (Corg) and nitrogen (N) of 3.0% and 0.31%, respectively, were observed in the northern Mentawai Basin and in the Savu and Lombok basins. Minimum d15N values of 3.7 per mil were measured in the northern Mentawai Basin, whereas they varied around 5.4 per mil at stations outside this region. Minimum bottom water oxygen concentrations of 1.1 mL L**1, corresponding to an oxygen saturation of 16.1%, indicate reduced ventilation of bottom water in the northern Mentawai Basin. This low bottom water oxygen reduces organic matter decomposition, which is demonstrated by the almost unaltered isotopic composition of nitrogen during early diagenesis. Maximum Corg accumulation rates (CARs) were measured in the Lombok (10.4 g C m**-2 yr**-1) and northern Mentawai basins (5.2 g C m**-2 yr**-1). Upwelling-induced high productivity is responsible for the high CAR off East Java, Lombok, and Savu Basins, while a better OM preservation caused by reduced ventilation contributes to the high CAR observed in the northern Mentawai Basin. The interplay between primary production, remineralisation, and organic carbon burial determines the regional heterogeneity. CAR in the Indian Ocean upwelling region off Indonesia is lower than in the Peru and Chile upwellings, but in the same order of magnitude as in the Arabian Sea, the Benguela, and Gulf of California upwellings, and corresponds to 0.1-7.1% of the global ocean carbon burial. This demonstrates the relevance of the Indian Ocean margin off Indonesia for the global OM burial.

Vertical distribution of relative numbers of the main groups of sound-scattering fishes in the Central Indian Ocean

Findings made in 31 catches with an Isaacs-Kidd midwater trawl in the light (09.00-16.00) and dark (21.00-04.00) periods of a day within a survey area of about 100 sq. miles with approximate center coordinates of 13°S and 78°E have been used to investigate vertical distribution of the main groups of sound-scattering fishes (35 species of the family Myctophidae and 16 species of other families). It has been shown that during daylight hours all fishes sink to depths deeper than 400 m. Data are presented concerning the fish population of night-time sound-scattering layers at depths of 70-150 m and about 400 m and of the daytime ones at depths of about 450 m.