A ‘core-complex-like structure’ formed by superimposed extension,folding and high-angle normal faulting. The Santi Petri dome (western Betics,Spain)

The Santi Petri dome (western Betics, southern Spain) shows a core-complex-like structure, where migmatitic gneisses and schists outcrop below low-grade slates and phyllites, all of which form the basement of the Neogene Málaga basin. The migmatites and schists suffered a coaxial-flattening event during isothermal decompression and were later exhumed by ductile ESE non-coaxial stretching. Further exhumation was achieved by W- to SW-transport brittle low-angle normal faulting. Subsequently these extensional structures were gently folded in the core of a NE/SW-oriented antiform during the Tortonian. Finally the Santi Petri domal geometry was accentuated by the interference of orthogonal high-angle faults with ENE–WSW and NNW–SSE orientation. This core-complex-like structure, formed by superposition of extensional and compressive tectonic events, does not represent a classical, purely extensional core complex, which shows that metamorphic structure and geometry are not decisive criteria to define a core-complex.     

Spatio-temporal variability of benthic macroalgae in a coral reef system highly influenced by fluvial discharge: Veracruz,Gulf of Mexico

The Veracruz Reef System, in the southern Gulf of Mexico, is a suitable area for the study of the temporal and spatial variability of macroalgae abundance, at reef settings influenced by the fluvial discharge of the Jamapa River, and by human activities in the city and port of Veracruz. With this purpose, the bottom cover of each morpho-functional group of benthic macroalgae (frondoses, turf, and crustose corallines), and hermatypic corals, was determined at ten selected coral reefs, on a seasonal basis (rainy and dry seasons), for the 2009–2015 period. The average cover of benthic macroalgae was high (53.1%), with turf as the dominant morpho-functional group (31.9%), as in several reef ecosystems in the tropical Atlantic, followed by crustose corallines (15.2%), and frondoses (6.1%). Although turf macroalgae is dominant, due to their high temporal and spatial variability, the Veracruz Reef System could not be considered to be in a stable state, but just in an intermediate unstable equilibrium state, which is highly influenced by a high sediment load. As expected, nearshore reefs presented higher macroalgae covers, and unexpectedly, the outer-shelf reefs presented the highest cover of frondoses. Despite fluvial discharge influence, no differences in cover were found between the rainy and dry seasons. There was a negative and significantly correlation between the cover of frondoses and turf, which suggests that the driver/s of the abundance of these macroalgae, act in opposite ways for each group. Three clusters of reefs, defined by community structure and conservation degree, were determined: nearshore or degraded, offshore or moderately conserved, and conserved; and the entire Veracruz Reef System is considered to be in a moderately state of conservation.     

From platform to basin to swell: orbital control on sedimentary sequences in the Oxfordian, Spain

Climatic, oceanographic and ecological changes that control the formation and deposition of sediment in shallow and deep depositional environments commonly occur with periodicities of a few 10 000 years. Consequently, in order to interpret sedimentary sequences in the geological past, high time resolution is required. This is best obtained by cyclostratigraphy. Three sections have been studied in the Oxfordian of north-eastern Spain: one represents a shallow, siliciclastic-carbonate platform with repetitive subaerial exposures, one an intraplatform basin with sponge bioherms, and one a swell where iron ooids and glauconite formed. The platform section displays a well-defined stacking pattern of depositional sequences; the deeper-water sections are well dated by ammonites. The correlation between the three sections is a best-fit solution integrating biostratigraphy, sequence stratigraphy and cyclostratigraphy. It is concluded that the small-scale depositional sequences formed in tune with the 100-ka orbital eccentricity cycle. An additional factor was differential subsidence that ruled basin morphology.     

Smoothing unit vector fields

Data, each consisting of a unit vector and a position, are modeled as a smooth unit vector field plus random directional errors; the smooth unit vector field is estimated by fitting a spline to observations. The estimate is nonlinear; it is obtained with an iterative procedure in which, at each step, a least-squares smoothing spline is fitted to a set of pseudovalues. The resulting estimate is exemplified with a data set extracted from the literature in geology. The proposed estimate is compared, using the same data set with that obtained through suggestions of other authors, and is found to have attractive qualities. The use of plots of the resulting estimate to find the appropriate value of the smoothing parameter is emphasized throughout.     

Incorporating variable source area hydrology into a curve-number-based watershed model

Many water quality models use some form of the curve number (CN) equation developed by the Soil Conservation Service (SCS; U.S. Depart of Agriculture) to predict storm runoff from watersheds based on an infiltration-excess response to rainfall. However, in humid, well-vegetated areas with shallow soils, such as in the northeastern USA, the predominant runoff generating mechanism is saturation-excess on variable source areas (VSAs). We reconceptualized the SCS–CN equation for VSAs, and incorporated it into the General Watershed Loading Function (GWLF) model. The new version of GWLF, named the Variable Source Loading Function (VSLF) model, simulates the watershed runoff response to rainfall using the standard SCS–CN equation, but spatially distributes the runoff response according to a soil wetness index. We spatially validated VSLF runoff predictions and compared VSLF to GWLF for a subwatershed of the New York City Water Supply System. The spatial distribution of runoff from VSLF is more physically realistic than the estimates from GWLF. This has important consequences for water quality modeling, and for the use of models to evaluate and guide watershed management, because correctly predicting the coincidence of runoff generation and pollutant sources is critical to simulating non-point source (NPS) pollution transported by runoff. Copyright © 2007 John Wiley & Sons, Ltd.     

Diffuse Emission of Hydrogen from Poás Volcano,Costa Rica,America Central

We report the results of four soil H₂ surveys carried out in 2000–2003 at Poás volcano, Costa Rica, to investigate the soil H₂ distribution and evaluate the diffuse H₂ emission as a potential surveillance tool for Poás volcano. Soil gas H₂ contents showed a wide range of concentration from 0.2 to 7,059 ppmV during the four surveys. Maps of soil gas H₂ based on Sequential Gaussian Simulation showed low H₂ concentration values in the soil atmosphere (<0.7 ppmV) for most of the study area, whereas high soil gas H₂ values were observed inside the active crater of Poás. A significant increase in soil gas H₂ concentration was observed inside the active crater during 2001 and 2002 with respect to year 2000, followed by a decrease in 2003. The observed spatial and temporal variations of soil H₂ concentration have been well correlated with seismicity, microgravimetry and fumarolic chemistry changes which occurred during this study. These observations evidence changes in the shallow magmatic-hydrothermal system of Poás, and it might be related to a potential magmatic intrusion during the period 1998–2004. Therefore, monitoring diffuse H₂ emission of Poás has become an important geochemical tool for the monitoring of its volcanic activity.     

Numerical solution of the dam failure problem by the random choice method

A numerical solution of the dam failure problem as described by the one-dimensional shallow water equations is presented. The construction of the solution is based on the random choice method consisting in the superposition of locally theoretical solutions and sampling techniques. The search of the optimal sampling is performed through the application of the random choice method to the scalar wave equation. The dam failure problem is then solved and a comparison with the theoretical solution is presented. It is shown that the random choice method computes the bore with almost infinite resolution, represents exactly the constant state behind it and calculates the depression wave with great accuracy.     

Hysteretic freezing characteristics of riparian peatlands in the Western Boreal Forest of Canada

Freezing characteristics were investigated for a sedge covered floating fen and spruce covered swamp located beside a shallow lake in the Western Boreal Forest of Canada. Thermal properties were measured in situ for one freeze‐thaw cycle, and for two freeze‐thaw cycles in laboratory columns. Thermal conductivity and liquid water content were related to a range of subsurface temperatures above and below the freezing thresholds, and clearly illustrate hysteresis between the freezing and thawing process. Thermal hysteresis occurs because of the large change in thermal conductivity between water and ice, high water content of the peat, and wide variation in pore sizes that govern ice formation. Field and laboratory results were combined to develop linear freezing functions, which were tested in a heat transfer model. For surface temperature boundary conditions, subsurface temperatures were simulated for the over‐winter period and compared with field measurements. Replication of the transient subsurface thermal regime required that freezing functions transition gradually from thawed to frozen state (spanning the ?0·25 to ?2 °C range) as opposed to a more abrupt step function. Subsurface temperatures indicate that the floating fen underwent complete phase change (from water to ice) and froze to approximately the same depth as lake ice thickness. Therefore, the floating fen peatland froze as a ‘shelf’ adjacent to the lake, whereas the spruce covered swamp had a higher capacity for thermal buffering, and subsurface freezing was both more gradual and limited in depth. These thermal properties, and the timing and duration of frozen state, are expected to control the interaction of water and nutrients between surface water and groundwater, which will be affected by changes in air temperature associated with global climate change. Copyright © 2009 John Wiley & Sons, Ltd.     

Characterisation of conduction phenomena in soils at the particle-scale: Finite element analyses in conjunction with synthetic 3D imaging

This paper describes an innovative method to characterise conduction parameters in geomaterials at the particle-scale. The technique is exemplified using 3D synthetic grain packing generated with discrete element approaches. This creates a geo-mechanically viable user-defined 3D granular image through which the particle skeleton and the corresponding pore network are constructed. Images are then imported into the finite element analyses to solve the governing equations of hydraulic and thermal conduction. Navier–Stokes equation is uniquely upscaled to Darcy’s law to assess hydraulic conductivity in soils, while a similar approach implements the Fourier equation to evaluate thermal conduction through grain chains and pore network. High performance computing is used to meet demanding numerical calculations of 3D meshed geometries. Packing density (i.e., porosity) and inter-particle contact areas are explored as variables to highlight the effects of pore volume and inter-particle contact condition in hydraulic and thermal conduction. This emerging technique allows not only characterising the macro-scale behaviour of conduction phenomena in soils but also quantifying and visualising the preferential and local conduction behaviour at the particle-scale. Laboratory measurements of hydraulic and thermal conductivities support numerically obtained results and validate the viability of the new methods used herein. This study introduces an alternative way to determine physical parameters of soils using emerging technology of rigorous numerical simulations in conjunction with 3D images, and to enable fundamental observation of particle-scale mechanisms of macro-scale manifestation.