Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Small-scale dust structures in Halley's coma: II. Disintegration of large dust bodies

Small-scale dust structures, SDSs, altogether ∼35 events with extent ∼30-220 km, have been recognized owing to electric field records, mostly near the closest approach of Vega-2 to Halley's nucleus. Several (8-9) morphological forms of SDS have been identified, and all they make one family. Among the family members, the key form (with respect to which, all other forms can be regarded as degenerate) is a sequence of 3-5 dust clouds. The morphological forms represent various Vega-2 passes through SDSs at different stages of development. SDSs observable as the key form consisted of several fairly regularly spaced dust subpopulations, whose plane of symmetry was parallel to the comet orbit plane. That regularity together with specific features of morphological forms strongly constrain disintegration scenarios and dynamics of fragments, and allow to draw a number of conclusions, the main of which are: SDS parent bodies were ice-free dust aggregates lifted from the nucleus near the comet perihelion, whose masses were in the range ∼0.1-1 of the biggest emitted mass (mass of a body accelerated to the escape velocity, i.e., ∼300-1500 kg); the disintegration scenario comprised a few steps, and the first-step disintegration consisted mainly in consecutive detachments of biggest first-step fragments (BF-SFs) from the parent body; a SDS observable as the key form included the dust minitail of parent body and a few BF-SF minitails, the former one being longer than the latter ones; SDS parent bodies had a fractal-like internal structure, and the BF-SF mass was a few percent of the parent body mass; the thermal conductivity of SDS parent body was less than ∼0.4 W m⁻¹ K⁻¹ or so, while the latent heat of gluing organics was roughly 80 kJ mol⁻¹; the disintegration mechanism was a combination of sintering and sublimation of organics. The multistep disintegration of SDS parent bodies can be reconciled with the basically one-step disintegration of aggregates responsible for the dust boundary (Oberc, P., Icarus 1996, 124, 195-208). The fractal-like structure and the relation between BF-SF mass and parent body mass are in agreement with predictions from the Weidenschilling model of comet formation. Large ice-free dust bodies, in particular SDS parent bodies, can be identified with refractory boulders postulated by some comet nucleus models.     

The organic composition of C/2001 A2 (LINEAR): II. Search for heterogeneity within a comet nucleus

The nucleus of Comet C/2001 A2 (LINEAR) split several times during its recent apparition, presenting an unusual opportunity to search for chemical differences in freshly exposed material. We conducted this search using NIRSPEC at the W.M. Keck Observatory on four dates in 2001: 9.5 and 10.5 July and 4.4 and 10.5 August. We detected the R0 and R1 lines of the ν₃ vibrational band of CH₄ near 3.3 μm on all dates. The R2 line was detected on 4.4 and 10.5 August. When we compare production rates of CH₄ to H₂O, we find evidence of a significant enhancement in August relative to that found in July. H₂CO was securely detected via its ν₁ and ν₅ bands on 9.5 July. On 10.5 July, H₂CO emission was much weaker, and its mixing ratio had dropped by a factor of about four. The mixing ratios for other detected volatile species did not change significantly over the course of the observations. We discuss the implications of this evidence for chemical heterogeneity in the nucleus of Comet C/2001 A2.     

The organic composition of Comet C/2001 A2 (LINEAR): I. Evidence for an unusual organic chemistry

We used the NIRSPEC instrument on the Keck-2 telescope atop Mauna Kea, HI to observe Comet C/2001 A2 (LINEAR) in a Target of Opportunity campaign on UT 2001 July 9.5, 10.5 August 4.4, 10.5. We measured seven organic parent volatiles (C₂H₆, C₂H₂, HCN, CH₄, CO, CH₃OH, H₂CO) simultaneously with H₂O. We obtained absolute production rates and relative abundances for parent volatiles, and also measured rotational temperatures for several of these species. The chemical composition of C/2001 A2 differs substantially from any comet we have observed to date. The abundances we measure (relative to H₂O) for C₂H₆, C₂H₂, HCN, and CH₃OH are enriched by a factor of ∼2 to 3 in C/2001 A2 compared with most comets in our database. Other molecular species were detected within the typical range of measured abundances. C/2001 A2 presented a unique opportunity to study the chemistry of a fragmenting comet where pristine areas are exposed to the Sun.     

The contribution of grains to the activity of comets: II. The brightness of the coma

We use the model of grain behavior in the coma developed by Beer et al. [Beer, E.H., Podolak, M., Prialnik, P., 2006. Icarus 180, 473-486] to compute the contribution of ice grains to the brightness of the coma. The motion of an ice grain along the comet-Sun axis is computed, taking into account gas drag, the gravity of the nucleus, and radiation pressure of sunlight. The sublimation of the grains is also included. We assume that the maximum distance that a grain travels along this axis is indicative of the size of the coma, and we compute the resultant brightness as a function of heliocentric distance. The results are then compared to observations.     

Temporal deconvolution of the hydrogen coma I. A hybrid model

This paper describes a semi-analytical model, tentatively called the time-resolved model (TRM), for calculating the proper temporally resolved column density and Lyman-alpha intensity profiles of the extended cometary hydrogen coma, with potentially arbitrary number of photodissociation sources with different initial velocity profiles. The non-analytical part consists of a true 3D particle simulator whose results are used to obtain the necessary corrective parameters for the analytical, spherically symmetric formulas for the column density and intensity profiles. The hybrid approach provides a physically accurate and computationally effective method to recreate the extended coma profiles and supports a new analysis method that can provide hydrogen lifetime estimates and a much enhanced temporal resolution for the water production curve, for which a method is described. The ability to recreate the water production profile by different analysis methods is studied with a combination of real and simulated data. This is the first part of a two part publication with the second part dedicated to real data analysis.     

Long- and Short-Term Photometric Behavior of Comet Hyakutake (1996 B2)

Narrowband filter photometry observations of Comet Hyakutake (1996 B2) were used to investigate this comet's short-term variability as well as its behavior for the apparition as a whole. Utilizing measurements obtained on a total of 13 nights between February 9, 1996, and April 14, 1996, we find that the heliocentric distance (r_H) dependence of the production rates of OH and NH were much shallower than those for either the carbon-bearing species or the visible dust. Based on the OH measurements, the derived water r_H-dependence was also significantly less steep than expected from a basic water vaporization model and required an effective active surface area of about 29 km² at r_H=1.8 AU, 16 km² at r_H=1 AU, and only 13 km² at r_H=0.6 AU. This decrease in active area may be due to seasonally induced variations of a heterogeneous surface, or due to a decreasing contribution of gas from icy grains in the innermost coma. The relative abundances of the minor gas species place Hyakutake into the “typical” category of comets in the A'Hearn et al. (1995, Icarus118, 223-270) taxonomic classification system. The spectrum is generally redder at shorter wavelengths throughout the apparition; however, the dust color progressively changes from being significantly reddened (37%/1000 Å) at large r_H to near-solar at small r_H. This change of color with distance implies a significant change in grain sizes or a changing proportion between two or more grain populations.A major outburst was initiated near March 19.9, just prior to the comet's close approach to Earth. The characteristic recovery from the outburst differed among the observed species, with OH recovering most rapidly, essentially returning to its baseline values by March 25. The spatial radial fall-off of OH throughout this interval was consistent with the expected nominal spatial distribution, while CN and C₂ displayed fall-offs consistent with a distributed source, and the dust fall-off was significantly less steep than 1/ρ, possibly due to fragmenting grains. Rotational lightcurve amplitudes were largest for the OH, NH, and dust, again consistent with the carbon-bearing species primarily originating from a distributed source. Significant variations were observed in the lightcurve amplitude and phase shifts as functions of aperture size. Finally, a refined value for the rotation period of 0.2614±0.0003 day was determined.     

Chemical and physical properties of gas jets in comets: II. Modeling OH, CN and C2 jets in Comet C/1995 O1 (Hale-Bopp) one month after perihelion

We present an analysis of OH, CN, and C₂ jets observed in Comet Hale-Bopp during April 22-26, 1997. We conclude that an extended source, which peaks in productivity after a certain amount of time has passed after being released from the nucleus (8.5, 2.5, and 42.6×¹⁰⁴ s, respectively) is responsible for the observed coma jet morphology in all three species. Sub-micron organic grains are the favored explanation for the extended source. Our models indicate that this extended source produces approximately 40% of the OH, 50% of the C₂, and 75% of the CN. The balance for each species is created by a diffuse nuclear gas source. Compared with the nuclear gas source and normalized to the CN abundance, the composition of the extended source is depleted in OH by a factor of ∼6, and depleted in C₂ by a factor of ∼2. The existence of anti-sunward jets do not require production of radicals throughout the cometary night. Instead, our model demonstrates that active areas exposed to near-twilight conditions throughout the comet's rotational period can produce the observed anti-sunward morphology.     

Navier-Stokes and direct Monte-Carlo simulations of the circumnuclear gas coma: III. Spherical, inhomogeneous sources

We pursue our program of comparative simulations of the cometary gas coma by the two most advanced techniques available: (1) numerical solution of Navier-Stokes equations coupled to the Boltzman equation in the surface boundary layer, and (2) direct Monte-Carlo simulation. Here, we consider two different spherical but compositionally inhomogeneous nuclei, at three very different levels of gas production. The results show the same excellent agreement between the two methods in a domain adjacent to the surface as found precedingly, practically down to free-molecular conditions. A wealth of coma density patterns with non-intuitive structure is obtained. Some of these structures appear even under free-molecular effusion from the surface. The physical origin of all structures is discussed, and their evolution with changing gas production is studied. The computed comae are compared to those computed by various authors precedingly. Intercomparison of the present results demonstrates that differing inhomogeneity patterns may lead to similar structures in the gas coma. Comparison between these structures and those created by homogeneous, aspherical surfaces shows that it is not possible to guess from empirical rules which one of the two processes is responsible for the creation of a given structure. The implications for the interpretation of future high resolution images, or of future in situ mass spectrometric samplings of the near-nucleus gas coma are discussed.     

Chemical and physical properties of gas jets in comets: I. Monte Carlo model of an inner cometary coma

We describe a 3-dimensional, time-dependent Monte Carlo model developed to analyze the chemical and physical nature of a cometary gas coma. Our model includes the necessary physics and chemistry to recreate the conditions applicable to Comet Hale-Bopp when the comet was near 1 AU from the Sun. Two base models were designed and are described here. The first is an isotropic model that emits particles (parents of the observed gases) from the entire nucleus; the second is a jet model that ejects parent particles solely from discrete active areas on the surface of the comet nucleus, resulting in coma jets. The two models are combined to produce the final model, which is compared with observations. The physical processes incorporated in both base models include: (1) isotropic ejection of daughter molecules (the observed gases) in the parent's frame of reference, (2) solar radiation pressure, (3) solar insolation effects, (4) collisions of daughter products with other molecules in the coma, and (5) acceleration of the gas in the coma. The observed daughter molecules are produced when a parent decays, which is represented by either an exponential decay distribution (photodissociation of the parent gas) or a triangular distribution (production from a grain extended source). Application of this model to the analysis the OH, C₂ and CN gas jets observed in the coma of Comet Hale-Bopp is the focus of the accompanying paper [Lederer, S.M., Campins, H., Osip, D.J., 2008. Icarus, in press (this issue)].     

Cassini UVIS observations of the Io plasma torus: I. Initial results

During the Cassini spacecraft's flyby of Jupiter (October, 2000-March, 2001), the Ultraviolet Imaging Spectrograph (UVIS) produced an extensive dataset consisting of 3349 spectrally dispersed images of the Io plasma torus. Here we present an example of the raw data and representative EUV spectra (561-1181 Å) of the torus, obtained on October 1, 2000 and November 14, 2000. For most of the flyby period, the entire Io torus fit within the UVIS field-of-view, enabling the measurement of the total power radiated from the torus in the extreme ultraviolet. A typical value for the total power radiated in the wavelength range of 580-1181 Å is 1.7×10¹² W, with observed variations of up to 25%. Several brightening events were observed. These events lasted for roughly 20 hours, during which time the emitted power increased rapidly by ∼20% before slowly returning to the pre-event level. Observed variations in the relative intensities of torus spectral features provide strong evidence for compositional changes in the torus plasma with time. Spatial profiles of the EUV emission show no evidence for a sharply peaked “ribbon” feature. The ratio of the brightness of the dusk ansa to the brightness of the dawn ansa is observed to be highly variable, with an average value of 1.30. Weak longitudinal variations in the brightness of the torus ansae were observed at the 2% level.     

Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

Analyses of dust coma morphology of Comet Hyakutake (1996 B2) near perigee: outburst behavior, jet motion, source region locations, and nucleus pole orientation

We present inner-coma dust imaging of Comet Hyakutake (1996 B2) obtained on 11 consecutive nights in late March 1996, an interval including a major outburst and the comet’s closest approach to Earth. The evolution of the outburst morphology is followed, along with the motion along the tail of several outburst fragments. Two spiral dust jets—a primary jet, along with a much weaker secondary jet—are visible throughout the interval and are produced by two source regions on a rotating nucleus. These are examined as a function of rotational phase and viewing geometry, with their appearance changing from a nearly face-on view on March 18 to side-on by March 28. The dust outflow velocity as a function of distance from the nucleus is derived, with the dust continuing to accelerate to a distance of 4000 km or more and reaching an average outflow velocity of 0.38 km s⁻¹ between 3000 and 8000 km. We present details of our Monte Carlo modeling of the jets and our methodology of fitting the model to the images. The modeling yields the pole orientation of the nucleus, with an obliquity of approximately 108°, corresponding to an RA of 13^h41^m and a Dec of −1.1°. For an assumed spherical nucleus, the primary active region is centered at approximately −66° latitude, has a radius of about 56°, and therefore covers about 22% of the surface. The source of the secondary jet is at a latitude of −28°, has a radius of about 16°, and is located at a longitude nearly 180° away from the primary source. Estimated uncertainties for the pole orientation and the source locations and sizes are each about 3°. This solution for the nucleus orientation and source locations explains the strong asymmetry in measured production rates before and after perihelion in radio observations (Biver et al., 1999, Astron. J. 118, 1850-1872). The modeling also tightly constrains the sidereal rotation period as 0.2618 ± 0.0001 day, completely consistent with the expected +0.0003 day difference from the observed solar rotation period of 0.2614 ± 0.0004 day determined by Schleicher and Osip (2002, Icarus 159, 210-233), given the pole orientation and position of the comet in its orbit.     

Cassini UVIS observations of the Io plasma torus: IV. Modeling temporal and azimuthal variability

In this fourth paper in a series, we present a model of the remarkable temporal and azimuthal variability of the Io plasma torus observed during the Cassini encounter with Jupiter. Over a period of three months, the Cassini Ultraviolet Imaging Spectrograph (UVIS) observed a dramatic variation in the average torus composition. Superimposed on this long-term variation, is a 10.07-h periodicity caused by an azimuthal variation in plasma composition subcorotating relative to System III longitude. Quite surprisingly, the amplitude of the azimuthal variation appears to be modulated at the beat frequency between the System III period and the observed 10.07-h period. Previously, we have successfully modeled the months-long compositional change by supposing a factor of three increase in the amount of material supplied to Io's extended neutral clouds. Here, we extend our torus chemistry model to include an azimuthal dimension. We postulate the existence of two azimuthal variations in the number of superthermal electrons in the torus: a primary variation that subcorotates with a period of 10.07 h and a secondary variation that remains fixed in System III longitude. Using these two hot electron variations, our model can reproduce the observed temporal and azimuthal variations observed by Cassini UVIS.     

Cassini UVIS observations of the Io plasma torus: III. Observations of temporal and azimuthal variability

In this third paper in a series presenting observations by the Cassini Ultraviolet Imaging Spectrometer (UVIS) of the Io plasma torus, we show remarkable, though subtle, spatio-temporal variations in torus properties. The Io torus is found to exhibit significant, near-sinusoidal variations in ion composition as a function of azimuthal position. The azimuthal variation in composition is such that the mixing ratio of S II is strongly correlated with the mixing ratio of S III and the equatorial electron density and strongly anti-correlated with the mixing ratios of both S IV and O II and the equatorial electron temperature. Surprisingly, the azimuthal variation in ion composition is observed to have a period of 10.07 h—1.5% longer than the System III rotation period of Jupiter, yet 1.3% shorter than the System IV period defined by [Brown, M.E., 1995. J. Geophys. Res. 100, 21683-21696]. Although the amplitude of the azimuthal variation of S III and O II remained in the range of 2-5%, the amplitude of the S II and S IV compositional variation ranged between 5 and 25% during the UVIS observations. Furthermore, the amplitude of the azimuthal variations of S II and S IV appears to be modulated by its location in System III longitude, such that when the region of maximum S II mixing ratio (minimum S IV mixing ratio) is aligned with a System III longitude of ∼200°±15°, the amplitude is a factor of ∼4 greater than when the variation is anti-aligned. This behavior can explain numerous, often apparently contradictory, observations of variations in the properties of the Io plasma torus with the System III and System IV coordinate systems.     

Direct Monte Carlo and multifluid modeling of the circumnuclear dust coma: Spherical grain dynamics revisited

This paper describes the first computations of dust distributions in the vicinity of an active cometary nucleus, using a multidimensional Direct Simulation Monte Carlo Method (DSMC). The physical model is simplistic: spherical grains of a broad range of sizes are liberated by H₂O sublimation from a selection of nonrotating sunlit spherical nuclei, and submitted to the nucleus gravity, the gas drag, and the solar radiation pressure. The results are compared to those obtained by the previously described Dust Multi-Fluid Method (DMF) and demonstrate an excellent agreement in the regions where the DMF is usable. Most importantly, the DSMC allows the discovery of hitherto unsuspected dust coma properties in those cases which cannot be treated by the DMF. This leads to a thorough reconsideration of the properties of the near-nucleus dust dynamics. In particular, the results show that (1) none of the three forces considered here can be neglected a priori, in particular not the radiation pressure; (2) hitherto unsuspected new families of grain trajectories exist, for instance trajectories leading from the nightside surface to the dayside coma; (3) a wealth of balistic-like trajectories leading from one point of the surface to another point exist; on the dayside, such trajectories lead to the formation of “mini-volcanoes.” The present model and results are discussed carefully. It is shown that (1) the neglected forces (inertia associated with a nucleus rotation, solar tidal force) are, in general, not negligible everywhere, and (2) when allowing for these additional forces, a time-dependent model will, in general, have to be used. The future steps of development of the model are outlined.     

Magnetic aggregation: II. Laboratory and microgravity experiments

In a previous publication (Dominik and Nübold, 2002, Icarus 157, 173-186), we presented analytical expressions and theoretical considerations concerning preplanetary dust aggregation with magnetized grains in the solar nebula. The present work is dedicated to the experimental study of magnetic aggregation in a ground-based laboratory as well as under microgravity conditions on parabolic flights. We conducted aggregation experiments with dust analogues in order to study the temporal evolution and the structural outcome of grain growth processes dominated by or comprising exclusively magnetic grains. Within aggregation times ranging from a couple of seconds to a few minutes only, formation of huge chain-like and/or web-like dust aggregates was observed. After aggregate retrieval we were able to study the sizes and structures of the aggregates in great detail. We established the fractal dimension of the aggregates as D_fs=1.20±0.05 for single chains and D_fc=1.50±0.21 for inter-connected web-like structures. This is considerably lower than for non-magnetic grain growth. The dynamic exponent z for the mass increase with time according to t^z was found to be z=2.7 from in-situ video images of the microgravity aggregation runs. The results are compared with the theoretical considerations presented earlier as well as with previous experimental work on the same and on related topics, respectively.     

Orbit determination with very short arcs: II. Identifications

When the observational data are not enough to compute a meaningful orbit for an asteroid/comet we can represent the data with an attributable, i.e., two angles and their time derivatives. The undetermined variables range and range rate span an admissible region of Solar System orbits, which can be sampled by a set of Virtual Asteroids (VAs) selected by means of an optimal triangulation [Milani, A., Gronchi, G.F., de' Michieli Vitturi, M., Kne?evi?, Z., 2004. Celest. Mech. Dyn. Astron. 90, 59-87]. The attributable 4 coordinates are the result of a fit and they have an uncertainty, represented by a covariance matrix. Two short arcs of observations, represented by two attributables, can be linked by considering for each VA (in the admissible region of the first arc) the covariance matrix for the prediction at the time of the second arc, and by comparing it with the attributable of the second arc with its own covariance. By defining an identification penalty we can select the VAs allowing to fit together both arcs and compute a preliminary orbit. Two attributables may not be enough to compute an orbit with convergent differential corrections. Thus the preliminary orbit is used in a constrained differential correction, providing solutions along the Line Of Variation which can be used as second generation VAs to further predict the observations at the time of a third arc. In general the identification with a third arc will ensure a well determined orbit, to which additional sets of observations can be attributed. To test these algorithms we use a large scale simulation and measure the completeness, the reliability and the efficiency of the overall procedure to build up orbits by accumulating identifications. Under the conditions expected for the next generation asteroid surveys, the methods developed in this and in the preceding papers are efficient enough to be used as primary identification methods, with very good results. One important property is that the completeness in finding the possible identifications is as good for comparatively rare orbits, such as the ones of Near-Earth Objects, as for main belt orbits.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: II. Results from a coupled dynamics-microphysics-photochemistry GCM

We present a 2D general circulation model of Titan's atmosphere, coupling axisymmetric dynamics with haze microphysics, a simplified photochemistry and eddy mixing. We develop a parameterization of latitudinal eddy mixing by barotropic waves based on a shallow-water, longitude-latitude model. The parameterization acts locally and in real time both on passive tracers and momentum. The mixing coefficient varies exponentially with a measure of the barotropic instability of the mean zonal flow. The coupled GCM approximately reproduces the Voyager temperature measurements and the latitudinal contrasts in the distributions of HCN and C₂H₂, as well as the main features of the zonal wind retrieved from the 1989 stellar occultation. Wind velocities are consistent with the observed reversal time of the North-South albedo asymmetry of 5 terrestrial years. Model results support the hypothesis of a non-uniform distribution of infrared opacity as the cause of the Voyager temperature asymmetry. Transport by the mean meridional circulation, combined with polar vortex isolation may be at the origin of the latitudinal contrasts of trace species, with eddy mixing remaining restricted to low latitudes most of the Titan year. We interpret the contrasts as a signature of non-axisymmetric motions.