Vent temperature trends at the Vulcano Fossa fumarole field: the role of permeability

Between 1994 and 2010, we completed 16 thermal surveys of Vulcano’s Fossa fumarole field (Aeolian Islands, Italy). In each survey, between 400 and 1,200 vent temperatures were collected using a thermal infrared thermometer from distances of ～1?m. The results show a general decrease in average vent temperature during 1994–2003, with the average for the entire field falling from ～220°C in 1994 to ～150°C by 2003. However, between 2004 and 2010, we witnessed heating, with the average increasing to ～190°C by 2010. Alongside these annual-scale field-wide trends, we record a spatial re-organisation of the fumarole field, characterised by shut down of vent zones towards the crater floor, matched by rejuvenation of zones located towards the crater rim. Heating may be expected to be associated with deflation because increased amounts of vaporisation will remove volume from the hydrothermal system Gambino and Guglielmino (J Geophys Res 113:B07402, 2008). However, over the 2004–2010 heating period, no ground deformation was observed. Instead, the number of seismic events increased from a typical rate of 37 events per month during 1994–2000 to 195 events per month during 2004–2010. As part of this increase, we noticed a much greater number of high-frequency events associated with rock fracturing. We thus suggest that the heating event of 2004–2010 was the result of changed permeability conditions, rather than change in the heat supply from the deeper magmatic source. Within this scenario, cooling causes shut down of lower sectors and re-establishment of pathways located towards the crater rim, causing fracturing, increased seismicity and heat flow in these regions. This is consistent with the zone of rejuvenation (which lies towards and at the rim) being the most favourable location for fracturing given the stress field of the Fossa cone Sch?pa et al. (J Volcanol Geotherm Res 203:133–145, 2011); it is also the most established zone, having been active at least since the early twentieth century. Our data show the value of deploying multi-disciplinary geophysical campaigns at degassing (fumarolic) hydrothermal systems. This allows more complete and constrained understanding of the true heat loss dynamics of the system. In the case study presented here, it allows us to distinguish true heating from apparent heating phases. While the former are triggered from the bottom-up, i.e. they are driven by increases in heat supply from the magmatic source, the latter are triggered from the top-down, i.e. by changing permeability conditions in the uppermost portion of the system to allow more efficient heat flow over zones predisposed to fracturing.     

Glacial isostatic adjustment in 3-D earth models: Implications for the analysis of tide gauge records along the U.S. east coast

Tide gauge records of recent sea-level change along the U.S. east coast have received significant attention within the literature of glacial isostatic adjustment (GIA). Geographic trends in these tide gauge rates are not reduced by a GIA correction based on a commonly adopted radial viscosity profile (characterized, in particular, by a lower mantle viscosity 1−2×10²¹ Pa s), and this has led to speculation that the residual trends reflect contributions from neotectonics or oceanographic processes. While the trends can be significantly reduced by adopting an Earth model with a stiffer lower mantle, such a model appears to be incompatible with independent constraints from post-glacial decay times in Hudson Bay. We use a finite-element model of the GIA process to investigate whether 3-D viscosity variations superimposed onto the “common” radial viscosity profile may provide a route to reconciling the east coast sea-level trends. We find that the specific 3-D structure we impose has little impact on the geographic trends in the GIA-corrected rates. However, we do find that the imposed lateral variations in lower mantle viscosity introduce a nearly uniform upward shift of 0.5 mm/yr in GIA-induced sea-level rates along the U.S. east coast. Thus, inferences of regional (U.S. east coast) sea-level rise due to modern melting of ice reservoirs, based on tide gauge rates corrected using 1-D GIA models, may be significantly biased by this simplifying assumption.     

Simulating the response of a closed‐basin lake to recent climate changes in tropical West Africa (Lake Bosumtwi,Ghana)

Historical changes in the level of Lake Bosumtwi, Ghana, have been simulated using a catchment‐scale hydrological model in order to assess the importance of changes in climate and land use on lake water balance on a monthly basis for the period 1939–2004. Several commonly used models for computing evaporation in data‐sparse regions are compared, including the Penman, the energy budget, and the Priestley–Taylor methods. Based on a comparison with recorded lake level variations, the model with the energy‐budget evaporation model subcomponent is most effective at reproducing observed lake level variations using regional climate records. A sensitivity analysis using this model indicates that Lake Bosumtwi is highly sensitive to changes in precipitation, cloudiness and temperature. However, the model is also sensitive to changes in runoff related to vegetation, and this factor needs to be considered in simulating lake level variations. Both interannual and longer‐term changes in lake level over the last 65 years appear to have been caused primarily by changes in precipitation, though the model also suggests that the drop in lake level over the last few decades has been moderated by changes in cloudiness and temperature over that time. Based on its effectiveness at simulating the magnitude and rate of lake level response to changing climate over the historical record, this model offers a potential future opportunity to examine the palaeoclimatic factors causing past lake level fluctuations preserved in the geological record at Lake Bosumtwi. Copyright © 2006 John Wiley & Sons, Ltd.     

A synthesis of morphological evolutions and Holocene stratigraphy of a wave-dominated estuary: The Arcachon lagoon,SW France

This work presents the first synthesis of secular to millenary morphological evolutions and stratigraphy of a wave-dominated estuary, the Arcachon lagoon, from a combination of unpublished bathymetric maps (1865 and 2001), core results and high-resolution seismic profiles recorded for the first time in this lagoon. The Arcachon lagoon is located on the Atlantic coast of France, facing the wave-dominated shelf of the Bay of Biscay. It is a mesotidal semi-enclosed environment of about 160 km².The sediment budget of the Arcachon lagoon was computed by subtracting the 1865 bathymetric map from that of 2001. The computed volume difference is low (?9.9±35×10⁶ m³ in 136 yrs) and is the result of the balance between erosion and accretion that occurs within tidal channels and tidal flats, respectively. This morphological evolution pattern is explained by low sediment supply and also by the tidal distortion resulting from the morphology of the lagoon. Deep channels connected to the inlet are dominated by ebb currents inducing erosion. Tidal flats and transverse channels display weak or flood-dominated tidal currents leading to the deposition of silts. The areas of tidal flat siltation locally correlate with the presence of oyster farms, suggesting the influence of Man on the lagoon sediment-fill. Transverse channel-infill is related to weak tidal currents resulting from the hydraulically inefficient orientation of these channels which served as an ancient drainage network.Evidence for tidal channel-infill and channel abandonment are also provided by seismic profiling and cores. The upper stratigraphic succession of the lagoon (about 10 m thick) includes four main stratigraphic units dominated by channel-fills. The two lower units (around 7500–2800 yrs BP) display tabular-shape sandy channels interpreted to be records of the open estuarine phase of the Arcachon lagoon. The two upper units (around 2800 yrs BP to present-day) display U-shaped mixed sand-and-mud channel-fills interpreted to be records of the closure of the lagoon. Given that the basal estuarine units are transgressive and the upper lagoonal units are regressive, the main stratigraphic change at around 2800 yrs BP is interpreted as being the maximum flooding surface (MFS). This late MFS is explained by the low sediment supply. It is proposed that the transition from the estuarine to the lagoonal phase is related to the development of the Cap-Ferret spit in response to an increase in the ratio between wave power to tide power. This change in wave-to-tide ratio may be triggered by wave power increase following the Subboreal/Subatlantic climate instability or a decrease in tide power following a decrease in tidal prism related to the lagoon sediment-fill.Thus, the evolution of the Arcachon lagoon over the last millenaries was mainly controlled by its spit development, leading to a wave-dominated estuary in terms of its geomorphology. Once it was partially closed, extensive mud flats developed in the lagoon which became ebb-dominated.     

Magnetohydrodynamic simulations of gamma-ray burst jets: Beyond the progenitor star

Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet’s Lorentz factor $\gamma$ and opening angle ${\Theta }_j$ satisfies $\gamma {\Theta }_j?1$ at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give $\gamma {\Theta }_j?1$, suggesting that MHD models may be inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet–envelope interface as a collimating rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to ${10}^{4.5}$ times that of the central compact object, the jet achieves an asymptotic Lorentz factor $\gamma ～500$ far outside the star and an asymptotic opening angle ${\Theta }_j?0.04\text{rad}?2°$, giving $\gamma {\Theta }_j～20$. These values are consistent with observations of typical long-duration GRBs, and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well.     

Characterization of uranium-contaminated sediments from beneath a nuclear waste storage tank from Hanford, Washington: Implications for contaminant transport and fate

The concentration and distribution of uranium (U) in sediment samples from three boreholes recovered near radioactive waste storage tanks at Hanford, Washington, USA, were determined in detail using bulk and micro-analytical techniques. The source of contamination was a plume that contained an estimated 7000 kg of dissolved U that seeped into the subsurface as a result of an accident that occurred during filling of tank BX-102. The desorption character and kinetics of U were also determined by experiment in order to assess the mobility of U in the vadose zone. Most samples contained too little moisture to obtain quantitative information on pore water compositions. Concentrations of U (and contaminant phosphate—P) in pore waters were therefore estimated by performing 1:1 sediment-to-water extractions and the data indicated concentrations of these elements were above that of uncontaminated “background” sediments. Further extraction of U by 8 N nitric acid indicated that a significant fraction of the total U is relatively immobile and may be sequestered in mobilization-resistant phases. Fine- and coarse-grained samples in sharp contact with one another were sub-sampled for further scrutiny and identification of U reservoirs. Segregation of the samples into their constituent size fractions coupled with microwave-assisted digestion of bulk samples showed that most of the U contamination was sequestered within the fine-grained fraction. Isotope exchange (²³³U) tests revealed that ∼51% to 63% of the U is labile, indicating that the remaining fund of U is locked up in mobilization-resistant phases. Analysis by Micro-X-ray Fluorescence and Micro-X-ray Absorption Near-Edge Spectroscopy (μ-XRF and μ-XANES) showed that U is primarily associated with Ca and is predominately U(VI). The spectra obtained on U-enriched “hot spots” using Time-Resolved Laser-Induced Fluorescence Spectroscopy (TRLIFS) provide strong evidence for uranophane-type [Ca(UO₂)₂(SiO₃OH)₂(H₂O)₅] and uranyl phosphate [Ca(UO₂)₂(PO₄)₂(H₂O)_10-12] phases. These data show that disseminated micro-precipitates can form in narrow pore spaces within the finer-grained matrix and that these objects are likely not restricted to lithic fragment environments. Uranium mobility may therefore be curtailed by precipitation of uranyl silicate and phosphate phases, with additional possible influence exerted by capillary barriers. Consequently, equilibrium-based desorption models that predict the concentrations and mobility of U in the subsurface matrix at Hanford are unnecessarily conservative.     

ESR dating of tooth enamel in Mousterian layer 20, El Castillo,Spain

Nine faunal teeth from layer 20 of El Castillo cave in Cantabrian Spain were dated using electron spin resonance (ESR). Two teeth were rejected due to inconsistent subsample ages, while the remaining teeth yielded a mean age that is consistent with the stratigraphic expectations: 42.7±3.5. Uncertainty in the external γ dose rate results in a potential systematic uncertainty of±6.4ka that should affect all samples equally. The results provide independent confirmation of previously reported ¹⁴C ages for layer 20. © 2010 Wiley Periodicals, Inc.