Delineation of geological facies from poorly differentiated data

The ability to delineate geologic facies and to estimate their properties from sparse data is essential for modeling physical and biochemical processes occurring in the subsurface. If such data are poorly differentiated, this challenging task is complicated further by the absence of a clear distinction between different hydrofacies at locations where data are   available. We consider three alternative approaches for analysis of poorly differentiated data: a k$k$-means clustering algorithm, an expectation–maximization algorithm, and a minimum-variance algorithm. Two distinct synthetically generated geological settings are used to analyze the ability of these algorithms to assign accurately the membership of such data in a given geologic facies. On average, the minimum-variance algorithm provides a more robust performance than its two counterparts, and when combined with a nearest neighbor algorithm, it also yields the most accurate reconstruction of the boundaries between the facies.     

Thermal conductive heating in fractured bedrock: Screening calculations to assess the effect of groundwater influx

A two-dimensional semi-analytical heat transfer solution is developed and a parameter sensitivity analysis performed to determine the relative importance of rock material properties (density, thermal conductivity and heat capacity) and hydrogeological properties (hydraulic gradient, fracture aperture, fracture spacing) on the ability to heat fractured rock using thermal conductive heating (TCH). The solution is developed using a Green’s function approach in which an integral equation is constructed for the temperature in the fracture. Subsurface temperature distributions are far more sensitive to hydrogeological properties than material properties. The bulk ground water influx (q) can provide a good estimate of the extent of influx cooling when influx is low to moderate, allowing the prediction of temperatures during heating without specific knowledge of the aperture and spacing of fractures. Target temperatures may not be reached or may be significantly delayed when the groundwater influx is large.     

A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.     

Effect of wave–bedform feedbacks on the formation of,and grain sorting over shoreface-connected sand ridges

The influence of wave–bedform feedbacks on both the initial formation of shoreface-connected sand ridges (sfcr) and on grain size sorting over these ridges on micro-tidal inner shelves is studied. Also, the effect of sediment sorting on the growth and the migration of sfcr is investigated. This is done by applying a linear stability analysis to an idealized process-based morphodynamic model, which simulates the initial growth of sfcr due to the positive coupling between waves, currents, and an erodible bed. The sediment consists of sand grains with two different sizes. New elements with respect to earlier studies on grain sorting over sfcr are that wave-topography interactions are explicitly accounted for, entrainment of sediment depends on bottom roughness, and transport of suspended sediment involves settling lag effects. The results of the model indicate that sediment sorting causes a reduction of the growth rate and migration speed of sfcr, whereas the wavelength is only slightly affected. In the case where the entrainment of suspended sediment depends on bottom roughness, the coarsest sediment is found in the troughs; otherwise, the finest sediment occurs in the troughs. Compared to previous work, modeled maximum variations in the mean grain size over the topography are in better agreement with field observations. Settling lag effects are important for the damping of high-wavenumber mode instabilities such that a preferred wavelength of the bedforms is obtained.     

Assessment of the hydraulic slope flow approach using a mesoscale model

The simplified hydraulic two-layer model for a katabatic flow is analysed using the outputs from a high-resolution mesoscale simulation. A stably stratified night is simulated for the Duero basin, a complex terrain area located in the northern Spanish plateau, with large vertical and horizontal spatial resolution. Well-defined katabatic flows on the basin slopes are generated by the simulation, that are relatively stationary and quasi-bidimensional for some areas in the central part of the night. The bulk quantities used in the two-layer approach as well as the different terms in the equations are computed from the three-dimensional information provided by the mesoscale simulation. This method allows to inspect how well the simplified approach represents the katabatic flow generated by the mesoscale model. The study shows that the hydraulic model allows for a comprehensive analysis of the basic mechanisms of the slope flows but is not able to close the budget equations, since the residuals are large.     

Stress sensitivity of stylolite morphology

Stylolites are rough surfaces that form by localized stress-induced dissolution. Using a set of limestone rock samples collected at different depths from a vertical section in Cirque de Navacelles (France), we study the influence of the lithostatic stress on the stylolites morphology on the basis of a recent morphogenesis model. We measured the roughness of a series of bedding-parallel stylolites and show that their morphology exhibits a scaling invariance with two self-affine scaling regimes separated by a crossover-length (L) at the millimeter scale consistent with previous studies. The importance of the present contribution is to estimate the stylolite formation stress σ from the sample position in the stratigraphic series and compare it to the crossover-length L using the expected relationship: L ～ σ ^?2. We obtained a successful prediction of the crossover behavior and reasonable absolute stress magnitude estimates using relevant parameters: depth of stylolite formation between 300 to 600 m with corresponding normal stress in the range of 10–18 MPa. Accordingly, the stylolite morphology contains a signature of the stress field during formation and we thus suggest that stylolites could be used as paleo-stress gauges of deformation processes in the upper crust.     

Short contribution: Buried Canopic channel identified near Egypt's Nile delta coast with radar (SRTM) imagery

Shuttle Radar Topography Mission (SRTM) data, with ground‐based geologic investigations, define the geographic position of the relict Canopic channel in Egypt's NW Nile delta. Two sinuous channel segments south of Abu Qir Bay are observed on a radar image: a more sinuous trace (36 km in length) west of Idku lagoon, and an eastern one (˜20 km), reaching the lagoon's southeast corner. Sediment cores recovered along channel traces show sand‐rich deposits (to > 4 m thick) beneath the surface, in contrast with silty muds in areas away from channels. Historic and archaeological findings indicate that the Ptolemaic Roman city of Schedia once lay directly along the Canopic that channeled water from the 3rd to 2nd centuries B.C. until ˜5th century A.D., after which Nile water was displaced to the east via Bolbitic and, later, Rosetta branches. Identification of buried, but well‐defined, Canopic channels provides baselines for renewed archaeological exploration in this delta sector. © 2006 Wiley Periodicals, Inc.     

Gas-Phase Formation of Water-Soluble Organic Compounds in the Atmosphere: A Retrosynthetic Analysis

Atmospheric particles contain a myriad of organiccompounds, including many multifunctional,water-soluble organic compounds. Many of thesecompounds are postulated to be secondary of origin. This work investigates the possible precursors ofseveral classes of multifunctional, water-solublesecondary organic compounds by analyzing the pathwaysleading to their formation, based on known gas-phasereactions. The analysis is termed `retrosynthetic'due to the backward direction of the analysis, fromproducts to precursors. Pathways for multi-functionalcompounds were generated combinatorially, consideringthe formation of one functional group at a time.Many multifunctional organic compounds with carboxylicacid, carbonyl, and hydroxy functional groups werefound to be first- or second-generation products ofcommon anthropogenic and biogenic volatile organiccompounds such as alkanes, alkenes, aromatics, andcyclic alkenes. The estimated yields of water-solubleorganic compounds from primary precursors ranged fromless than 1% to over 10%, based on stoichiometricconsiderations. The SOA formation index, whichcombines the concepts of yields and rates, was used tocompare the feasibility of the retrosyntheticpathways. Many of the candidate pathways involve theisomerization reaction of alkoxy radicals andoxygenated intermediate products such asmonocarboxylic acids and hydroxyaldehydes.     

Community forest management in Mexico: carbon mitigation and biodiversity conservation through rural development

Forest management is an important carbon mitigation strategy for developing countries. As demonstrated by the case of Mexico, community forest management is especially effective because it offers tangible local benefits while conserving forests and sequestering carbon. Community forestry receives minimal government support now, but the clean development mechanism (CDM) of the Kyoto Protocol could leverage additional resources to promote the approach in Mexico and elsewhere. We argue that adequately designed and implemented, community forestry management projects can avoid deforestation and restore forest cover and forest density. They comprise promising options for providing both carbon mitigation and sustainable rural development. These kinds of projects should be included in the CDM.