A locally conservative variational multiscale method for the simulation of porous media flow with multiscale source terms

We present a variational multiscale mixed finite element method for the solution of Darcy flow in porous media, in which both the permeability field and the source term display a multiscale character. The formulation is based on a multiscale split of the solution into coarse and subgrid scales. This decomposition is invoked in a variational setting that leads to a rigorous definition of a (global) coarse problem and a set of (local) subgrid problems. One of the key issues for the success of the method is the proper definition of the boundary conditions for the localization of the subgrid problems. We identify a weak compatibility condition that allows for subgrid communication across element interfaces, a feature that turns out to be essential for obtaining high-quality solutions. We also remove the singularities due to concentrated sources from the coarse-scale problem by introducing additional multiscale basis functions, based on a decomposition of fine-scale source terms into coarse and deviatoric components. The method is locally conservative and employs a low-order approximation of pressure and velocity at both scales. We illustrate the performance of the method on several synthetic cases and conclude that the method is able to capture the global and local flow patterns accurately.     

Quantitative reconstruction of Late Cenozoic landscapes: a case study in the Sierra de Atapuerca (Burgos,Spain)

We have developed a method to reconstruct palaeorelief by means of detailed geomorphological and geological studies, geostatistical tools, GIS and a DEM. This method has been applied to the Sierra de Atapuerca (NE Duero Basin, Burgos, Spain), allowing us to model a three‐dimensional reconstruction of the relief evolution from the Middle Miocene to the present. The modelling procedure is based on geostatistical recovery of the palaeosurfaces characteristic of each geomorphological evolution stage, using polynomial regressions, trend surfaces and kriging. The modelling of morphology trends has been useful in establishing new geological and geomorphological relationships in the geodynamic evolution of this basin, such as uplift quantification, correlation of erosion surfaces and sedimentary units, and the evolution of fluvial base levels. The palaeosurface reconstruction together with an analysis of the slope retreat have allowed us to reconstruct the palaeoreliefs that define the Late Cenozoic landscape evolution of this area, where the Lower and Middle Pleistocene archaeopalaeontological sites of the Sierra de Atapuerca are located. Copyright © 2007 John Wiley & Sons, Ltd.     

Origins and genesis of loess deposits in central Spain,as indicated by heavy mineral compositions and grain‐size variability

The formation of loess deposits strongly depends on the availability of deflatable dust‐sized material. Identification of source areas and mechanisms of dust production is essential for an appropriate characterization of related palaeoenvironmental conditions. So far, little research has been done on loess deposits in the Mediterranean region and information about mechanisms of dust production is very rare. In this context, it is not clear from where the loess deposits from the upper Tagus Basin in central Spain originated. The main objective of this study was to find out whether these loess deposits were formed primarily by mountain processes and fluvial comminution, or whether they have their origin in the weathering of Tertiary marls in the centre of the Madrid Basin. A further concern was to link the determined source areas with loess formation mechanisms and concomitant environmental contexts. Following a comprehensive approach using heavy minerals and grain‐size data, together with information on stratigraphic features and geomorphic positions, different local loess sources could be identified. This study shows that during the last 35 kyr a major proportion of the Tagus loess deposits was deflated from river floodplains, while the mountain region of the Iberian Range was identified as a significant source of respective floodplain sediments. Based on heavy mineral compositions and dating results, it was found that during Heinrich Stadial 3, sediment supply from the Iberian Range strongly increased, suggesting an environmental shift from initially warmer towards colder temperatures in late Marine Isotope Stage 3. Furthermore, it was found that grain‐size patterns clearly indicate maximum wind strengths during Heinrich Stadial 3, followed by Heinrich Stadial 2 and Heinrich Stadial 1. These results demonstrate that the formation of dust‐sized particles in this part of the Mediterranean operated similar to many other places in temperate and continental regions, and that weathering processes linked to hot and dry environments were only of minor importance.     

Paleomagnetism of early cretaceous arapey formation (Northern Uruguay)

We report detailed rock-magnetic and paleomagnetic directional data from 35 lava flows (302 standard paleomagnetic cores) sampled in the Central-Northern region of Uruguay in order to contribute to the study of the paleosecular variation of the Earth’s magnetic field during early Cretaceous and to obtain precise Cretaceous paleomagnetic pole positions for stable South America. The average unit direction is rather precisely determined from 29 out of 35 sites. All A₉₅ confidence angles are less than 8°, which points to small within-site dispersion and high directional stability. Normal polarity magnetizations are revealed for 19 sites and 10 are reversely magnetized. Two other sites yield well defined intermediate polarities. The mean direction, supported by a positive reversal test is in reasonably good agreement with the expected paleodirection for Early Cretaceous stable South America and in disagreement with a 10° clockwise rotation found in the previous studies. On the other hand, paleomagnetic poles are significantly different from the pole position suggested by hotspot reconstructions, which may be due to true polar wander or the hotspot motion. Our data suggest a different style of secular variation during (and just before) the Cretaceous Normal Superchron and the last 5 Ma, supporting a link between paleosecular variation and reversal frequency.     

Deriving Wind Force Terms from Nautical Reports Through Content Analysis. The Spanish and French Cases

Records of wind strength taken onboard Spanish and French ships during the 1750–1850 period have been digitized and examined using content analysis techniques to derive the equivalent wind strength in terms of the current Beaufort scale, this conversion being a key step in any attempt to compare ancient records with modern climatologies. During the analysis it was evident that Spanish and French officers used a great number of different terms to describe the wind force. However, when the records are analyzed and homogenized, a broadly common and essentially narrow vocabulary was identified, indicating that, at this period, a great effort had been made to regulate the maritime language. Using contemporary dictionaries and navigation manuals, an equivalence was established between the original wind force terms and the Beaufort scale.     

Origin of the Rubian carbonate-hosted magnesite deposit, Galicia, NW Spain: mineralogical, REE, fluid inclusion and isotope evidence

The Rubian magnesite deposit (West Asturian—Leonese Zone, Iberian Variscan belt) is hosted by a 100-m-thick folded and metamorphosed Lower Cambrian carbonate/siliciclastic metasedimentary sequence—the Cándana Limestone Formation. It comprises upper (20-m thickness) and lower (17-m thickness) lens-shaped ore bodies separated by 55 m of slates and micaceous schists. The main (lower) magnesite ore body comprises a package of magnesite beds with dolomite-rich intercalations, sandwiched between slates and micaceous schists. In the upper ore body, the magnesite beds are thinner (centimetre scale mainly) and occur between slate beds. Mafic dolerite dykes intrude the mineralisation. The mineralisation passes eastwards into sequence of bedded dolostone (Buxan) and laminated to banded calcitic marble (Mao). These show significant Variscan extensional shearing or fold-related deformation, whereas neither Rubian dolomite nor magnesite show evidence of tectonic disturbance. This suggests that the dolomitisation and magnesite formation postdate the main Variscan deformation. In addition, the morphology of magnesite crystals and primary fluid inclusions indicate that magnesite is a neoformed hydrothermal mineral. Magnesite contains irregularly distributed dolomite inclusions (<50 μm) and these are interpreted as relics of a metasomatically replaced dolostone precursor. The total rare earth element (REE) contents of magnesite are very similar to those of Buxan dolostone but are depleted in light rare earth elements (LREE); heavy rare earth element concentrations are comparable. However, magnesite REE chondrite normalised profiles lack any characteristic anomaly indicative of marine environment. Compared with Mao calcite, magnesite is distinct in terms of both REE concentrations and patterns. Fluid inclusion studies show that the mineralising fluids were MgCl₂–NaCl–CaCl₂–H₂O aqueous brines exhibiting highly variable salinities (3.3 to 29.5 wt.% salts). This may be the result of a combination of fluid mixing, migration of pulses of variable-salinity brines and/or local dissolution and replacement processes of the host dolostone. Fluid inclusion data and comparison with other N Iberian dolostone-hosted metasomatic deposits suggest that Rubian magnesite probably formed at temperatures between 160 and 200°C. This corresponds, at hydrostatic pressure (500 bar), to a depth of formation of ~~5 km. Mineralisation-related Rubian dolomite yields δ ¹⁸O values (δ ¹⁸O: 12.0–15.4‰, mean: 14.4±1.1‰) depleted by around 5‰ compared with barren Buxan dolomite (δ ¹⁸O: 17.1–20.2‰, mean: 19.4±1.0‰). This was interpreted to reflect an influx of ¹⁸O-depleted waters accompanied by a temperature increase in a fluid-dominated system. Overlapping calculated δ ¹⁸O_fluid values (~+5‰ at 200°C) for fluids in equilibrium with Rubian dolomite and magnesite show that they were formed by the same hydrothermal system at different temperatures. In terms of δ ¹³C values, Rubian dolomite (δ ¹³C: −1.4 to 1.9‰, mean: 0.4±1.3‰) and magnesite (δ ¹³C: −2.3 to 2.4‰, mean: 0.60±1.0‰) generally exhibit more negative δ ¹³C values compared with Buxan dolomite (δ ¹³C: −0.2 to 1.9‰, mean: 0.8±0.6‰) and Mao calcite (δ ¹³C: −0.3 to 1.5‰, mean: 0.6±0.6‰), indicating progressive modification to lower δ ¹³C values through interaction with hydrothermal fluids. ⁸⁷Sr/⁸⁶Sr ratios, calculated at 290 Ma, vary from 0.70849 to 0.70976 for the Mao calcite and from 0.70538 to 0.70880 for the Buxan dolostone. The ⁸⁷Sr/⁸⁶Sr ratios in Rubian magnesite are more radiogenic and range from 0.71123 to 0.71494. The combined δ ¹⁸O–δ ¹³C and ⁸⁷Sr/⁸⁶Sr data indicate that the magnesite-related fluids were modified basinal brines that have reacted and equilibrated with intercalated siliciclastic rocks. Magnesite formation is genetically linked to regional hydrothermal dolomitisation associated with lithospheric delamination, late-Variscan high heat flow and extensional tectonics in the NW Iberian Belt. A comparison with genetic models for the Puebla de Lillo talc deposits suggests that the formation of hydrothermal replacive magnesite at Rubian resulted from a metasomatic column with magnesite forming at higher fluid/rock ratios than dolomite. In this study, magnesite generation took place via the local reaction of hydrothermal dolostone with the same hydrothermal fluids in very high permeability zones at high fluid/rock ratios (e.g. faults). It was also possibly aided by additional heat from intrusive dykes or sub-cropping igneous bodies. This would locally raise isotherms enabling a transition from the dolomite stability field to that of magnesite.Editorial handling: F. Tornos     

Post-processing HAZUS earthquake damage and loss assessments for individual buildings

Natural Hazards - Building damage probabilities are invaluable for assessing short-term losses from natural hazards. In many countries however, the individual building level data required for...     

The onset of the Messinian salinity crisis in the deep Eastern Mediterranean basin

Astronomical tuning of the Messinian pre‐salt succession in the Levant Basin allows for the first time the reconstruction of a detailed chronology of the Messinian salinity crisis (MSC) events in deep setting and their correlation with marginal records that supports the CIESM ( 2008 ) 3‐stage model. Our main conclusions are (1) MSC events were synchronous across marginal and deep basins, (2) MSC onset in deep basins occurred at 5.97 Ma, (3) only foraminifera‐barren, evaporite‐free shales accumulated in deep settings between 5.97 and 5.60 Ma, (4) deep evaporites (anhydrite and halite) deposition started later, at 5.60 Ma and (5) new and published ⁸⁷Sr/⁸⁶Sr data indicate that during all stages, evaporites precipitated from the same water body in all the Mediterranean sub‐basins. The wide synchrony of events and ⁸⁷Sr/⁸⁶Sr homogeneity implies inter‐sub‐basin connection during the whole MSC and is not compatible with large sea‐level fall and desiccation of the Mediterranean.     

Extra-tropical atmospheric response to ENSO in the CMIP5 models

The seasonal mean extra-tropical atmospheric response to El Niño/Southern Oscillation (ENSO) is assessed in the historical and pre-industrial control CMIP5 simulations. This analysis considers two types of El Niño events, characterized by positive sea surface temperature (SST) anomalies in either the central equatorial Pacific (CP) or eastern equatorial Pacific (EP), as well as EP and CP La Niña events, characterized by negative SST anomalies in the same two regions. Seasonal mean geopotential height anomalies in key regions typify the magnitude and structure of the disruption of the Walker circulation cell in the tropical Pacific, upper tropospheric ENSO teleconnections and the polar stratospheric response. In the CMIP5 ensembles, the magnitude of the Walker cell disruption is correlated with the strength of the mid-latitude responses in the upper troposphere i.e., the North Pacific and South Pacific lows strengthen during El Niño events. The simulated responses to El Niño and La Niña have opposite sign. The seasonal mean extra-tropical, upper tropospheric responses to EP and CP events are indistinguishable. The ENSO responses in the MERRA reanalysis lie within the model scatter of the historical simulations. Similar responses are simulated in the pre-industrial and historical CMIP5 simulations. Overall, there is a weak correlation between the strength of the tropical response to ENSO and the strength of the polar stratospheric response. ENSO-related polar stratospheric variability is best simulated in the “high-top” subset of models with a well-resolved stratosphere.