Assessing the accuracy of phase equilibrium software in reproducing the liquidus multiple saturation conditions of lunar and Martian basalt compositions

Accurate computational modeling allows the use of software as a first approach to some petrological problems that typically require experimentation, but most programs have not yet been fully tested for accuracy with lunar or Martian melt compositions. The programs pMELTS, MAGPOX, and Perple_X stand out for phase equilibrium modeling, as their calibrations include experiments of lunar compositions or have precise thermodynamic constraints for similar compositions. A set of lunar mare basalts, picritic glasses, and basaltic Martian compositions with known experimentally determined multiple saturation point (MSP) conditions were used here for phase equilibrium modeling. The accuracy of each program was tested through the determination of MSPs on the liquidus of the selected compositions. This point in pressure–temperature space can be considered as a direct proxy of the stable phases and the equilibrium conditions during partial melting of mantle sources. We identify a trend in experimental data between MSP temperature and MgO, CaO, and SiO₂ concentrations, and similar trends are found in model results. However, only Perple_X is able to closely match the experimental data, despite the fact it does not accurately model ilmenite saturation for high-Ti lunar basalts. We find that pMELTS miscalculates olivine saturation for MgO-rich compositions and MAGPOX systematically underestimates MSP pressure and temperatures and can only be used when olivine is the liquidus phase. For modeling lunar or Martian basalt compositions, Perple_X can be used for optimal results, although no software is yet capable of bypassing the need to constrain MSP conditions through experimentation.     

Investigating potential future changes in surface water flooding hazard and impact

Surface water flooding (SWF) is a recurrent hazard that affects lives and livelihoods. Climate change is projected to change the frequency of extreme rainfall events that can lead to SWF. Increasingly, data from Regional Climate Models (RCMs) are being used to investigate the potential water-related impacts of climate change; such assessments often focus on broad-scale fluvial flooding and the use of coarse resolution (>12 km) RCMs. However, high-resolution (<4 km) convection-permitting RCMs are now becoming available that allow impact assessments of more localised SWF to be made. At the same time, there has been an increasing demand for more robust and timely real-time forecast and alert information on SWF. In the UK, a real-time SWF Hazard Impact Model framework has been developed. The system uses 1-km gridded surface runoff estimates from a hydrological model to simulate the SWF hazard. These are linked to detailed inundation model outputs through an Impact Library to assess impacts on property, people, transport, and infrastructure for four severity levels. Here, a set of high-resolution (1.5 km and 12 km) RCM data has been used as input to a grid-based hydrological model over southern Britain to simulate Current (1996–2009) and Future (~2100s; RCP8.5) surface runoff. Counts of threshold-exceedance for surface runoff and precipitation (at 1-, 3- and 6-hr durations) are analysed. Results show that the percentage increases in surface runoff extremes, are less than those of precipitation extremes. The higher-resolution RCM simulates the largest percentage increases, which occur in winter, and the winter exceedance counts are greater than summer exceedance counts. For property impacts, the largest percentage increases are also in winter; however, it is the 12-km RCM output that leads to the largest percentage increase in impacts. The added-value of high-resolution climate model data for hydrological modelling is from capturing the more intense convective storms in surface runoff estimates.     

Cascades of sub-decadal,channel-floodplain changes in low-gradient,non-vegetated reaches near a dryland river terminus: Salar de Uyuni,Bolivia

The terminus of the ephemeral Río Colorado is located at the margins of Salar de Uyuni, Bolivia, the world's largest salt lake. The low-gradient (<0.0006 m m^-1), non-vegetated reaches approaching the terminus provide an excellent natural laboratory for investigating cascades of channel-floodplain changes that occur in response to quasi-regular flows (at least once annually) and fine-grained sediment supply (dominantly silt and clay). High-resolution satellite imagery (<0.65 m, various dates from 2004 onwards) and field data reveal widespread, pronounced and rapid morphodynamics on sub-decadal timescales, including channel erosion and chute cutoff formation, and development of crevasse channels and splays, floodouts (unchannelled surfaces at channel termini), and erosion cells (floodplain scour-transport-fill features). In particular, following high annual precipitation (>400 mm) in 2004–2005 and two subsequent high magnitude daily precipitation events (~40 mm), all of which led to widespread flooding, numerous crevasse splays formed between 2004 and 2016, avulsions occurred at nearby floodouts, and erosion cells downstream of the splays and floodouts underwent striking morphological changes. High-precision GPS data reveal two preferential localities for erosion cell development: partially or fully abandoned channels with crevasse splay remnants, and topographic lows between channels. In this overall low-gradient setting, comparatively high gradients (up to ~0.0006 m m^-1) at the edge of splay deposits and topography created by crevasses and abandoned channels may initiate knickpoint retreat and thereafter erosion cell development. Abandoned channels with splays tend to give rise to narrow, deep erosion cells, while topographic lows promote relatively shallow, wide erosion cells. In both situations, erosion cells may extend upslope and downslope, and eventually connect to form straight channels. The channel-floodplain morphodynamics near the Río Colorado terminus extend previous analyses of low-gradient, dryland river systems, particularly because the lack of vegetation and quasi-regular floods drive cascades of rapid changes on sub-decadal timescales. © 2018 John Wiley & Sons, Ltd.     

Application of HVSR to estimating thickness of laterite weathering profiles in basalt

Lateritic weathering profiles (LWPs) are widespread in the tropics and comprise an important component of the Critical Zone (CZ). The Hawaiian Islands make an excellent natural laboratory for examining the tropical CZ, where the bedrock composition (basalt) is nearly uniform and rainfall varies greatly. This natural laboratory is employed to assess the utility of the HVSR (horizontal/vertical spectral ratio) method to characterize the shear-wave velocity (V_s) structure of LWPs, particularly the depth to the contact between saprolite and basalt bedrock. LWP thicknesses determined from HVSR provide good agreement with multi-channel analysis of surface waves (MASW) profiles, well logs and outcrop. LWP thicknesses may be estimated from the fundamental mode equation or through forward models. Prior knowledge about the subsurface from well, outcrop, and MASW profiles may greatly aid modeling in some cases. For the 3.2 to 1.8 Ma Koolau Volcano on Oahu, the downward rate of advance of the weathering front varies from 0.004 to 0.041 m/ka. For the 0.44 to 0.10 Ma Kohala Volcano (Big Island of Hawaii) rates vary from 0.013 to 0.047 m/ka. Simple H/V spectra develop in areas where the combined effects of time and elevated rainfall produce thick LWPs with a flat base and a general absence of core stones with an ideal layered geometry. Abundant buried core stones violate the assumption of simple layered geometries and scatter acoustic energy, leading to uninterpretable results. This is common where low rainfall and a young basaltic substrate leave abundant core stones as well as an undulating contact between saprolite and bedrock. Velocity inversions (high V_s intervals within low V_s saprolite) may also be present and originate from relatively intact bedrock horizons or mineralogical changes within saprolite. At Kohala, a gibbsite-rich horizon produces such a velocity inversion due to enhanced weathering and subsequent collapse of saprolite in a discrete horizon. © 2019 John Wiley & Sons, Ltd.     

Detection of subsurface archaeological architecture by computer assisted airphoto interpretation

Computer-assisted enhancement of scanned historical airphotos demonstrates successful detection and mapping of subsurface archaeological features. Results show a locally adaptive contrast enhancement filter produced the best resolution of lineaments. Such enhancements can enable reconstruction of protohistoric architecture despite modern construction. The city-state of Idalion is used here as a case study to demonstrate enhancement. Results are confirmed by excavation of Idalion architecture. Idalion was a thriving metropolis during Cypro-Archaic to Roman times, surrounded by a 9.6 km defensive wall. The area within may still contain significant amounts of unmapped architecture. Some of this architecture was mapped previously from aerial photographs using manual interpretation techniques, but never excavated as proof of detection. Guided by enhanced imagery, our excavations confirmed the location of five towers/bastions along the city wall. Enhanced airphotos also revealed a region disturbed by contemporary construction of a water pipeline. Such information is important to archaeologists planning digs, and can be acquired using standard digital image processing techniques. © 1999 John Wiley & Sons, Inc.     

Analysis of noble gases in water using a quadrupole mass spectrometer in static mode

A quadrupole-based mass spectrometer used in static mode has been employed to determine noble gas concentrations of water samples by isotope dilution. Water samples are degassed and spiked with minor isotopes of each noble gas. After separation of the various gas components, the isotopic ratios are measured in the mass spectrometer. Because a spike is added to each sample, the recovery of each gas component does not have to be quantitative. Equilibration temperatures are calculated from the gas concentrations, based on noble gas solubility data. An accuracy of ±1°C, from the air equilibration temperature, is attainable for air-saturated water samples prepared under controlled conditions in the laboratory. The method has been applied to groundwaters from the London Basin Chalk aquifer. The noble-gas-derived temperatures are consistent with climatic conditions prevailing at the periods of recharge.     

‘Bond cycles’ recorded in terrestrial Pleistocene sediments of southwestern British Columbia,Canada

Recent data from exposures of terrestrial Pleistocene sediments in the Fraser Lowland of southwestern British Columbia reveal at least two ‘Bond cycles’ within Oxygen Isotope Stage 2. The maximum of the Coquitlam Stade coincides with the timing of Heinrich event H2, the Port Moody Interstade with Dansgaard–Oeschger (D–O) interstade 2, the maximum of the Vashon Stade with H1, and the Fort Langley interval with D–O interstade 1. The Sumas Stade apparently preceded H0 (Younger Dryas) but could have been in response to the same climatic signal. The timing of Sumas advances may be explained by a combination of glacio-isostatic rebound, destabilisation of the ice margin, and rapid movement over a short distance on soft muddy beds of a rising sea floor, thereby leading the timing of North Atlantic events by hundreds of years. In contrast, Coquitlam and Vashon advances were mainly over permeable glaciofluvial sediments and because of this their maxima probably did not precede the timing of H2 and H1. The Port Moody Interstade coincided with the global Last Glacial Maximum, due in part to the moderating effect of moist summer storms in a southward-shifted jet stream that influenced the Fraser Lowland at that time. Copyright © 1999 John Wiley & Sons, Ltd.     

Bed material texture in low order streams on the Queen Charlotte Islands,British Columbia

Low order channels comprise a large proportion of the links of every drainage basin, and are often at the centre of land management concerns. They exhibit hydrological and geomorphological characteristics atypical of higher order links. This paper examines the nature and causes of variations in the bed material texture of two streams on the Queen Charlotte Islands, British Columbia. The extant, functional exponential model is found to be inadequate for explaining observed changes in grain size parameters with distance downstream. Recurrent disruption of sediment transport by large organic debris jams, and the sporadic contamination of the fluvial sediment population by colluvial inputs, preclude the development of longitudinal structure. Rather, grain size varies erratically over short distances. A stochastic model best describes the observed variations, and should be adopted as an alternative to the exponential model in low order links. Characteristic variances are controlled by the degree of hillslope-channel coupling, and the extent and characteristics of non-alluvial storage mechanisms.     

Spatial and temporal variation in subtidal molluscan diversity amongst temperate estuarine habitats

Effective management of marine ecosystems is enhanced when detailed information on biodiversity is available. Key information to underpin management actions and conservation planning includes relationships between species assemblages and environmental gradients, and information on species distributions. We conducted a subtidal biodiversity assessment of surface‐dwelling subtidal molluscs in eight a priori defined habitat types using underwater visual censuses to quantitatively explore relationships between molluscan assemblages, and their correlation with benthic habitats and abiotic variables. In addition, variations in diversity were examined for two key habitat types (areas dominated by Dendronephthya australis and by filter feeders) over a period of 15 months to examine temporal change. We found that molluscs form distinct assemblages within subtidal habitats, but that assemblages within key habitats show inherent temporal variability. Regional (gamma) diversity of molluscs was found to result from a combination of: (i) within habitat alpha diversity, which increased with habitat complexity; (ii) between habitat beta diversity, with significant differences in molluscan assemblages amongst habitats with differing benthic growth, substrate type, and depth; and (iii) temporal beta diversity, with significant changes detected in molluscan assemblages over time. The results demonstrate how habitats and abiotic variables (principally depth and substrate type) combine to contribute to molluscan biodiversity in temperate estuaries, and illustrate the value of these factors as surrogates for surface‐dwelling subtidal molluscs in conservation planning.