Cloud computing and virtualization within the regional climate model and evaluation system

The Regional Climate Model Evaluation System (RCMES) facilitates the rapid, flexible inclusion of NASA observations into climate model evaluations. RCMES provides two fundamental components. A database (RCMED) is a scalable point-oriented cloud database used to elastically store remote sensing observations and to make them available using a space time query interface. The analysis toolkit (RCMET) is a Python-based toolkit that can be delivered as a cloud virtual machine, or as an installer package deployed using Python Buildout to users in order to allow for temporal and spatial regridding, metrics calculation (RMSE, bias, PDFs, etc.) and end-user visualization. RCMET is available to users in an “offline”, lone scientist mode based on a virtual machine dynamically constructed with model outputs and observations to evaluate; or on an institution’s computational cluster seated close to the observations and model outputs. We have leveraged RCMES within the content of the Coordinated Regional Downscaling Experiment (CORDEX) project, working with the University of Cape Town and other institutions to compare the model output to NASA remote sensing data; in addition we are also working with the North American Regional Climate Change Assessment Program (NARCCAP). In this paper we explain the contribution of cloud computing to RCMES’s specifically describing studies of various cloud databases we evaluated for RCMED, and virtualization toolkits for RCMET, and their potential strengths in delivering user-created dynamic regional climate model evaluation virtual machines for our users.     

Snow distribution from SSM/I and its relationships to the hydroclimatology of the Mackenzie River Basin,Canada

The spatial and temporal distribution of snow cover extent (SCE) and snow water equivalent (SWE) play vital roles in the hydrology of northern watersheds. We apply remotely sensed Special Sensor Microwave Imager (SSM/I) data from 1988 to 2007 to explore the relationships between snow distribution and the hydroclimatology of the Mackenzie River Basin (MRB) of Canada and its major sub-basins. The Environment Canada (EC) algorithm is adopted to retrieve the SWE from SSM/I data. Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day maximum snow cover extent products (MOD10A2) are used to estimate the different thresholds of retrieved SWE from SSM/I to classify the land cover as snow or no snow for various sub-basins in the MRB. The sub-basins have varying topography and hence different thresholds that range from 10 mm to 30 mm SWE. The accuracy of snow cover mapping based on the combination of several thresholds for the different sub-basins reaches ≈ 90%. The northern basins are found to have stronger linear relationships between the date on which snow cover fraction (SCF) reaches 50% or when SWE reaches 50% and mean air temperatures, than the southern basins. Correlation coefficients between SCF, SWE, and hydroclimatological variables show the new SCF products from SSM/I perform better than SWE from SSM/I to analyze the relationships with the regional hydroclimatology. Statistical models relating SCF and SWE to runoff indicate that the SCF and SWE from EC algorithms are able to predict the discharge in the early snow ablation seasons in northern watersheds.     

Autosuspension and the energetics of two-phase flows: Reply to comments on ‘experimental test of autosuspension’ by J. B. Southard and M. E. Mackintosh

The experimental programme of Southard and Mackintosh (1981) provides a definitive refutation of the autosuspension criterion as it was originally defined. A simple energy budget for dilute two-phase flow is derived from continuum theory of mixtures; it shows that Bagnold's criterion is based on an arbitrary comparison of two elements of the complete two-phase system. It is recommended that facile, ab hoc criteria for the behaviour of such complex flows be avoided in favour of analyses based on fundamental equations of motion.     

Characterization of magnetite particles in shocked quartz by means of electron- and magnetic force microscopy: Vredefort, South Africa

Submicroscopic opaque particles from highly shocked granite-gneisses close to the core of the Vredefort impact structure have been investigated by means of micro-analytical techniques with high spatial resolution such as electron diffraction, orientation contrast imagery and magnetic force microscopy. The opaque particles have been identified as nano- to micro-sized magnetite that occur in several distinct modes. In one sample magnetite occurs along relict planar deformation features (PDFs) in quartz, generally accepted as typical shock lamellae. The magnetite particles along shock lamellae in quartz grains virtually all show uniform crystallographic orientations. In most instances, the groups of magnetite within different quartz grains are systematically misorientated such that they share a subparallel <101> direction. The magnetite groups of all measured quartz grains thus appear to have a crystallographic preferred orientation in space. In a second sample, orientations of magnetite particles have been measured in microfractures (non-diagnostic of shock) of quartz, albite and in the alteration halos, (e.g. biotite grains breaking down to chlorite). The crystallographic orientations of magnetite particles are diverse, with only a minor portion having a preferred orientation. Scanning electron microscopy shows that magnetite along the relict PDFs is invariably associated with other microcrystalline phases such as quartz, K-feldspar and biotite. Petrographic observations suggest that these microcrystalline phases crystallized from locally formed micro-melts that intruded zones of weakness such as microfractures and PDFs shortly after the shock event. The extremely narrow widths of the PDFs suggest that heat may have dissipated rapidly resulting in melts crystallizing relatively close to where they were generated. Magnetic force microscopy confirms the presence of magnetic particles along PDFs. The smallest particles, <5 μm with high aspect ratios 15:1 usually exhibit intense, uniform magnetic signals characteristic of single-domain magnetite. Consistent offsets between attractive and repulsive magnetic signals of individual single-domain particles suggest consistent directions of magnetization for a large proportion of particles. Received: 16 November 1998 / Accepted: 17 May 1999     

Optical and near-infrared colours as a discriminant of the age and metallicity of stellar populations

We present a comprehensive analysis of the ability of current stellar population models to reproduce the optical ( ugriz ) and near-infrared ( JHK ) colours of a small sample of well-studied nearby elliptical and S0 galaxies. We find broad agreement between the ages and metallicities derived using different population models, although different models show different systematic deviations from the measured broad-band fluxes. Although it is possible to constrain simple stellar population models to a well-defined area in age–metallicity space, there is a clear degeneracy between these parameters even with such a full range of precise colours. The precision to which age and metallicity can be determined independently, using only broad-band photometry with realistic errors, is  Δ[Fe/H]≃ 0.18  and  Δlog Age ≃ 0.25  . To constrain the populations and therefore the star formation history further, it will be necessary to combine broad-band optical–IR photometry with either spectral line indices, or else photometry at wavelengths outside this range.     

The information system of the French Peatland Observation Service: Service National d'Observation Tourbières – A valuable tool to assess the impact of global changes on the hydrology and biogeochemistry of temperate peatlands through long term monitoring

Mitigating and adapting to global changes requires a better understanding of the response of the Biosphere to these environmental variations. Human disturbances and their effects act in the long term (decades to centuries) and consequently, a similar time frame is needed to fully understand the hydrological and biogeochemical functioning of a natural system. To this end, the ‘Centre National de la Recherche Scientifique’ (CNRS) promotes and certifies long-term monitoring tools called national observation services or ‘Service National d'Observation’ (SNO) in a large range of hydrological and biogeochemical systems (e.g., cryosphere, catchments, aquifers). The SNO investigating peatlands, the SNO ‘Tourbières’, was certified in 2011 ( https://www.sno-tourbieres.cnrs.fr/ ). Peatlands are mostly found in the high latitudes of the northern hemisphere and French peatlands are located in the southern part of this area. Thus, they are located in environmental conditions that will occur in northern peatlands in coming decades or centuries and can be considered as sentinels. The SNO Tourbières is composed of four peatlands: La Guette (lowland central France), Landemarais (lowland oceanic western France), Frasne (upland continental eastern France) and Bernadouze (upland southern France). Thirty target variables are monitored to study the hydrological and biogeochemical functioning of the sites. They are grouped into four datasets: hydrology, fluvial export of organic matter, greenhouse gas fluxes and meteorology/soil physics. The data from all sites follow a common processing chain from the sensors to the public repository. The raw data are stored on an FTP server. After operator or automatic processing, data are stored in a database, from which a web application extracts the data to make them available ( https://data-snot.cnrs.fr/data-access/ ). Each year at least, an archive of each dataset is stored in Zenodo, with a digital object identifier (DOI) attribution ( https://zenodo.org/communities/sno_tourbieres_data/ ).     

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.     

Meso-scale testing and development of test procedures to maintain mass balance

The Conrad Blucher Institute for Surveying and Science (Texas A&M University––Corpus Christi) has conducted numerous petroleum experiments at the Shoreline Environmental Research Facility (Corpus Christi, Texas, USA). The meso-scale facility has multiple wave tanks, permitting some control in experimental design of the investigations, but allowing for real-world conditions. This paper outlines the evolution of a materials balance approach in conducting petroleum experiments at the facility. The first attempt at a materials balance was during a 1998 study on the fate/effects of dispersant use on crude oil. Both water column and beach sediment samples were collected. For the materials balance, the defined environmental compartments for oil accumulation were sediments, water column, and the water surface, while the discharge from the tanks was presumed to be the primary sink. The “lessons learned” included a need to quantify oil adhesion to the tank surfaces. This was resolved by adhering strips of the polymer tank lining to the tank sides that could be later removed and extracted for oil. Also, a protocol was needed to quantify any floating oil on the water surface. A water surface (oil slick) quantification protocol was developed, involving the use of solid-phase extraction disks. This protocol was first tested during a shoreline cleaner experiment, and later refined in subsequent dispersant effectiveness studies. The effectiveness tests were designed to simulate shallow embayments which created the need for additional adjustments in the tanks. Since dispersant efficacy is largely affected by hydrodynamics, it was necessary to scale the hydrodynamic conditions of the tanks to those expected in our prototype system (Corpus Christi Bay, Texas). The use of a scaled model permits the experiment to be reproduced and/or evaluated under different conditions. To minimize wave reflection in the tank, a parabolic wave dissipater was built. In terms of materials balance, this design reduced available surface area as a sink for oil adsorption.     

Displacement-length scaling relations for faults on the terrestrial planets

Displacement-length (D/L)scaling relations for normal and thrust faults from Mars, and thrust faults from Mercury, for which sufficiently accurate measurements are available, are consistently smaller than terrestrial D/L ratios by a factor of about 5, regardless of fault type (i.e. normal or thrust). We demonstrate that D/L ratios for faults scale, to first order, with planetary gravity. In particular, confining pressure modulates: (1) the magnitude of shear driving stress on the fault; (2) the shear yield strength of near-tip rock; and (3) the Young's (or shear) modulus of crustal rock. In general, all three factors decrease with gravity for the same rock type and pore-pressure state (e.g. wet conditions). Faults on planets with lower surface gravities, such as Mars and Mercury, demonstrate systematically smaller D/L ratios than faults on larger planets, such as Earth. Smaller D/L ratios of faults on Venus and the Moon are predicted by this approach, and we infer still smaller values of D/L ratio for faults on icy satellites in the outer solar system. Collection of additional displacement-length and down-dip height data from terrestrial normal, strike-slip, and thrust faults, located within fold-and-thrust belts, plate margins, and continental interiors, is required to evaluate the influence of fault shape and progressive deformation on the scaling relations for faults from Earth and elsewhere.