Numerical Method for Conditional Simulation of Levy Random Fields

Stochastic simulations of subsurface heterogeneity require accurate statistical models for spatial fluctuations. Incremental values in subsurface properties were shown previously to be approximated accurately by Levy distributions in the center and in the start of the tails of the distribution. New simulation methods utilizing these observations have been developed. Multivariate Levy distributions are used to model the multipoint joint probability density. Explicit bounds on the simulated variables prevent nonphysical extreme values and introduce a cutoff in the tails of the distribution of increments. Long-range spatial dependence is introduced through off-diagonal terms in the Levy association matrix, which is decomposed to yield a maximum likelihood type estimate at unobserved locations. This procedure reduces to a known interpolation formula developed for Gaussian fractal fields in the situation of two control points. The conditional density is not univariate Levy and is not available in closed form, but can be constructed numerically. Sequential simulation algorithms utilizing the numerically constructed conditional density successfully reproduce the desired statistical properties in simulations.     

“Extreme boiling” model for variable salinity of the Hokko low-sulfidation epithermal Au prospect, southwestern Hokkaido, Japan

The Hokko prospect is located in the Minamikayabe area southwestern Hokkaido, Japan, where gold-bearing quartz veins of Pliocene age are exposed at the surface. The alteration mineral assemblage is typical of low-sulfidation epithermal systems, with the quartz veins associated with adularia alteration overprinted on Late Miocene propylitic alteration. Fluid inclusion studies of the vein quartz reveal mean homogenization temperatures of approximately 220 °C, and the co-existence of low-salinity (<2 wt.% NaCl equivalent) and moderate salinity (2 to 12 wt.% NaCl equivalent) fluid inclusions within the same veins. The moderate salinity fluid inclusions (2–12 wt.% NaCl equivalent) typically have relatively low homogenization temperatures between 150° to 200 °C. The results obtained from stable isotope analysis of  δ¹⁸O in quartz vein material showed a gradual decrease in  δ¹⁸O signatures with increasing depth. The majority of the samples have calculated fluid source signatures (δ¹⁸O_H2O) between −8.0 and −10.0‰, but there is a significant change in the composition above 185 m drill depth. The shallower samples in particular show a wide range of oxygen isotope signatures that are associated with the moderate salinity fluid inclusions. We interpret that low-salinity inclusions within the Hokko system represent the composition of the liquid phase of the fluid, before boiling, and that the moderate-salinity inclusions are representative of the residual liquid phase, after extensive non-adiabatic boiling and vapor loss in an open system. This mechanism resulted in the entrapment of fluids with variable salinities at the same time, and in close proximity to each other. This is also reflected in the  δ¹⁸O_H2O values which become more variable and heavier where the moderate-salinity inclusions occur. Deposition of ore minerals within the Hokko vein system also occurred at this time as a result of boiling and gas loss. Received: 30 May 1997 / Accepted: 6 January 1998     

Linear response of galactic haloes to adiabatic gravitational perturbations

We determine the response of a self-similar isothermal stellar system to small adiabatic gravitational perturbations. For odd spherical harmonics the response is identical to the response of the analogous isothermal fluid system. For even spherical harmonics the response can be regarded as an infinite series of wavetrains in log r , implying alternating compression and rarefaction in equal logarithmic radius intervals. Partly because of the oscillatory nature of the solutions, tidal fields from external sources are not strongly amplified by an intervening isothermal stellar system, except at radii ≲10^−3.5 times the satellite radius; at some radii the stellar system can even screen the external tidal field in a manner analogous to Debye screening. As Weinberg has pointed out, individual resonances in a stellar system can strongly amplify external tidal fields over a limited radial range, but we cannot address this possibility because we examine only adiabatic perturbations. We also discuss the application of our method to the halo response caused by the slow growth of an embedded thin disc.     

The kinematics of the ionized gas in the Circinus galaxy

The TAURUS-2 Fabry–Perot interferometer, mounted on the 3.9-m Anglo-Australian Telescope, has been used to observe the Circinus galaxy. We have mapped the intensity and velocity distribution of the ionized hydrogen in the galaxy using the Balmer series Hα spectral line.
The semiresolved core (observed with a seeing disc of 30 pc) appears amorphous in shape, which is commonly observed in Seyfert 2 galaxies. Its peak coincides with the core position measured in the radio continuum, suggesting that ionized gas surrounds a non-thermal source.
A circumnuclear ring or spiral of radius 220 pc and a rotational velocity of 350 km s⁻¹ (assuming circular motions) surrounds the core. The inclination angle of this feature, i =40°±10°, is less than that of the previously observed radio continuum disc. The velocity channel maps obtained for the Hα ring show that the kinematics resemble those of a rotating ring and the intensity displays a complex structure indicative of several, unresolved, H II regions. We believe the ring to represent a circumnuclear starburst.
Our Hα data also show the presence of the previously detected [O III ] ionization cone to the north-west of the core, measuring more than 400 pc in length. We suggest that the ionization cone lies in a different plane from that of the starburst ring and is directed away from us. Several kinematic components of the core are derived and we calculate an outflow velocity in excess of 150–200 km s⁻¹ for gas above the core of Circinus. We also present evidence for inflowing ionized gas at the centre of Circinus.
The correlation of the Hα and radio continuum features is discussed, as well as the possible presence of a starburst-driven superwind in the Circinus galaxy.     

Oscillations in radiation-driven outflows

In this paper we investigate the dynamical behaviour of radiation-driven winds, specifically winds that arise when Compton scattering transfers momentum from the radiation field to the gas flow. Such winds occur during strong X-ray bursts from slowly accreting neutron stars, and also may be driven from the inner regions of a black hole or neutron star accretion disc when the mass transfer rate is very high. By linearizing the radiation hydrodynamic equations around steady spherical outflow, we evaluate the time-dependent response of these winds to perturbations introduced at their inner boundaries. We find that although radiation-driven winds are generally stable, they act as mechanical filters that should produce quasi-periodic oscillations or peaked noise in their radiation output when perturbations force them stochastically. This behaviour may underlie the photospheric oscillations observed in some strong Type I X-ray bursts.     

Mean crustal velocity: a critical parameter for interpreting crustal structure and crustal growth

Mean crustal velocity is a critical parameter for genesis of continental crystalline crust because it is a function of mean crustal composition and therefore may be used to resolve continental crustal growth in space and time. Although the best values of mean crustal velocity are determined from wide-angle reflection measurements, most studied here necessarily come from vertical averages in crustal refraction determinations. The mode of 158 values of mean crustal velocity is 6.3 km/s, a velocity which corresponds to a mean crustal composition of granodiorite to felsic quartz diorite; Archean crust may be slightly more mafic. Mean crustal velocities range from 5.8 to 7.0 km/s. The lowest values invariably are found in thermally disturbed rift zones and the highest values correspond to velocities in gabbro. Velocities in island arcs may be as low as 6.0 km/s but are typically 6.5–6.9 km/s which corresponds to andesitic composition; estimates of island arc composition are andesitic. If values of mean crustal velocity are not biased, this observation suggests that continental crust did not grow simply by addition of island arc material. Possibilities are that crust formed from fusion of island arcs and was later changed to more felsic composition by addition of material from the mantle or that the late Archean episode of major crustal growth did not involve processes similar to younger island arcs. Some crustal blocks might be changed in composition and thickness by such processes as underplating, interthrusting, necking and sub-crustal erosion. Specially designed experiments are suggested to determine this parameter so critical for understanding genesis of continental crust.     

Tungsten in iron meteorites

Tungsten concentrations have been determined by instrumental neutron activation in 104 iron meteorites, and range from 0.07 to 5 μg/g. In individual groups, concentrations vary by factors of between 1.5 and 8, but there are negative W-Ni correlations in 8 groups: IAB, IC, IIAB, IID, IIE, IIIAB, IIICD, and IIIF. The lowest W concentrations are found in groups IAB and IIICD, which also have the smallest slopes on a W-Ni plot. Eighteen anomalous irons have W concentrations between 5 μg/g (Butler) and 0.11 μg/g (Rafrüti). The distribution of W in irons shows similarities to that of other refractory siderophilic elements (except Mo), but is closest to the distribution of Ru and Pt.Assuming that chemical trends in group IIIAB were produced by fractional crystallization, a value of 1.6 can be deduced for the distribution coefficient of W between solid and liquid metal, cf. 0.89 for Mo. Experimental evidence in support of these values is tenuous.