Capture effects in carbonaceous material: A Stardust analogue study

Abstract— It is reasonable to expect that cometary samples returned to Earth by the Stardust space probe have been altered to some degree during capture in aerogel at 6.1 km/s. In order to help interpret the measured structure of these particles with respect to their original cometary nature, a series of coal samples of known structure and chemical composition was fired into aerogel at Stardust capture velocity. This portion of the study analyzed the surfaces of aerogel‐embedded particles using Raman spectroscopy. Results show that particle surfaces are largely homogenized during capture regardless of metamorphic grade or chemical composition, apparently to include a devolatilization step during capture processing. This provides a possible mechanism for alteration of some aliphatic compound‐rich phases through devolatilization of cometary carbonaceous material followed by re‐condensation within the particle. Results also show that the possibility of alteration must be considered for any particular Stardust grain, as examples of both graphitization and amorphization are found in the coal samples. It is evident that Raman G band (～1580 cm^?1) parameters provide a means of characterizing Stardust carbonaceous material to include identifying those grains which have been subjected to significant capture alteration.     

Stardust impact analogs: Resolving pre‐ and postimpact mineralogy in Stardust Al foils

Abstract– The grains returned by NASA’s Stardust mission from comet 81P/Wild 2 represent a valuable sample set that is significantly advancing our understanding of small solar system bodies. However, the grains were captured via impact at ～6.1 km s^?1 and have experienced pressures and temperatures that caused alteration. To ensure correct interpretations of comet 81P/Wild 2 mineralogy, and therefore preaccretional or parent body processes, an understanding of the effects of capture is required. Using a two‐stage light‐gas gun, we recreated Stardust encounter conditions and generated a series of impact analogs for a range of minerals of cometary relevance into flight spare Al foils. Through analyses of both preimpact projectiles and postimpact analogs by transmission electron microscopy, we explore the impact processes occurring during capture and distinguish between those materials inherent to the impactor and those that are the product of capture. We review existing and present additional data on olivine, diopside, pyrrhotite, and pentlandite. We find that surviving crystalline material is observed in most single grain impactor residues. However, none is found in that of a relatively monodisperse aggregate. A variety of impact‐generated components are observed in all samples. Al incorporation into melt‐derived phases allows differentiation between melt and shock‐induced phases. In single grain impactor residues, impact‐generated phases largely retain original (nonvolatile) major element ratios. We conclude that both surviving and impact‐generated phases in residues of single grain impactors provide valuable information regarding the mineralogy of the impacting grain whilst further studies are required to fully understand aggregate impacts and the role of subgrain interactions during impact.     

Subarctic catchment water storage and carbon cycling – Leading the way for future studies using integrated datasets at Pallas,Finland

Subarctic ecohydrological processes are changing rapidly, but detailed and integrated ecohydrological investigations are not as widespread as necessary. We introduce an integrated research catchment site (Pallas) for atmosphere, ecosystems, and ecohydrology studies in subarctic conditions in Finland that can be used for a new set of comparative catchment investigations. The Pallas site provides unique observational data and high-intensity field measurement datasets over long periods. The infrastructure for atmosphere- to landscape-scale research in ecosystem processes in a subarctic landscape has recently been complemented with detailed ecohydrological measurements. We identify three dominant processes in subarctic ecohydrology: (a) strong seasonality drives ecohydrological regimes, (b) limited dynamic storage causes rapid stream response to water inputs (snowmelt and intensive storms), and (c) hydrological state of the system regulates catchment-scale dissolved carbon dynamics and greenhouse (GHG) fluxes. Surface water and groundwater interactions play an important role in regulating catchment-scale carbon balances and ecosystem respiration within subarctic peatlands, particularly their spatial variability in the landscape. Based on our observations from Pallas, we highlight key research gaps in subarctic ecohydrology and propose several ways forward. We also demonstrate that the Pallas catchment meets the need for sustaining and pushing the boundaries of critical long-term integrated ecohydrological research in high-latitude environments.     

Spatially distributed modelling of soil erosion and sediment yield at regional scales in Spain

Initiated by the need to quantify erosion rates and the impacts of global changes on erosion, several attempts have been made to apply erosion models at regional scales. However, these models have often been directed towards on-site soil erosion estimates, emphasising sheet and rill erosion processes, and disregarding gully erosion, channel erosion and sediment transport. These models are therefore of limited use for the assessment of sediment yield, off-site impacts of erosion, and for the development of environmental management to control these impacts at regional scale. This study analyses and compares three spatially distributed models for the prediction of soil erosion and/or sediment yield at regional scales: the WATEM-SEDEM model that is based on the empirical Revised Universal Soil Loss Equation (RUSLE) in combination with a sediment transport equation, the physics-based Pan European Soil Erosion Risk Assessment model (PESERA), and a newly developed Spatially Distributed Scoring model (SPADS). The three models were applied to 61 Spanish drainage basins and model predictions were evaluated against data on measured reservoir sedimentation rates. Global data sets on land use, climate, elevation and soil characteristics were used as model input for WATEM-SEDEM and SPADS, whereas published soil erosion estimates of PESERA at 1 km² resolution were used directly. SPADS and WATEM-SEDEM provided best results after separate calibration for basins with a Sediment Delivery Ratio (SDR) higher than 5% and those with an SDR lower than 5%. In this way, SPADS explained 67% of variation in sediment yield, while WATEM-SEDEM explained 48% of the variation. PESERA represents a promising alternative to the use of empirical models at the regional scale as it can be applied to very diverse environments with little calibration. However, PESERA provides soil erosion rates and not sediment yield estimates. For most basins PESERA soil erosion rates vary between fifty and close to zero percent of total sediment yield. Two major factors may explain this discrepancy between modelled soil erosion rates and measured sediment yield. First, it may be that PESERA underestimates soil erosion under Mediterranean conditions, although PESERA soil erosion rates are of the same order of magnitude as erosion rates measured in erosion plot studies. Second, gully-, river channel erosion and sediment transport processes may be much more important than sheet- and rill erosion for regional scale sediment yield in these environments. These issues therefore require further attention in future model development. Although spatially lumped models provide better predictions of sediment yield at the basin scale, and while validation of the predicted spatial patterns of sources and sinks of sediment requires further research, spatially distributed models are expected to be of value to support management decisions regarding the assessment of on-site and off-site impacts of erosion at the regional scale.     

Experiments on the effect of inflow and outflow sequences on suspended sediment exchange rates

In laboratory experiments, the influence of inflow and outflow sequences on the behavior of fine sedi-ment was investigated. The experimental set-up consisted of two interconnected rectangular basins, between which water was moved back and forth. Suspended sediment concentration in the main basin as well as the sediment exchange rates were derived from turbidity measurements.The suspended sediment ratio, SSR, and sediment exchange rates (influx sediment rate, ISR, and evacuated sediment rate, ESR) were measured. In twenty test runs, a parametric study on the magnitude and frequency of inflow and outflow cycles, the relative duration between inflow and outflow sequences, the initial sediment concentration, and the intake position was done. An initial test with stagnant water described the set-tling behavior of fine sediment and served as a reference scenario.The test results show that settling of fine particles near the intake/outlet structure can be considerably reduced by the nature of the inflow and outflow sequences. High cycle magnitude and frequency lead to maximum suspended sediment ratio in the system. For low discharges, the evolution of suspended sediment concentration cannot be directly correlated to the inflow and outflow cycles. However, compared to"no operation"conditions, the suspended sediment ratio could be increased by 10%to 40%locally. For high discharge, the evolution of suspended sediment concentration correlated with discharge cycles and suspended sediment ratios between 50%and 80%higher than for stagnant water could be achieved. Similar ratios could be obtained when the intake is located closer to the bottom or to the free water surface.Meanwhile, the overall sediment balance remained in equilibrium over the test period, indicating that the influx and evacuated sediment rates are not significantly influenced by the inflow and outflow cycles.     

Approximate traveltime inversion in downhole microseismic monitoring

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.     

Saltwater Upconing Due to Cyclic Pumping by Horizontal Wells in Freshwater Lenses

This article deals with the quantification of saltwater upconing below horizontal wells in freshwater lenses using analytical solutions as a computationally fast alternative to numerical simulations. Comparisons between analytical calculations and numerical simulations are presented regarding three aspects: (1) cyclic pumping; (2) dispersion; and (3) finite horizontal wells in a finite domain (a freshwater lens). Various hydrogeological conditions and pumping regimes within a dry half year are considered. The results show that the influence of elastic and phreatic storage (which are not taken into account in the analytical solutions) on the upconing of the interface is minimal. Furthermore, the analytical calculations based on the interface approach compare well with numerical simulations as long as the dimensionless interface upconing is below 1/3, which is in line with previous studies on steady pumping. Superimposing an analytical solution for mixing by dispersion below the well over an analytical solution based on the interface approach is appropriate in case the vertical flow velocity around the interface is nearly constant but should not be used for estimating the salinity of the pumped groundwater. The analytical calculations of interface upconing below a finite horizontal well compare well with the numerical simulations in case the distance between the horizontal well and the initial interface does not vary significantly along the well and in case the natural fluctuation of the freshwater lens is small. In order to maintain a low level of salinity in the well during a dry half year, the dimensionless analytically calculated interface upconing should stay below 0.25.     

Isotopic Composition of Waters in the Dem?novská dolina Valley and Its Underground Hydrologic System During Winter and Spring of 2010/2011

In the Dem?novská dolina Cave system (Slovakia) and its vicinity, 32 sampling places for regular observation (in 2-months interval) of δ18O and δ2H in water were established. This monitoring included precipitation waters, waters in the surface streams, waters of the underground hydrological system as well as the dripping seepage waters of the cave system. Altitudinal extent of the area was from 800 m a.s.l. (lowermost cave entrance) to 2024 m a.s.l. (Chopok Mt. on the top of the crystalline range). Initial results show some similarities but also differences within the analyzed water types. For precipitation, a high variability of isotopic composition was confirmed, from quite depleted up to more enriched waters (δ18O from -16.8‰ up to -5.7‰; δ2H from -121.6‰ to -32.7‰). During the recharge process and groundwater/surface water formation, precipitation water is homogenized, what is reflected in much more stable isotope content. The most depleted (δ18O ≈ -11.7‰ to -10.8‰; δ2H ≈ -78.9‰ to -73.4‰) were the waters of surface streams, running from the northern slopes of the Nízke Tatry Mts., formed by crystalline rocks, alochtonous to the under?ground hydrological system. Smaller autochtonous surface water streams (formed in the side valleys of the main karstic canyon) are slightly enriched (heavier, as δ18O ≈ -11.4‰ to -10.6‰; δ2H ≈ -78.3‰ to -71.5‰), what reflects lower altitudes of their watersheds. Interesting is the distribution of the isotope content of the underground streams in the cave system. The most depleted are the underground streams directly (visibly) communicating with surface waters (δ18O≈-11.33±0.13‰; δ2H≈-76.88±1.68‰). Extent of the relationship of underground streams to the autochtonous seepage waters (slow circulation through the fissures) is manifested by respective portion of iso?topically enriched waters–as the underground streams show different isotope composition. The combination of the alochtonous water components (from surface streams reaching the karstic area from the adjacent crystalline via swallow holes) and autochtonous water components (recharged directly in karstified limestones) is visible especially on the subsurface stream of Dem?novka. The most isotopically enriched (heaviest) of all water types are dripping seepage waters (δ18O ≈ -10.4‰ to -9.4‰; δ2H ≈ -71.6‰ to -65.0‰).     

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.