Photometric modeling of Saturn's rings: I. Monte Carlo method and the effect of nonzero volume filling factor     

Heikki Salo  Raine Karjalainen 《Icarus》2003,164(2):428-460

The scattering properties of particulate rings with volume filling factors in the interval D=0.001-0.3 are studied, with photometric Monte Carlo ray tracing simulations combining the advantages of direct (photons followed from the source) and indirect methods (brightness as seen from the observing direction). Besides vertically homogeneous models, ranging from monolayers to classical many-particle thick rings, particle distributions obtained from dynamical simulations are studied, possessing a nonuniform vertical profile and a power law distribution of particle sizes. Self-gravity is not included to assure homogeneity in planar directions. Our main goal is to check whether the moderately flattened ring models predicted by dynamical simulations (with central plane D>0.1) are consistent with the basic photometric properties of Saturn's rings seen in ground-based observations, including the brightening near zero phase angle (opposition effect), and the brightening of the B-ring with increasing elevation angle (tilt effect). Our photometric simulations indicate that dense rings are typically brighter in reflected light than those with D→0, due to enhanced single scattering. For a vertically illuminated layer of identical particles this enhancement amounts at intermediate viewing elevations to roughly 1+2D. Increased single scattering is also obtained for low elevation illumination, further augmented at low phase angles α by the opposition brightening when D increases: the simulated opposition effect agrees very well with the Lumme and Bowell (1981, Astron. J. 86, 1694-1704) theoretical formula. For large α the total intensity may also decrease, due to reduced amount of multiple scattering. For the low (α=13°) and high (α=155°) phase angle geometries analyzed in Dones et al. (1993, Icarus 105, 184-215) the brightness change for D=0.1 amounts to 20% and −17%, respectively. In the case of an extended size distribution, dynamical simulations indicate that the smallest particles typically occupy a layer several times thicker than the largest particles. Even if the large particles form a dynamically dense system, a narrow opposition peak can arise due to mutual shadowing among the small particles: for example, a size distribution extending about two decades can account for the observed about 1° wide opposition peak, solely in terms of mutual shadowing. The reduced width of the opposition peak for extended size distribution is in accordance with Hapke's (1986, Icarus 67, 264-280) treatment for semi-infinite layers. Due to vertical profile and particle size distribution, the photometric behavior is sensitive to the viewing elevation: this can account for the tilt-effect of the B-ring, as dense and thus bright central parts of the ring become better visible for larger elevation, whereas in the case of smaller elevation, mainly low volume density upper layers are visible. Since multiple scattering is not involved, the explanation works also for albedo well below unity. Inclusion of nonzero volume density helps also to model some of the Voyager observations. For example, the discrepancy between predicted and observed brightness at large phase angles for much of the A-ring (Dones et al., 1993, Icarus 105, 184-215) is removed when the enhanced low α single scattering and reduced large α multiple scattering is allowed for. Also, a model with vertical thickness increasing with saturnocentric distance offers at least a qualitative explanation for the observed contrast reversal between the inner and outer A-ring in low and high phase Voyager images. Differences in local size distribution and thus on the effective D may also account for the contrast reversal in resonance sites.  相似文献   

HST observations of spokes in Saturn's B ring     

Colleen A. McGhee  Richard G. French  Jeffrey N. Cuzzi  Rebecca Danos 《Icarus》2005,173(2):508-521

As part of a long-term study of Saturn's rings, we have used the Hubble Space Telescope's (HST) Wide Field and Planetary Camera (WFPC2) to obtain several hundred high resolution images from 1996 to 2004, spanning the full range of ring tilt and solar phase angles accessible from the Earth. Using these multiwavelength observations and HST archival data, we have measured the photometric properties of spokes in the B ring, visible in a substantial number of images. We determined the spoke particle size distribution by fitting the wavelength-dependent extinction efficiency of a prominent, isolated spoke, using a Mie scattering model. Following Doyle and Grün (1990, Icarus 85, 168-190), we assumed that the spoke particles were sub-micron size spheres of pure water ice, with a Hansen-Hovenier size distribution (Hansen and Hovenier, 1974, J. Atmos. Sci. 31, 1137-1160). The WFPC2 wavelength coverage is broader than that of the Voyager data, resulting in tighter constraints on the nature of spoke particles. The effective particle size was reff=0.57±0.05 μm, and the size distribution was quite narrow with a variance of b=0.09±0.03, very similar to the results of Doyle and Grün (1990, Icarus 85, 168-190), and consistent with predictions of plasma cloud models for spoke production from meteoritic impacts (Goertz and Morfill, 1983, Icarus 53, 219-229; Goertz, 1984, Adv. Space Res. 4, 137-141). In all, we identified 36 spokes or spoke complexes, predominantly on the morning (east) ansa. The photometric contrast of the spokes is strongly dependent on effective ring opening angle, Beff. Spokes were clearly visible on the north face of the rings in 1994, just prior to the most recent ring plane crossing (RPX) epoch, and on the south face shortly after RPX. However, spokes were both less abundant and fainter as the rings opened up, and no spokes were detected after 18 October 1998 (Beff=−15.43°), when a single faint spoke was seen on the morning ansa. The high resolution and photometric quality of the WFPC2 images enabled us to set a detection limit of ?1% in fractional brightness contrast for spokes for the post-1998 observations. We compare the observed trend of spoke contrast with Beff to radiative transfer calculations based on three models of the distribution of spoke material. In the first, the spoke “haze” is uniformly mixed with macroscopic B ring particles. No variation in spoke contrast is predicted for single-scattering, in this case, and only a modest decrease in contrast with Beff is predicted when multiple scattering is taken into account. In the second model, the spoke dust occupies an extended layer that is thicker than the B ring, which gives virtually identical results to a third case, when the haze layer lies exclusively above the ring. Multiple-scattering Monte Carlo calculations for these two extended haze models match the trend of spoke contrast exceptionally well. We compute the predicted spoke contrast for a wide variety of viewing geometries, including forward- and backscattering. Based on these results, spokes should be easily detectable during the Cassini mission when the rings are viewed at relatively small (|B|?10°) ring opening angles.  相似文献   

Mass transfer and preferred orientation development during extensional microcracking in slate-belt folds, Elura Mine, Australia     

J. A. DE ROO 《Journal of Metamorphic Geology》1989,7(3):311-322

Microstructures in slate belt rocks at the Elura Mine, near Cobar, south-eastern Australia, indicate that volume loss by syntectonic dissolution is coupled with mass accretion by reprecipitation of the dissolved material in dilational sites. The mass accretion is sustained primarily by repetitive tensile microfracturing at high pore-fluid pressures. Oriented growth in the inter- and intragranular microcracks is locally host-controlled, creating lattice- and shape-preferred orientations. The grain-scale crack-seal features throughout the rock reflect rhythmic fluid pressure fluctuations; a balance is achieved between the fracture-induced permeability (and consequent flushing rates), and the rate of fluid build-up in a relatively sealed environment.
Instability in the balancing factors can lead to localization and intensification of tensile failure (and hence, tension vein formation) in the grain aggregate. Growth of veins by crack-seal also reflects a steady state, but with more localized fluctuations of fluid flow on the aggregate scale. Still larger imbalances between flushing and fluid accumulation (i.e. pressure variations) induce breccia veining. The larger pressure gradients over greater distances, associated with dilation localization (from pervasive microfracturing to spaced breccia domains), allow fluid channelling with an increased potential for chemical fluid/rock disequilibrium. Therefore, large breccia vein systems tend to be sites of extensive fluid/rock interaction and replacement, as spectacularly illustrated by the syntectonic sulphide orebodies at Elura. The huge amounts of silicate, carbonate and sulphide accumulated during folding at Elura illustrate the large scale of source and sink couples possible in solute mass transfer.  相似文献   

Modeling the Lyα radiation of high-redshift galaxies     

Sabine Richling 《Astrophysics and Space Science》2003,284(2):361-364

The Lyα line emission of high-redshift galaxies depends on the density and temperature distribution of the gas, the kinematics and the dust content. We use a finite element method to model theLyα radiation of different 3D configurations considering complete frequency redistribution and the influence of velocity fields. Our results show that the central absorption feature of the double-peaked Lyα line profile observed in many radio galaxies with z=2-4 is probably the consequence of frequency redistribution rather than foreground absorption. The blue peak of the profile is enhanced for models with in fall motion and the red peak for models with outflow motion. In particular, we attempt to model the extendedLyα emission of high-redshift radio galaxies, where we consider results of corresponding hydrodynamical simulations to select possible model configurations. We find that Lyα photons scattered outside a jet-influenced low-density region are able to produce an extended Lyα halo. This revised version was published online in August 2006 with corrections to the Cover Date.  相似文献   

Smoothed particle hydrodynamics with radiation transfer   总被引：1，自引：0，他引：1  

S. Oxley  M. M. Woolfson 《Monthly notices of the Royal Astronomical Society》2003,343(3):900-912

  相似文献   

Methane absorption in the atmosphere of Jupiter from 1800 to 9500 cm and implications for vertical cloud structure     

P.G.J. Irwin  K. Sihra  F.W. Taylor 《Icarus》2005,176(2):255-271

New measurements of the low-temperature near-infrared absorption of methane (Sihra, 1998, Laboratory measurements of near-infrared methane bands for remote sensing of the jovian atmosphere, Ph.D. thesis, University of Oxford) have been combined with existing, longer path-length, higher-temperature data of Strong et al. (1993, Spectral parameters of self- and hydrogen-broadened methane from 2000 to 9500 cm⁻¹ for remote sounding of the atmosphere of Jupiter, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) and fitted with band models. The combined data set is found to be more consistent with previous low-temperature methane absorption measurements than that of Strong et al. (1993, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) but covers the same wider wavelength range and accounts for both self- and hydrogen-broadening conditions. These data have been fitted with k-coefficients in the manner described by Irwin et al. (1996, Calculated k-distribution coefficients for hydrogen- and self-broadened methane in the range 2000-9500 cm⁻¹ from exponential sum fitting to band modelled spectra, J. Geophys. Res. 101, 26,137-26,154) and have been used in multiple-scattering radiative transfer models to assess their impact on our previous estimates of the jovian cloud structure obtained from Galileo Near-Infrared Mapping Spectrometer (NIMS) observations (Irwin et al., 1998, Cloud structure and atmospheric composition of Jupiter retrieved from Galileo NIMS real-time spectra, J. Geophys. Res. 103, 23,001-23,021; Irwin et al., 2001, The origin of belt/zone contrasts in the atmosphere of Jupiter and their correlation with 5-μm opacity, Icarus 149, 397-415; Irwin and Dyudina, 2002, The retrieval of cloud structure maps in the equatorial region of Jupiter using a principal component analysis of Galileo/NIMS data, Icarus 156, 52-63). Although significant differences in methane opacity are found at cooler temperatures, the difference in the optical depth of the atmosphere due to methane is found to diminish rapidly with increasing pressure and temperature and thus has negligible effect on the cloud structure inferred at deeper levels. Hence the main cloud opacity variation is still found to peak at around 1-2 bar using our previous analytical approach, and is thus still in disagreement with Galileo Solid State Imager (SSI) determinations (Banfield et al., 1998, Jupiter's cloud structure from Galileo imaging data, Icarus 135, 230-250; Simon-Miller et al., 2001, Color and the vertical structure in Jupiter's belts, zones and weather systems, Icarus 154, 459-474) which place the main cloud deck near 0.9 bar. Further analysis of our retrievals reveals that this discrepancy is probably due to the different assumptions of the two analyses. Our retrievals use a smooth vertically extended cloud profile while the SSI determinations assume a thin NH₃ cloud below an extended haze. When the main opacity in our model is similarly assumed to be due to a thin cloud below an extended haze, we find the main level of cloud opacity variation to be near the 1 bar level—close to that determined by SSI and moderately close to the expected condensation level of ammonia ice of 0.85 bar, assuming that the abundance of ammonia on Jupiter is (7±1)×¹⁰⁻⁴ (Folkner et al., 1998, Ammonia abundance in Jupiter's atmosphere derived from the attenuation of the Galileo probe's radio signal, J. Geophys. Res. 103, 22,847-22,855; Atreya et al., 1999, A comparison of the atmospheres of Jupiter and Saturn: deep atmospheric composition, cloud structure, vertical mixing, and origin, Planet. Space Sci. 47, 1243-1262). However our data in the 1-2.5 μm range have good height discrimination and our lowest estimate of the cloud base pressure of 1 bar is still too great to be consistent with the most recent estimates of the ammonia abundance of 3.5 × solar. Furthermore the observed limited spatial distribution of ammonia ice absorption features on Jupiter suggests that pure ammonia ice is only present in regions of localised vigorous uplift (Baines et al., 2002, Fresh ammonia ice clouds in Jupiter: spectroscopic identification, spatial distribution, and dynamical implications, Icarus 159, 74-94) and is subsequently rapidly modified in some way which masks its pure absorption features. Hence we conclude that the main cloud deck on Jupiter is unlikely to be composed of pure ammonia ice and instead find that it must be composed of either NH₄SH or some other unknown combination of ammonia, water, and hydrogen sulphide and exists at pressures of between 1 and 2 bar.  相似文献   

Updated Opacity Project radiative accelerations     

M. J. Seaton 《Monthly notices of the Royal Astronomical Society》2007,382(1):245-250

Updated data for the calculation of radiative accelerations, g _rad, referred to as OP 2005, have been made generally available on the Web. They differ from the earlier 1997 OP data mainly in the inclusion of contributions from inner-shell processes. The frequency resolution used for OP 2005 is shown to be satisfactory except for extreme cases of low concentrations and low densities. The 2005 data are in reasonably good agreement with results from OPAL, given by Richer et al. in graphical form, for a model with   T _eff= 10⁴  and  log ( R ) =−3  . They are in less good agreement with OPAL results of Turcotte et al. for the solar radiative interior.  相似文献   

The effects of clumping on derived abundances in H ii regions     

John S. Mathis  Kenneth Wood 《Monthly notices of the Royal Astronomical Society》2005,360(1):227-235

  相似文献   

Colliding Blast Waves Driven by the Interaction of a Short-Pulse Laser with a Gas of Atomic Clusters     

Roland A. Smith  James Lazarus  Matthias Hohenberger  Alastair S. Moore  Joseph S. Robinson  Edward T. Gumbrell  Mike Dunne 《Astrophysics and Space Science》2007,307(1-3):131-137

Collisions between shocks are commonly found in many astrophysical objects, however robust numerical models or laboratory analogues of these complex systems remain challenging to implement. We report on the development of scaled laboratory experiments which employ new techniques for launching and diagnosing colliding shocks and high Mach number blast waves, scalable to a limited subset of astrophysically-relevant regimes. Use of an extended medium of atomic clusters enables efficient (>80%) coupling of 700 fs, 1 J, 1054 nm laser pulses to a “cluster” gas with an average density of ≈10¹⁹ particles cm⁻³, producing an initial energy density >10⁵ J cm⁻³, equivalent to ≈5×10⁹ J/g. Multiple laser foci are used to tailor the spatial profile of energy deposition, or to launch pairs of counter-propagating cylindrical shocks which then collide. By probing the collision interferometrically at multiple view angles in 5^{^} increments and applying an inverse Radon transform to the resulting phase projections we have been able to tomographicall reconstruct the full three-dimensional, time-framed electron density profile of the system.  相似文献   

Radiative Transfer in Inhomogeneous Atmospheres. II     

A. G. Nikoghossian 《Astrophysics》2004,47(2):248-259

This paper is a continuation of a study of radiative transfer in one-dimensional inhomogeneous atmospheres. Two of the most important characteristics of multiple scattering in these media are calculated: the photon escape probability and the average number of scattering events. The latter is determined separately for photons leaving the medium and for photons that have undergone thermalization in the medium. The problem of finding the radiation field in an inhomogeneous atmosphere containing energy sources is also examined. It is assumed that the power of these sources, as well as the scattering coefficient, can vary arbitrarily with depth. It is shown that knowledge of the reflection and transmission coefficients of the atmosphere makes it possible to reduce all these problems to solving some first order linear differential equations with specified initial conditions. A series of new analytic results are obtained. Numerical calculations are done for two types of atmosphere with different depth dependences for the scattering coefficient. These are interpreted physically.  相似文献