Insights into the complexity of crustal differentiation: K2O‐poor leucosomes within metasedimentary migmatites from the Southern Marginal Zone of the Limpopo Belt,South Africa

Differentiation of the continental crust is the result of complex interactions between a large number of processes, which govern partial melting of the deep crust, magma formation and segregation, and magma ascent to significantly higher crustal levels. The anatectic metasedimentary rocks exposed in the Southern Marginal Zone of the Limpopo Belt represent an unusually well‐exposed natural laboratory where the portion of these processes that operate in the deep crust can be directly investigated in the field. The formation of these migmatites occurred via absent incongruent melting reactions involving biotite, which produced cm‐ to m‐scale, K₂O‐poor garnet‐bearing stromatic leucosomes, with high Ca/Na ratios relative to their source rocks. Field investigation combined with geochemical analyses, and phase equilibrium modelling designed to investigate some aspects of disequilibrium partial melting show that the outcrop features and compositions of the leucosomes suggest several steps in their evolution: (1) Melting of a portion of the source, with restricted plagioclase availability due to kinetic controls, to produce a magma (melt + entrained peritectic minerals in variable proportions relative to melt); (2) Segregation of the magma at near peak metamorphic conditions into melt accumulation sites (MAS), also known as future leucosome; (3a) Re‐equilibration of the magma with a portion of the bounding mafic residuum via chemical diffusion (H₂O, K₂O), which triggers the co‐precipitation of quartz and plagioclase in the MAS; (3b) Extraction of melt‐dominated magma to higher crustal levels, leaving peritectic minerals entrained from the site of the melting reaction, and the minerals precipitated in the MASs to form the leucosome in the source. The key mechanism controlling this behaviour is the kinetically induced restriction of the amount of plagioclase available to the melting reaction. This results in elevated melt H₂O and K₂O and chemical potential gradient for these components across the leucosome/mafic residuum contact. The combination of all of these processes accurately explains the composition of the K₂O‐poor leucosomes. These findings have important implications for our understanding of melt segregation in the lower crust and minimum melt residency time which, according to the chemical modelling, is <5 years. We demonstrate that in some migmatitic granulites, the leucosomes constitute a type of felsic refractory residuum, rather than evidence of failed magma extraction. This provides a new insight into the ways that source heterogeneity may control anatexis.     

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.     

Mid‐Weichselian interstadial in Kolari,western Finnish Lapland

Finnish Lapland is known as an area where numerous sites with sediments from Pleistocene glacial and interglacial periods occur. Recent sedimentological observations and dating call for reinterpretation of the record, which shows a complicated Mid‐Weichselian ice‐sheet evolution within the ice‐divide zone. Here, a large, previously unstudied section from a former Hannukainen iron mine was investigated sedimentologically and dated with optically stimulated luminescence (OSL). Ten sedimentary units were identified displaying a variety of depositional environments (glacial, glaciolacustrine, fluvial and aeolian). They are all – except for the lowermost, deeply weathered till – interpreted to be of Mid‐ or Late Weichselian/Holocene age. Five OSL samples from fluvial sediments give ages ranging from 55 to 35 ka, indicating two MIS 3 ice‐free intervals of unknown duration. The Mid‐Weichselian interstadial was interrupted by a re‐advance event, which occurred later than 35 ka and caused glaciotectonic deformation, folding and stacking of older sediments. This new evidence emphasizes the importance of the Kolari area when unravelling the complex Late Pleistocene glacial history of northern Finland and adjacent regions.     

Zircon fission‐track ages from Newfoundland—A proxy for high geothermal gradients and exhumation before opening of the Central Atlantic Ocean

Following Appalachian orogenesis, metamorphic rocks in central Newfoundland were exhumed and reburied under Tournaisian strata. New zircon fission‐track (ZFT) ages of metamorphic rocks below the Tournaisian unconformity yield post‐depositionally reset ages of 212–235 Ma indicating regional fluid‐absent reheating to at least ≥220°C. Post‐Tournaisian sedimentary thicknesses in surrounding basins show that burial alone cannot explain such temperatures, thus requiring that palaeo‐geothermal gradients increased to ≥30–40°C/km before final late Triassic accelerated cooling. We attribute these elevated palaeo‐geothermal gradients to localized thermal blanketing by insulating sediments overlying radiogenic high‐heat‐producing granitoids. Late Triassic rifting and magmatism before break up of Pangaea likely also contributed to elevated heat flow, as well as uplift, triggering late Triassic accelerated cooling and exhumation. Thermochronological ages of 240–200 Ma are seen throughout Atlantic Canada, and record rifting and basaltic magmatism on the conjugate margins of the Central Atlantic Ocean preceding the onset of oceanic spreading at ~190 Ma.     

Modeling the 2D behavior of dry‐stone retaining walls by a fully discrete element method

The dry‐stone retaining walls (DSRW) have been tipped as a promising solution for sustainable development. However, before recently, their behavior is relatively obscure. In this study, discrete element method (DEM) approach was applied to simulate the plane strain failure of these walls. A commercial DEM package (PFC2D™) was used throughout this study. The authors used a fully discrete approach; thus, both the wall and the backfill were modeled as discrete elements. The methodology for obtaining the micromechanical parameters was discussed in detail; this includes the three mechanical sub‐systems of DSRWs: wall, backfill and interface. The models were loaded progressively until failure, and then the results were compared with the full‐scale experimental results where the walls were loaded, respectively, with hydrostatic load and backfill. Despite its complexity and its intensive calculation time, DEM model can then be used to validate a more simplified approach. Copyright © 2015 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.     

Evaluating soil moisture estimation from ground-penetrating radar hyperbola fitting with respect to a systematic time-domain reflectometry data collection in a boreal podzolic agricultural field

The upper 30 cm of the soil profile, which hosts the majority of the root biomass, can be considered as the shallow agricultural root zone of most temperate crops. The electromagnetic wave velocity in the soil obtained from reflection hyperbolas in ground-penetrating radar (GPR) data can be used to estimate soil moisture (SM). Finding shallow hyperbolas in a radargram and minimizing the subjective error associated with the hyperbola fitting are the main challenges in this approach. Nevertheless, we were motivated by the recent improvements of hyperbola fitting algorithms, which can reduce the subjective error and processing time. To overcome the difficulty of finding very shallow hyperbolas, we applied the hyperbola fitting method to reflections ranging from 27- to 50-cm depth using a 500-MHz centre-frequency GPR and compared the estimated moisture with vertically installed, 30-cm-long time-domain reflectometry (TDR) sensors. We also compared TDR and GPR sample areas in a 2-D plane using different GPR survey types and different hyperbola depths. SM measured with TDR and GPR were not significantly different according to Mann–Whitney's test. Our analyses showed that a root mean square error of 0.03 m³ m⁻³ was found between the two methods. In conclusion, the proposed method might be suitable to estimate SM with an acceptable accuracy within the root zone if the soil profile is fairly uniform within the application depth range.     

The June 2013 Alberta catastrophic flooding event – part 2: fine‐scale precipitation and associated features

Data obtained from a variety of sources including the Canadian Lightning Detection Network, weather radars, weather stations and operational numerical weather model analyses were used to address the evolution of precipitation during the June 2013 southern Alberta flood. The event was linked to a mid‐level closed low pressure system to the west of the region and a surface low pressure region initially to its south. This configuration brought warm, moist unstable air into the region that led to dramatic, organized convection with an abundance of lightning and some hail. Such conditions occurred in the southern parts of the region whereas the northern parts were devoid of lightning. Initially, precipitation rates were high (extreme 15‐min rainfall rates up to 102 mm h^?1 were measured) but decreased to lower values as the precipitation shifted to long‐lived stratiform conditions. Both the convective and stratiform precipitation components were affected by the topography. Similar flooding events, such as June 2002, have occurred over this region although the 2002 event was colder and precipitation was not associated with substantial convection over southwest Alberta. Copyright © 2016 Her Majesty the Queen in Right of Canada. Hydrological Processes. © John Wiley & Sons, Ltd.