Effects of clumping on temperature – I. Externally heated clouds

We present a study of radiative transfer in dusty, clumpy star-forming regions. A series of self-consistent, 3D, continuum radiative transfer models are constructed for a grid of models parametrized by central luminosity, filling factor, clump radius and face-averaged optical depth. The temperature distribution within the clouds is studied as a function of this parametrization. Among our results, we find that: (i) the effective optical depth in clumpy regions is less than in equivalent homogeneous regions of the same average optical depth, leading to a deeper penetration of heating radiation in clumpy clouds, and temperatures higher by over 60 per cent; (ii) penetration of radiation is driven by the fraction of open sky (FOS) – which is a measure of the fraction of solid angle along which no clumps exist; (iii) FOS increases as clump radius increases and as filling factor decreases; (iv) for values of   FOS >0.6–0.8  the sky is sufficiently 'open' that the temperature distribution is relatively insensitive to FOS; (v) the physical process by which radiation penetrates is preferentially through streaming of radiation between clumps as opposed to diffusion through clumps; (vi) filling factor always dominates the determination of the temperature distribution for large optical depths, and for small clump radii at smaller optical depths; (vii) at lower face-averaged optical depths, the temperature distribution is most sensitive to filling factors of 1–10 per cent, in accordance with many observations; (viii) direct shadowing by clumps can be important for distances approximately one clump radius behind a clump.     

Spatially intermittent fields in photospheric magnetoconvection

Motivated by recent high-resolution observations of the solar surface, we investigate the problem of non-linear magnetoconvection in a three-dimensional compressible layer. We present results from a set of numerical simulations which model the situation in which there is a weak imposed magnetic field. This weak-field regime is characterized by vigorous granular convection and spatially intermittent magnetic field structures. When the imposed field is very weak, magnetic flux tends to accumulate at the edges of the convective cells, where it forms compact, almost 'point-like' structures which are reminiscent of those observed in the quiet Sun. If the imposed field is slightly stronger, there is a tendency for magnetic flux to become concentrated into 'ribbon-like' structures which are comparable to those observed in solar plages. The dependence of these simulations upon the strength of the imposed magnetic field is analysed in detail, and the concept of the fractal dimension is used to make a further, more quantitative comparison between these simulations and photospheric observations.     

Shock Propagation Through Multiphase Media*

This paper presents two and three dimensional simulations of the interaction of shocks with media with large numbers of dense inclusions. An approximate model of the interaction of a starburst wind with the surrounding galactic ISM illustrates issues which must be addressed in global models of ISM dynamics. As a step towards developing the sub-grid model of multiphase turbulence, we define and study a form of ‘multiphase Riemann problem’. This allows us to develop macroscopic characteristics of the flows which may be compared to such subgrid models.     

Deep-ocean vertical noise directionality

The structure of beam noise measured at the output of a vertical array in a range dependent ocean basin was investigated using the modified wide-angle parabolic equation (PE). Noise sources were distributed throughout the basin, and the field due to each noise source at an array located in the midbasin was calculated. The response of the array to the superposition of the noise sources was found by beamforming. An efficient and direct approach that superimposes the noise sources on the PE field as the field is marched toward the array was developed. Downslope calculations of the midbasin vertical directionality were made between 50 and 400 Hz with this technique. Use of a geoacoustic model shows that the bottom behaves as a low-pass filter     

Inverted echo sounder measurement of dynamic height through an ENSOcycle in the central equatorial Pacific

A four-year record from an inverted echo sounder deployed near Palmyra Island at 6°N in the central Pacific Ocean is compared with a simultaneous record of subsurface pressure from this island lagoon. A factor m, converting round-trip acoustic travel time to surface dynamic height relative to a deep pressure level, was estimated from the ratio of the spectra of the two records in the energetic synoptic oscillation band. Year-to-year variation in m was not statistically significant. For the overall record, m was found to be -70±8 dynamic m/s, where the error bounds represent a 90% confidence interval. This is consistent with first-baroclinic-mode excitation     

The “Shiva Hypothesis”: Impacts, mass extinctions, and the galaxy

The “Shiva Hypothesis”, in which recurrent, cyclical mass extinctions of life on Earth result from impacts of comets or asteroids, provides a possible unification of important processes in astrophysics, planetary geology, and the history of life. Collisions with Earth-crossing asteroids and comets ≥ a few km in diameter are calculated to produce widespread environmental disasters (dust clouds, wildfires), and occur with the proper frequency to account for the record of five major mass extinctions (from ≥ 10⁸ Mt TNT impacts) and ~ 20 minor mass extinctions (from 10⁷–10⁸ Mt impacts) recorded in the past 540 million years. Recent studies of a number of extinctions show evidence of severe environmental disturbances and mass mortality consistent with the expected after-effects (dust clouds, wildfires) of catastrophic impacts. At least six cases of features generally considered diagnostic of large impacts (e.g., large impact craters, layers with high platinum-group elements, shock-related minerals, and/or microtektites) are known at or close to extinction-event boundaries. Six additional cases of elevated iridium levels at or near extinction boundaries are of the amplitude that might be expected from collision of relatively low-Ir objects such as comets. The records of cratering and mass extinction show a correlation, and might be explained by a combination of periodic and stochastic impactors. The mass extinction record shows evidence for a periodic component of about 26 to 30 Myr, and an ~ 30 Myr periodic component has been detected in impact craters by some workers, with recent pulses of impacts in the last 2–3 million years, and at ~ 35, 65, and 95 million years ago. A cyclical astronomical pacemaker for such pulses of impacts may involve the motions of the Earth through the Milky Way Galaxy. As the Solar System revolves around the galactic center, it also oscillates up and down through the plane of the disk-shaped galaxy with a half-cycle ~ 30±3 Myr. This cycle should lead to quasi-periodic encounters with interstellar clouds, and periodic variations in the galactic tidal force with maxima at times of plane crossing. This “galactic carrousel” effect may provide a viable perturber of the Oort Cloud comets, producing periodic showers of comets in the inner Solar System. These impact pulses, along with stochastic impactors, may represent the major punctuations in earth history.     

Three-dimensional linear solution for wave propagation with slopingbottom

This paper presents a three-dimensional analytic linear wave solution for surface gravity wave propagation over a sloping bottom that is valid for small, but realistic, slopes. The sloping-bottom linear model is compared to published laboratory data and to predictions of two-dimensional, constant-bottom nonlinear theories. The model is shown to describe the measured wave-height growth in the wave transformation region up to a limiting local Ursell number U_r of 0.35-1.0, depending on the wave type, although, as a linear model, it does not predict the harmonics observed in that range. For U_r<0.35, the harmonics can generally be neglected and the sloping-bottom linear theory agrees closely with both the published wave-height data and third-order Stokes nonlinear theory. As a three-dimensional linear model, superposition can be invoked to synthesize and relate wave structure in the transformation region to complex incident ocean spectra with both wind wave and swell components that arrive with a range of incidence angles. As such, the sloping-bottom linear model presented here should be a convenient useful tool for ocean modeling through a significant portion of the wave transformation region     

Ambipolar diffusion in the solar atmosphere

The process of non-linear ambipolar diffusion in the region overlying the solar surface can be an effective mechanism for producing sharp magnetic structures and current sheets. These may be the sites responsible for the occurrence of connectivity of magnetic field lines, and the subsequent explosive input of energy for heating of some of the features in the atmosphere of the Sun..