Solid-state astrochemistry in star-forming regions

Surface-catalysed chemical reactions, the molecules they form, and the implications for star formation were key themes in an international workshop on "Solid-state chemistry in star-forming regions" at the Lorentz Centre, Leiden University, the Netherlands from 14–17 April 2003. Helen Fraser, David Williams, Ian Sims, Anita Richards and Jeremy Yates report.     

The Phase A study of the ESA M4 mission candidate ARIEL

ARIEL, the Atmospheric Remote sensing Infrared Exoplanet Large survey, is one of the three M-class mission candidates competing for the M4 launch slot within the Cosmic Vision science programme of the European Space Agency (ESA). As such, ARIEL has been the subject of a Phase A study that involved European industry, research institutes and universities from ESA member states. This study is now completed and the M4 down-selection is expected to be concluded in November 2017. ARIEL is a concept for a dedicated mission to measure the chemical composition and structure of hundreds of exoplanet atmospheres using the technique of transit spectroscopy. ARIEL targets extend from gas giants (Jupiter or Neptune-like) to super-Earths in the very hot to warm zones of F to M-type host stars, opening up the way to large-scale, comparative planetology that would place our own Solar System in the context of other planetary systems in the Milky Way. A technical and programmatic review of the ARIEL mission was performed between February and May 2017, with the objective of assessing the readiness of the mission to progress to the Phase B1 study. No critical issues were identified and the mission was deemed technically feasible within the M4 programmatic boundary conditions. In this paper we give an overview of the final mission concept for ARIEL as of the end of the Phase A study, from scientific, technical and operational perspectives.     

Impact history of the Apollo 17 landing site revealed by U‐Pb SIMS ages

Secondary ion mass spectrometry (SIMS) U‐Pb ages of Ca‐phosphates from four texturally distinct breccia samples (72255, 76055, 76015, 76215) collected at the Apollo 17 landing site were obtained in an attempt to identify whether they represent a single or several impact event(s). The determined ages, combined with inferences from petrologic relationships, may indicate two or possibly three different impact events at 3920 ± 3 Ma, 3922 ± 5 Ma, and 3930 ± 5 Ma (all errors 2σ). Searching for possible sources of the breccias by calculating the continuous ejecta radii of impact basins and large craters as well as their expected ejecta thicknesses, we conclude that Nectaris, Crisium, Serenitatis, and Imbrium are likely candidates. If the previous interpretation that the micropoikilitic breccias collected at the North Massif represent Serenitatis ejecta is correct, then the average ²⁰⁷Pb/²⁰⁶Pb age of 3930 ± 5 Ma (2σ) dates the formation of the Serenitatis basin. The occurrence of zircon in the breccias sampled at the South Massif, which contain Ca‐phosphates yielding an age of 3922 ± 5 Ma (2σ), may indicate that the breccia originated from within the Procellarum KREEP terrane (PKT) and the Imbrium basin appears to be the only basin that could have sourced them. However, this interpretation implies that all basins suggested to fall stratigraphically between Serenitatis and Imbrium formed within a short (<11 Ma) time interval, highlighting serious contradictions between global stratigraphic constraints, sample interpretation, and chronological data. Alternatively, the slightly older age of the two micropoikilitic breccias may be a result of incomplete resetting of the U‐Pb system preserved in some phosphate grains. Based on the currently available data set this possibility cannot be excluded.     

Regolith breccia Northwest Africa 7533: Mineralogy and petrology with implications for early Mars

Northwest Africa 7533, a polymict Martian breccia, consists of fine‐grained clast‐laden melt particles and microcrystalline matrix. While both melt and matrix contain medium‐grained noritic‐monzonitic material and crystal clasts, the matrix also contains lithic clasts with zoned pigeonite and augite plus two feldspars, microbasaltic clasts, vitrophyric and microcrystalline spherules, and shards. The clast‐laden melt rocks contain clump‐like aggregates of orthopyroxene surrounded by aureoles of plagioclase. Some shards of vesicular melt rocks resemble the pyroxene‐plagioclase clump‐aureole structures. Submicron size matrix grains show some triple junctions, but most are irregular with high intergranular porosity. The noritic‐monzonitic rocks contain exsolved pyroxenes and perthitic intergrowths, and cooled more slowly than rocks with zoned‐pyroxene or fine grain size. Noritic material contains orthopyroxene or inverted pigeonite, augite, calcic to intermediate plagioclase, and chromite to Cr‐bearing magnetite; monzonitic clasts contain augite, sodic plagioclase, K feldspar, Ti‐bearing magnetite, ilmenite, chlorapatite, and zircon. These feldspathic rocks show similarities to some rocks at Gale Crater like Black Trout, Mara, and Jake M. The most magnesian orthopyroxene clasts are close to ALH 84001 orthopyroxene in composition. All these materials are enriched in siderophile elements, indicating impact melting and incorporation of a projectile component, except for Ni‐poor pyroxene clasts which are from pristine rocks. Clast‐laden melt rocks, spherules, shards, and siderophile element contents indicate formation of NWA 7533 as a regolith breccia. The zircons, mainly derived from monzonitic (melt) rocks, crystallized at 4.43 ± 0.03 Ga (Humayun et al. 2013 ) and a ¹⁴⁷Sm‐¹⁴³Nd isochron for NWA 7034 yielding 4.42 ± 0.07 Ga (Nyquist et al. 2016 ) defines the crystallization age of all its igneous portions. The zircon from the monzonitic rocks has a higher Δ¹⁷O than other Martian meteorites explained in part by assimilation of regolith materials enriched during surface alteration (Nemchin et al. 2014 ). This record of protolith interaction with atmosphere‐hydrosphere during regolith formation before melting demonstrates a thin atmosphere, a wet early surface environment on Mars, and an evolved crust likely to have contaminated younger extrusive rocks. The latest events recorded when the breccia was on Mars are resetting of apatite, much feldspar and some zircons at 1.35–1.4 Ga (Bellucci et al. 2015 ), and formation of Ni‐bearing pyrite veins during or shortly after this disturbance (Lorand et al. 2015 ).     

3D joint inversion of magnetotelluric and airborne tipper data: a case study from the Morrison porphyry Cu–Au–Mo deposit,British Columbia,Canada

Z‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.     

Swash saturation: an assessment of available models

An extensive previously published (Hughes et al. Mar Geol 355, 88–97, 2014) field data set representing the full range of micro-tidal beach states (reflective, intermediate and dissipative) is used to investigate swash saturation. Two models that predict the behavior of saturated swash are tested: one driven by standing waves and the other driven by bores. Despite being based on entirely different premises, they predict similar trends in the limiting (saturated) swash height with respect to dependency on frequency and beach gradient. For a given frequency and beach gradient, however, the bore-driven model predicts a larger saturated swash height by a factor 2.5. Both models broadly predict the general behavior of swash saturation evident in the data, but neither model is accurate in detail. While swash saturation in the short-wave frequency band is common on some beach types, it does not always occur across all beach types. Further work is required on wave reflection/breaking and the role of wave-wave and wave-swash interactions to determine limiting swash heights on natural beaches.     

Triangulation hand‐held laser‐scanning (TriHaLaS) for micro‐ and meso‐habitat surveys in streams

A novel hand‐held laser‐based stream bed survey system is presented. The system facilitates the capture of detailed 3D mapping of shallow (< 0.7 m) riverbed topography in sections approximately 4 m by 2 m. The system includes a stationary reference system, which projects three laser sheets (two at offset angles), within which a hand‐held monitoring pole is moved. The unique configuration of the light sheets intercepts with the pole as it moves within the survey area providing an exact horizontal location. Pole tilt is compensated for by an inertial measurement unit on the pole, and the height above the bed of the pole and submerged scanning laser are monitored relative to the horizontal laser sheet. Verification and application measurements demonstrate high resolution and accuracy in the horizontal (~5 mm) and vertical (~1 mm) direction. The system can be applied at sites where a free view is blocked and other optical through‐water methods fail. It is appropriate for studies on riverbed statistics and dynamics, which necessitate non‐invasive in‐situ surveys. Copyright © 2017 John Wiley & Sons, Ltd.     

Revisiting “An Exercise in Groundwater Model Calibration and Prediction” After 30 Years: Insights and New Directions

In 1988, an important publication moved model calibration and forecasting beyond case studies and theoretical analysis. It reported on a somewhat idyllic graduate student modeling exercise where many of the system properties were known; the primary forecasts of interest were heads in pumping wells after a river was modified. The model was calibrated using manual trial-and-error approaches where a model's forecast quality was not related to how well it was calibrated. Here, we investigate whether tools widely available today obviate the shortcomings identified 30 years ago. A reconstructed version of the 1988 true model was tested using increasing parameter estimation sophistication. The parameter estimation demonstrated the inverse problem was non-unique because only head data were available for calibration. When a flux observation was included, current parameter estimation approaches were able to overcome all calibration and forecast issues noted in 1988. The best forecasts were obtained from a highly parameterized model that used pilot points for hydraulic conductivity and was constrained with soft knowledge. Like the 1988 results, however, the best calibrated model did not produce the best forecasts due to parameter overfitting. Finally, a computationally frugal linear uncertainty analysis demonstrated that the single-zone model was oversimplified, with only half of the forecasts falling within the calculated uncertainty bounds. Uncertainties from the highly parameterized models had all six forecasts within the calculated uncertainty. The current results outperformed those of the 1988 effort, demonstrating the value of quantitative parameter estimation and uncertainty analysis methods.