Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.     

Glacier inventory of the Gran Campo Nevado Ice Cap in the Southern Andes and glacier changes observed during recent decades

The Gran Campo Nevado (GCN) forms an isolated ice cap on the Península Muñoz Gamero (PMG) located 200 km to the south of the Southern Patagonia Icefield (SPI). We present a glacier inventory of the GCN made up by 27 drainage basins (in total 199.5 km²) and other small cirque and valley glaciers of the southern part of PMG (in total 53 km²). The glacier inventory is based on a digital elevation model (DEM) and ortho-photos. Contour lines from maps, relief information derived from Landsat TM satellite imagery from 1986 and 2002 and stereoscopic data from aerial photos were combined in a knowledge-based scheme to obtain a DEM of the area. A digital ortho-photo map based on aerial photos from 1998 and several ortho-photos based on aerial photos from 1942 and 1984 could be produced from the initial DEM. A geographical information system (GIS) served to outline the extent of the present glaciation. All major glaciers of the GCN show a significant glacier retreat during the last 60 yr. Some of the outlet glaciers lost more than 20% of their total area during this period. Overall glacier retreat amounts to 2.8% of glacier length per decade and the glacier area loss is 2.4% per decade in the period from 1942 to 2002. We hypothesise that GCN glaciers may have reacted faster and more synchronously with the observed warming trend during recent decades when compared with the SPI.     

Oxygen and carbon isotope ratios in the martian atmosphere

Oxygen and carbon isotope ratios in the martian CO₂ are key values to study evolution of volatiles on Mars. The major problems in spectroscopic determinations of these ratios on Mars are uncertainties associated with: (1) equivalent widths of the observed absorption lines, (2) line strengths in spectroscopic databases, and (3) thermal structure of the martian atmosphere during the observation. We have made special efforts to reduce all these uncertainties. We observed Mars using the Fourier Transform Spectrometer at the Canada–France–Hawaii Telescope. While the oxygen and carbon isotope ratios on Mars were byproducts in the previous observations, our observation was specifically aimed at these isotope ratios. We covered a range of 6022 to 6308 cm⁻¹ with the highest resolving power of ν/δν=3.5×10⁵ and a signal-to-noise ratio of 180 in the middle of the spectrum. The chosen spectral range involves 475 lines of the main isotope, 184 lines of ¹³CO₂, 181 lines of CO¹⁸O, and 119 lines of CO¹⁷O. (Lines with strengths exceeding 10⁻²⁷ cm at 218 K are considered here.) Due to the high spectral resolution, most of the lines are not blended. Uncertainties of retrieved isotope abundances are in inverse proportion to resolving power, signal-to-noise ratio, and square root of the number of lines. Laboratory studies of the CO₂ isotope spectra in the range of our observation achieved an accuracy of 1% in the line strengths. Detailed observations of temperature profiles using MGS/TES and data on temperature variations with local time from two GCMs are used to simulate each absorption line at various heights in each part of the instrument field of view and then sum up the results. Thermal radiation of Mars' surface and atmosphere is negligible in the chosen spectral range, and this reduces errors associated with uncertainties in the thermal structure on Mars. Using a combination of all these factors, the highest accuracy has been achieved in measuring the CO₂ isotope ratios: ¹³C/¹²C = 0.978 ± 0.020 and ¹⁸O/¹⁶O = 1.018 ± 0.018 times the terrestrial standards. Heavy isotopes in the atmosphere are enriched by nonthermal escape and sputtering, and depleted by fractionation with solid-phase reservoirs. The retrieved ratios show that isotope fractionation between CO₂ and oxygen and carbon reservoirs in the solid phase is almost balanced by nonthermal escape and sputtering of O and C from Mars.     

The formation of Ganymede's grooved terrain: Numerical modeling of extensional necking instabilities

Ganymede's grooved terrain likely formed during an epoch of global expansion, when unstable extension of the lithosphere resulted in the development of periodic necking instabilities. Linear, infinitesimal-strain models of extensional necking support this model of groove formation, finding that the fastest growing modes of an instability have wavelengths and growth rates consistent with Ganymede's grooves. However, several questions remain unanswered, including how nonlinearities affect instability growth at large strains, and what role instabilities play in tectonically resurfacing preexisting terrain. To address these questions we numerically model the extension of an icy lithosphere to examine the growth of periodic necking instabilities over a broad range of strain rates and temperature gradients. We explored thermal gradients up to 45 K km⁻¹ and found that, at infinitesimal strain, maximum growth rates occur at high temperature gradients (45 K km⁻¹) and moderate strain rates (10⁻¹³ s⁻¹). Dominant wavelengths range from 1.8 to 16.4 km (post extension). Our infinitesimal growth rates are qualitatively consistent with, but an order of magnitude lower than, previous linearized calculations. When strain exceeds ∼10% growth rates decrease, limiting the total amount of amplification that can result from unstable extension. This fall-off in growth occurs at lower groove amplitudes for high-temperature-gradient, thin-lithosphere simulations than for low-temperature-gradient, thick-lithosphere simulations. At large strains, this shifts the ideal conditions for producing large amplitude grooves from high temperature gradients to more moderate temperature gradients (15 K km⁻¹). We find that the formation of periodic necking instabilities can modify preexisting terrain, replacing semi-random topography up to 100 m in amplitude with periodic ridges and troughs, assisting the tectonic resurfacing process. Despite this success, the small topographic amplification produced by our model presents a formidable challenge to the necking instability mechanism for groove formation. Success of the necking instability mechanism may require rheological weakening or strain localization by faulting, effects not included in our analysis.     

Comparison of Five Numerical Codes for Automated Tracing of Coronal Loops

The three-dimensional (3D) modeling of coronal loops and filaments requires algorithms that automatically trace curvilinear features in solar EUV or soft X-ray images. We compare five existing algorithms that have been developed and customized to trace curvilinear features in solar images: i) the oriented-connectivity method (OCM), which is an extension of the Strous pixel-labeling algorithm (developed by Lee, Newman, and Gary); ii) the dynamic aperture-based loop-segmentation method (developed by Lee, Newman, and Gary); iii) unbiased detection of curvilinear structures (developed by Steger, Raghupathy, and Smith); iv) the oriented-direction method (developed by Aschwanden); and v) ridge detection by automated scaling (developed by Inhester). We test the five existing numerical codes with a TRACE image that shows a bipolar active region and contains over 100 discernable loops. We evaluate the performance of the five codes by comparing the cumulative distribution of loop lengths, the median and maximum loop length, the completeness or detection efficiency, the accuracy, and flux sensitivity. These algorithms are useful for the reconstruction of the 3D geometry of coronal loops from stereoscopic observations with the STEREO spacecraft, or for quantitative comparisons of observed EUV loop geometries with (nonlinear force-free) magnetic field extrapolation models.     

The Maria asteroid family: Genetic relationships and a plausible source of mesosiderites near the 3:1 Kirkwood Gap

We present a mineralogical assessment of 12 Maria family asteroids, using near-infrared spectral data obtained over the years 2000-2009 combined with visible spectral data (when available) to cover the spectral interval of 0.4-2.5 μm. Our analysis indicates the Maria asteroid family, which is located adjacent to the chaotic region of the 3:1 Kirkwood Gap, appears to be a true genetic family composed of assemblages analogous to mesosiderite-type meteorites. Dynamical models by Farinella et al. (Farinella, P., Gunczi, R., Froeschlé, Ch., Froeschlé, C., [1993]. Icarus 101, 174-187) predict this region should supply meteoroids into Earth-crossing orbits. Thus, the Maria family is a plausible source of some or all of the mesosiderites in our meteorite collections. These individual asteroids were most likely once part of a larger parent object that was broken apart and dispersed. One of the Maria dynamical family members investigated, ((695) Bella), was found to be unrelated to the genetic Maria family members. The parameters of (695) Bella indicate an H-chondrite assemblage, and that Bella may be a sister or daughter of Asteroid (6) Hebe.     

Mineralogical characterization of potential targets for the ASTEX mission scenario

In-situ investigation of asteroids is the next logical step in understanding their exact surface mineralogy, petrology, elemental abundances, particle size distribution, internal structure, and collisional evolution. Near-Earth asteroids (NEAs) provide us with ample opportunities for in-situ scientific exploration with lower Δv requirements and subsequently lower costs than their main belt counterparts. The ASTEX mission concept aims at surface characterization of two compositionally diverse NEAs, one with primitive and the other with a strong thermally evolved surface mineralogy. Here we present the first results of our ground-based characterization of potential ASTEX mission targets using the SpeX instrument on the NASA IRTF. Of the four potential targets we characterized, two (1991 JW and 1998 PA) have compositions similar to ordinary chondrite mineralogy. The other two targets (1994 CC and 1999 TA10) are thermally evolved objects with igneous formation histories. While 1994 CC is a triplet system and thus very challenging to orbit the V-type NEA, 1999 TA10 is the most interesting scientific ASTEX target identified so far.     

Three‐dimensional shapes and Fe contents of Stardust impact tracks: A track formation model and estimation of comet Wild 2 coma dust particle densities

Abstract– We investigated three‐dimensional structures of comet Wild 2 coma particle impact tracks using synchrotron radiation (SR) X‐ray microtomography at SPring‐8 to elucidate the nature of comet Wild 2 coma dust particles captured in aerogel by understanding the capture process. All tracks have a similar entrance morphology, indicating a common track formation process near the entrance by impact shock propagation irrespective of impactor materials. Distributions of elements along the tracks were simultaneously measured using SR‐XRF. Iron is distributed throughout the tracks, but it tends to concentrate in the terminal grains and at the bottoms of bulbs. Based on these results, we propose an impact track formation process. We estimate the densities of cometary dust particles based on the hypothesis that the kinetic energy of impacting dust particles is proportional to the track volume. The density of 148 cometary dust particles we investigated ranges from 0.80 to 5.96 g cm^?3 with an average of 1.01 (±0.25) g cm^?3. Moreover, we suggest that less fragile crystalline particles account for approximately 5 vol% (20 wt%) of impacting particles. This value of crystalline particles corresponds to that of chondrules and CAIs, which were transported from the inner region of the solar system to the outer comet‐forming region. Our results also suggest the presence of volatile components, such as organic material and perhaps ice, in some bulbous tracks (type‐C).     

Constraining albedo,diameter and composition of near-Earth asteroids via near-infrared spectroscopy

We present a method to constrain the albedo and diameters of near-Earth asteroids (NEAs) based on thermal flux in their near-infrared spectra (0.7–2.5 μm) using the Standard Thermal Model. Near-infrared spectra obtained with the SpeX instrument on NASA Infrared Telescope Facility are used to estimate the albedo and diameters of 12 NEAs (1992 JE, 1992 UY4, 1999 JD6, 2004 XP14, 2005 YY93, 2007 DS84, 2005 AD13, 2005 WJ56, 1999 JM8, 2005 RC34, 2003 YE45, and 2008 QS11). Albedo estimates were compared with average albedo for various taxonomic classes outlined by Thomas et al. (Thomas, C.A. et al. [2011]. Astron. J. 142(3)) and are consistent with their results. Spectral band parameters, like band centers, are derived and compared to spectra of laboratory mineral mixtures and meteorites to constrain their composition and possible meteorite analogs. Based on our study we estimate the albedos and diameters of these NEAs and compare them with those obtained by other techniques such as ground-based mid-infrared, Spitzer thermal infrared and Arecibo radar. Our results are broadly consistent with the results from other direct methods like radar. Determining the compositions of low albedo asteroids is a challenge due to the lack of deep silicate absorption features. However, based on weak absorption features and albedo, we suggest possible meteorite analogs for these NEAs, which include black chondrites, CM2 carbonaceous chondrites and enstatite achondrites. We did not find any specific trends in albedo and composition among the asteroids we observed.