Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Rainfall partitioning through a mixed cedar swamp and associated C and N fluxes in Southern Ontario,Canada

Partitioning of rainfall through a forest canopy into throughfall, stemflow, and canopy interception is a critical process in the water cycle, and the contact of precipitation with vegetated surfaces leads to increased delivery of solutes to the forest floor. This study investigates the rainfall partitioning over a growing season through a temperate, riparian, mixed coniferous‐deciduous cedar swamp, an ecosystem not well studied with respect to this process. Seasonal throughfall, stemflow, and interception were 69.2%, 1.5%, and 29.3% of recorded above‐canopy precipitation, respectively. Event throughfall ranged from a low of 31.5 ± 6.8% for a small 0.8‐mm event to a high of 82.9 ± 2.4% for a large 42.7‐mm event. Rain fluxes of at least 8 mm were needed to generate stemflow from all instrumented trees. Most trees had funnelling ratios <1.0, with an exponential decrease in funnelling ratio with increasing tree size. Despite this, stand‐scale funnelling ratios averaged 2.81 ± 1.73, indicating equivalent depth of water delivered across the swamp floor by stemflow was greater than incident precipitation. Throughfall dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) averaged 26.60 ± 2.96 and 2.02 ± 0.16 mg L^?1, respectively, which were ~11 and three times above‐canopy rain levels. Stemflow DOC averaged 73.33 ± 7.43 mg L^?1, 35 times higher than precipitation, and TDN was 4.45 ± 0.56 mg L^?1, 7.5 times higher than rain. Stemflow DOC concentration was highest from Populus balsamifera and TDN greatest from Thuja occidentalis trees. Although total below‐canopy flux of TDN increased with increasing event size, DOC flux was greatest for events 20–30 mm, suggesting a canopy storage threshold of DOC was readily diluted. In addition to documenting rainfall partitioning in a novel ecosystem, this study demonstrates the excess carbon and nitrogen delivered to riparian swamps, suggesting the assimilative capacity of these zones may be underestimated.     

Calibration of paired watersheds: Utility of moving sums in presence of externalities

Historically, paired watershed studies have been used to quantify the hydrological effects of land use and management practices by concurrently monitoring 2 similar watersheds during calibration (pretreatment) and post‐treatment periods. This study characterizes seasonal water table and flow response to rainfall during the calibration period and tests a change detection technique of moving sums of recursive residuals (MOSUM) to select calibration periods for each control–treatment watershed pair when the regression coefficients for daily water table elevation were most stable to minimize regression model uncertainty. The control and treatment watersheds were 1 watershed of 3–4‐year‐old intensely managed loblolly pine (Pinus taeda L.) with natural understory, 1 watershed of 3–4‐year‐old loblolly pine intercropped with switchgrass (Panicum virgatum), 1 watershed of 14–15‐year‐old thinned loblolly pine with natural understory (control), and 1 watershed of switchgrass only. The study period spanned from 2009 to 2012. Silvicultural operational practices during this period acted as external factors, potentially shifting hydrologic calibration relationships between control and treatment watersheds. MOSUM results indicated significant changes in regression parameters due to silvicultural operations and were used to identify stable relationships for water table elevation. None of the calibration relationships developed using this method were significantly different from the classical calibration relationship based on published historical data. We attribute that to the similarity of historical and 2010–2012 leaf area index on control and treatment watersheds as moderated by the emergent vegetation. Although the MOSUM approach does not eliminate the need for true calibration data or replace the classic paired watershed approach, our results show that it may be an effective alternative approach when true data are unavailable, as it minimizes the impacts of external disturbances other than the treatment of interest.     

Factors controlling the interannual variations of MJO intensity

The interannual variations of intensity of the Madden-Julian Oscillation (MJO) during boreal winter are investigated by using the observed outgoing longwave radiation (OLR) and the reanalysis data of ECMWF and NCEP. The standard deviation of 20-80-day filtered OLR anomaly is used to measure the MJO intensity. The dominant spatial structure of the interannual variability is revealed by an EOF analysis of the MJO intensity field. It is found that the leading mode is associated with eastern Pacific type ENSO, whereas the second mode is related to central Pacific type ENSO. A simple atmospheric model is used to investigate the relative roles of background moisture and wind changes in affecting the overall strength of MJO. The numerical experiments indicate that the background moisture effect is dominant while the background wind change has a minor effect.     

Evaluation of coseismic landslide hazard on the proposed Haast-Hollyford Highway,South Island,New Zealand

Coseismic landsliding presents a major hazard to infrastructure in mountains during large earthquakes. This is particularly true for road networks, as historically coseismic landsliding has resulted in road losses larger than those due to ground shaking. Assessing the exposure of current and planned highway links to coseismic landsliding for future earthquake scenarios is therefore vital for disaster risk reduction. This study presents a method to evaluate the exposure of critical infrastructure to landsliding from scenario earthquakes from an underlying quantitative landslide hazard assessment. The method is applied to a proposed new highway link in South Island, New Zealand, for a scenario Alpine Fault earthquake and compared to the current network. Exposure (the likelihood of a network being affected by one or more landslides) is evaluated from a regional-scale coseismic landslide hazard model and assessed on a relative basis from 0 to 1. The results show that the proposed Haast-Hollyford Highway (HHH) would be highly exposed to coseismic landsliding with at least 30–40?km likely to be badly affected (the Simonin Pass route being the worse affected of the two routes). In the current South Island State Highway network, the HHH would be the link most exposed to landsliding and would increase the total network exposure by 50–70% despite increasing the total road length by just 3%. The present work is intended to provide an effective method to assess coseismic landslide hazard of infrastructure in mountains with seismic hazard, and potentially identify mitigation options and critical network segments.     

The Anglo-Brabant Massif: Persistent but enigmatic palaeo-relief at the heart of western Europe

The surface geology of central England and Belgium obscures a large ‘basement’ massif with a complex history and stronger crust and lithosphere than surrounding regions. The nucleus was forged by subduction-related magmatism at the Gondwana margin in Ediacaran time. Partitioning into a platform, in the English Midlands, and a basin stretching to Belgium, in the east, was already evident in Cambrian/earliest Ordovician time. The accretion of the Monian Composite Terrane during the Penobscotian deformation phase preceded late Tremadocian rifting, and Floian separation, of the Avalonia Terrane from the Gondwana margin. Late Ordovician magmatism in a belt from the Lake District to Belgium records subduction beneath Avalonia of part of the Tornquist Sea. This ‘Western Pacific-style’ oceanic basin closed in latest Ordovician time, uniting Avalonia and Baltica. Closure of the Iapetus Ocean in early Silurian time was soon followed by closure of the Rheic Ocean, recorded by subduction along the southern margin of the massif. The causes of late Caledonian deformation are poorly understood and controversial. Partitioned behaviour of the massif persisted into late Palaeozoic time. Late Devonian and Carboniferous sequences show strong onlap onto the massif, which was little affected by crustal extension. Compressional deformation during the Variscan Orogeny also appears slight, and was focussed in the west where a wedge-shaped mountain foreland uplift was driven by orogenic indentation, splitting the massif from the Welsh Massif along the reactivated Malvern Line. Permian to Mesozoic sequences exhibit persistent but variable degrees of onlap onto the massif.     

The role of hydrothermal fluids in sedimentation in saline alkaline lakes: Evidence from Nasikie Engida,Kenya Rift Valley

Saline alkaline lakes that precipitate sodium carbonate evaporites are most common in volcanic terrains in semi‐arid environments. Processes that lead to trona precipitation are poorly understood compared to those in sulphate‐dominated and chloride‐dominated lake brines. Nasikie Engida (Little Magadi) in the southern Kenya Rift shows the initial stages of soda evaporite formation. This small shallow (<2 m deep; 7 km long) lake is recharged by alkaline hot springs and seasonal runoff but unlike neighbouring Lake Magadi is perennial. This study aims to understand modern sedimentary and geochemical processes in Nasikie Engida and to assess the importance of geothermal fluids in evaporite formation. Perennial hot‐spring inflow waters along the northern shoreline evaporate and become saturated with respect to nahcolite and trona, which precipitate in the southern part of the lake, up to 6 km from the hot springs. Nahcolite (NaHCO₃) forms bladed crystals that nucleate on the lake floor. Trona (Na₂CO₃·NaHCO₃·2H₂O) precipitates from more concentrated brines as rafts and as bottom‐nucleated shrubs of acicular crystals that coalesce laterally to form bedded trona. Many processes modify the fluid composition as it evolves. Silica is removed as gels and by early diagenetic reactions and diatoms. Sulphate is depleted by bacterial reduction. Potassium and chloride, of moderate concentration, remain conservative in the brine. Clastic sedimentation is relatively minor because of the predominant hydrothermal inflow. Nahcolite precipitates when and where pCO₂ is high, notably near sublacustrine spring discharge. Results from Nasikie Engida show that hot spring discharge has maintained the lake for at least 2 kyr, and that the evaporite formation is strongly influenced by local discharge of carbon dioxide. Brine evolution and evaporite deposition at Nasikie Engida help to explain conditions under which ancient sodium carbonate evaporites formed, including those in other East African rift basins, the Eocene Green River Formation (western USA), and elsewhere.     

Detailing the impact of the Storegga Tsunami at Montrose,Scotland

The Storegga tsunami, dated in Norway to 8150±30 cal. years BP, hit many countries bordering the North Sea. Run-ups of >30 m occurred and 1000s of kilometres of coast were impacted. Whilst recent modelling successfully generated a tsunami wave train, the wave heights and velocities, it under-estimated wave run-ups. Work presented here used luminescence to directly date the Storegga tsunami deposits at the type site of Maryton, Aberdeenshire in Scotland. It also undertook sedimentological characterization to establish provenance, and number and relative power of the tsunami waves. Tsunami model refinement used this to better understand coastal inundation. Luminescence ages successfully date Scottish Storegga tsunami deposits to 8100±250 years. Sedimentology showed that at Montrose, three tsunami waves came from the northeast or east, over-ran pre-existing marine sands and weathered igneous bedrock on the coastal plain. Incorporation of an inundation model predicts well a tsunami impacting on the Montrose Basin in terms of replicate direction and sediment size. However, under-estimation of run-up persisted requiring further consideration of palaeotopography and palaeo-near-shore bathymetry for it to agree with sedimentary evidence. Future model evolution incorporating this will be better able to inform on the hazard risk and potential impacts for future high-magnitude submarine generated tsunami events.     

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground-penetrating radar

Across the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large-scale storms, while providing important habitats and recreational opportunities. Understanding their long-term evolution is vital in managing their condition in a changing climate. Recently, ground-penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre-scale resolution images of near-surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo-environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground-truth control. Nonetheless, GPR data interpretation can be non-intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi-automated GPR feature-extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo-landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty-organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR-interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.