Structure of the Earth’s circumsolar dust ring

Interplanetary dust particles from comets and asteroids pervade the Solar System and become temporarily trapped into orbital resonances with Earth, leading to a circumsolar dust ring. Using the unique vantage point of the Spitzer Space Telescope from its Earth-trailing solar orbit, we have measured for the first time the azimuthal structure of the Earth’s resonant dust ring. There is a relative paucity of particles within 0.1 AU of the Earth, followed by an enhancement in a cloud that is centered 0.2 AU behind Earth with a width of 0.08 AU along the Earth’s orbit. The North ecliptic pole is ∼3% brighter at 8 μm wavelength when viewed from inside the enhancement. The presence of azimuthal asymmetries in debris disks around other stars is considered strong evidence for planets. By measuring the properties of the Earth’s resonant ring, we can provide “ground truth” to models for interactions of planets and debris disks, possibly leading to improved predictions for detectability of life-bearing planets. The low amplitude of the azimuthal asymmetry in the Earth’s circumsolar ring suggests significant contributions to the zodiacal light from particles that are large (>30 μm) or have large orbital eccentricity that makes capture into mean motion resonances inefficient.     

The effects and characteristics of atmospheric dust during martian global dust storm 2001A

We present retrieved trends in dust optical depth, dust effective radius and surface temperature from our analysis of Mars Global Surveyor Thermal Emission Spectrometer daytime data from global dust storm 2001A, and describe their significance for the martian dust cycle. The dust optical depth becomes correlated with surface pressure during southern spring and summer in years both with and without a global dust storm, indicating that global dust mixing processes are important at those seasons. The correlation is low at other times of the year. We found that the observed decay of optical depths at the later stages of the dust storm match, to first-order, theoretical values of clearing from Stokes–Cunningham fallout of the dust. Zonally averaged effective radius is constant within standard deviation of results (between 1.2 and 2.0 μm, with a global mean for all seasons of 1.7 μm), at all latitudes and seasons except at southern latitudes of 35° and higher around equinoxes in both martian years, where it is larger than average (2–3 μm). The emergence and disappearance of these larger particles correlates with observations of polar cap edge storms at those latitudes. Northern latitude observations under similar conditions did not yield a similar trend of larger average effective radii during the equinoxes. We also report on a linear correlation between daytime surface temperature drop and rise in optical depth during the global dust storm. Global dust storm 2001A produced a significant optical depth and surface temperature change.     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

The opposition and tilt effects of Saturn’s rings from HST observations

The two major factors contributing to the opposition brightening of Saturn’s rings are (i) the intrinsic brightening of particles due to coherent backscattering and/or shadow hiding on their surfaces, and (ii) the reduced interparticle shadowing when the solar phase angle α → 0°. We utilize the extensive set of Hubble Space Telescope observations (Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 199–223) for different elevation angles B and wavelengths λ to disentangle these contributions. We assume that the intrinsic contribution is independent of B, so that any B dependence of the phase curves is due to interparticle shadowing, which must also act similarly for all λ’s. Our study complements that of Poulet et al. (Poulet, F., Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 224), who used a subset of data for a single B ∼ 10°, and the French et al. (French, R.G., Verbiscer, A., Salo, H., McGhee, C.A., Dones, L. [2007b] PASP 119, 623–642) study for the B ∼ 23° data set that included exact opposition. We construct a grid of dynamical/photometric simulation models, with the method of Salo and Karjalainen (Salo and Karjalainen [2003]. Icarus 164, 428–460), and use these simulations to fit the elevation-dependent part of opposition brightening. Eliminating the modeled interparticle component yields the intrinsic contribution to the opposition effect: for the B and A rings it is almost entirely due to coherent backscattering; for the C ring, an intraparticle shadow hiding contribution may also be present.Based on our simulations, the width of the interparticle shadowing effect is roughly proportional to B. This follows from the observation that as B decreases, the scattering is primarily from the rarefied low filling factor upper ring layers, whereas at larger B’s the dense inner parts are visible. Vertical segregation of particle sizes further enhances this effect. The elevation angle dependence of interparticle shadowing also explains most of the B ring tilt effect (the increase of brightness with elevation). From comparison of the magnitude of the tilt effect at different filters, we show that multiple scattering can account for at most a 10% brightness increase as B → 26°, whereas the remaining 20% brightening is due to a variable degree of interparticle shadowing. The negative tilt effect of the middle A ring is well explained by the the same self-gravity wake models that account for the observed A ring azimuthal brightness asymmetry (Salo, H., Karjalainen, R., French, R.G. [2004]. Icarus 170, 70–90; French, R.G., Salo, H., McGhee, C.A., Dones, L. [2007]. Icarus 189, 493–522).     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

Correlations between VIMS and RADAR data over the surface of Titan: Implications for Titan’s surface properties

We apply a multivariate statistical method to Titan data acquired by different instruments onboard the Cassini spacecraft. We have searched through Cassini/VIMS hyperspectral cubes, selecting those data with convenient viewing geometry and that overlap with Cassini/RADAR scatterometry footprints with a comparable spatial resolution. We look for correlations between the infrared and microwave ranges the two instruments cover. Where found, the normalized backscatter cross-section obtained from the scatterometer measurement, corrected for incidence angle, and the calibrated antenna temperature measured along with the scatterometry echoes, are combined with the infrared reflectances, with estimated errors, to produce an aggregate data set, that we process using a multivariate classification method to identify homogeneous taxonomic units in the multivariate space of the samples.In medium resolution data (from 20 to 100 km/pixel), sampling relatively large portions of the satellite’s surface, we find regional geophysical units matching both the major dark and bright features seen in the optical mosaic. Given the VIMS cubes and RADAR scatterometer passes considered in this work, the largest homogeneous type is associated with the dark equatorial basins, showing similar characteristics as each other on the basis of all the considered parameters.On the other hand, the major bright features seen in these data generally do not show the same characteristics as each other. Xanadu, the largest continental feature, is as bright as the other equatorial bright features, while showing the highest backscattering coefficient of the entire satellite. Tsegihi is very bright at 5 μm but it shows a low backscattering coefficient, so it could have a low roughness on a regional scale and/or a different composition. Another well-defined region, located southwest of Xanadu beyond the Tui Regio, seems to be detached from the surrounding terrains, being bright at 2.69, 2.78 and 5 μm but having a low radar brightness. In this way, other units can be found that show correlations or anti-correlations between the scatterometric response and the spectrophotometric behavior, not evident from the optical remote sensing data.     

A meteoroid stream survey using the Canadian Meteor Orbit Radar: II: Identification of minor showers using a 3D wavelet transform

A 7 year survey using the Canadian Meteor Orbit Radar (CMOR), a specular backscattering orbital radar, has produced three million individually measured meteoroid orbits for particles with mean mass near 10⁻⁷ kg. We apply a 3D wavelet transform to our measured velocity vectors, partitioning them into 1° solar longitude bins while stacking all 7 years of data into a single “virtual” year to search for showers which show annual activity and last for at least 3 days. Our automated stream search algorithm has identified 117 meteor showers. We have recovered 42 of the 45 previously described streams from our first reconnaissance survey (Brown, P., Weryk, R.J., Wong, D.K., Jones, J. [2008]. Icarus 195, 317-339). Removing possible duplicate showers from the automated results leaves 109 total streams. These include 42 identified in survey I and at least 62 newly identified streams. Our large data sample and the enhanced sensitivity of the 3D wavelet search compared to our earlier survey have allowed us to extend the period of activity for several major showers. This includes detection of the Geminid shower from early November to late December and the Quadrantids from early November to mid-January. Among our newly identified streams are the Theta Serpentids which appears to be derived from 2008 KP and the Canum Venaticids which have a similar orbit to C/1975 X1 (Sato). We also find evidence that nearly 60% of all our streams are part of seven major stream complexes, linked via secular invariants.     

Mid-infrared spectra of the shocked Murchison CM chondrite: Comparison with astronomical observations of dust in debris disks

We present laboratory mid-infrared transmission/absorption spectra obtained from matrix of the hydrated Murchison CM meteorite experimentally shocked at peak pressures of 10-49 GPa, and compare them to astronomical observations of circumstellar dust in different stages of the formation of planetary systems. The laboratory spectra of the Murchison samples exhibit characteristic changes in the infrared features. A weakly shocked sample (shocked at 10 GPa) shows almost no changes from the unshocked sample dominated by hydrous silicate (serpentine). Moderately shocked samples (21-34 GPa) have typical serpentine features gradually replaced by bands of amorphous material and olivine with increasing shock pressure. A strongly shocked sample (36 GPa) shows major changes due to decomposition of the serpentine and due to devolatilization. A shock melted sample (49 GPa) shows features of olivine recrystallized from melted material.The infrared spectra of the shocked Murchison samples show similarities to astronomical spectra of dust in various young stellar objects and debris disks. The spectra of highly shocked Murchison samples (36 and 49 GPa) are similar to those of dust in the debris disks of HD113766 and HD69830, and the transitional disk of HD100546. The moderately shocked samples (21-34 GPa) exhibit spectra similar to those of dust in the debris disks of Beta Pictoris and BD+20307, and the transitional disk of GM Aur. An average of the spectra of all Murchison samples (0-49 GPa) has a similarity to the spectrum of the older protoplanetary disk of SU Auriga. In the gas-rich transitional and protoplanetary disks, the abundances of amorphous silicates and gases have widely been considered to be a primary property. However, our study suggests that impact processing may play a significant role in generating secondary amorphous silicates and gases in those disks. Infrared spectra of the shocked Murchison samples also show similarities to the dust from comets (C/2002 V1, C/2001 RX14, 9P/Tempel 1, and Hale Bopp), suggesting that the comets also contain shocked Murchison-like material.     

Mercury radar speckle dynamics

Current data reveal that Mercury is a dynamic system with a core which has not yet solidified completely and is at least partially decoupled from the mantle. Radar speckle displacement experiments have demonstrated that the accuracy in spin-dynamics determination for Earth-like planets can approach 10⁻⁵. The extended analysis of space-time correlation properties of radar echoes shows that the behavior of speckles does not prevent estimation of Mercury’s instantaneous spin-vector components to accuracy of a few parts in 10⁷. This limit can be reached with more powerful radar facilities and leads to constraining the interior in more detail from effects of spin dynamics, e.g., from observation of the core-mantle interplay through high precision monitoring of the 88-day spin-variation of Mercury’s crust.