Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Lake bank filtration at Nainital,India: water-quality evaluation

There are different water-supply schemes in Uttarakhand, India to tap the water from streams, rivers and lakes. At Nainital, seven tube-wells (depths 22.6–36.7 m), located at a distance of <100 m from the lake, are being used to abstract (1) lake water after passage through the soil and (2) subsurface water/groundwater flowing towards the lake. Water samples from the lake and tube-wells were analyzed in monsoon and non-monsoon periods from 1997 to 2006. Total dissolved solids, EC, alkalinity and hardness were found to be marginally greater in tube-well waters. The difference in hydrochemistry of tube-well water was mainly due to variation in flow regimes during monsoon and non-monsoon periods. Results clearly indicate that lake water as such is not potable as it contains unacceptable levels of organic matter in terms of COD (~44 mg/L), coliforms (~15.6 × 10⁴ MPN/100 mL) and nutrients. Coliform bacteria and COD have not been detected in any of the tube-well water samples over the years. Lake water, treated by sand filters did not conform to drinking water standards. These investigations have led to the closure of the treatment facility and installation of two tube-wells in addition to the existing five tube-wells.

---

Résumé  Il existe divers projets d’alimentation en eau dans l’état d’Uttarakhand, Inde, afin de capter l’eau de ruisseaux, de rivières et de lacs. A Nainital, sept puits tubés (profondeur de 22.6–36.7 m), situés à une distance < à 100 m du lac, sont utilisés pour prélever (1) de l’eau du lac après transit à travers le sol et (2) de l’eau de sub-surface/eau souterraine s’écoulant vers le lac. Des échantillons d’eau du lac et des puits tubés ont été analysés en périodes de mousson et de celles sans mousson 1997 à 2006. Résidu sec, C.E., alcalinité et dureté ont été trouvés marginalement supérieurs dans l’eau des puits tubés. La différence d’hydrochimie de l’eau des puits tubés était surtout due à la variation des régimes d’écoulement pendant les périodes de mousson et de celles sans mousson. Les résultats indiquent clairement que l’eau du lac en tant que telle n’est pas potable car elle contient des teneurs inacceptables de matière organique en termes de COD (~44 mg/L), de coliformes (~15.6 × 10⁴ MPN/100 mL) et d’éléments nutritifs. Des bactéries coliformes et du COD n’ont été détectés dans aucun des échantillons d’eau de puits tubés au fil des années. L’eau du lac traitée par des filtres à sable ne se conformait pas aux normes de l’eau potable. Ces recherches ont conduit à la fermeture de l’installation de traitement et à l’implantation de deux puits tubés en plus des cinq puits existants.

---

Resumen  En el estado de Uttarakhand, India, existen diferentes esquemas de abastecimiento de agua que explotan agua de arroyos, ríos y lagos. En Nainital, siete pozos (profundidades entre 22.6–36.7 m), ubicados a una distancia de <100 m del lago, se usan para extraer (1) agua del lago luego de su pasaje a través del suelo y (2) agua superficial y subterránea que fluye hacia el lago. En períodos de monzón y de no monzón de 1997 a 2006 se han analizado muestras de agua del lago y de las captaciones. Se halló que el agua de los pozos es ligeramente mayor en términos del total de sólidos disueltos, la conductividad eléctrica, la alcalinidad y la dureza. La diferencia en la hidroquímica del agua de las perforaciones se debe principalmente a la variación de los regímenes de flujo durante los períodos de monzón y de no monzón. Los resultados claramente indican que el agua del lago no es potable por su contenido inaceptable de materia orgánica medida como demanda de carbono orgánico/oxígeno –DCO- (~44 mg/L), coliformes (~15.6 × 10⁴ NMP/100 mL) y nutrientes. En el período, no se han detectado bacterias coliformes ni DCO en las muestras de agua de las captaciones. El agua del lago, tratada con filtros de arena, no conformó los estándares de agua para bebida. Estas investigaciones han demostrado la necesidad de clausurar las instalaciones de tratamiento y la adición de dos captaciones a las cinco ya existentes.

     

Hydrogeologic assessment of escalating groundwater exploitation in the Indus Basin, Pakistan

Groundwater development has contributed significantly to food security and reduction in poverty in Pakistan. Due to rapid population growth there has been a dramatic increase in the intensity of groundwater exploitation leading to declining water tables and deteriorating groundwater quality. In such prevailing conditions, the hydrogeological appraisal of escalating groundwater exploitation has become of paramount importance. Keeping this in view, a surface water–groundwater quantity and quality model was developed to assess future groundwater trends in the Rechna Doab (RD), a sub-catchment of the Indus River Basin. Scenario analysis shows that if dry conditions persist, there will be an overall decline in groundwater levels of around 10 m for the whole of RD during the next 25 years. The lower parts of RD with limited surface water supplies will undergo the highest decline in groundwater levels (10 to 20 m), which will make groundwater pumping very expensive for farmers. There is a high risk of groundwater salinization due to vertical upconing and lateral movement of highly saline groundwater into the fresh shallow aquifers in the upper parts of RD. If groundwater pumping is allowed to increase at the current rate, there will be an overall decline in groundwater salinity for the lower and middle parts of RD because of enhanced river leakage.     

Geochemistry of magmatic rocks of the Koksharovka alkaline-ultrabasic massif,primorye, and results of trace-element modeling

New geochemical data are presented on the magmatic rocks of the Late Jurassic Koksharovka alkaline-ultrabasic massif, which is associated with deposits of vermiculite, apatite, V-bearing titanomagnetite, and placer isoferroplatinum. The REE geochemistry and strontium, oxygen, and carbon isotopic composition of carbonatites and related ijolites and pyroxenites, together with geological observations, point to the magmatic origin of the Koksharovka carbonatites. The origin of associated magmatic rocks is discussed. Trace element modeling of partial melting of mantle sources was conducted to decipher the genesis of the melts of the Koksharovka carbonatites and host titanite-kaersutite pyroxenites.     

Developing a SPASE Query Language

The advent of the Virtual Observatory has begun an evolution in the space physics data environment. A number of nascent and discipline specific Virtual Observatories have started to emerge with an emphasis on data search and retrieval. As this new data environment takes shape an emphasis will be placed on interdisciplinary communication in attempts to address large scale and global problems. To this end we formulate the development of a query language to facilitate Virtual Observatory to Virtual Observatory communication. Furthermore, we outline the goals of such a language, how it would work and how existing community efforts can be leveraged to speed the development of this query language.

Spatial orientation of the spin vectors of the disk galaxies in the Virgo cluster I. A catalogue of the disk galaxies in the Virgo cluster

Based on four UK Schmidt Telescope (UKST) IIIa-j direct plates which had been digitized by the Automated Plate Measuring System (APM), a detailed identification of the disk galaxies in the Virgo cluster has been carried out. The certain and possible member galaxies are mainly selected from the updated Catalogue of 2096 Galaxies in the Virgo Cluster Area (VCC). The area of this identification covers a 85 square degree region centered on this cluster. As a result, we compiled a new catalogue of 340 certain and possible member disk galaxies, including 85 lenticulars (SO), 136 spirals (S) and 119 irregulars (Irr), where the APM-measured values of the position angle (P.A.) and diameters of 299 disk galaxies are provided. For at least 149 galaxies, the values of the P.A. and diameters are published for the first time.Comparison between the APM-measured P.A. and diameters of galaxies and those given by the Uppsala General Catalog of Galaxies (UGC) demonstrates the reliability of our data. The statistical properties and completeness of this database are also discussed. Compared with the UGC-given P.A. and diameters which were determined by eyes, our data has no selection effect owing to the instrument measurements. Within our searching region, this might be the largest and relatively complete photometric database published so far for investigating the spatial orientation of the spin vectors of galaxies in the Virgo cluster.     

Evolutionary status of the pre-protostellar core L1498

L1498 is a classic example of a dense cold pre-protostellar core. To study the evolutionary status, the structure, dynamics, and chemical properties of this core we have obtained high spatial and high spectral resolution observations of molecules tracing densities of 10(3)-10(5) cm-3. We observed CCS, NH3, C3H2, and HC7N with NASA's DSN 70 m antennas. We also present large-scale maps of C18O and 13CO observed with the AT&T 7 m antenna. For the high spatial resolution maps of selected regions within the core we used the VLA for CCS at 22 GHz, and the Owens Valley Radio Observatory (OVRO) MMA for CCS at 94 GHz and CS (2-1). The 22 GHz CCS emission marks a high-density [n(H2) > 10(4) cm -3] core, which is elongated with a major axis along the SE-NW direction. NH3 and C3H2 emissions are located inside the boundary of the CCS emission. C18O emission traces a lower density gas extending beyond the CCS boundary. Along the major axis of the dense core, CCS, NH3 and C3H2 emission show evidence of limb brightening. The observations are consistent with a chemically differentiated onion-shell structure for the L1498 core, with NH3 in the inner and CCS in the outer parts of the core. The high angular resolution (9"-12") spectral line maps obtained by combining NASA Goldstone 70 m and VLA data resolve the CCS 22 GHz emission in the southeast and northwest boundaries into arclike enhancements, supporting the picture that CCS emission originates in a shell outside the NH3 emitting region. Interferometric maps of CCS at 94 GHz and CS at 98 GHz show that their emitting regions contain several small-scale dense condensations. We suggest that the differences between the CCS, CS, C3H2, and NH3 emission are caused by a time-dependent effect as the core evolves slowly. We interpret the chemical and physical properties of L1498 in terms of a quasi-static (or slowly contracting) dense core in which the outer envelope is still growing. The growth rate of the core is determined by the density increase in the CCS shell resulting from the accretion of the outer low-density gas traced by C18O. We conclude that L1498 could become unstable to rapid collapse to form a protostar in less than 5 x 10(6) yr.     

Detection of a new interstellar molecular ion, H2COH+ (protonated formaldehyde)

A new interstellar molecular ion, H2COH+ (protonated formaldehyde), has been detected toward Sgr B2, Orion KL, W51, and possibly in NGC 7538 and DR21(OH). Six transitions were detected in Sgr B2(M). The 1(1,0)-1(0,1) transition was detected in all sources listed above. Searches were also made toward the cold, dark clouds TMC-1 and L134N, Orion (3N, 1E), and a red giant, IRC + 10216, without success. The excitation temperatures of H2COH+ are calculated to be 60-110 K, and the column densities are on the order of 10(12)-10(14) cm-2 in Sgr B2, Orion KL, and W51. The fractional abundance of H2COH+ is on the order of 10(-11) to 10-(9), and the ratio of H2COH+ to H2CO is in the range 0.001-0.5 in these objects. The values in Orion KL seem to be consistent with the "early time" values of recent model calculations by Lee, Bettens, & Herbst, but they appear to be higher than the model values in Sgr B2 and W51 even if we take the large uncertainties of column densities of H2CO into account. We suggest production routes starting from CH3OH may play an important role in the formation of H2COH+.     

Trace Element and Isotope Geochemistry of the Volcanic Rocks of Bequia, Grenadine Islands, Lesser Antilles Arc: a Study of Subduction Enrichment and Intra-crustal Contamination

Pliocene volcanics on the island of Bequia comprise two interbeddedsuites of basalts and andesites. The isotopically homogeneoussuite (IHS) has a limited range of Sr—Nd—Pb isotopes(⁸⁷Sr/⁸⁶Sr 0.7040–0.7046, ¹⁴³ Nd/¹⁴⁴ Nd 0.5130 and ²⁰⁶Pb/²⁰⁴Pb 19.36–19.51), and mantle-like ¹⁸O values (5.5in clinopyroxene). The isotopically diverse suite (IDS) is characterizedby much wider ranges of radiogenic isotopes (⁸⁷ Sr/⁸⁶Sr 0.7048–0.7077,¹⁴³ Nd/¹⁴⁴ Nd 0.5128–0.5123 and ²⁰⁶ Pb/²⁰⁴ Pb 19.7–20.2),in which all of the Sr and Pb ratios are higher and Nd ratiosare lower than those of the IHS. The IDS is also characterizedby high ¹⁸ O values, up to 7.6 in clinopyroxene. The Sr andPb isotope ratios are too high, and the Nd isotope ratios aretoo low in the IDS for any of these lavas to be derived fromunmodified depleted mantle. Both suites are petrologically very similar and their majorelement compositions and phenocryst contents suggest that theywere formed largely by fractional crystallization of a hydroustholeiitic melt at pressures <3 kbar. The isotopic ratiosand enrichments in large ion lithophile elements (LILE), andto some extent light rare earth elements (LREE), as comparedwith mid-ocean ridge basalts (MORB), of the IHS lavas suggestthat they were derived from a depleted mantle source which hadbeen re-enriched by the addition of 1–4% of a subductioncomponent. This component probably comprised a mixture of dehydrationfluids, and perhaps minor siliceous melts, released from subductingsediments and mafic crust. The extreme isotopic ranges, largeenrichments in incompatible elements, more fractionated LREEpatterns and higher ¹⁸ O values of the IDS lavas are interpretedas resulting from 10–55% assimilation—fractionalcrystallization of sediments, derived from the Guyana Shield,which are present in the arc crust, by IHS type melts. KEY WORDS: trace elements; radiogenic isotopes; arc lavas; Lesser Antilles ^*Corresponding author.     

Design of land-air parameterization scheme (LAPS) for modelling boundary layer surface processes

Summary A land-air parametrization scheme (LAPS) describes mass, energy and momentum transfer between the land surface and the atmosphere. The scheme is designed as a software package which can be run as part of an atmospheric model or a stand-alone scheme. A single layer approach is chosen for the physical and biophysical scheme background. The scheme has six prognostic variables: two temperatures (one for the canopy vegetation and one for soil surface), one interception storage, and three soil moisture storage variables. The scheme's upper boundary conditions are: air temperature, water vapour pressure, wind speed, radiation and precipitation at some reference level within the atmospheric boundary layer. The sensible and latent heat are calculated using resistance representation. The evaporation from the bare soil is parametrized using the scheme. The soil part is designed as a three-layer model which is used to describe the vertical transfer of water in the soil.The performances of the LAPS scheme were tested using the results of meteorological measurements over a maize field at the experimental site De Sinderhoeve (The Netherlands). The predicted partitioning of the absorbed radiation into sensible and latent heat fluxes is in good agreement with observations. Also, the predicted leaf temperature agrees quite well with the observed values.With 9 Figures