Properties of natural supraglacial debris in relation to modelling sub‐debris ice ablation

As debris‐covered glaciers become a more prominent feature of a shrinking mountain cryosphere, there is increasing need to successfully model the surface energy and mass balance of debris‐covered glaciers, yet measurements of the processes operating in natural supraglacial debris covers are sparse. We report measurements of vertical temperature profiles in debris on the Ngozumpa glacier in Nepal, that show: (i) conductive processes dominate during the ablation season in matrix‐supported diamict; (ii) ventilation may be possible in coarse surface layers; (iii) phase changes associated with seasonal change have a marked effect on the effective thermal diffusivity of the debris. Effective thermal conductivity determined from vertical temperature profiles in the debris is generally ~30% higher in summer than in winter, but values depend on the volume and phase of water in the debris. Surface albedo can vary widely over small spatial scales, as does the debris thickness. Measurements indicate that debris thickness is best represented as a probability density function with the peak debris thickness increasing down‐glacier. The findings from Ngozumpa glacier indicate that the probability distribution of debris thickness changes from positively skewed in the upper glacier towards a more normal distribution nearer the terminus. Although many of these effects remain to be quantified, our observations highlight aspects of spatial and temporal variability in supraglacial debris that may require consideration in annual or multi‐annual distributed modelling of debris‐covered glacier surface energy and mass balance. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of drought over the past two millennia using near-shore sediment cores from a Canadian boreal lake

Lack of long-term studies on drought in the boreal region of northwest Ontario limits our ability to assess the vulnerability of this region to climate change. We investigated changes in diatoms, scaled chrysophytes, and sedimentary pigments in two near-shore cores from Gall Lake to infer limnological and water-level changes over the last two millennia. The two coring locations, at modern water depths of 7.5 and 11.5 m, were selected to contrast inferences for past lake level based on distance from the modern water-depth boundary between predominantly benthic and planktonic diatom assemblages in surface sediments (6.0 m). Diatom-inferred depth inferences were more variable in the core from 7.5-m water depth, consistent with the hypothesis that the greatest changes occurred at water depths closest to the benthic:planktonic boundary. Both sites revealed a pronounced drought from ~AD 950 to 1300, synchronous with the medieval climate anomaly (MCA). This finding suggests a northeast expansion of the arid MCA into northwest Ontario, extending the known spatial extent of this megadrought. Scaled chrysophytes and sedimentary pigments suggest a recent increase in thermal stratification. Our findings indicate this region is more susceptible to changes in moisture than was previously suspected based on the instrumental record for the past century.     

Sixty thousand years of magmatic volatile history before the caldera-forming eruption of Mount Mazama, Crater Lake, Oregon

The well-documented eruptive history of Mount Mazama, Oregon, provides an excellent opportunity to use pre-eruptive volatile concentrations to study the growth of an explosive silicic magmatic system. Melt inclusions (MI) hosted in pyroxene and plagioclase crystals from eight dacitic–rhyodacitic eruptive deposits (71–7.7?ka) were analyzed to determine variations in volatile-element concentrations and changes in magma storage conditions leading up to and including the climactic eruption of Crater Lake caldera. Temperatures (Fe–Ti oxides) increased through the series of dacites, then decreased, and increased again through the rhyodacites (918–968 to ~950 to 845–895?°C). Oxygen fugacity began at nickel–nickel-oxide buffer (NNO) +0.8 (71?ka), dropped slightly to NNO +0.3, and then climbed to its highest value with the climactic eruption (7.7?ka) at NNO +1.1 log units. In parallel with oxidation state, maximum MI sulfur concentrations were high early in the eruptive sequence (~500?ppm), decreased (to ~200?ppm), and then increased again with the climactic eruption (~500?ppm). Maximum MI sulfur correlates with the Sr content (as a proxy for LREE, Ba, Rb, P₂O₅) of recharge magmas, represented by basaltic andesitic to andesitic enclaves and similar-aged lavas. These results suggest that oxidized Sr-rich recharge magmas dominated early and late in the development of the pre-climactic dacite–rhyodacite system. Dissolved H₂O concentrations in MI do not, however, correlate with these changes in dominant recharge magma, instead recording vapor solubility relations in the developing shallow magma storage and conduit region. Dissolved H₂O concentrations form two populations through time: the first at 3–4.6 wt% (with a few extreme values up to 6.1 wt%) and the second at ≤2.4 wt%. CO₂ concentrations measured in a subset of these inclusions reach up to 240?ppm in early-erupted deposits (71?ka) and are below detection in climactic deposits (7.7?ka). Combined H₂O and CO₂ concentrations and solubility models indicate a dominant storage region at 4–7?km (up to 12?km), with drier inclusions that diffusively re-equilibrated and/or were trapped at shallower depths. Boron and Cl (except in the climactic deposit) largely remained in the melt, suggesting vapor–melt partition coefficients and gas fractions were low. Modeled Li, F, and S vapor–melt partition coefficients are higher than those of B and Cl. The decrease in maximum MI CO₂ concentration following the earliest dacitic eruptions is interpreted to result from a broadening of the shallow storage region to greater than the diameter of subjacent feeders, so that greater proportions of reservoir magma were to the side of CO₂-bearing vapor bubbles ascending vertically from the locus of recharge magma injection, thereby escaping recarbonation by streaming vapor bubbles. The Mazama melt inclusions provide a picture of a growing magma storage region, where chemical variations in melt and magma occur due to changes in the nature and supply rate of magma recharge, the timing of degassing, and the possible degree of equilibration with gases from below.     

Seismic interferometry and ambient noise tomography in the British Isles

Traditional methods of imaging the Earth's subsurface using seismic waves require an identifiable, impulsive source of seismic energy, for example an earthquake or explosive source. Naturally occurring, ambient seismic waves form an ever-present source of energy that is conventionally regarded as unusable since it is not impulsive. As such it is generally removed from seismic data and subsequent analysis. A new method known as seismic interferometry can be used to extract useful information about the Earth's subsurface from the ambient noise wavefield. Consequently, seismic interferometry is an important new tool for exploring areas which are otherwise seismically quiescent, such as the British Isles in which there are relatively few strong earthquakes. One of the possible applications of seismic interferometry is ambient noise tomography (ANT). ANT is a way of using interferometry to image subsurface seismic velocity variations using seismic (surface) waves extracted from the background ambient vibrations of the Earth. To date, ANT has been used successfully to image the Earth's crust and upper-mantle on regional and continental scales in many locations and has the power to resolve major geological features such as sedimentary basins and igneous and metamorphic cores. Here we provide a review of seismic interferometry and ANT, and show that the seismic interferometry method works well within the British Isles. We illustrate the usefulness of the method in seismically quiescent areas by presenting the first surface wave group velocity maps of the Scottish Highlands using only ambient seismic noise. These maps show low velocity anomalies in sedimentary basins such as the Moray Firth, and high velocity anomalies in igneous and metamorphic centres such as the Lewisian complex. They also suggest that the Moho shallows from south to north across Scotland which agrees with previous geophysical studies in the region.     

The potential for convection and implications for geothermal energy in the Perth Basin, Western Australia

Convection of groundwater in aquifers can create areas of anomalously high temperature at shallow depths which could be exploited for geothermal energy. Temperature measurements in the Perth Basin (Western Australia) reveal thermal patterns that are consistent with convection in the Yarragadee Aquifer. This observation is supported by Rayleigh number calculations, which show that convection is possible within the range of aquifer thickness, geothermal gradient, salinity gradient and permeability encountered in the Yarragadee Aquifer, assuming that the aquifer can be treated as a homogeneous anisotropic layer. Numerical simulations of convection in a simplified model of the Yarragadee Aquifer show that: (1) the spacing of convective upwellings can be predicted from aquifer thickness and permeability anisotropy; (2) convective upwellings may be circular or elongate in plan view; (3) convective upwellings create significant temperature enhancements relative to the conductive profile; (4) convective flow rates are similar to regional groundwater flow rates; and (5) convection homogenises salinity within the aquifer. Further work is required to constrain the average horizontal and vertical permeability of the Yarragadee Aquifer, to assess the validity of treating the aquifer as a homogeneous anisotropic layer, and to determine the impact of realistic aquifer geometry and advection on convection.     

Polycyclic aromatic hydrocarbons in commercial fish and lobsters from the coastal waters of Madagascar following an oil spill in August 2009

Concentrations of polycyclic aromatic hydrocarbons were determined in species of commercial fish and lobsters following an oil-spill just off the protected Madagascan coastline. Samples were collected along the coastline within and outside the affected area. Summed PAH concentrations ranged from 1.9 μg kg⁻¹ to 63 μg kg⁻¹ wet weight, but with no higher molecular weight PAHs (>202 Da) being detected. All concentrations of benzo[a]pyrene, benz[a]anthracene and dibenz[a,h]anthracene were <0.1 μg kg⁻¹ wet weight, well within the EU and UK set limits for the protection of human health. Additionally, samples were calculated as the benzo[a]pyrene toxic equivalency quotient (TEQ) and found to be well below the level of concern in relation to health of human consumers. Evaluation of the biota PAH data indicated the origin of PAH was predominantly petrogenic with >80% arising from oil sources. Profile studies indicate a low-level multisource petrogenic contamination probably representing a pre-spill background for the area.     

High-temperature emplacement of the Cerro Galán and Toconquis Group ignimbrites (Puna plateau,NW Argentina) determined by TRM analyses

Estimates of pyroclastic flow emplacement temperatures in the Cerro Galán ignimbrite and Toconquis Group ignimbrites were determined using thermal remanent magnetization of lithic clasts embedded within the deposits. These ignimbrites belong to the Cerro Galán volcanic system, one of the largest calderas in the world, in the Puna plateau, NW Argentina. Temperature estimates for the 2.08-Ma Cerro Galán ignimbrite are retrieved from 40 sites in 14 localities (176 measured clasts), distributed at different distances from the caldera and different stratigraphic heights. Additionally, temperature estimates were obtained from 27 sample sites (125 measured clasts) from seven ignimbrite units forming the older Toconquis Group (5.60–4.51 Ma), mainly outcropping along a type section at Rio Las Pitas, Vega Real Grande. The paleomagnetic data obtained by progressive thermal demagnetization show that the clasts of the Cerro Galán ignimbrite have one single magnetic component, oriented close to the expected geomagnetic field at the time of emplacement. Results show therefore that most of the clasts acquired a new magnetization oriented parallel to the magnetic field at the moment of the ignimbrite deposition, suggesting that the clasts were heated up to or above the highest blocking temperature (T _b) of the magnetic minerals (T _b = 580°C for magnetite; T _b = 600–630°C for hematite). We obtained similar emplacement temperature estimations for six out of the seven volcanic units belonging to the Toconquis Group, with the exception of one unit (Lower Merihuaca), where we found two distinct magnetic components. The estimation of emplacement temperatures in this latter case is constrained at 580–610°C, which are lower than the other ignimbrites. These estimations are also in agreement with the lowest pre-eruptive magma temperatures calculated for the same unit (i.e., 790°C; hornblende–plagioclase thermometer; Folkes et al. 2011b). We conclude that the Cerro Galán ignimbrite and Toconquis Group ignimbrites were emplaced at temperatures equal to or higher than 620°C, except for Lower Merihuaca unit emplaced at lower temperatures. The homogeneity of high temperatures from proximal to distal facies in the Cerro Galán ignimbrite provides constraints for the emplacement model, marked by a relatively low eruption column, low levels of turbulence, air entrainment, surface–water interaction, and a high level of topographic confinement, all ensuring minimal heat loss.     

Biocalcification in the Eastern Oyster (<Emphasis Type="Italic">Crassostrea virginica)</Emphasis> in Relation to Long-term Trends in Chesapeake Bay pH

Anthropogenic carbon dioxide (CO₂) emissions reduce pH of marine waters due to the absorption of atmospheric CO₂ and formation of carbonic acid. Estuarine waters are more susceptible to acidification because they are subject to multiple acid sources and are less buffered than marine waters. Consequently, estuarine shell forming species may experience acidification sooner than marine species although the tolerance of estuarine calcifiers to pH changes is poorly understood. We analyzed 23 years of Chesapeake Bay water quality monitoring data and found that daytime average pH significantly decreased across polyhaline waters although pH has not significantly changed across mesohaline waters. In some tributaries that once supported large oyster populations, pH is increasing. Current average conditions within some tributaries however correspond to values that we found in laboratory studies to reduce oyster biocalcification rates or resulted in net shell dissolution. Calcification rates of juvenile eastern oysters, Crassostrea virginica, were measured in laboratory studies in a three-way factorial design with 3 pH levels, two salinities, and two temperatures. Biocalcification declined significantly with a reduction of ∼0.5 pH units and higher temperature and salinity mitigated the decrease in biocalcification.     

Cainozoic Stratigraphy of the lake Eyre basin and part of the arid region lying to the South

The Cainozoic history of the Lake Eyre region opened with a period of deep weathering during which many of the older rocks were extensively kaolinized. Following erosion and later deposition of a thin sheet of Tertiary fluviatile deposits, a period of weathering resulted in the widespread formation of silcrete. Another period of erosion and deposition was followed by soil formation and minor silici‐fication.

An important period of erosion followed during which some of the main elements of the present landscape were outlined. Warping during this interval gave rise to shallow basins in which lacustrine sediments accumulated. At about the same time, a system of mound springs developed near the western margin of the Great Artesian Basin. Another major period of erosion followed, by which time the main topographic features of the present landscape had evolved. This last event probably took place near the close of the Tertiary. Throughout the Tertiary, drainage was external and ancestral Lake Eyre remained fresh.

The Quaternary was characterized by four periods of aeolian and, to a lesser extent, water erosion and deposition alternating with periods of landscape stability, when weathering and soil formation took place.

Throughout the Cainozoic there was an alternation of relatively humid and dry periods, but true aridity and internal drainage did not appear until the Quaternary. Sand ridges were not formed until the late Quaternary. Intervals of gentle warping occurred from time to time during the Tertiary, but the Quaternary has for the most part been a period of stability. These events have given rise to a sequence of distinctive rock and soil‐stratigraphic units whose characteristics are considered in some detail.     

Hydrothermal models of the Perth metropolitan area, Western Australia: implications for geothermal energy

Hydrothermal simulations are used to provide insight into the subsurface thermal regime of the Perth metropolitan area (PMA) in Western Australia. High average permeabilities and estimated fluid flow rates in shallow aquifers of the PMA suggest that advection and convection may occur in these aquifers. These processes are simulated, using a new geological model of the PMA to constrain the geometry of aquifers, aquitards and faults. The results show that advection has a strong influence on subsurface temperature, especially in the north of the PMA, where aquifer recharge creates an area of anomalously low temperature. Convection may be important, depending on the permeability of the Yarragadee Aquifer. If convection occurs, it creates thermal highs and lows with a spacing of approximately 5 km. Some of these thermal anomalies migrate over geological time due to coupling between advection and convection, but they are stationary on human timescales. Fault permeability influences the pattern of convection. Advection and convection cause variations in the geothermal gradient which cannot be predicted by conductive models; therefore, these processes should be considered in any model that is used for assessment of geothermal resources in the PMA.