Fresnel zone inversion for lateral heterogeneities in the earth

We propose a different kind of seismic inversion from travel-time or waveform inversion for lateral heterogeneities in the earth: Fresnel zone inversion. Amplitude and phase delay of data in several frequency ranges are inverted for model space around ray paths with a width corresponding to the considered frequency so that primary effect of finiteness of wavelength be included. For vertically heterogeneous media, Fréchet derivatives for inversion are obtained very efficiently using the paraxial ray approximation, with nearly similar amounts of computation compared to travel-time inversion. As an example, Fréchet derivatives are computed for a teleseismic observation system for a three-dimensional structure in the lithosphere beneath an array of seismic stations. Even if the used frequency is around 2 Hz, the width of Fréchet derivatives cannot be neglected, particularly near the bottom of the lithosphere. Sensitivity of model parameters to observations is, moreover, different in our approach from conventional travel-time inversion: it is zero along ray paths but large slightly away from them. Some model calculations show that travel-time inversion, particularly with models divided into very fine meshes or blocks, might give misleading results. An example of inversion for a simple Camembert model, in the event that travel-time inversion gives no reliable results, shows how this technique works with much smaller data sets and computation than waveform inversions.     

Strontium concentrations and strontium-chlorinity ratios in sea water of the North Pacific and the adjacent seas of Japan

Strontium contents of 232 sea water samples collected at various stations in the North Pacific and adjacent seas of Japan were measured by the atomic absorption spectrophotometry and strontium-chlorinity ratios were determined. Mean Sr concentration is 8.08 mg/kg and mean Sr/Cl ratio is 0.425 mg/kg/Contrary to some recent reports, regional and vertical variations of Sr/Cl ratios were statistically insignificant, and presence of particulate strontium was not confirmed.     

Internal characteristics and numerical analysis of plunging breaker on a slope

Many investigations about the direct measurements of velocities to clarify the internal mechanism of the breaker have been carried out as a result of recent progress in the measuring techniques.This research attempts to clarify the breaking wave transformation system on a slope by an experiment and numerical analysis. In an experiment, the velocities in the surf zone were measured directly using an electromagnetic current meter, and the space distribution characteristic of the vorticity ω = (∂u/∂y − ∂u/∂x) and the skewness γ = (∂u/∂y + ∂u/∂x) were examined. Also, occurrence situations of the vortices at the time of water mass inrush were measured by video tape recorder (VTR) image processing. However, because the breaker is a violent phenomenon that is entrained with plentiful bubbles, the extent to which we can clarify breaker transformation in experiments is limited. Numerical simulations are substituted for experiments as a method to clarify breaker transformation.In numerical analysis, finite amplitude wave analysis based on the potential theory (non-viscous fluid) is possible before wave breaking; however, the analysis must take into account the viscous fluid after breaking. So, we use the Reynolds equations to develop a numerical simulation system of the breaker transformation on a sloping bottom. The numerical energy dissipation model of the breaker was compared to the experimental results, and a modified Simplified Marker and Cell (SMAC) method is presented. The internal characteristics of the breaker transformation are described using application examples.     

Deep convection in a lake triggered by wind: Two-dimensional numerical experiments with a nonhydrostatic model

We have investigated the fundamental processes of deep convection in a lake at high latitudes triggered by wind during spring or autumn and the associated deep water formation, executing vertically two-dimensional numerical experiments with a nonhydrostatic model. The water column in which a relatively cold mixed layer overlies a relatively warm layer becomes unstable, when the Ekman convergence on the shore due to along-shore wind deepens the mixed layer below the compensation depth, where water densities in both layers becomes equal to each other because of the thermobaric effect. At the onset of deep convection, the critical Rayleigh number agrees with that predicted by the linear theory. The onset time of deep convection is inversely proportional to the magnitude of wind stress. On the other hand, the onset time is minimal when water temperature in the mixed layer _m is 3.1°C because a change of _m has two effects oppositely acting on the stability of the water column. After the first onset, deep convection occurs intermittently for a few days. The sinking of the mixed layer water occurs in a thermal-like shape, and its amount is 4184% of the time-integrated Ekman transport when _m 3°C while it decreases to less than 10% for _m lower than 1.5°C. The present process can explain 30% of the amount of deep water renewal which is expected from the observation in Lake Baikal.     

239, 240Pu,137Cs and90Sr in the central North Pacific

^{239, 240}Pu,¹³⁷Cs and⁹⁰Sr concentrations were determined in sea waters from the central and western North Pacific in 1980 and 1982. The results are consistent with those reported earlier for North Pacific waters. The profiles of⁹⁰Sr and¹³⁷Cs show a monotonic decrease with depth, whilst^{239, 240}Pu shows a distinct subsurface maximum at a depth between 400 and 1,000 m. The calculated inventories of these nuclides significantly exceed the global mean fallout inputs for these latitudes. This may be due to local fallout input to the ocean at times of large-scale nuclear weapon tests in the equatorial North Pacific. The existence of measurable amounts of¹³⁷Cs and^{239, 240}Pu in deep waters suggests that these nuclides are transported by sinking particulate matter from the surface to the deep ocean.     

Prediction of subsurface water level change from satellite data

This study explores the potential for predicting the spatial variation in subsurface water level change with crop growth stage from satellite data in Thabua Irrigation Project, situated in the northern central region of Thailand. The relationship between subsurface water level change from pumping water to irrigate rice in the dry season and the age of the rice was analysed. The spatial model of subsurface water level change was developed from the classification using greenness or (normalized difference vegetation index NDVI) derived from Landsat 5 Thematic Mapper data. The NDVI of 52 rice fields was employed to assess its relationship to the age of the rice. It was found that NDVI and rice age have a good correlation (R² = 0·73). The low NDVI values (−0·059 to 0·082) in these fields were related to the young rice stage (0–30 days). NDVI and subsurface water level change were also correlated in this study and found to have a high correlation (Water level change (m day⁻¹) = 0·3442 × NDVI − 0·0372; R² = 0·96). From this model, the water level change caused by rice at different growth stages was derived. This was used to show the spatial variation of water level change in the project during the 1998–99 dry‐season cropping. This simple method of using NDVI relationships with water level change and crop growth stages proves to be useful in determining the areas prone to excessive lowering of the subsurface water level during the dry season. This could assist in the appropriate planning of the use of subsurface water resources in dry‐season cropping. Copyright © 2004 John Wiley & Sons, Ltd.     

Ridge collision and in situ greenstones in accretionary complexes: An example from the Late Cretaceous Ryukyu Islands and southwest Japan margin

Abstract Mélange units containing greenstones are common throughout the Cretaceous-Miocene Shimanto Supergroup in the Ryukyu Is and southwest Japan. Most greenstones in the accretionary complex originated in oceanic spreading ridges and seamounts, and they formed far from the convergent margin. Some mélange-like units in the supergroup, however, contain greenstones that were extruded upon and intruded into unconsolidated fine-grained terrigenous clastic sediments. It is inferred that eruption of the in situ greenstones resulted from igneous activity in the trench area. Geochemical signatures indicate that the greenstone protoliths were similar to mafic lavas generated at spreading ridges. Fossil ages of the strata containing in situ greenstones become younger over a distance of 1300 km eastward from Amami-Oshima (Cenomanian-Turonian) in the Ryukyu Is to central Japan (Late Maestrichtian-earliest Paleocene), implying that a site of igneous activity in the trench area migrated eastward along the Ryukyu Is and southwest Japan margin. Plate reconstructions of the northwest Pacific Ocean suggest the presence of the Kula-Pacific ridge near Late Cretaceous to early Paleogene Japan. In this context, it is suggested that the greenstones formed in response to Kula-Pacific ridge-forearc collision.
Ancient ridge-forearc collisions are best recognized by the presence of mid-ocean ridge basalt (MORB) extruded on sediments inferred to have accumulated in the trench area. Diachronous occurrences of the strata associated with these MORB in an orogenic belt are useful for documenting the ridge collision through time.     

Urban albedo as a function of the urban structure — A two-dimensional numerical simulation

Urban albedo change as a function of urban geometrical structure has been examined by using a two-dimensional urban block-canyon array model. The complex multiple reflections of incident photons in the urban canyon are simulated by using a Monte-Carlo method. The photons are tracked until they leave the canyon or are completely absorbed. In the model, the direct and diffuse components of incident solar radiation are introduced and the specular and isotropic reflection characteristics are considered for the relevant urban surfaces. The result shows that the urban albedo decreases as the urban irregularity increases as indicated by the model experiment of Aida (1982). The dependence of albedo on the incident solar zenith angle observed in the experiment is also confirmed for various urban models.As an application, some actual urban structures in the Marunouchi area in Tokyo are examined. Urban planning for absorption and reflection of solar radiation in urban areas is also discussed based on the analysis of the change in albedo with canyon dimensions and solar zenith angle.     

A new fracture criterion for brittle rocks

A general form of a “fracture function” for isotropic brittle materials is expressed in terms of the three invariants of the stress tensor. The coefficients in the function are determined by use of the small number of experimental data under specific conditions. This function is applicable to an estimate of the fracture condition of brittle rocks under a general stress state i.e., σ₁ ≠ σ₂ ≠ σ₃. The application of this function is attempted for the data of three brittle rocks i.e., Dunham dolomite, Mizuho trachyte, and Westerly granite, reported by previous workers. For the first two, this criterion gives a good estimation of the effect of the intermediate principal stress σ₂ on failure. For the last, the fracture strength at high confining pressure is estimated by use of the several data obtained under very low confining pressures, and the agreement with experimental data is also satisfactory.