Dynamics of Jupiter’s equatorial region at cloud top level from Cassini and HST images

We present a study of the equatorial region of Jupiter, between latitudes ∼15°S and ∼15°N, based on Cassini ISS images obtained during the Jupiter flyby at the end of 2000, and HST images acquired in May and July 2008. We examine the structure of the zonal wind profile and report the detection of significant longitudinal variations in the intensity of the 6°N eastward jet, up to 60 m s⁻¹ in Cassini and HST observations. These longitudinal variations are, in the HST case, associated with different cloud morphology. Photometric and radiative transfer analysis of the cloud features used as tracers in HST images show that at most there is only a small height difference, no larger than ∼0.5-1 scale heights, between the slow (∼100 m s⁻¹) and fast (∼150 m s⁻¹) moving features. This suggests that speed variability at 6°N is not dominated by vertical wind shears but instead we propose that Rossby wave activity is the responsible for the zonal variability. Removing this variability, we find that Jupiter’s equatorial jet is actually symmetric relative to equator with two peaks of ∼140-150 m s⁻¹ located at latitudes 6°N and 6°S and at a similar pressure level. We also study the local dynamics of particular equatorial features such as several dark projections associated with 5 μm hot spots and a large, long-lived feature called the White Spot (WS) located at 6°S. Convergent flow at the dark projections appears to be a characteristic which depends on the particular morphology and has only been detected in some cases. The internal flow field in the White Spot indicates that it is a weakly rotating quasi-equatorial anticyclone relative to the ambient meridionally sheared flow.     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

Experimental determination of the coefficient of restitution for meter-scale granite spheres

We present results of a series of large-scale experiments to measure the coefficient of restitution for 1-m-diameter rocky bodies in impacts with collision speeds up to ∼1.5 m s⁻¹. The experiments were conducted in an outdoor setting, with two 40-ton cranes used to suspend the ∼1300-kg granite spheres pendulum-style in mutual contact at the bottoms of their respective paths of motion. The spheres were displaced up to ∼1 m from their rest positions and allowed to impact each other in normal-incidence collisions at relative speeds up to ∼1.5 m s⁻¹. Video data from 66 normal-incidence impacts suggest a value for the coefficient of restitution of 0.83 ± 0.06 for collisions between ∼1-m-scale spheres at speeds of order 1 m s⁻¹. No clear trend of coefficient of restitution with impact speed is discernable in the data.     

Deep Impact, Stardust-NExT and the behavior of Comet 9P/Tempel 1 from 1997 to 2010

We present observational data for Comet 9P/Tempel 1 taken from 1997 through 2010 in an international collaboration in support of the Deep Impact and Stardust-NExT missions. The data were obtained to characterize the nucleus prior to the Deep Impact 2005 encounter, and to enable us to understand the rotation state in order to make a time of arrival adjustment in February 2010 that would allow us to image at least 25% of the nucleus seen by the Deep Impact spacecraft to better than 80 m/pixel, and to image the crater made during the encounter, if possible. In total, ∼500 whole or partial nights were allocated to this project at 14 observatories worldwide, utilizing 25 telescopes. Seventy percent of these nights yielded useful data. The data were used to determine the linear phase coefficient for the comet in the R-band to be 0.045 ± 0.001 mag deg⁻¹ from 1° to 16°. Cometary activity was observed to begin inbound near r ∼ 4.0 AU and the activity ended near r ∼ 4.6 AU as seen from the heliocentric secular light curves, water-sublimation models and from dust dynamical modeling. The light curve exhibits a significant pre- and post-perihelion brightness and activity asymmetry. There was a secular decrease in activity between the 2000 and 2005 perihelion passages of ∼20%. The post-perihelion light curve cannot be easily explained by a simple decrease in solar insolation or observing geometry. CN emission was detected in the comet at 2.43 AU pre-perihelion, and by r = 2.24 AU emission from C₂ and C₃ were evident. In December 2004 the production rate of CN increased from 1.8 × 10²³ mol s⁻¹ to Q_CN = 2.75 × 10²³ mol s⁻¹ in early January 2005 and 9.3 × 10²⁴ mol s⁻¹ on June 6, 2005 at r = 1.53 AU.     

Maintenance of equatorial superrotation in the atmospheres of Venus and Titan

This paper extends Leovy's theory on Venus’ equatorial superrotation by analytically examining additional terms in the mean zonal momentum equation that stably balances the momentum source of pumping by thermal tides. The general analytical solution is applied to the atmospheres of both Venus and Saturn's moon Titan. The main results are: (i) Venus’ equatorial superrotation of 118 m s⁻¹ results primarily from a balance between the momentum source of pumping by thermal tides and the momentum sink of meridional advection of wind shear by horizontal branches of the Hadley circulation; (ii) no solution is found for Titan's stratospheric equatorial superrotation centered at the 1-hPa level; (iii) however, if the main solar radiation absorption layer in Titan's stratosphere is lifted from 1 hPa (∼185 km) to 0.1 hPa (∼288 km), an equatorial superrotation of ∼110 m s⁻¹ centered at 0.1-hPa could be maintained. Titan's equatorial superrotation results mainly from a balance between the momentum source of tidal pumping and the momentum sink of frictional drag.     

Generation and emplacement of fine-grained ejecta in planetary impacts

We report here on a survey of distal fine-grained ejecta deposits on the Moon, Mars, and Venus. On all three planets, fine-grained ejecta form circular haloes that extend beyond the continuous ejecta and other types of distal deposits such as run-out lobes or ramparts. Using Earth-based radar images, we find that lunar fine-grained ejecta haloes represent meters-thick deposits with abrupt margins, and are depleted in rocks ?1 cm in diameter. Martian haloes show low nighttime thermal IR temperatures and thermal inertia, indicating the presence of fine particles estimated to range from ∼10 μm to 10 mm. Using the large sample sizes afforded by global datasets for Venus and Mars, and a complete nearside radar map for the Moon, we establish statistically robust scaling relationships between crater radius R and fine-grained ejecta run-out r^* for all three planets. On the Moon, r^* ∼ R^−0.18 for craters 5-640 km in diameter. For Venus, radar-dark haloes are larger than those on the Moon, but scale as r^* ∼ R^−0.49, consistent with ejecta entrainment in Venus’ dense atmosphere. On Mars, fine-ejecta haloes are larger than lunar haloes for a given crater size, indicating entrainment of ejecta by the atmosphere or vaporized subsurface volatiles, but scale as R^−0.13, similar to the ballistic lunar scaling. Ejecta suspension in vortices generated by passage of the ejecta curtain is predicted to result in ejecta run-out that scales with crater size as R^1/2, and the wind speeds so generated may be insufficient to transport particles at the larger end of the calculated range. The observed scaling and morphology of the low-temperature haloes leads us rather to favor winds generated by early-stage vapor plume expansion as the emplacement mechanism for low-temperature halo materials.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

Laboratory simulations using the Arizona State University Vortex Generator (ASUVG) were run to simulate sediment flux in dust devils in terrestrial ambient and Mars-analog conditions. The objective of this study was to measure vortex sediment flux in the laboratory to yield estimations of natural dust devils on Earth and Mars, where all parameters may not be measured. These tests used particles ranging from 2 to 2000 μm in diameter and 1300 to 4800 kg m⁻³ in density, and the results were compared with data from natural dust devils on Earth and Mars. Typically, the cores of dust devils (regardless of planetary environment) have a pressure decrease of ∼0.1-1.5% of ambient atmospheric pressure, which enhances the lifting of particles from the surface. Core pressure decreases in our experiments ranged from ∼0.01% to 5.00% of ambient pressure (10 mbar Mars cases and 1000 mbar for Earth cases) corresponding to a few tenths of a millibar for Mars cases and a few millibars for Earth cases. Sediment flux experiments were run at vortex tangential wind velocities of 1-45 m s⁻¹, which typically correspond to ∼30-70% above vortex threshold values for the test particle sizes and densities. Sediment flux was determined by time-averaged measurements of mass loss for a given vortex size. Sediment fluxes of ∼10⁻⁶-10⁰ kg m⁻² s⁻¹ were obtained, similar to estimates and measurements for fluxes in dust devils on Earth and Mars. Sediment flux is closely related to the vortex intensity, which depends on the strength of the pressure decrease in the core (ΔP). This study found vortex size is less important for lifting materials because many different diameters can have the same ΔP. This finding is critical in scaling the laboratory results to natural dust devils that can be several orders of magnitude larger than the laboratory counterparts.     

New laboratory measurements of CH4 in Titan's conditions and a reanalysis of the DISR near-surface spectra at the Huygens landing site

Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm⁻¹ spectral resolution, adequate to resolve the lines at high pressure (p_N2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν₂+2ν₃ band (7515-7620 cm⁻¹), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm⁻¹ sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH₄ mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.     

Interplanetary survival probability of Aspergillus terreus spores under simulated solar vacuum ultraviolet irradiation

This work is a part of ESA/EU SURE project aiming to quantify the survival probability of fungal spores in space under solar irradiation in the vacuum ultraviolet (VUV) (110-180 nm) spectral region. The contribution and impact of VUV photons, vacuum, low temperature and their synergies on the survival probability of Aspergillus terreus spores is measured at simulated space conditions on Earth. To simulate the solar VUV irradiation, the spores are irradiated with a continuous discharge VUV hydrogen photon source and a molecular fluorine laser, at low and high photon intensities at 10¹⁵ photon m⁻² s⁻¹ and 3.9×10²⁷ photons pulse⁻¹ m⁻² s⁻¹, respectively. The survival probability of spores is independent from the intensity and the fluence of photons, within certain limits, in agreement with previous studies. The spores are shielded from a thin carbon layer, which is formed quickly on the external surface of the proteinaceous membrane at higher photon intensities at the start of the VUV irradiation. Extrapolating the results in space conditions, for an interplanetary direct transfer orbit from Mars to Earth, the spores will be irradiated with 3.3×10²¹ solar VUV photons m⁻². This photon fluence is equivalent to the irradiation of spores on Earth with 54 laser pulses with an experimental ∼92% survival probability, disregarding the contribution of space vacuum and low temperature, or to continuous solar VUV irradiation for 38 days in space near the Earth with an extrapolated ∼61% survival probability. The experimental results indicate that the damage of spores is mainly from the dehydration stress in vacuum. The high survival probability after 4 days in vacuum (∼34%) is due to the exudation of proteins on the external membrane, thus preventing further dehydration of spores. In addition, the survival probability is increasing to ∼54% at 10 K with 0.12 K/s cooling and heating rates.     

The Unusual Absorption Line Spectrum of Quasar SDSS J2220+0109

We report the Balmer broad absorption lines (BALs) in the quasar SDSS J2220 + 0109 discovered from the SDSS data, and present a detailed analysis of the peculiar absorption line spectrum, including the He I* multiplet at λλ3189, 3889 arising from the metastable 2³s-state helium and the Balmer Hα and Hβ lines from the excited hydrogen H I of n = 2 level, which are rarely seen in quasar spectra, as well as many absorption lines arising from the excited Fe II* of the levels 7 955 cm⁻¹, 13 474 cm⁻¹ and 13 673 cm⁻¹ in the wavelength range 3100∼3300 Å. Ca II H, K absorption line doublets also clearly appear in the SDSS spectrum. All absorption lines show a similar blueshifted velocity structure of Δv ≈ − 1500 ∼ 0 km·s⁻¹ relative to the quasar's systematic redshift determined from the emission lines. Detailed analysis suggests that the Balmer absorption lines should arise from the partially ionized region with a column density of N_HI ≈ 10²¹ cm⁻² for an electron density of ne ∼ 10⁶ cm⁻³; and that the hydrogen n = 2 level may be populated via collisional excitation with Lyα pumping.     

The polar contribution to the heat flow of Io

Polar brightness temperatures on Io are higher than expected for any passive surface heated by absorbed sunlight. This discrepancy implies large scale volcanic activity from which we derive a new component of Io's heat flow. We present a ‘Three Component’ thermal background, infrared emission model for Io that includes active polar regions. The widespread polar activity contributes an additional ∼0.6 W m⁻² to Io's heat flow budget above the ∼2.5 W m⁻² from thermal anomalies. Our estimate for Io's global average heat flow increases to ∼3±1 W m⁻² and ∼1.3±0.4×10¹⁴ watts total.     

The evolving flow of Jupiter’s White Ovals and adjacent cyclones

We present results regarding the dynamical meteorology of Jupiter’s White Ovals at different points in their evolution. Starting from the era with three White Ovals FA, BC, and DE (Galileo), continuing to the post-merger epoch with only one Oval BA (Cassini), and finally to Oval BA’s current reddened state (New Horizons), we demonstrate that the dynamics of their flow have similarly evolved along with their appearance. In the Galileo epoch, Oval DE had an elliptical shape with peak zonal wind speeds of ∼90 m s⁻¹ in both its northern and southern peripheries. During the post-merger epoch, Oval BA’s shape was more triangular and less elliptical than Oval DE; in addition to widening in the north-south direction, its northern periphery was 20 m s⁻¹ slower, and its southern periphery was 20 m s⁻¹ faster than Oval DE’s flow during the Galileo era. Finally, in the New Horizons era, the reddened Oval BA had evolved back to a classical elliptical form. The northern periphery of Oval BA increased in speed by 20 m s⁻¹ from Cassini to New Horizons, ending up at a speed nearly identical to that of the northern periphery of Oval DE during Galileo. However, the peak speeds along the southern rim of the newly formed Oval BA were consistently faster than the corresponding speeds in Oval DE, and they increased still further between Cassini and New Horizons, ending up at ∼140-150 m s⁻¹. Relative vorticity maps of Oval BA reveal a cyclonic ring surrounding its outer periphery, similar to the ring present around the Great Red Spot. The cyclonic ring around Oval BA in 2007 appears to be moderately stronger than observed in 1997 and 2001, suggesting that this may be associated with the coloration of the vortex. The modest strengthening of the winds in Oval BA, the appearance of red aerosols, and the appearance of a turbulent, cyclonic feature to Oval BA’s northwest create a strong resemblance with the Great Red Spot from both a dynamical and morphological perspective.In addition to the White Ovals, we also measure the winds within two compact cyclonic regions, one in the Galileo data set and one in the Cassini data set. In the images, these cyclonic features appear turbulent and filamentary, but our wind field reveals that the flow manifests as a coherent high-speed collar surrounding relatively quiescent interiors. Our relative vorticity maps show that the vorticity likewise concentrates in a collar near the outermost periphery, unlike the White Ovals which have peak relative vorticity magnitudes near the center of the vortex. The cyclones contain several localized bright regions consistent with the characteristics of thunderstorms identified in other studies. Although less studied than their anticyclonic cousins, these cyclones may offer crucial insights into the planet’s cloud-level energetics and dynamical meteorology.     

Venus Doppler winds at cloud tops observed with ESPaDOnS at CFHT

We present new wind measurements in Venus’ lower mesosphere from visible spectroscopy during the 2007 worldwide coordinated ground campaign in support of ESA's Venus Express mission. These observations consisted of high-resolution spectra of Fraunhofer lines in the entire visible range (0.37-1.05 μm) to measure the winds near 68 km using the Doppler shift of solar radiation scattered by clouds toward the observer's direction. The observations included various points of the dayside hemisphere at a phase angle of ∼109°. We took advantage of two symmetrical elongations in July and September 2007 at Canada-France-Hawaii's 3.6-m telescope. Kinematical fits to the Doppler winds provide a mean equatorial velocity of (104±10) m s⁻¹ for the zonal retrograde flow. This velocity agrees quite well with the mean value obtained by tracking the UV markings from several spacecraft.     

Vertical structure of Jupiter's troposphere from nonlinear simulations of long-lived vortices

Numerical simulations of jovian vortices at tropical and temperate latitudes, under different atmospheric conditions, have been performed using the EPIC code [Dowling, T.E., Fisher, A.S., Gierasch, P.J., Harrington, J., LeBeau, R.P., Santori, C.M., 1998. Icarus 132, 221-238] to simulate the high-resolution observations of motions and of the lifetimes presented in a previous work [Legarreta, J., Sánchez-Lavega, A., 2005. Icarus 174, 178-191] and infer the vertical structure of Jupiter's troposphere. We first find that in order to reproduce the longevity and drift rate of the vortices, the Brunt-Väisälä frequency of the atmosphere in the upper troposphere (pressures P∼1 to 7 bar) should have a lower limit value of 5×¹⁰⁻³ s⁻¹, increasing upward up to 1.25×¹⁰⁻² s⁻¹ at pressures P∼0.5 bar (latitudes between 15° and 45° in both hemispheres). Second, the vortices drift also depend on the vertical structure of the zonal wind speed in the same range of altitudes. Simulations of the slowly drifting Southern hemisphere vortices (GRS, White Ovals and anticyclones at 40° S) require a vertically-constant zonal-wind with depth, but Northern hemisphere vortices (cyclonic “barges” and anticyclones at 19, 41 and 45° N) require decreasing winds at a rate of ∼5 m s⁻¹ per scale height. However vortices drifting at a high speed, close to or in the peak of East or West jets and in both hemispheres, require the wind speed slightly increasing with depth, as is the case for the anticyclones at 20° S and at 34° N. We deduce that the maximum absolute vertical shear of the zonal wind from P∼1 bar up to P∼7 bar in these jets is ∼15 m s⁻¹ per scale height. Intense vortices with tangential velocity at their periphery ∼100 m s⁻¹ tend to decay asymptotically to velocities ∼40 to 60 m s⁻¹ with a characteristic time that depends on the vortex intensity and static stability of the atmosphere. The vortices adjust their tangential velocity to the averaged peak to peak velocity of the opposed eastward and westward jets at their boundary. We show through our simulations that large-scale and long-lived vortices whose maximum tangential velocity is ∼100 m s⁻¹ can survive by absorbing smaller intense vortices.     

Martian post-impact hydrothermal systems incorporating freezing

We simulate the evolution of post-impact hydrothermal systems within 45 km and 90 km diameter craters on Mars. We focus on the effects of freezing, which alters the permeability structure and fluid flow compared with unfrozen cases. Discharge rates, total discharge and water-rock ratios increase with permeability. Systems with permeabilities of 10⁻¹⁰ m² or higher exhibit convection in the hydrosphere, allowing them to derive heat from greater depths. Surface discharges persist for ∼10³-10⁵ years under freezing surface conditions, with higher permeabilities permitting longer lifetimes. Maximum discharge rates and total discharges range from 0.1 to 10 m³ s⁻¹ and 10⁹ to 10¹² m³, respectively, for systems with permeabilities between 10⁻¹⁴ and 10⁻¹² m². Near-surface water-rock ratios range from <1 for low permeability, frozen cases to ∼10³ for high permeabilities and/or unfrozen cases. Propagation of the freezing front radially inwards focuses flow towards the center of the crater resulting in a diagnostic increase in water-rock ratios there. This process may explain the phyllosilicate assemblages observed at some crater central peaks.     

The thermal conductivity of meteorites: New measurements and analysis

We have measured the thermal conductivity at low temperatures (5-300 K) of six meteorites representing a range of compositions, including the ordinary chondrites Cronstad (H5) and Lumpkin (L6), the enstatite chondrite Abee (E4), the carbonaceous chondrites NWA 5515 (CK4 find) and Cold Bokkeveld (CM2), and the iron meteorite Campo del Cielo (IAB find). All measurements were made using a Quantum Design Physical Properties Measurement System, Thermal Transport Option (TTO) on samples cut into regular parallelepipeds of ∼2-6 mm dimension. The iron meteorite conductivity increases roughly linearly from 15 W m⁻¹ K⁻¹ at 100 K to 27 W m⁻¹ K⁻¹ at 300 K, comparable to typical values for metallic iron. By contrast, the conductivities of all the stony samples except Abee appear to be controlled by the inhomogeneous nature of the meteorite fabric, resulting in values that are much lower than those of pure minerals and which vary only slightly with temperature above 100 K. The L and CK sample conductivities above 100 K are both about 1.5 W m⁻¹ K⁻¹, that of the H is 1.9 W m⁻¹ K⁻¹, and that of the CM sample is 0.5 W m⁻¹ K⁻¹; by contrast the literature value at 300 K for serpentine is 2.5 W m⁻¹ K⁻¹ and those of enstatite and olivine range from 4.5 to 5 W m⁻¹ K⁻¹ (which is comparable to the Abee value). These measurements are among the first direct measurements of thermal conductivity for meteorites. The results compare well with previous estimates for meteorites, where conductivity was derived from diffusivity measurements and modeled heat capacities; our new values are of a higher precision and cover a wider range of temperatures and meteorite types. If the rocky material that makes up asteroids and provides the dust to comets, Kuiper Belt objects, and icy satellites has the same low thermal conductivities as the ordinary and carbonaceous chondrites measured here, this would significantly change models of their thermal evolution. These values would also lower their thermal inertia, thus affecting the Yarkovsky and YORP evolution of orbits and spin for solid objects; however, in this case the effect would not be as great, as thermal inertia only varies as the square root of the conductivity and, for most asteroids, is controlled by the dusty nature of asteroidal surfaces rather than the conductivity of the material itself.     

Global contraction of planetary bodies due to despinning: Application to Mercury and Iapetus

We extend previous work on the global tectonic patterns generated by despinning with a self-consistent treatment of the isotropic despinning contraction that has been ignored. We provide simple analytic approximations that quantify the effect of the isotropic despinning contraction on the global shape and tectonic pattern. The isotropic despinning contraction of Mercury is ∼93 m (T/1 day)⁻², where T is the initial rotation period. If we take into account both the isotropic contraction and the degree-2 deformations associated with despinning, the preponderance of compressional tectonic features on Mercury’s surface requires an additional isotropic contraction ?1 km (T/1 day)⁻², presumably due to cooling of the interior and growth of the solid inner core. The isotropic despinning contraction of Iapetus is ∼9 m (T/16 h)⁻², and it is not sensitive to the presence of a core or the thickness of the elastic lithosphere. The tectonic pattern expected for despinning, including the isotropic contraction, does not explain Iapetus’ ridge. Furthermore, the ridge remains unexplained with the addition of any isotropic compressional stresses, including those generating by cooling.