The stiff upper LIP: investigating the High Arctic Large Igneous Province

The Canadian Arctic Islands expose a complex network of dykes and sills that belong to the High Arctic Large Igneous Province (HALIP), which intruded volatile‐rich sedimentary rocks of the Sverdrup Basin (shale, limestone, sandstone and evaporite) some 130 to 120 million years ago. There is thus great potential in studying the HALIP to learn how volatile‐rich sedimentary rocks respond to magmatic heating events during LIP emplacement. The HALIP remains, however, one of the least well known LIPs on the planet due to its remote location, short field season, and harsh climate. A Canadian–Swedish team of geologists set out in summer 2015 to further explore HALIP sills and their sedimentary host rocks, including the sampling of igneous and meta‐sedimentary rocks for subsequent geochemical analysis, and high pressure‐temperature petrological experiments to help define the actual processes and time‐scales of magma–sediment interaction. The research results will advance our understanding of how climate‐active volatiles such as CO_2, SO₂ and CH₄ are mobilised during the magma–sediment interaction related to LIP events, a process which is hypothesised to have drastically affected Earth's carbon and sulphur cycles. In addition, assimilation of sulphate evaporites, for example, is anticipated to trigger sulphide immiscibility in the magma bodies and in so doing could promote the formation of Ni‐PGE ore bodies. Here we document the joys and challenges of ‘frontier arctic fieldwork’ and discuss some of our initial observations from the High Arctic Large Igneous Province.     

Middle Paleolithic occupation on a Marine Isotope Stage 5 lakeshore in the Nefud Desert,Saudi Arabia

Major hydrological variations associated with glacial and interglacial climates in North Africa and the Levant have been related to Middle Paleolithic occupations and dispersals, but suitable archaeological sites to explore such relationships are rare on the Arabian Peninsula. Here we report the discovery of Middle Paleolithic assemblages in the Nefud Desert of northern Arabia associated with stratified deposits dated to 75,000 years ago. The site is located in close proximity to a substantial relict lake and indicates that Middle Paleolithic hominins penetrated deeply into the Arabian Peninsula to inhabit landscapes vegetated by grasses and some trees. Our discovery supports the hypothesis of range expansion by Middle Paleolithic populations into Arabia during the final humid phase of Marine Isotope Stage 5, when environmental conditions were still favorable.     

Response of suspended sediment particle size distributions to changes in water chemistry at an Andean mountain stream confluence receiving arsenic rich acid drainage

Acid drainage is an important water quality issue in Andean watersheds, affecting the sustainability of urban, agricultural, and industrial activities. Mixing zones receiving acid drainage are critical sites where changes in pH and chemical environment promote the formation and dissolution of iron and aluminum oxy/hydroxides. These particles can significantly change the speciation of toxic metals and metalloids throughout drainage networks via sorption, desorption, and settling processes. However, little is known about the behavior of particle size distributions (PSDs) in streams affected by acid drainage and their relationship to metal speciation. This work studied: (a) the PSDs for a wide range of mixing ratios found at a fluvial confluence affected by acid drainage, and (b) the response of PSDs and arsenic speciation to environmental changes found when the particles approach complete mixing conditions. The confluence between the Azufre River (pH ~ 2, high concentration of dissolved metals) and Caracarani River (pH = 8.6, low concentration of dissolved metals) was used as a representative model for study. Field measurements show a bimodal PSD with modal diameters of ~50 and 300 μm. At shorter distances from the junction, the smaller modes with smaller particle volumes were dominant across the stream cross‐sections. A systematic shift towards larger particle sizes and larger particle volumes occurred downstream. The analysis of laboratory PSDs for Azufre/Caracarani mixing ratios between 0.01 and 0.5 (pH from 6.2 to 2.3) showed a bimodal trend with ~15 and 50 μm characteristic diameters; larger particles formed at pH>4. When particle suspensions were transferred in laboratory experiments from very low pH to full mixing conditions (pH ~ 2.8 and mixing ratio ~ 0.25) particle sizes varied, and the dissolved arsenic concentration decreased. The observed reaction kinetics were slow compared to the time scale of advective transport, creating opportunities for engineered controls for arsenic. This work contributes to a better understanding of the chemical‐hydrodynamic interactions in watersheds affected by mining, and identifying opportunities to improve water quality at points of use.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Exploration of remotely sensed forest structure and ultrasonic range sensor metrics to improve empirical snow models

Current methods to estimate snow accumulation and ablation at the plot and watershed levels can be improved as new technologies offer alternative approaches to more accurately monitor snow dynamics and their drivers. Here we conduct a meta‐analysis of snow and vegetation data collected in British Columbia to explore the relationships between a wide range of forest structure variables – obtained from Light Detection and Ranging (LiDAR), hemispherical photography (HP) and Landsat Thematic Mapper – and several indicators of snow accumulation and ablation estimated from manual snow surveys and ultrasonic range sensors. By merging and standardizing all the ground plot information available in the study area, we demonstrate how LiDAR‐derived forest cover above 0.5 m was the variable explaining the highest percentage of absolute peak snow water equivalent (SWE) (33%), while HP‐derived leaf area index and gap fraction (45° angle of view) were the best potential predictors of snow ablation rate (explaining 57% of variance). This study reveals how continuous SWE data from ultrasonic sensors are fundamental to obtain statistically significant relationships between snow indicators and structural metrics by increasing mean r² by 20% when compared to manual surveys. The relationships between vegetation and spectral indices from Landsat and snow indicators, not explored before, were almost as high as those shown by LiDAR or HP and thus point towards a new line of research with important practical implications. While the use of different data sources from two snow seasons prevented us from developing models with predictive capacity, a large sample size helped to identify outliers that weakened the relationships and suggest improvements for future research. A concise overview of the limitations of this and previous studies is provided along with propositions to consistently improve experimental designs to take advantage of remote sensing technologies, and better represent spatial and temporal variations of snow. Copyright © 2013 John Wiley & Sons, Ltd.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

A GLUE‐based uncertainty assessment framework for tritium‐inferred transit time estimations under baseflow conditions

The last decade has seen major technical and scientific improvements in the study of water transfer time through catchments. Nevertheless, it has been argued that most of these developments used conservative tracers that may disregard the oldest component of water transfer, which often has transit times greater than 5 years. Indeed, although the analytical reproducibility of tracers limits the detection of the older flow components associated with the most dampened seasonal fluctuations, this is very rarely taken into account in modelling applications. Tritium is the only environmental tracer at hand to investigate transfer times in the 5‐ to 50‐year range in surface waters, as dissolved gases are not suitable due to the degassing process. Water dating with tritium has often been difficult because of the complex history of its atmospheric concentration, but its current stabilization together with recent analytical improvements open promising perspectives. In this context, the innovative contribution of this study lies in the development of a generalized likelihood uncertainty estimation‐based approach for analysing the uncertainties associated with the modelling of transit time due to both parameter identification and tracer analytical precision issues. A coupled resampling procedure allows assessment of the statistical significance of the transfer time differences found in diverse waters. This approach was developed for tritium and the exponential‐piston model but can be implemented for virtually any tracer and model. Stream baseflow, spring and shallow aquifer waters from the Vallcebre research catchments, analysed for tritium in different years with different analytical precisions, were investigated by using this approach and taking into account other sources of uncertainty. The results showed three groups of waters of different mean transit times, with all the stream baseflow and spring waters older than the 5‐year threshold needing tritium. Low sensitivity of the results to the model structure was also demonstrated. Dual solutions were found for the waters sampled in 2013, but these results may be disambiguated when additional analyses will be made in a few years. Copyright © 2016 John Wiley & Sons, Ltd.