Characterization of slip zone associated with the 1999 Taiwan Chi-Chi earthquake: X-ray CT image analyses and microstructural observations of the Taiwan Chelungpu fault

To characterize the fault-related rocks within the Chelungpu fault, we performed X-ray computed tomography (CT) image analyses and microstructural observations of Hole B core samples from the Taiwan Chelungpu-fault Drilling Project. We identified the slip zone associated with the 1999 Chi-Chi earthquake, within the black gouge zone in the shallowest major fault zone, by comparison with previous reports. The slip zone was characterized by low CT number, cataclastic (or ultracataclastic) texture, and high possibility to have experienced a mechanically fluidized state. Taking these characteristics and previous reports of frictional heating in the slip zone into consideration, we suggested that thermal pressurization was the most likely dynamic weakening mechanism during the earthquake.     

Decision-theoretic sensitivity analysis for reservoir development under uncertainty using multilevel quasi-Monte Carlo methods

At various stages of petroleum reservoir development, we encounter a large degree of geological uncertainty under which a rational decision has to be made. In order to identify which parameter or group of parameters significantly affects the output of a decision model, we investigate decision-theoretic sensitivity analysis and its computational issues in this paper. In particular, we employ the so-called expected value of partial perfect information (EVPPI) as a sensitivity index and apply multilevel Monte Carlo (MLMC) methods to efficient estimation of EVPPI. In a recent paper by Giles and Goda, an antithetic MLMC estimator for EVPPI is proposed and its variance analysis is conducted under some assumptions on a decision model. In this paper, for an improvement on the performance of the MLMC estimator, we incorporate randomized quasi-Monte Carlo methods within the inner sampling, which results in an multilevel quasi-Monte Carlo (MLQMC) estimator. We apply both the antithetic MLMC and MLQMC estimators to a simple waterflooding decision problem under uncertainty on absolute permeability and relative permeability curves. Through numerical experiments, we compare the performances of the MLMC and MLQMC estimators and confirm a significant advantage of the MLQMC estimator.     

The formation of floor-fractured craters in Xanthe Terra

Floor-fractured craters (FFC) are a peculiar form of degradation of impact craters defined by the presence of crevice networks and mesas affecting crater floors. They are preferentially distributed near chaotic terrains and outflow channels. The scope of this paper is to present a detailed systematic analysis of FFC at Xanthe Terra. FFC morphologies in this region are classified into five types making a picture of different stages of the same degradation process. FFC are geographically intermixed with un-fractured normal craters (non-FFC). Young craters are less prone to show this type of degradation, as suggested by fresh ejecta layer with preserved crater floor. Size distributions of FFC and non-FFC indicate that larger craters are preferentially fractured. Careful examinations of the crater floor elevations reveal that the crevices often extend deeper than the original crater cavity. Furthermore, an onset depth for the formation of FFC is evidenced from the difference of spatial distributions between FFC and non-FFC. Roof-collapsed depressions observed in the same region have been also documented and their characteristics suggest the removal of subsurface material at depth from about 1200 to 4000 m. These observations taken together suggest a subsurface zone of volume deficit at depth from 1 to 2 km down to several kilometers responsible for FFC formation. Then a scenario of FFC formations is presented in the context of groundwater discharge events at the late Hesperian. This scenario involves two key processes, Earth fissuring and piping erosion, known to occur with rapid groundwater migrations on Earth.     

Numerical Investigations of Mean Winds Within Canopies of Regularly Arrayed Cubical Buildings Under Neutral Stability Conditions

Recently, several attempts have been made to model the wind velocity in an urban canopy in order to accurately predict the mixing and transport of momentum, heat, and pollutants within and above the canopy on an urban scale. For this purpose, unverified assumptions made by Macdonald (Boundary-Layer Meteorol 97:25–45, 2000) to develop a model for the profile of the mean wind velocity within an urban canopy have been used. In the present study, in order to provide foundations for improving the urban canopy models, the properties of the spatially-averaged mean quantities used to make these assumptions have been investigated by performing large-eddy simulations (LES) of the airflow around square and staggered arrays of cubical blocks with the following plan area densities: λ _p = 0.05, 0.11, 0.16, 0.20, 0.25, and 0.33. The LES results confirm that the discrepancy between the spatial average of wind velocity and Macdonald’s five-point average of wind velocity can be large in both types of arrays for large λ _p . It is also confirmed that Prandtl’s mixing length varies significantly with height within the canopy, contrary to Macdonald’s assumption for both types of arrays and for both small and large λ _p . On the other hand, in accordance with Macdonald’s assumption, the sectional drag coefficient is found to be almost constant with height except in the case of staggered arrays with high λ _p .     

The origin of chaotic behaviour in the Miranda-Umbriel 3 : 1 resonances

We have investigated the pericentric resonances through which Miranda and Umbriel are believed to have passed when, due to tidal evolution, their orbital mean motions reached a 3 : 1 commensurability. Our investigation is based upon a perturbative treatment. The predictions of this theory are in good agreement with the results of numerical integrations concerning both the extend of the chaotic layers generated by the separatrices of the primary resonances and the location of the secondary resonances. The effect of tidal evolution is discussed on the bases of the adiatatic invariant theory and its extension to separatrix crossing. We recover qualitatively the mean features of the numerical experiments of Tittermore and Wisdom (1988–1989), Dermott et al (1988) and Malhotra and Dermott (1989).     

The quasi-stationary coronal magnetic field and electron density as determined from a Faraday rotation experiment

Pioneer VI was launched into a circumsolar orbit on December 16, 1965, and was occulted by the sun in the latter half of November, 1968. During the occultation period, the 2292-MHz S-band telemetry carrier underwent Faraday rotation due to the interaction of this signal with the plasma and magnetic field in the solar corona. The NASA/JPL 210-ft diameter antenna of the Deep Space Network near Barstow, California, was used for the measurement. The antenna feed was modified for automatic polarization tracking for this experiment. The measurement results are interpreted with a theoretical model of the solar corona. This model consists of a modified Allen-Baumbach electron density and a coronal magnetic field calculated both from Mount Wilson magnetograph observations using a source surface model and field extrapolations from the Explorer 33 satellite magnetometer. The observations and the calculated rotation show general agreement with respect to magnitude, sense, and timing, suggesting the source-surface model and field extrapolations from 1 AU are a valid technique to obtain the magnetic field in the corona from 4 to 12 solar radii. Variations present can easily be ascribed to density enhancements known to be present in the corona. Longitudinal variations of the density in the corona cannot be obtained from coronagraph observations, and thus a purely radial variation was assumed. An improved fit to the Faraday rotation data is obtained with an equatorial electron density $$N = 10^8 \left( {\frac{{6000}}{{R^{10} }} + \frac{{0.002}}{{R^2 }}} \right)...{\text{ cm}}^{{\text{ - 3}}} {\text{ (4 < }}R < 12){\text{ }}...$$ where R is in solar radii. The work of W. V. T. Rusch and J. E. Ohlson was supported in part by research sponsored by the Joint Services Electronics Program through the Air Force Office of Scientific Research under Grant AF-AFOSR 69-1622A at the University of Southern California. The work done by K. H. Schatten was in part supported by the National Academy of Science on a National Research Council postdoctoral fellowship. The work of J. M. Wilcox was supported in part by the Office of Naval Research under Contract Nonr 3656(26), by the National Aeronautics and Space Administration under Grant NGR 05-003-230, and by the National Science Foundation under Grant GA-1319 at the University of California at Berkeley.     

Laboratory analogy of crystalline enstatite grain formation in plasma field

Crystalline enstatite (MgSiO₃) grains were produced by the simultaneous evaporation of SiO grains and Mg vapor in a plasma field. The MgSiO₃ grains were spherical or needlelike. The necessity of a plasma field in astromineralogy is suggested in the present study.     

A spectroscopic study of Epsilon Aurigae

Epsilon Aurigae has been observed during ingress and totality between 1982 and 1983 at Okayama. Analyses of profiles of H line and of radial velocities of neutral hydrogen and metals show that the secondary component consists of at least three parts in structure.Paper presented at the Lembang-Bamberg IAU Colloquium No. 80 on Double Stars: Physical Properties and Generic Relations, held at Bandung, Indonesia, 3–7 June, 1983.