Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.     

Evaluating suspended sediment and turbidity reduction projects in a glacially conditioned catchment through dynamic regression and fluvial process-based modelling

Elevated turbidity (T_n) and suspended sediment concentrations (SSC) during and following flood events can degrade water supply quality and aquatic ecosystem integrity. Streams draining glacially conditioned mountainous terrain, such as those in the Catskill Mountains of New York State, are particularly susceptible to high levels of T_n and SSC sourced from erosional contact with glacial-related sediment. This study forwards a novel approach to evaluate the effectiveness of stream restoration best management practices (BMPs) meant to reduce stream T_n and SSC, and demonstrates the approach within the Stony Clove sub-basin of the Catskills, a water supply source for New York City. The proposed approach is designed to isolate BMP effects from natural trends in T_n and SSC caused by trends in discharge and shifts in average T_n or SSC per unit discharge (Q) following large flood events. We develop Dynamic Linear Models (DLMs) to quantify how T_n-Q and SSC-Q relationships change over time at monitoring stations upstream and downstream of BMPs within the Stony Clove and in three other sub-basins without BMPs, providing observational evidence of BMP effectiveness. A process-based model, the River Erosion Model, is then developed to simulate natural, hydrology-driven SSC-Q dynamics in the Stony Clove sub-basin (absent of BMP effects). We use DLMs to compare the modelled and observed SSC-Q dynamics and isolate the influence of the BMPs. Results suggest that observed reductions in SSC and T_n in the Stony Clove sub-basin have been driven by a combination of declining streamflow and the installed BMPs, confirming the utility of the BMPs for the monitored hydrologic conditions.     

Some observations and characteristics of lightning ground discharges in southern Ontario

Abstract

Lightning ground flash and stroke observations were made with a single‐station gated, wideband magnetic direction‐finding system with a nominal range of 180 km located in Southern Ontario during the May‐September lightning‐active seasons of 1982 and 1983. Information was recorded on the azimuth of arrival, time, amplitude, stroke multiplicity and order, and polarity. The local climatology and seasonal statistics of lightning are analysed and summarized, and compared with standard observations of thunderstorm days and hours. Regional daily flashing rates and extremes for periods of 5 to 60 min were found to have a good empirical relationship. About 15% of the flashes had multiple strokes, generally less than 10 but with as many as 14 strokes. About 8% of the flashes were positive discharges; 3% of these were multistroke with no more than 2 strokes.

The lightning activity exhibits well defined diurnal peaks in the afternoon and at night 1–3 h before sunrise. The time interval between strokes was found to have a lognormal distribution with modal and median values of 60 and 75 ms, respectively, and no significant dependence on the order of stroke. The stroke‐to‐stroke amplitude changes within the same flash show that subsequent stroke amplitudes are often greater than the first. Subsequent strokes follow many patterns of change, the most common being an amplitude oscillating with ascending stroke order. The multistroke flash duration median values rose from about 80 ms for 2‐stroke flashes to about 650 ms for 8‐stroke flashes.

Under certain assumptions of system detection efficiency and range limits a regional ground flash density of 1.62 and 2.44 km^?2 a^?1 was estimated for the two years.     

Response of low‐angle dunes to variable flow

Current understanding of bedform dynamics is largely based on field and laboratory observations of bedforms in steady flow environments. There are relatively few investigations of bedforms in flows dominated by unsteadiness associated with rapidly changing flows or tides. As a consequence, the ability to predict bedform response to variable flow is rudimentary. Using high‐resolution multibeam bathymetric data, this study explores the dynamics of a dune field developed by tidally modulated, fluvially dominated flow in the Fraser River Estuary, British Columbia, Canada. The dunes were dominantly low lee angle features characteristic of large, deep river channels. Data were collected over a field ca 1·0 km long and 0·5 km wide through a complete diurnal tidal cycle during the rising limb of the hydrograph immediately prior to peak freshet, yielding the most comprehensive characterization of low‐angle dunes ever reported. The data show that bedform height and lee angle slope respond to variable flow by declining as the tide ebbs, then increasing as the tide rises and the flow velocities decrease. Bedform lengths do not appear to respond to the changes in velocity caused by the tides. Changes in the bedform height and lee angle have a counterclockwise hysteresis with mean flow velocity, indicating that changes in the bedform geometry lag changes in the flow. The data reveal that lee angle slope responds directly to suspended sediment concentration, supporting previous speculation that low‐angle dune morphology is maintained by erosion of the dune stoss and crest at high flow, and deposition of that material in the dune trough.     

Significance of shallow core transects for reef models and sea‐level curves,Heron Reef,Great Barrier Reef

A sequence of shallow reef cores from Heron Reef, Great Barrier Reef, provides new insights into Holocene reef growth models. Isochron analysis of a leeward core transect suggests that the north‐western end of Heron Reef reached current sea‐level by ca 6·5 kyr bp and then prograded leeward at a rate of ca 19·6 m/kyr between 5·1 kyr and 4·1 kyr bp (pre‐1950) to the present reef margin. A single short core on the opposing margin of the reef is consistent with greater and more recent progradation there. Further to the east, one windward core reached modern sea‐level by ca 6·3 kyr bp , suggesting near ‘keep‐up’ behaviour at that location, but the opposing leeward margin behind the lagoon reached sea‐level much more recently. Hence, Heron Reef exhibited significantly different reef growth behaviour on different parts of the same margin. Mean reef accretion rates calculated from within 20 m of one another in the leeward core transect varied between ca 2·9 m and 4·7 m/kyr depending on relative position in the prograding wedge. These cores serve as a warning regarding the use of isolated cores to inform reef growth rates because apparent aggradation at any given location on a reef varies depending on its location relative to a prograding margin. Only transects of closely spaced cores can document reef behaviour adequately so as to inform reef growth models and sea‐level curves. The cores also emphasize potential problems in U‐series dates for corals within a shallow (ca 1·5 m) zone beneath the reef flat. Apparent age inversions restricted to that active diagenetic zone may reflect remobilization and concentration of Th in irregularly distributed microbialites or biofilms that were missed during sample vetting. Importantly, the Th‐containing contaminant causes ages to appear too old, rather than too young, as would be expected from younger cement.     

Anisotropic mechanical properties of zircon and the effect of radiation damage

This study provides new insights into the relationship between radiation-dose-dependent structural damage due to natural U and Th impurities and the anisotropic mechanical properties (Poisson’s ratio, elastic modulus and hardness) of zircon. Natural zircon samples from Sri Lanka (see Muarakami et al. in Am Mineral 76:1510–1532, 1991) and synthetic samples, covering a dose range of zero up to 6.8 × 10¹⁸ α-decays/g, have been studied by nanoindentation. Measurements along the [100] crystallographic direction and calculations, based on elastic stiffness constants determined by Özkan (J Appl Phys 47:4772–4779, 1976), revealed a general radiation-induced decrease in stiffness (~54 %) and hardness (~48 %) and an increase in the Poisson’s ratio (~54 %) with increasing dose. Additional indentations on selected samples along the [001] allowed one to follow the amorphization process to the point that the mechanical properties are isotropic. This work shows that the radiation-dose-dependent changes of the mechanical properties of zircon can be directly correlated with the amorphous fraction as determined by previous investigations with local and global probes (Ríos et al. in J Phys Condens Matter 12:2401–2412, 2000a; Farnan and Salje in J Appl Phys 89:2084–2090, 2001; Zhang and Salje in J Phys Condens Matter 13:3057–3071, 2001). The excellent agreement, revealed by the different methods, indicates a large influence of structural and even local phenomena on the macroscopic mechanical properties. Therefore, this study indicates the importance of acquiring better knowledge about the mechanical long-term stability of radiation-damaged materials.