We measured maps of atmospheric water (H2O) and its deuterated form (HDO) across the martian globe, showing strong isotopic anomalies and a significant high deuterium/hydrogen (D/H) enrichment ...indicative of great water loss. The maps sample the evolution of sublimation from the north polar cap, revealing that the released water has a representative D/H value enriched by a factor of about 7 relative to Earth's ocean Vienna standard mean ocean water (VSMOW). Certain basins and orographic depressions show even higher enrichment, whereas high-altitude regions show much lower values (1 to 3 VSMOW). Our atmospheric maps indicate that water ice in the polar reservoirs is enriched in deuterium to at least 8 VSMOW, which would mean that early Mars (4.5 billion years ago) had a global equivalent water layer at least 137 meters deep.
ABSTRACT
The observed rate of water production of comet 67P/Churyumov–Gerasimenko near its perihelion can be approximated by a very steep power function of the heliocentric distance. Widely used ...thermophysical models based on a static dust layer on top of the icy/refractory matrix are poorly consistent with these observations. We analyse published model results and demonstrate that thermophysical models with a uniform and static ice free layer do not reproduce the observed steep water production rates of 67P near perihelion. Based on transient thermal modeling we conclude that the accelerated gas activity can be explained assuming that the active area fraction near the south pole is increased. The deeper penetration of the heat wave during polar day (no sunset) can activate sublimation through thicker inert dust layers. This can also lead to removal of thicker dust layers and consequently to an expansion of the active area.
We present the analysis of ≈100 molecular maps of the coma of comet 67P/Churyumov-Gerasimenko that were obtained with the MIRO submillimeter radiotelescope on board the Rosetta spacecraft. From the ...spectral line mapping of H 216 $_2^{16}$ 216 O, H 218 $_2^{18}$218 O, H 217 $_2^{17}$217 O, CH3OH, NH3, and CO and some fixed nadir pointings, we retrieved the outgassing pattern and total production rates for these species. The analysis covers the period from July 2014, inbound to perihelion, to June 2016, outbound, and heliocentric distances rh = 1.24–3.65 AU. A steep evolution of the outgassing rates with heliocentric distance is observed, typically in rh−16 $r_{\textrm{h}}^{-6}$rh−6 , with significant differences between molecules (e.g. steeper variation for H2O post-perihelion than for methanol). As a consequence, the abundances relative to water in the coma vary. The CH3OH and CO abundances increase after perihelion, while the NH3 abundance peaks around perihelion and then decreases. Outgassing patterns have been modeled as 2D Gaussian jets. The width of these jets is maximum around the equinoxes when the bulk of the outgassing is located near the equator. From July 2014 to February 2015, the outgassing is mostly restricted to a narrower jet (full width at half-maximum ≈80°) originating from high northern latitudes, while around perihelion, most of the gaseous production comes from the southernmost regions ( − 80 ± 5° cometocentric latitude) and forms a 100°–130° (full width at half-maximum) wide fan. We find a peak production of water of 0.8 × 1028 molec. s−1, 2.5 times lower than measured by the ROSINA experiment, and place an upper limit to a 50% additional production that could come from the sublimation of icy grains. We estimate the total loss of ices during this perihelion passage to be 4.18 ± 0.18 × 109 kg. We derive a dust-to-gas ratio in the lost material of 0.7–2.3 (including all sources of errors) based on the nucleus mass loss of 10.5 ± 3.4 × 109 kg estimated by the RSI experiment. We also obtain an estimate of the H 218 $_2^{18}$218 O/H 217 $_2^{17}$217 O ratio of 5.6 ± 0.8.
Aims. We investigate the influence of three basic factors on water production rate as a function of heliocentric distance: nucleus shape, the spin axis orientation, and the distribution of activity ...on a comet’s surface. Methods. We used a basic water sublimation model driven by solar insolation to derive total production rates for different nuclei shapes and spin axis orientations using the orbital parameters of 67P/Churyumov-Gerasimenko. We used known shape models derived from prior missions to the Jupiter Family and short period comets. The slopes of production rates versus heliocentric distance were calculated for the different model setups. Results. The standard (homogeneous) outgassing model confirms the well-known result regarding the heliocentric dependence of water production rate that remains invariant for different nuclei shapes as long as the rotation axis is perpendicular to the orbital plane. When the rotation axis is not perpendicular, the nucleus shape becomes a critically important factor in determining the water production curves as the illuminated cross section of the nucleus changes with heliocentric distance. Shape and obliquity can produce changes in the illuminated cross section of up to 50% over an orbit. In addition, different spin axis orientations for a given shape can dramatically alter the pre- and post-perihelion production curves, as do assumptions about the activity distribution on the surface. If, however, the illuminated cross section of the nucleus is invariant, then the dependence on the above parameters is weak, as demonstrated here with the 67P/Churyumov-Gerasimenko shape. The comets Hartley 2 and Wild 2 are shown to yield significantly different production curve shapes for the same orbit and orientation as 67P/CG, varying by as much as a factor of three as a result of only changing the nucleus shape. Finally, we show that varying just three basic parameters, shape, spin axis orientation, and active spots distribution on the surface can lead to arbitrary deviations from the expected inverse square law dependence of water production rates near 1 au. Conclusions. With the results obtained, we cannot avoid the conclusion that, without prior knowledge of basic parameters (shape, spin axis orientation, activity locations), it is difficult to reveal the nature of cometary outgassing from the heliocentric water production rates. Similarly, the inter-comparison of water production curves of two such comets may not be meaningful.
Aims. Using spectroscopic and continuum data measured by the MIRO instrument on board Rosetta of comet 67P/Churyumov-Gerasimenko, it is possible to derive and track the change in the water production ...rate, to learn how the outgassing evolves with heliocentric distance. The MIRO data are well suited to investigate the evolution of 67P, in unprecedented spatial and temporal detail. Methods. To obtain estimates of the local effective Haser production rates we developed an efficient and reliable retrieval approach with precalculated lookup tables. We employed line area ratios (H216O/H218O) from pure nadir observations as the key variable, along with the Doppler shift velocity, and continuum temperature. This method was applied to the MIRO data from August 2014 until April 2016. Perihelion occurred on August 13, 2015 when the comet was 1.24 AU from the Sun. Results. During the perihelion approach, the water production rates increased by an order of magnitude, and from the observations, the derived maximum for a single observation on August 29, 2015 is (1.42 ± 0.51) ×1028. Modeling the data indicates that there is an offset in the peak outgassing, occurring 34 ± 10 days after perihelion. During the pre-perihelion phase, the production rate changes with heliocentric distance as rh−3.8±0.2; during post-perihelion, the dependence is rh−4.3±0.2. The comet is calculated to have lost 0.12 ± 0.06 % of its mass during the perihelion passage, considering only water ice sublimation. Additionally, this method provides well sampled data to determine the spatial distribution of outgassing versus heliocentric distance. The time evolution is definitely not uniform across the surface. Pre- and post-perihelion, the surface temperature on the southern hemisphere changes rapidly, as does the sublimation rate with an exponent of ~−6. There is a strong latitudinal dependence on the rh exponent with significant variation between northern and southern hemispheres, and so the average over the comet surface may only be of limited importance. We present more detailed regional variation in the outgassing, demonstrating that the highest derived production rates originate from the Wosret, Neith and Bes regions during perihelion.
We describe the design and construction of GREAT (German REceiver for Astronomy at Terahertz frequencies) operated on the Stratospheric Observatory For Infrared Astronomy (SOFIA). GREAT is a modular ...dual-color heterodyne instrument for high-resolution far-infrared (FIR) spectroscopy. Selected for SOFIA’s Early Science demonstration, the instrument has successfully performed three Short and more than a dozen Basic Science flights since first light was recorded on its April 1, 2011 commissioning flight. We report on the in-flight performance and operation of the receiver that – in various flight configurations, with three different detector channels – observed in several science-defined frequency windows between 1.25 and 2.5 THz. The receiver optics was verified to be diffraction-limited as designed, with nominal efficiencies; receiver sensitivities are state-of-the-art, with excellent system stability. The modular design allows for the continuous integration of latest technologies; we briefly discuss additional channels under development and ongoing improvements for Cycle 1 observations. GREAT is a principal investigator instrument, developed by a consortium of four German research institutes, available to the SOFIA users on a collaborative basis.
ABSTRACT
A hierarchical model of microstructure of cometary dust seems to accurately capture the morphological complexity of these particles as observed by the Rosetta mission to the ...67P/Churyumov–Gerasimenko. The main aim of this work is to investigate how uncertainties in the knowledge of the microstructure of the surface layer affect our estimates of gas production. New models that incorporate scale of inhomogenieties in a sound theoretical framework are used for hierarchical dust layers. The studied layers are constructed in two steps: first we design ballistic aggregates as the building units and then, using these porous blocks, assemble the layers constrained by the known porosity range. The mean pore size and permeability are studied. Modelling is performed for various values of porosity, grain size, and layer thickness. The simulation results are embedded in the thermal model, explicitly including a radiation thermal conductivity and a resistance of the dust layer for the gas flow. It is shown that the average pore size is satisfactorily approximated by a linear function of the ratio of the effective porosity of the layer to the degree of filling. Simple fittings were obtained for the layer permeability. Our results indicate that in the expected range of nucleus porosity, the gas production rate is weakly dependent on the detailed layer microstructure, and appropriate effective values of homogeneous dust layers can be used to evaluate the gas activity. We also note that adding complex elements into the model yields unavoidable statistical uncertainties within several tens of per cent.
Aims.
The aim of this study is to measure the vertical distribution of HCN on Titan’s stratosphere using ground-based submillimetre observations acquired quasi-simultaneously with the
Herschel
ones. ...This allows us to perform a consistency check between space and ground-based observations and to build a reference mean HCN vertical profile in Titan’s stratosphere.
Methods.
Using APEX and IRAM 30-m, we obtained the spectral emission of HCN (4-3) and (3-2) lines. Observations were reduced with GILDAS-CLASS. We applied a line-by-line radiative transfer code to calculate the synthetic spectra of HCN, and a retrieval algorithm based on optimal estimation to retrieve the temperature and HCN vertical distributions. We used the standard deviation-based metric to quantify the dispersion between the ground-based and
Herschel
HCN profiles and the mean one.
Results.
Our derived HCN abundance profiles are consistent with an increase from 40 ppb at ~100 km to 4 ppm at ~200 km, which is an altitude region where the HCN signatures are sensitive. We also demonstrate that the retrieved HCN distribution is sensitive to the data information and is restricted to Titan’s stratosphere. The HCN obtained from APEX data is less accurate than the one from IRAM data because of the poorer data quality, and covers a narrower altitude range. Comparisons between our results and the values from
Herschel
show similar abundance distributions, with maximum differences of 2.5 ppm ranging between 100 and 300 km in the vertical range. These comparisons also allow us to inter-validate both data sets and indicate reliable and consistent measurements. The inferred abundances are also consistent with the vertical distribution in previous observational studies, with the profiles from ALMA, Cassini/CIRS, and SMA (the latest ones below ~230 km). Our HCN profile is also comparable to photochemical models by Krasnopolsky (2014) and Vuitton et al. (2019) below 230 km and consistent with that of Loison et al. (2015) above 250 km. However, it appears to show large differences with respect to the estimates by Loison et al. (2015), Dobrijevic & Loison (2018), and Lora et al. (2018) below 170 km, and by Dobrijevic & Loison (2018) and Lora et al. (2018) above 400 km, although they are similar in shape. We conclude that these particular photochemical models need improvement.
ABSTRACT
Explanation of the observed gas activity based on the parameters of the comet’s nucleus is not an obvious matter. Solutions based on certain thermal models can be obtained, but the use of ...arbitrary values for poorly known model parameters is always required. In this work, we carry out a study of the dependence of gas activity on these parameters and assess the importance of each of them. We consider model porous dust layers of diverse structures. Solid state and gas thermal conductivities are examined to estimate the possible range of effective thermal conductivity. The simulation results are embedded in a thermal model, explicitly including a radiative thermal conductivity and a resistance of the dust layer against the gas flow. Sublimation of water ice and supervolatiles (CO2 and CO) is tested at different heliocentric distances. It is shown that when sublimation is the main energy sink, the role of uncertainties in the structure of the layer is small. As the relative contribution of sublimation decreases, the scatter of solutions reaches tens of per cent. The expected large uncertainties in the effective thermal conductivity can also significantly change the gas production. The analysis performed shows that, in spite of the parameter range having been narrowed down by results from the Rosetta mission, the unavoidable uncertainty in the values of some model parameters (e.g. thermal conductivity) blurs the theoretical simulation estimates. Instead of presenting a narrow set of specific solutions, it is desirable to analyse the entire range of possible solutions.
Aims. Using data from the Rosetta mission to comet 67P/Churyumov–Gerasimenko, we evaluate the physical properties of the surface and subsurface of the nucleus and derive estimates for the thermal ...inertia (TI) and roughness in several regions on the largest lobe of the nucleus. Methods. We have developed a thermal model to compute the temperature on the surface and in the uppermost subsurface layers of the nucleus. The model takes heat conduction, self-heating, and shadowing effects into account. To reproduce the brightness temperatures measured by the MIRO instrument, the thermal model is coupled to a radiative transfer model to derive the TI. To reproduce the spatially resolved infrared measurements of the VIRTIS instrument, the thermal model is coupled to a radiance model to derive the TI and surface roughness. These methods are applied to Rosetta data from September 2014. Results. The resulting TI values from both instruments are broadly consistent with each other. From the millimetre channel on MIRO, we determine the TI in the subsurface to be <80 JK−1 m−2 s−0.5 for the Seth, Ash, and Aten regions. The submillimetre channel implies similar results but also suggests that higher values could be possible. A low TI is consistent with other MIRO measurements and in situ data from the MUPUS instrument at the final landing site of Philae. The VIRTIS results give a best-fitting value of 80 JK−1 m−2 s−0.5 and values in the range 40–160 JK−1 m−2 s−0.5 in the same areas. These observations also allow the subpixel scale surface roughness to be estimated and compared to images from the OSIRIS camera. The VIRTIS data imply that there is significant roughness on the infrared scale below the resolution of the available shape model and that, counter-intuitively, visually smooth terrain (centimetre scale) can be rough at small (micrometre–millimetre) scales, and visually rough terrain can be smooth at small scales.