Trends in middle- and upper-level tropospheric humidity from NCEP reanalysis data

The National Centers for Environmental Prediction (NCEP) reanalysis data on tropospheric humidity are examined for the period 1973 to 2007. It is accepted that radiosonde-derived humidity data must be treated with great caution, particularly at altitudes above the 500 hPa pressure level. With that caveat, the face-value 35-year trend in zonal-average annual-average specific humidity q is significantly negative at all altitudes above 850 hPa (roughly the top of the convective boundary layer) in the tropics and southern midlatitudes and at altitudes above 600 hPa in the northern midlatitudes. It is significantly positive below 850 hPa in all three zones, as might be expected in a mixed layer with rising temperatures over a moist surface. The results are qualitatively consistent with trends in NCEP atmospheric temperatures (which must also be treated with great caution) that show an increase in the stability of the convective boundary layer as the global temperature has risen over the period. The upper-level negative trends in q are inconsistent with climate-model calculations and are largely (but not completely) inconsistent with satellite data. Water vapor feedback in climate models is positive mainly because of their roughly constant relative humidity (i.e., increasing q) in the mid-to-upper troposphere as the planet warms. Negative trends in q as found in the NCEP data would imply that long-term water vapor feedback is negative—that it would reduce rather than amplify the response of the climate system to external forcing such as that from increasing atmospheric CO₂. In this context, it is important to establish what (if any) aspects of the observed trends survive detailed examination of the impact of past changes of radiosonde instrumentation and protocol within the various international networks.     

Three-dimensional finite element modelling of stress relaxation tests in anisotropic clayey medium: direct problem and back analysis

This paper deals with a new strategy for initial stress identification by stress relaxation methods, coupled with finite element calculation, and applied to the overcoring test. The back-analysis of such a test uses an inversion method which consists in the minimisation, with a gradient-based algorithm, of a cost functional of least-squares type, which quantizes the difference between measured and computed strains. The computed strains are assessed by three-dimensional finite element modellings of the overcoring test. The inversion methodology is applied to a recent in situ overcoring test performed at Mont Terri laboratory, Switzerland. The inversion gave good results and allows us to validate the inversion methodology. The constitutive law considered for this application is transverse isotropic elasticity but the inversion method developed is applicable to most constitutive laws and every kind of in situ test.     

Evaluation of Methane Sources in Groundwater in Northeastern Pennsylvania

Testing of 1701 water wells in northeastern Pennsylvania shows that methane is ubiquitous in groundwater, with higher concentrations observed in valleys vs. upland areas and in association with calcium‐sodium‐bicarbonate, sodium‐bicarbonate, and sodium‐chloride rich waters—indicating that, on a regional scale, methane concentrations are best correlated to topographic and hydrogeologic features, rather than shale‐gas extraction. In addition, our assessment of isotopic and molecular analyses of hydrocarbon gases in the Dimock Township suggest that gases present in local water wells are most consistent with Middle and Upper Devonian gases sampled in the annular spaces of local gas wells, as opposed to Marcellus Production gas. Combined, these findings suggest that the methane concentrations in Susquehanna County water wells can be explained without the migration of Marcellus shale gas through fractures, an observation that has important implications for understanding the nature of risks associated with shale‐gas extraction.     

Secondary turbidity maximum in a partially mixed microtidal estuary

Data from a two-year period of monthly slackwater surveys reveal that in addition to the classical estuary turbidity maximum (ETM), another peak of bottom total suspended sediment (TSS) concentration, or a so-called secondary turbidity maximum (STM), often exists in the middle part of the York River estuary, Virginia. This STM, observed in most (but not all) of the slackwater surveys, moves back and forth in the region of about 20 to 40 km from the York River mouth where the mud percentage of bottom sediment is very high. The distribution of the potential energy anomaly, which was calculated using salinity data, indicates that the STM usually resides in the transition zone between the upstream well mixed and the downstream more stratified water columns. An analysis using the conservation equation of suspended sediment concentration in the water column reveals that four processes may contribute to the formation of the STM: convergence of bottom residual flow, tidal asymmetry, inhibition of turbulent diffusion by stratification, and bottom resuspension. The along-channel variations of the strength of bottom residual flow, the effect of tidal asymmetry, and the stratification patterns are probably due to the geometric features of the York River estuary.     

Assimilation-fractional crystallization of Polvadera Group rocks in the northwestern Jemez Volcanic Field,New Mexico

Pliocene Polvadera Group rocks in the northwestern Rio Grande rift-marginal portion of the Jemez Volcanic Field record the rapid transition from weakly alkaline Lobato Basalt magmatism (48–52% SiO₂; 7.9 Ma) through calc-alkaline Lobato andesite and dacite (53–64% SiO₂) and Tschicoma dacite-rhyodacite magmatism (63–69% SiO₂; 7.4 Ma). Petrologically, Lobato andesite and dacite and Tschicoma dacite-rhyodacite represent a cogenetic suite of differentiates (the La Grulla Plateau or LGP suite) distinctive from the bulk of Polvadera Group rocks including Tschicoma andesite. Increasing (⁸⁷Sr/⁸⁶Sr)_O ratios with differentiation within the LGP suite from 0.7051 (54% SiO₂) to 0.7064 (68% SiO₂), trace element variations, and disequilibrium mineral assemblages suggest open system differentiation involving ⁸⁷Sr-enriched upper crust. A likely parental magma is the voluminous Lobato Basalt ((⁸⁷Sr/⁸⁶)_O= 0.7043–0.7050) which was erupted predominantly earlier and to the east toward the rift axis. The best model for petrogenesis involves bulk assimilation of locally wide-spread Proterozoic (1.4–1.6 Ga) upper crustal granite by fractionally crystallizing Lobato Basalt. Assimilation-fractional crystallization (AFC) modeling of Sr-isotope and trace element variation (DePaolo 1981) indicates that 40% crystallization of Lobato Basalt accompanied by 10% addition of granite reproduces the observed geochemical and Sr-isotopic diversity. Neither magma mixing, nor mafic recharge have complicated the magmatic evolution of the LGP suite. Crustal thickness and/or retarded tectonism could have facilitated conditions necessary for evolution by AFC to occur within the upper crust.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

FRACTIONAL TRANSPORT OF SEDIMENT MIXTURES

A new method based on the Transport Capacity Fraction (TCF) concept is proposed to compute the fractional transport rates for nonuniform sediment mixtures in sand-bed channels. The TCF concept is derived from the understanding that the measurements and predictions of bed-material load are more accurate and reliable than the measurements and predictions of fractional loads. First the bed-material load is computed using an appropriate equation, then the fractional transport rates are determined by distributing the bed-material load into size groups through a transport capacity distribution function. For the computation of bed-material loads, the Ackers and White, Engelund and Hansen, and Yang equations are used in this study. Two new transport capacity distribution functions are developed for flows in sand-bed channels. The new expressions presented in this paper account for the sheltering and exposure effects that exist in mixtures. Comparisons with measured data show that the proposed method can signifi     

Conduction models of the temperature distribution in the East Rift Zone of Kilauea Volcano

Temperature variations in the 1966-meter Hawaii Geothermal Project well HGP-A are simulated by model studies using a finite element code for conductive heat flow. Three models were generated: a constant temperature source from a vertical dike; a constant heat-generating magma chamber; and a transient heat source from a tapered vertical dike. Fair correlation is obtained between the HGP-A well temperature and the tapered dike 125 years after it is injected with an initial (transient) 1200° C temperature. Results provide background information from which to evaluate the importance of water convection in maintaining the temperature distribution in the East Rift Zone.     

Characterizing European urban settlements from population data: a cartographic approach

ABSTRACT

The cartographic identification and characterization of urban settlements is problematic, particularly in varied areas or in large geographical extensions. This is due to the divergence between morphological area structures and data attribution, which is generally available at the level of administrative units. Most existing experiences focus on small regions or use indirect data, from which only indirect information can be derived. This paper proposes an identification of urban settlements in the European Union space, which combines urban population and shape through geographic information system (GIS) techniques. For that purpose, 1 km² cells with associated population data are used. A procedure is developed establishing three different types of urban settlements, according to population and population density values: high-density urban clusters, small and medium-sized towns, and very small towns. A validation of this procedure is performed using direct survey from local experts in 10 selected study areas. The final results show that 78,307 urban settlements can be identified: 88.17% of them are classified as “very small towns,” 10.75% as “small and medium-sized towns,” and 1.08% as high-density urban clusters. This exercise generates new information for the European Union and surrounding countries, and yields the basis for further research. All GIS operations are performed using vector layers, which is not usual in large-area regional studies.