Mass transport around comets and its impact on the seasonal differences in water production rates

Comets are surrounded by a thin expanding atmosphere, and although the nucleus' gravity is small, some molecules and grains, possibly with the inclusion of ices, can get transported around the nucleus through scattering (atoms/molecules) and gravitational pull (grains). Based on the obliquity of the comet, it is also possible that volatile material and icy grains get trapped in regions, which are in shadow until the comet passes its equinox. When the Sun rises above the horizon and the surface starts to heat up, this condensed material starts to desorb and icy grains will sublimate off the surface, possibly increasing the comet's neutral gas production rate on the outbound path. In this paper we investigate the mass transport around the nucleus, and based on a simplified model, we derive the possible contribution to the asymmetry in the seasonal gas production rate that could arise from trapped material released from cold areas once they come into sunlight. We conclude that the total amount of volatiles retained by this effect can only contribute up to a few percent of the asymmetry observed in some comets.

Cosmogenic production rates and recoil loss effects in micrometeorites and interplanetary dust particles

We present a purely physical model to determine cosmogenic production rates for noble gases and radionuclides in micrometeorites (MMs) and interplanetary dust particles (IDPs) by solar cosmic-rays (SCR) and galactic cosmic-rays (GCR) fully considering recoil loss effects. Our model is based on various nuclear model codes to calculate recoil cross sections, recoil ranges, and finally the percentages of the cosmogenic nuclides that are lost as a function of grain size, chemical composition of the grain, and the spectral distribution of the projectiles. The main advantage of our new model compared with earlier approaches is that we consider the entire SCR particle spectrum up to 240 MeV and not only single energy points. Recoil losses for GCR-produced nuclides are assumed to be equal to recoil losses for SCR-produced nuclides. Combining the model predictions with Poynting-Robertson orbital lifetimes, we calculate cosmic-ray exposure ages for recently studied MMs, cosmic spherules, and IDPs. The ages for MMs and the cosmic-spherule are in the range <2.2–233 Ma, which corresponds, according to the Poynting-Robertson drag, to orbital distances in the range 4.0–34 AU. For two IDPs, we determine exposure ages of longer than 900 Ma, which corresponds to orbital distances larger than 150 AU. The orbital distance in the range 4–6 AU for one MM and the cosmic spherule indicate an origin either in the asteroid belt or release from comets coming either from the Kuiper Belt or the Oort Cloud. Three of the studied MMs have orbital distances in the range 23–34 AU, clearly indicating a cometary origin, either from short-period comets from the Kuiper Belt or from the Oort Cloud. The two IDPs have orbital distances of more than 150 AU, indicating an origin from Oort Cloud comets.

Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko

Context. The complex shape of comet 67P and its oblique rotation axis cause pronounced seasonal effects. Irradiation and hence activity vary strongly. Aims. We investigate the insolation of the cometary surface in order to predict the sublimation of water ice. The strongly varying erosion levels are correlated with the topography and morphology of the present cometary surface and its evolution. Methods. The insolation as a function of heliocentric distance and diurnal (spin dependent) variation is calculated using >10(5) facets of a detailed digital terrain model. Shading, but also illumination and thermal radiation by facets in the field of view of a specific facet are iteratively taken into account. We use a two-layer model of a thin porous dust cover above an icy surface to calculate the water sublimation, presuming steady state and a uniform surface. Our second model, which includes the history of warming and cooling due to thermal inertia, is restricted to a much simpler shape model but allows us to test various distributions of active areas. Results. Sublimation from a dirty ice surface yields maximum erosion. A thin dust cover of 50 pm yields similar rates at perihelion. Only about 6% of the surface needs to be active to match the observed water production rates at perihelion. A dust layer of 1 mm thickness suppresses the activity by a factor of 4 to 5. Erosion on the south side can reach more than 10 m per orbit at active spots. The energy input to the concave neck area (Hapi) during northern summer is enhanced by about 50% owing to self-illumination. Here surface temperatures reach maximum values along the foot of the Hathor wall. Integrated over the whole orbit this area receives the least energy input. Based on the detailed shape model, the simulations identify "hot spots" in depressions and larger pits in good correlation with observed dust activity. Three-quarters of the total sublimation is produced while the sub-solar latitude is south, resulting in a distinct dichotomy in activity and morphology. Conclusions. The northern areas display a much rougher morphology than what is seen on Imhotep, an area at the equator that will be fully illuminated when 67P is closer to the Sun. Self-illumination in concave regions enhance the energy input and hence erosion. This explains the early activity observed at Hapi. Cliffs are more prone to erosion than horizontal, often dust covered, areas, which leads to surface planation. Local activity can only persist if the forming cliff walls are eroding. Comet 67P has two lobes and also two distinct sides. Transport of material from the south to the north is probable. The morphology of the Imhotep plain should be typical for the terrains of the yet unseen southern hemisphere.

The dust environment of comet 67P/Churyumov-Gerasimenko from Rosetta OSIRIS and VLT observations in the 4.5 to 2.9 AU heliocentric distance range inbound

Context. The ESA Rosetta spacecraft, currently orbiting around cornet 67P/Churyumov-Gerasimenko, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. Aims. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. Methods. A Monte Carlo dust tail code, which has already been used to characterise the dust environments of several comets and active asteroids, has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from. Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. Results. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is alpha = -3, instead of alpha = -2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 AU, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 AU. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

Modeling the marine aragonite cycle: changes under rising carbon dioxide and its role in shallow water CaCO3 dissolution

The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by 29% and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.

Calibration of cosmogenic noble gas production in ordinary chondrites based on 36 Cl- 36 Ar ages. Part 1: Refined produced rates for cosmogenic 21 Ne and 38 Ar

We measured the concentrations and isotopic compositions of He, Ne, and Ar in bulk samples and metal separates of 14 ordinary chondrite falls with long exposure ages and high metamorphic grades. In addition, we measured concentrations of the cosmogenic radionuclides 10Be, 26Al, and 36Cl in metal separates and in the nonmagnetic fractions of the selected meteorites. Using cosmogenic 36Cl and 36Ar measured in the metal separates, we determined 36Cl-36Ar cosmic-ray exposure (CRE) ages, which are shielding-independent and therefore particularly reliable. Using the cosmogenic noble gases and radionuclides, we are able to decipher the CRE history for the studied objects. Based on the correlation 3He/21Ne versus 22Ne/21Ne, we demonstrate that, among the meteorites studied, only one suffered significant diffusive losses (about 35%). The data confirm that the linear correlation 3He/21Ne versus 22Ne/21Ne breaks down at high shielding. Using 36Cl-36Ar exposure ages and measured noble gas concentrations, we determine 21Ne and 38Ar production rates as a function of 22Ne/21Ne. The new data agree with recent model calculations for the relationship between 21Ne and 38Ar production rates and the 22Ne/21Ne ratio, which does not always provide unique shielding information. Based on the model calculations, we determine a new correlation line for 21Ne and 38Ar production rates as a function of the shielding indicator 22Ne/21Ne for H, L, and LL chondrites with preatmospheric radii less than about 65 cm. We also calculated the 10Be/21Ne and 26Al/21Ne production rate ratios for the investigated samples, which show good agreement with recent model calculations.

Helium transfer from water into quartz crystals: A new approach for porewater dating

Several important fundamental and applied problems require a quantification of slow rates of groundwater flow. To resolve these problems helium appears to be a promising tracer. In this contribution we discuss a new approach, which gives the helium inventory in a rock – pore water system by using the relevant mineral record, i.e., without extraction and investigation of the porewater samples. Some U- and Th-poor minerals such as quartz (quartz separates from Permo-Carboniferous Formation, sandstone–shale interlayering, Molasses Basin, Northern Switzerland, hereafter PCF, are used in this study) contain excessive helium having migrated into their internal helium-accessible volume (HAV) from the surrounding porewater [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. These volumes are estimated by using helium as a nano-size penetrating tool, i.e., by saturation of the minerals with helium under controlled pressure–temperature conditions and subsequent measurements of the helium-saturated concentrations. In the quartz separates HAV/total volume ratios vary from 0.017% to 0.16%; along with the measured initial (unsaturated) He concentration the HAV gives the internal helium pressure, the mean value obtained for 7 samples (25 sample aliquots) is P=0.45F0.15 atm (1 r). The product of helium pressure and solubility (7.35_10_3 cc STP He/cc H2O for the temperature and salinity of PCF aquifers reported in [F.J. Pearson, W. Balderer, H.H. Loosli, B.E. Lehmann, A. Matter, T. Peters, H. Schmassmann, A. Gautschi, Applied Isotope Hydrogeology–A Case Study in Northern Switzerland, Elsevier Amsterdam, 1991, 439 pp.]) is the mineral-derived He concentration in the respective porewater, CPW=0.0035F0.0017 cc He/cc H2O. This value is in full accord with measured He concentrations in PCF aquifers, CPCF, varying from 0.0045 to 0.0016 cc He/cc H2O. This agreement validates the proposed approach and also shows that the mineral–porewater helium–concentration equilibrium has been established. Indeed, estimates of the He-migration rates through our quartz samples show that in ~6000 years the internal pressure should equilibrate with He-concentration in related porewater of PCF, and this time interval is short compared to independent estimates [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. The helium inventory in the rock–porewater assemblage shows that helium abundance in pore waters is indeed important. In shale samples (with ~15% porosity) porewaters contain more helium than the host minerals altogether. Porewater heliumconcentration profiles, available from the mineral record, along with helium production rates are input parameters allowing model(s) of helium migration through a hydrological structure to be developed. Quite high helium concentrations in PCF porewaters imply slow removal mechanisms, which will be discussed elsewhere.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

A relation between screening masses and real-time rates

Thermal screening masses related to the conserved vector current are determined for the case that the current carries a non-zero Matsubara frequency, both in a weak-coupling approach and through lattice QCD. We point out that such screening masses are sensitive to the same infrared physics as light-cone real-time rates. In particular, on the perturbative side, the inhomogeneous Schrödinger equation determining screening correlators is shown to have the same general form as the equation implementing LPM resummation for the soft-dilepton and photon production rates from a hot QCD plasma. The static potential appearing in the equation is identical to that whose soft part has been determined up to NLO and on the lattice in the context of jet quenching. Numerical results based on this potential suggest that screening masses overshoot the free results (multiples of 2πT) more strongly than at zero Matsubara frequency. Four-dimensional lattice simulations in two-flavour QCD at temperatures of 250 and 340 MeV confirm the non-static screening masses at the 10% level. Overall our results lend support to studies of jet quenching based on the same potential at T ≳ 250 MeV.