Laboratory measurements of ‘porosity‐free’ intrinsic Vp and Vs in an olivine gabbro of the Oman ophiolite: Implication for interpretation of the seismic structure of lower oceanic crust

In order to determine ‘porosity‐free’ intrinsic ultrasonic compressional (Vp) and shear wave (Vs) velocities and Vp/Vs of an olivine gabbro from the Oman ophiolite, we developed a new experimental system using a piston‐cylinder type high‐pressure apparatus. The new system allowed us to measure velocities at pressures ranging from 0.20 to 1.00 GPa and at temperatures up to 300°C for Vp and 400°C for Vs. At room temperature, the Vp and Vp/Vs increase rapidly with pressure up to 0.40 GPa, while between 0.45 and 1.00 GPa the increase is more gradual. The change in increasing rate is attributed to closure of porosity at pressures above 0.45 GPa. Based on the linear regression of data obtained at higher pressures (0.45–1.00 GPa) and extrapolation to the lower pressures, combined with temperature derivatives of velocities of the sample measured at 1.00 GPa, we determined the intrinsic Vp and Vs of the sample as a function of pressure (P, in GPa) and temperature (T, in °C). The intrinsic velocities can be expressed as Vp (km/s) = 7.004 + 0.096 × P ? 0.00015 × T, and Vs (km/s) = 3.827 + 0.007 × P ? 0.00008 × T. We evaluated the intrinsic Vp and Vs of the olivine gabbro at oceanic crustal conditions and compared them with a velocity depth‐profile of the borehole seismic observatory WP‐2 area in the northwestern Pacific Basin. Although the intrinsic Vp (～7.0 km/s) and Vs (～3.8 km/s) for the olivine gabbro studied are comparable to those of seismic layer 3 in the WP‐2 area, the estimated vertical gradients of intrinsic velocities are significantly smaller than those reported from layer 3. These results suggest that velocity profiles of layer 3 in the WP‐2 area may reflect the presence of a minor porosity in lower oceanic crust, which closes with increasing depth and/or continuous changes in mineralogy of layer 3 rocks.     

Subhorizontal boundary between ultrahigh-pressure and low-pressure metamorphic units in the Sulu–Tjube area of the Kokchetav Massif, Kazakhstan

Abstract The Kokchetav Massif of Kazakhstan includes high to ultrahigh-pressure (HP–UHP) metamorphic rocks (some of which were recrystallized at depths in excess of 150 km), juxtaposed against much lower pressure metamorphic components. We investigated the relationship between the HP–UHP metamorphic unit and the low pressure (LP) unit (Daulet Suite) in the Sulu–Tjube area, where the metamorphic rocks have previously been interpreted as constituting a megamelange with subvertical structural attitudes. Analyses of fold structures suggest that the HP–UHP metamorphic unit overlies the LP unit across a west-dipping subhorizontal boundary. In addition, kinematic indicators display top-to-the-north senses of shear along the tectonic contact between the two units, indicating that the HP–UHP unit has been extruded northward onto the LP unit. Following the juxtaposition of the two units, upright folds developed in both units, and these are associated with the previously reported steeply dipping metamorphic foliations. These data have important implications for the mode of exhumation of the UHP rocks from upper mantle to shallow crustal depths.     

Content and behavior of fluorine in Japanese quaternary volcanic rocks and petrogenetic application

Fluorine contents in about 160 representative Quaternary volcanic rocks and 15 hornblende and biotite phenocrysts in a calc-alkali series in Japan have been determined by a selective ion-electrode method. Tholeiites have the lowest contents and the narrowest range (58–145 ppm), while alkali basalts have the highest contentws and the widest range (301–666 ppm), high-alumina basalts have intermediate values (188–292 ppm). F contents in basalts clearly increase from east to west across the Japanese Islands, as do alkalies, P₂O₅ REE, U, Th and H₂O.The volcanic rocks studied are divided into two groups on the basis of F: (1) witt, increasing % SiO₂ or advancing fractionation, F contents show either progressive enrichment; or (2) with increasing fractionation, F contents show rather constant values. The former is produced by fractionation of anhydrous phases from basalt to mafic andesite magmas; the tholeiite series of Nasu volcanic zone (outer zone), northeastern, Japan is a typical example. The latter group is derived through separation of amphibole-bearing phases from basaltic magmas at various depths from upper mantle (about 30 km) to upper crust; the alkali series in southwestern Japan and the calc-alkali series of Chokai volcanic zone (inner zone), northeastern Japan, are examples.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

Radiative characteristics in a Japanese forested drainage basin during snowmelt

Radiative characteristics in a forested drainage basin during the snowmelt season were examined in order to better understand and predict snowmelt runoff in the basin. A method for estimating net radiation in a forest (R_nf) was presented using the total sky view factor (P) and the sun path sky view factor (Q). Solar radiation, albedo, atmospheric radiation and air temperature observed at an open site were also required. The total and the sun path sky view factors were determined from all‐sky photographs. Q was expressed as a linear function of P for 0·15<P<0·86 regardless of forest type. For P<0·15, Q was set to zero, and for P>0·86, Q was equal to unity. The short‐wave radiation budget at the forest floor (S_nf) increased with P, whereas the long‐wave radiation budget (L_nf) decreased with P. R_nf increased with P for 0·15<P<0·86, and changed little with P for P<0·15 and P>0·86, as the increase in S_nf was offset by the decrease in L_nf . The forest effect on R_nf was diminished under cloudy or high albedo conditions, because S_nf was easily offset by L_nf . This estimation method was extended to the whole basin, and R_nf was obtained over a watershed covered by trees. At the beginning of the snowmelt season when the albedo remained high, the forest effect became null because the decrease in S_nf was balanced by the increase in L_nf . As the albedo gradually lowered with the advance of the snowmelt season, the decrease in S_nf owing to forest covers exceeded the increase in L_nf , and the forest effect to decrease R_nf became evident. Copyright © 1999 John Wiley & Sons, Ltd.     

Damage analysis of simulated I-880 structures under the Loma Prieta Earthquake

This paper analyzes the collapse of the Cypress portion of the Nimitz freeway in Oakland, California, during the Loma Prieta Earthquake of 17 October 1989. The most important bent configuration is chosen to carry out a nonlinear dynamic analysis and damage evaluation using the program SARCF-III. The acceleration history recorded at a nearby location with similar soil properties is used as the input ground motion. In order to cover the uncertainty about the structural properties and reinforcing details six different structures models are utilized. The maximum dynamic internal forces are compared to the structural capacity. It is investigated whether the collapse can be explained using a simplified structural model and one of the available computer codes.     

Blended-acquisition design of irregular geometries towards faster,cheaper, safer and better seismic surveying

The application of blended acquisition has drawn considerable attention owing to its ability to improve the operational efficiency as well as the data quality and health, safety and environment performance. Furthermore, the acquisition of less data contributes to the business aspect, while the desired data density is still realizable via subsequent data reconstruction. The use of fewer detectors and sources also minimizes operational risks in the field. Therefore, a combined implementation of these technologies potentially enhances the value of a seismic survey further. One way to encourage this is to minimize any imperfection in deblending and data reconstruction during processing. In addition, one may derive survey parameters that enable a further improvement in these processes as introduced in this study. The proposed survey design workflow iteratively performs the following steps to derive the survey parameters responsible for source blending as well as the spatial sampling of detectors and sources. The first step is the application of blending and sampling operators to unblended and well-sampled data. We then apply closed-loop deblending and data reconstruction. The residue for a given design from this step is evaluated and subsequently used by genetic algorithms to simultaneously update the survey parameters related to both blending and spatial sampling. The updated parameters are fed into the next iteration until they satisfy the given termination criteria. We also propose a repeated encoding sequence to form a parameter sequence in genetic algorithms, making the size of problem space manageable. The results of the proposed workflow are outlined using blended dispersed source array data incorporating different scenarios that represent acquisition in marine, transition zone and land environments. Clear differences attributed solely to the parameter design are easily recognizable. Additionally, a comparison among different optimization schemes illustrates the ability of genetic algorithms along with a repeated encoding sequence to find better solutions within a computationally affordable time. The optimized parameters yield a notable enhancement in the deblending and data reconstruction quality and consequently provide optimal acquisition scenarios.     

Evaluation of SEBS for estimation of actual evapotranspiration using ASTER satellite data for irrigation areas of Australia

Spatial knowledge of land surface evapotranspiration (ET) is of prime interest for environmental applications, such as optimizing irrigation water use, irrigation system performance, crop water deficit, drought mitigation strategies, and accurate initialization of climate prediction models especially in arid and semiarid catchments where water shortage is a critical problem. The recent drought in Australia and concerns about climate change have highlighted the need to manage water resources more sustainably especially in the Murrumbidgee catchment which utilizes bulk water for food production. This study deals with the application of a Surface Energy Balance System (SEBS) algorithm based on Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) data and field observations has been proposed and tested for deriving ET over Coleambally Irrigation Area, located in the southwest of NSW, Australia. We have used 12 ASTER scenes covering the time period of 2002, 2003, 2004, 2005, 2006, and 2009 for estimating the actual ET over the study area. To validate the proposed methodology, the ground-measured ET was compared to the ASTER-derived actual ET values for the study area. The derived ET value over the study area is much closer to the field measurement. From the remote sensing results and observations, the root mean square error is 0.89 and the mean absolute percentage difference is 2.87 %, which demonstrate the reasonability of SEBS ET estimation for the study area.     

Spectral similarity between scalars at very low frequencies in the unstable atmospheric surface layer over the Tibetan plateau

Similarity between temperature and water vapour was investigated in the Fourier domain by using their coherency spectra, defined as the correlation coefficient between their Fourier modes, and the relative efficiency of their vertical transport. The class-averaged values of these indices were computed from the turbulence measurements over sparse grasslands on the Tibetan plateau during the intensive observations of GEWEX (Global Energy and Water Experiment) Asian Monsoon Experiment (GAME). It was found that the energy-containing eddies at scales up to 100z (z being height) are characterised by well-maintained similarity between these scalars. The scalars are highly correlated with each other, and their transport efficiencies are almost equal within this scale range. In contrast, similarity was not always maintained at scales larger than 1000z. Detailed analyses showed that this breakdown of similarity occurs occasionally or sporadically, suggesting that it is caused by events whose average return period is not much smaller than the current averaging time, i.e. 30 min. We speculate that entrainment of drier and warmer air at the top of the atmospheric boundary layer caused the scalar dissimilarity at this scale range.     

Crustal thermal regime inferred from magnetic anomaly data and its relationship to seismogenic layer thickness: The Japanese islands case study

The seismogenic layer thickness correlates with surface heat flow beneath the Japanese islands. However, this correlation is shown at restricted area, where seismic activity is high. In order to overcome this spatial limitation, we used another approach to estimate the regional thermal structure in the crust beneath the Japanese islands with more uniform coverage. The bottom depths of the magnetized crust determined from the spectral analysis of residual magnetic anomalies is generally interpreted as the level of the Curie point isotherm. We applied this method to estimate the crustal thermal structure in square windows of 2.125° × 2.125°. The obtained depths ranging from 11 to 30 km with average value of 18 km, correlate with the seismogenic layer thickness. It suggests that the Curie point depth is a useful indicator of the crustal thermal structure in these regions.