Zonal mean properties of Jupiter's upper troposphere from voyager infrared observations

The highest spatial resolution Voyager IRIS spectra are used to produce zonal averages of the temperature at the 150- and 270-mb pressure levels, of the para-hydrogen fraction at 270 mb, of the ammonia abundance near the 680-mb level, and of two infrared cloud optical depths, one near 5 μm and one near 45 μm wavelength. There are two cloud components, one uniformly distributed and only apparent at 5 μm, and another that correlates strongly with the ammonia abundance and that is apparent at both 5 and 45 μm. From the ratio of optical depths at the two wavelengths, the particles in the variable cloud are between 3 and 10 μm in radius. This cloud is located near the ammonia condensation level. The other particles are either smaller or deeper. The cloud and ammonia distribution is consistent with concentration by upward vertical motion at the equatorward edges of prograde atmospheric jets. The temperature field is also consistent with such vertical motion, with radiative heating balancing adiabatic expansional cooling. The para-hydrogen distribution also appears consistent, but noise levels are high. The thermal wind shear indicates decay of the jets with height within the upper troposphere, with a vertical scale of two or three scale heights. The entire set of upper troposphere data is consistent with a simple axisymmetric dynamical model with Coriolis acceleration of the zonal wind balanced by a linear drag. The meridional residual mean circulation in the model, if interpreted also as a Lagrangian mean circulation, would explain nicely the distribution of ammonia and para-hydrogen. The circulation is a response to a deeper tropospheric flow of unknown origin. However, the horizontal scale of jets is on the order of the deformation radius based on a scale height at the base of the upper troposphere. It is conjectured that the physics of the flow may require this to be true, and may also require that the relative vorticity gradient be of the same order as the planetary vorticity gradient, thereby fixing both the dimensions and amplitudes of the jets.     

A co-accretional model of satellite formation

Numerical calculation of a simple accretion model including the effects of tidal friction indicate that coformation is tenable only if the planet's Q is less than about 10³. The parameter which most strongly affects the final mass ratio of the pair is the time at which the secondary embryo is introduced. Our model yields the proper Moon-Earth mass ratio if the Moon embryo is introduced when the Earth is only about $\text{1}\text{10}$ of its final mass. The lunar orbit remains at about 10 Earth radii throughout most of the growth.This model of satellite formation overcomes two difficulties of the “circumterrestrial cloud” model of Ruskol (1960, 1963, 1972): (1) The difficulty of accumulating a mass as great as the entire Moon before gravitational instability reduces the cloud to a small number of moonlets is removed. (2) The differences between terrestrial and outer planet satellite systems is easily understood in terms of the differences in Q between these planets. The high Q of the outer planets does not allow a satellite embryo to survive a significant portion of the accretion process, thus only small bodies which formed very late in the accumulation of the planet remain as satellites. The low Q of the terrestrial planets allows satellite embryos of these planets to survive during accretion, thus massive satellites such as the Earth's Moon are expected. The present lack of such satellites of the other terrestrial planets may be the result of tidal evolution, either infall following primary despinning (Burns, 1973) or escape due to increase in orbit eccentricity.     

The effect of ammonia ice on the outgoing thermal radiance from the atmosphere of Jupiter

We examine the effects of NH₃ ice particle clouds in the atmosphere of Jupiter on outgoing thermal radiances. The cloud models are characterized by a number density at the cloud base, by the ratio of the scale height of the vertical distribution of particles (H_p) to the gas scale height (H_g), and by an effective particle radius. NH₃ ice particle-scattering properties are scaled from laboratory measurements. The number density for the various particle radius and scale height models is inferred from the observed disk average radiance at 246 cm^?1, and preliminary lower limits on particle sizes are inferred from the lack of apparent NH₃ absorption features in the observed spectral radiances as well as the observed minimum flux near 2100 cm^?1. We find lower limits on the particle size of 3 μm if H_p/H_g = 0.15, or 10μmif H_p/H_g = 0.50 or 0.05. NH₃ ice particles are relatively dark near the far-infrared and 8.5-μm atmospheric windows, and the outgoing thermal radiances are not very sensitive to various assumptions about the particle-scattering function as opposed to radiances at 5 μm, where particles are relatively brighter. We examined observations in these three different spectral window regions which provide, in principle, complementary constraints on cloud parameters. Characterization of the cloud scale height is difficult, but a promising approach is the examination of radiances and their center-to-limb variation in spectral regions where there is significant opacity provided by gases of known vertical distribution. A blackbody cloud top model can reduce systematic errors due to clouds in temperature sounding to the level of 1K or less. The NH₃ clouds provide a substantial influence on the internal infrared flux field near the 600-mbar level.     

Response of an optically thin,isothermal atmosphere to a convective overshoot

Radiation is believed to be hostile to the generation of gravity waves by granulation at the base of photosphere where the radiation is effective. A convective overshoot from subphotosphere seems able to penetrate to a height where the solar temperature is minimum and to excite the gravity waves in a stable region there.The response of the solar atmosphere to a Gaussian disturbance characterizing such a convective overshoot is studied in an unbounded isothermal atmosphere. Radiative effects are included, but only in regions which are optically thin. The response is measured in terms of mean vertical kinetic energy density (E _z) and mean vertical external energy flux (Q _z). E _z and Q _z were calculated for a wide range of frequencies centered at the observed 5-min velocity oscillation period. The computed sharp and broad power spectra at the lower chromosphere and the upper photosphere, respectively, are attributed to the combined effects of space damping and source function. Low-frequency waves (2000 s or longer) are found to be not responsible for depositing energy in the upper solar atmosphere.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Lunar viscosity as obtained from the selenotherms

Based on the selenothermsT(z) (= temperature-depth functions) and melting point-depth functionsT _m(z) viscosity valuesη(z) are calculated. According to two different creep laws used, two sets of viscosity values are obtained. Viscosities in the outer part of the Moon are found to be larger than those anywhere on Earth. These high values ofη explain the large elasticityQ found in lunar seismograms. Viscosities below about 500 km in depth are so small that, at present, some kind of convection or a flow of matter is possible. Tidegenerated moonquakes at depths of around 1000 km seem to be connected with some viscous process. From considerations of viscosities at the time period of mare filling, some selection of ancient selenotherms may be performed.     

Existence and amount of intergalactic dust

The existence of intergalactic dust has been proved by the following observational facts: the decrease of the numbers of distant galaxies and clusters of galaxies behind the central regions of near clusters of galaxies; the different distributions of RR Lyrae stars and galaxies near ι Microscopii (Hoffmeister's cloud); the dependence of colour excesses of galaxies on supergalactic coordinates as well as on the surface density of bright galaxies; the colour index vs redshift correlation of quasistellar objects. The densities of intergalactic dust are estimated to be between 5×10^?30 g cm^?3 (near the centers of clusters of galaxies) and 2×10^?34 g cm^?3 (in general intergalactic space). The grains may be formed either in the early phases of the Universe (25相似文献   

Accretion of dust by chondrules in a MHD-turbulent solar nebula

Numerical magnetohydrodynamic (MHD) simulations of a turbulent solar nebula are used to study the growth of dust mantles swept up by chondrules. A small neighborhood of the solar nebula is represented by an orbiting patch of gas at a radius of 3 AU, and includes vertical stratification of the gas density. The differential rotation of the nebular gas is replaced by a shear flow. Turbulence is driven by destabilization of the flow as a result of the magnetorotational instability (MRI), whereby magnetic field lines anchored to the gas are continuously stretched by the shearing motion. A passive contaminant mimics small dust grains that are aerodynamically well coupled to the gas, and chondrules are modeled by Lagrangian particles that interact with the gas through drag. Whenever a chondrule enters a region permeated by dust, its radius grows at a rate that depends on the local dust density and the relative velocity between itself and the dust. The local dust abundance decreases accordingly. Compaction and fragmentation of dust aggregates are not included. Different chondrule volume densities ρ_c lead to varying depletion and rimmed-chondrule size growth times. Most of the dust sweep-up occurs within ～1 gas scale-height of the nebula midplane. Chondrules can reach their asymptotic radius in 10–800 years, although short growth times due to very high ρ_c may not be altogether realistic. If the sticking efficiency Q of dust to chondrules depends on their relative speed δv, such that Q < 10^?2 whenever δv > v_stick ≈ 34 cm/s (with v_stick a critical sticking velocity), then longer growth times result due to the prevalence of high MRI-turbulent relative velocities. The vertical variation of nebula turbulent intensity results in a moderate dependence of mean rimmed-chondrule size with nebula height, and in a ～20% dispersion in radius values at every height bin. The technique used here could be combined with Monte Carlo (MC) methods that include the physics of dust compaction, in a self-consistent MHD-MC model of dust rim growth around chondrules in the solar nebula.     

Stellar occultations by turbulent planetary atmospheres: A heuristic scattering model

A model of ray refraction by an isothermal atmosphere with a scattering screen at the center of bending is used to generate analytic results which simulate the effects of real atmospheric turbulence on occultations. Calculations are carried through for scattering which is constant with height and for exponential height dependence. The effect of the scattering is to bias the the mean intensity of the occulted source, and hence systematically to distort bending angles and height differences obtained from inversion of the intensity data. However, the effect is of order 〈δ?²〉/?² for either model, where 〈δ?²〉 is the mean square scattering angle and ? is the average bending angle. The effect turns out to be small for plausible turbulence, since 〈δ?²〉/gfe² is of approximately the same order as the relative mean square density fluctuation. Thus the random effects of turbulence are unlikely to be a source of large systematic error in occultations, provided that the data can be meaningfully averaged either temporally or over a number of occultation events.     

The analytic properties of the flapping current sheets in the earth magnetotail

The current sheet in Earth’s magnetotail often flaps, and the flapping waves could be induced propagating towards the dawn and dusk flanks, which could make the current sheet dynamic. To explore the dynamic characteristics of current sheet associated with the flapping motion holistically and provide reasonable physical interpretations, detailed direct calculation and analysis have been applied to one approximate analytic model of magnetic field in the flapping current sheet. The main results from the model demonstrate: (1) the magnetic fluctuation amplitude is attenuated from the center of current sheet to the lobe regions; The larger wave amplitude would induce the larger magnetic amplitude; (2) the curvature of magnetic field lines (MFLs), with maximum at the center of current sheet, is only dependent on the displacement Z along the south-north direction from the center of current sheet, regardless of the tilt of current sheet; (3) the half-thickness of neutral sheet, h, the minimum curvature radius of MFLs, R_cmin, and the tilt angle of current sheet, δ, satisfies h=R_cmin cos δ; (4) the gradient of magnetic strength forms a double-peak profile, and the peak value would be more intense if the local current sheet is more tilted; (5) current density j and its j_y, j_z components reach the extremum at the center of CS. j and j_z would be more intense if the local current sheet is more tilted, but it is not the case for j_y; and (6) the field-aligned component of current density mainly appears in the neutral sheet, and the sign of it would change alternatively as the flapping waves passing by. To check the validity of the model, one simulation on the virtual measurements has been made, and the results are in well consistence with actual observations of Cluster.     

Delta spots and great flares

Using eighteen years of observations at Big Bear, we summarize the development of δ spots and the great flares they produce. We find δ groups to develop in three ways: eruption of a single complex active region formed below the surface, eruption of large satellite spots near (particularly in front of) a large older spot, or collision of spots of opposite polarity from different dipoles. Our sample of twenty-one δ spots shows that once they lock together, they never separate, although rarely an umbra is ejected. The δ spots are already disposed to their final form when they emerge. The driving force for the shear is spot motion, either flux emergence or the forward motion of p spots in an inverted magnetic configuration. We observe the following phenomena preceding great flares:

1. δ spots, preferentially Types 1 and 2.
2. Umbrae obscured by Hα emission.
3. Bright Hα emission marking flux emergence and reconnection.
4. Greatly sheared magnetic configurations, marked by penumbral and Hα fibrils parallel to the inversion line.

We assert that with adequate spatial resolution one may predict the occurrence of great flares with these indicators.     

Mechanisms of the vertical secular heating of a stellar disk

We investigate the nonlinear growth stages of the bending instability in stellar disks with exponential radial density profiles. We found that the unstable modes are global (the wavelengths are larger than the disk scale lengths) and that the instability saturation level is much higher than that following from a linear criterion. The instability saturation time scales are of the order of one billion years or more. For this reason, the bending instability can play an important role in the secular heating of a stellar disk in the z direction. In an extensive series of numerical N-body simulations with a high spatial resolution, we were able to scan in detail the space of key parameters (the initial disk thickness z₀, the Toomre parameter Q, and the ratio of dark halo mass to disk mass M_h/M_d). We revealed three distinct mechanisms of disk heating in the z direction: bending instability of the entire disk, bending instability of the bar, and heating on vertical inhomogeneities in the distribution of stellar matter.     

Local time asymmetry in the characteristics of Pc5 magnetic pulsations

The wave characteristics of Pc5 magnetic pulsations are analyzed with data of OGO-5, ISEE-1 and -2 satellites. The toroidal modes (δB_D >δB_H) of Pc5 pulsations are observed at a higher magnetic latitude in the dawnside outer magnetosphere. The compressional and poloidal modes (δB_z.dfnc; ～ δB_H >δB_D) of Pc5 pulsations are mostly observed near the magnetic equator in the duskside outer magnetosphere. This L.T. asymmetry in the occurrence of dominant modes of Pc5's in space can be explained by the velocity shear instability (Yumoto and Saito, 1980) in the magnetospheric boundary layer, where Alfvénic signals in the IMF medium are assumed to penetrate into the magnetospheric boundary layer along the Archimedean spiral. The asymmetrical behaviour of Pc5 pulsation activity on the ground across the noon meridian can be also explained by the ionospheric screening effect on the compressional Pc5 magnetic pulsations. The compressional modes with a large horizontal wave number in the duskside magnetosphere are expected to be suppressed on the ground throughout the ionosphere and atmosphere.     

On photospheric and chromospheric penumbral waves

The full magnetoatmospheric wave equation is numerically solved for the particular penumbral model of Nye and Thomas (1974). For chromospheric running penumbral waves it is found that the maximum vertical velocity occurs at the base of the Hα region (z ≈ 300 km) much lower than previous WKB estimates have suggested. The maximum ‘vertical’ kinetic energy occurs at z ≈ - 130 km. Very significant horizontal velocities are also found for these waves, and, in the absence of shear flow it appears that previous estimates of photospheric vertical velocities of order 10^?1 km s^?1 could be substantial underestimates. For the photospheric events of Musman et al. (1976), a high vertical velocity maximum is found in the corona, and the modes appear highly dispersive for periods ≥ 220 s. The effects of a sinusoidal shear flow profile on running penumbral waves are examined, and it is found that the eigenvalues (horizontal wavenumber k) are changed little, but the eigenmodes become significantly distorted; the position of the vertical velocity peak rises compared to the zero flow case, and the velocity below that peak drops significantly. This effect may well cancel the increased estimates based on zero flow.     

The stability of vertical motion in the N-body circular Sitnikov problem

We present results about the stability of vertical motion and its bifurcations into families of 3-dimensional (3D) periodic orbits in the Sitnikov restricted N-body problem. In particular, we consider ν = N ? 1 equal mass primary bodies which rotate on a circle, while the Nth body (of negligible mass) moves perpendicularly to the plane of the primaries. Thus, we extend previous work on the 4-body Sitnikov problem to the N-body case, with N = 5, 9, 15, 25 and beyond. We find, for all cases we have considered with N ≥ 4, that the Sitnikov family has only one stability interval (on the z-axis), unlike the N = 3 case where there is an infinity of such intervals. We also show that for N = 5, 9, 15, 25 there are, respectively, 14, 16, 18, 20 critical Sitnikov periodic orbits from which 3D families (no longer rectilinear) bifurcate. We have also studied the physically interesting question of the extent of bounded dynamics away from the z-axis, taking initial conditions on x, y planes, at constant z(0) = z ₀ values, where z ₀ lies within the interval of stable rectilinear motions. We performed a similar study of the dynamics near some members of 3D families of periodic solutions and found, on suitably chosen Poincaré surfaces of section, “islands” of ordered motion, while away from them most orbits become chaotic and eventually escape to infinity. Finally, we solve the equations of motion of a small mass in the presence of a uniform rotating ring. Studying the stability of the vertical orbits in that case, we again discover a single stability interval, which, as N grows, tends to coincide with the stability interval of the N-body problem, when the values of the density and radius of the ring equal those of the corresponding system of N ? 1 primary masses.     

Microscale simulations of Venus’ convective adjustment and mixing near the surface: Thermal and material transport processes

Heat and material transport processes caused by convective adjustment and mixing are important in modeling of Venus’ atmosphere. In the present study, microscale atmospheric simulations near the venusian surface were conducted using a Weather Research and Forecasting model to elucidate the thermal and material transport processes of convective adjustment and mixing. When convective adjustment occurs, the heat and passive tracer are rapidly mixed into the upper stable layer with convective penetration. The convective adjustment produces large eddy diffusions of heat and passive tracer, which may explain the large eddy diffusions estimated in the radiative-convective equilibrium model.For values of surface heat flux Q greater than a threshold (=0.064 K m s⁻¹ in the present study), the convectively mixed layer with high eddy diffusion coefficients grows with time. In contrast, the mixed layer decays with time for Q values smaller than the threshold. The thermal structure near the surface is controlled not only by extremely long-term radiative processes, but also by microscale dynamics with time scales of several hours. A mixed layer with high eddy diffusion coefficients may be maintained or grow with time if the surface heat flux is high in the volcanic hotspot and adjacent areas.     

The abundance of acetylene in the atmosphere of Jupiter

The Q and R branches of the C₂H₂ ν₅ fundamental, observed in emission in an aircraft spectrum of Jupiter near 750 cm^?1, have been analyzed with the help of an improved line listing for this band. The line parameters have been certified in the laboratory with the same interferometer used in the Jovian observations. The maximum mixing ratio of C₂H₂ is found to be between 5 × 10^?8 and 6 × 10^?9, depending on the form of its vertical distribution and the temperature structure assumed for the lower stratosphere. Most consistent with observations of both Q and R branches are: (1) distributions of C₂H₂ with a constant mixing ratio in the stratosphere and a cutoff at a total pressure of 100 mbar or less, and (2) the assumption of a temperature at 10^?2 bar which is near 155°K.     

Jupiter's 24° N highest speed jet: Vertical structure deduced from nonlinear simulations of a large-amplitude natural disturbance

The evolution of a large-amplitude disturbance at cloud level in Jupiter's 24° N jet stream in 1990 is used to constrain the vertical structure of a realistic atmospheric model down to the 6 bar pressure level. We use the EPIC model (Dowling et al., 1998, The explicit planetary isentropic-coordinate (EPIC) atmospheric model, Icarus 132, 221-238) to perform long-term, three-dimensional, nonlinear simulations with a series of systematic variations in vertical structure and find that the details of the 1990 disturbance combine with the characteristics of the 24° N jet, the fastest on Jupiter, to yield a tight constraint on the solution space. The most important free parameters are the vertical dependence of the zonal-wind profile, and the thermal structure, below the cloud tops (p>0.7 bar) at the jet's central latitude. The temporal evolution of the disturbed cloud patterns, which spans more than 2 years, can be reproduced if the jet peak reaches ∼180 ms⁻¹ at the cloud level and increases to ∼210 ms⁻¹ at 1 bar and up to ∼240 ms⁻¹ at 6 bar; the observations were not reproduced for other configurations investigated. This trend is consistent with that measured by the Galileo Probe at 7° N; the implication is that this jovian jet extends well below the solar radiation penetration level situated near the 2 bar level.     

Correlation between the Lyman alpha absorption line density and the emission line redshift of QSOs

A statistical analysis of the absorption line density N(z_abs) in the Lαforest of 16 high-redshift quasars shows that the line density is clearly dependent on the emission redshift of the quasar, z_em. Not only is the mean density of all the Lα absorption lines larger for larger z_em, but more importantly, for each given absorption redshift interval, N(z_abs) is also positively correlated with z_em. A two-factor variance analysis, with ejection velocity as the second factor, shows that the line density in the Lα forest depends sensitively on z_em but not on the ejection velocity. These results are difficult to interpret by either the intervening or the ejecting hypothesis and are unlikely to be due to the selection effects suggested by Carswell et al. [13].     

Morphology and variability of the Titan ringlet and Huygens ringlet edges

We present a forward modeling approach for determining, in part, the ring particle spatial distribution in the vicinity of sharp ring or ringlet edges. Synthetic edge occultation profiles are computed based on a two-parameter particle spatial distribution model. One parameter, h, characterizes the vertical extent of the ring and the other, δ, characterizes the radial scale over which the ring optical depth transitions from the background ring value to zero. We compare our synthetic occultation profiles to high resolution stellar occultation light curves observed by the Cassini Ultraviolet Imaging Spectrograph (UVIS) High Speed Photometer (HSP) for occultations by the Titan ringlet and Huygens ringlet edges.More than 100 stellar occultations of the Huygens ringlet and Titan ringlet edges were studied, comprising 343 independent occultation cuts of the edges of these two ringlets. In 237 of these profiles the measured light-curve was fit well with our two-parameter edge model. Of the remaining edge occultations, 69 contained structure that could only be fit with extremely large values of the ring-plane vertical thickness (h > 1 km) or by adopting a different model for the radial profile of the ring optical depth. An additional 37 could not be fit by our two-parameter model.Certain occultations at low ring-plane incidence angles as well as occultations nearly tangent to the ring edge allow the direct measurement of the radial scale over which the particle packing varies at the edge of the ringlet. In 24 occultations with these particular viewing geometries, we find a wide variation in the radial scale of the edge. We are able to constrain the vertical extent of the rings at the edge to less than ∼300 m in the 70% of the occultations with appropriate viewing geometry, however tighter constraints could not be placed on h due to the weaker sensitivity of the occultation profile to vertical thickness compared to its sensitivity to δ.Many occultations of a single edge could not be fit to a single value of δ, indicating large temporal or azimuthal variability, although the azimuthal variation in δ with respect to the longitudes of various moons in the system did not show any discernible pattern.