Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation of the thermal regime and subsurface properties of a tidally affected,variably saturated streambed

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.     

Exponential time integrators for stochastic partial differential equations in 3D reservoir simulation

The transport of chemically reactive solutes (e.g. surfactants, CO₂ or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system, we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds of thousands unknowns.     

An ecosystem-based composite spatial model for floodplain vulnerability assessment: a case study of Artigas,Uruguay

Floods are natural processes that constitute a hazard to society when associated to improper land use. Anthropic activities in floodplains are a factor of vulnerability that converts a natural hazard into a threat factor, eventually leading to disaster. Nowadays, natural and social complex processes demand integrated assessments in order to improve their understanding, helping decision making over sustainable use of territory, as well as integrating society’s activity in ecosystems and potentials, restrictions and benefits that society obtain from them. In this context, the objective of this work was to build a composite vulnerability model for a floodplain under urban influence, using an integrated assessment approach. This model was based on three dimensions; threat, fragility and an ecosystem services provision. These dimensions were calculated using both primary and secondary information, and weights by specialists. Main results show that the area presents high vulnerability with an increasing gradient towards high and urbanized areas, associated with an important number and relevant ecosystem services. Also, a spatial heterogeneity of the three dimensions emerged, making evident this area’s complexity and the need of integrated assessments to approach it. The composite vulnerability model proposed presents an elevated potential for natural and social processes analysis in floodplains, which is crucial for these territory management. Moreover, these integrated dimensions could contribute to decision making in different levels, as well as generating important supplies for environmental management and land planning.

     

http://www.sciencedirect.com/science/article/pii/S1674987111000405

This study has focused on the processes of soil degradation and chemical element concentration in tea-growing regions of Rwanda, Africa. Soil degradation accelerated by erosion is caused not only by topography but also by human activities. This soil degradation involves both the physical loss and reduction in the amount of topsoil associated with nutrient decline. Soil samples were collected from eleven tropical zones in Rwanda and from variable depth within each collecting site. Of these, Samples from three locations in each zone were analyzed in the laboratory, with the result that the pH of all soil samples is shown to be less than 5 (pH < 5) with a general average of 4.4. The elements such as iron (Fe), copper (Cu), manganese (Mn), and zinc (Zn) are present in high concentration levels. In contrast calcium (Ca) and sodium (Na) are present at low-level concentrations and carbon (C) was found in minimal concentrations. In addition, elements derived from fertilizers, such as nitrogen (N), phosphorous (P), and potassium (K) which is also from minerals such as feldspar, are also present in low-level concentrations. The results indicate that the soil in certain Rwandan tea plantations is acidic and that this level of pH may help explain, in addition to natural factors, the deficiency of some elements such as Ca, Mg, P and N. The use of chemical fertilizers, land use system and the location of fields relative to household plots are also considered to help explain why tea plantation soils are typically degraded.     

A powder diffraction quantification method (ALJOR) for atmospheric particulate matter

The percentage calculation of crystalline phases in atmospheric particulate matter samples by means of an adiabatic approach of the matrix-flushing method with preferred preparations that avoids preferential orientations, should never be interpreted as absolute values. On the contrary, it should be seen as an indicator of the significant differences between each and every analyzed sample. Factors such as the application of constants, which are only useful for preparations of polycrystalline samples free of preferential orientation, and/or like the deficient resolution shown by the diffraction diagrams for certain phyllosilicate phases, bring a high degree of uncertainty to these percentage calculations. An alternative method of crystalline and amorphous phase quantification is presented below. Once the majority phases in the samples were identified, appropriate pure phases corresponding to each of them were selected. Then, the corresponding calibration curves were built using corundum (number 676 NIST) as a reference pattern phase. After validating the obtained calibration curves, the constants corresponding to each phase and selected spacing were calculated. The ALJOR method has been used for the quantification of five samples.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Investigation of a large collapse sinkhole affecting a multi-storey building by means of geophysics and the trenching technique (Zaragoza city,NE Spain)

An active sinkhole around 100 m long has been investigated in the city of Zaragoza (NE Spain). Subsidence activity on this depression, including the sudden occurrence of a collapse sinkhole 5 m across, led to the abandonment of a factory in the 1990s. At the present time, a building with 100 flats and shallow pad foundations partially built on the sinkhole, is affected by rapid differential settlement. The development of the sinkhole results from the karstification of the halite- and glauberite- bearing bedrock and the sagging and collapse of the overlying bedrock and alluvium, more than 30 m thick. GPR and electrical resistivity profiles have provided information on the distribution and geometry of the subsidence structure. The application of the trenching technique and geochronological methods (AMS and OSL dating) has allowed us to infer objective and practical data on the sinkhole including (1) Limits of the subsidence structure, (2) subsidence mechanisms, (3) cumulative subsidence (>408 cm), (4) subsidence rates on specific failure planes (>1.8 cm/year), (5) episodic displacement regime of some fault planes. The available information indicates that the progressive deformation recorded in the building will continue and might be punctuated by events of more rapid displacement. This work illustrates the practicality of the trenching technique for the study of sinkholes in mantled karst areas.