Magnetic characterisation and correlation of a Younger Dryas tephra in North Atlantic marine sediments

A technique for identifying non‐visible basaltic tephra‐rich horizons of Younger Dryas (YD)/Greenland Stadial (GS) 1 age in northeast Atlantic sediments using rapid, non‐destructive magnetic measurements is presented. Three high‐resolution marine sediment cores have been studied in an E–W transect across the Hebridean margin: St Kilda Basin (MD95‐2007), Barra Fan (MD95‐2006) and Rockall Trough (MD04‐2822). Magnetic susceptibilities and remanent magnetisations were measured at contiguous 1 cm resolution on bulk sediments. In all three cores, an interval with higher proportions of hard magnetic minerals coincides with a clearly defined peak in basaltic tephra shard (>250 µm) counts, which can be constrained to the early part of the YD/GS1 based on faunal climate proxies. Electron microprobe analyses of the magnetically distinct basaltic tephra interval, in all three cores, displays the same major element geochemistry as published for the Vedde basaltic (I Tab. 1), i.e. sourced from the Icelandic volcano Katla. The identification of transitional alkalic basaltic tephras within marine sediments could potentially be facilitated by magnetic analysis as a useful chronostratigraphic screening tool. Copyright © 2009 John Wiley & Sons, Ltd.     

Wind‐driven depth‐averaged circulation in Queen Charlotte Sound and Hecate Strait

Abstract

We develop a wind‐driven depth‐averaged model of the circulation on the continental shelf around the Queen Charlotte Islands. The model captures a major feature of the winter current‐meter observations: a flow in Moresby Trough against the direction of the prevailing winds. Moresby Trough is a steep submarine canyon cutting across the shelf from the Pacific Ocean to the mainland. The flow patterns revealed by simulated drifters lead to four generalizations about the depth‐averaged, wind‐driven flow: (1) the flow is subject to strong topographic steering, (2) the exchange between Queen Charlotte Sound and the Pacific Ocean is limited to small regions near Cape St James and Cape Scott, (3) the exchange between Queen Charlotte Sound and Hecate Strait is controlled by Moresby Trough, and (4) the observed outflows past Cape St James are not explained by the dynamics of this model.     

A goodness of fit measure related to r2 for model performance assessment

Checking the predictive worth of an environmental model inevitably includes a goodness of fit metric to quantify the degree of matching to recorded data, thereby giving a measure of model performance. Considerable analysis and discussion have taken place over fit indices in hydrology, but a neglected aspect is the degree of communicability to other disciplines. It is suggested that a fit index is best communicated to colleagues via reference to models giving unbiased predictions, because unbiased environmental models are a desirable goal across disciplines. That is, broad recognition of a fit index is aided if it simplifies in the unbiased case to a familiar and logical expression. This does not hold for the Nash–Sutcliffe Efficiency E which reduces to the somewhat awkward unbiased expression E = 2 – 1/r², where r² is the coefficient of determination. A new goodness of fit index V is proposed for model validation as V = r²/(2‐E), which simplifies to the easily communicated V = r⁴ in the unbiased case. The index is defined over the range 0 ≤ V ≤ 1, and it happens that V < E for larger values of E. Some synthetic and recorded data sets are used to illustrate characteristics of V in comparison to E. Copyright © 2013 John Wiley & Sons, Ltd.     

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.     

Evaluating digital elevation models for glaciologic applications: An example from Nevado Coropuna, Peruvian Andes

This paper evaluates the suitability of readily available elevation data derived from recent sensors – the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and the Shuttle Radar Topography Mission (SRTM) – for glaciological applications. The study area is Nevado Coropuna (6426 m), situated in Cordillera Ampato of Southern Peru. The glaciated area was 82.6 km² in 1962, based on aerial photography. We estimate the glacier area to be ca. 60.8 km² in 2000, based on analysis of the ASTER L1B scene.We used two 1:50,000 topographic maps constructed from 1955 aerial photography to create a digital elevation model with 30 m resolution, which we used as a reference dataset. Of the various interpolation techniques examined, the TOPOGRID algorithm was found to be superior to other techniques, and yielded a DEM with a vertical accuracy of ± 14.7 m. The 1955 DEM was compared to the SRTM DEM (2000) and ASTER DEM (2001) on a cell-by-cell basis. Steps included: validating the DEM's against field GPS survey points on rock areas; visualization techniques such as shaded relief and contour maps; quantifying errors (bias) in each DEM; correlating vertical differences between various DEM's with topographic characteristics (elevation, slope and aspect) and subtracting DEM elevations on a cell-by-cell basis.The RMS error of the SRTM DEM with respect to GPS points on non-glaciated areas was 23 m. The ASTER DEM had a RMS error of 61 m with respect to GPS points and displayed 200–300 m horizontal offsets and elevation ‘spikes’ on the glaciated area when compared to the DEM from topographic data.Cell-by-cell comparison of SRTM and ASTER-derived elevations with topographic data showed ablation at the toes of the glaciers (− 25 m to − 75 m surface lowering) and an apparent thickening at the summits. The mean altitude difference on glaciated area (SRTM minus topographic DEM) was − 5 m, pointing towards a lowering of the glacier surface during the period 1955–2000. Spurious values on the glacier surface in the ASTER DEM affected the analysis and thus prevented us from quantifying the glacier changes based on the ASTER data.     

Nature and timing of large landslides in the Himalaya and Transhimalaya of northern India

Four large landslides, each with a debris volume >10⁶ m³, in the Himalaya and Transhimalaya of northern India were examined, mapped, and dated using ¹⁰Be terrestrial cosmogenic radionuclide surface exposure dating. The landslides date to 7.7±1.0 ka (Darcha), 7.9±0.8 ka (Patseo), 6.6±0.4 ka (Kelang Serai), and 8.5±0.5 ka (Chilam). Comparison of slip surface dips and physically reasonable angles of internal friction suggests that the landslides may have been triggered by increased pore water pressure, seismic shaking, or a combination of these two processes. However, the steepness of discontinuities in the Darcha rock-slope, suggests that it was more likely to have started as a consequence of gravitationally-induced buckling of planar slabs. Deglaciation of the region occurred more than 2000 years before the Darcha, Patseo, and Kelang Serai landslides; it is unlikely that glacial debuttressing was responsible for triggering the landslides. The four landslides, their causes, potential triggers and mechanisms, and their ages are compared to 12 previously dated large landslides in the region. Fourteen of the 16 dated landslides occurred during periods of intensified monsoons. Seismic shaking, however, cannot be ruled out as a mechanism for landslide initiation, because the Himalaya has experienced great earthquakes on centennial to millennial timescales. The average Holocene landscape lowering due to large landslides for the Lahul region, which contains the Darcha, Patseo, and Kelang Serai landslides, is ～0.12 mm/yr. Previously published large-landslide landscape-lowering rates for the Himalaya differ significantly. Furthermore, regional glacial and fluvial denudation rates for the Himalaya are more than an order of magnitude greater. This difference highlights the lack of large-landslide data, lack of chronology, problems associated with single catchment/large landslide-based calculations, and the need for regional landscape-lowering determinations over a standardized time period.