Aerodynamic roughness over an inhomogeneous ground surface

The aerodynamic roughness parameter z ₀ over inhomogeneous ground surfaces, such as cities, rural towns and so on, is determined by analyzing the wind data at AMeDAS observatories in the Tohoku and Kanto districts of Japan, by making use of Rossby number similarity theory. It is found that the aerodynamic roughness parameter is proportional to the average size of the roughness elements.A practical method of estimating the aerodynamic roughness parameter over an extensive area with various inhomogeneities is developed. In this method, the Digital National Land Information data bank is employed. As an example, the roughness parameter distribution around Tsukuba Academic City is presented.     

The effect of the Mt St Helens eruption on tropospheric and stratospheric ions

Balloon-borne observations of electrical conductivity in the troposphere and stratosphere were performed using conductivity sondes at Garmisch-Partenkirchen, West Germany, from June to November, 1980, after the Mt St Helens eruption. A significant decrease of atmospheric ions in the altitudes from the troposphere to lower stratosphere has been detected until several months after the eruption in comparison with the observational results obtained before the eruption. Simulteneous ruby lidar observation a month after the eruption indicates an increased amount of aerosol at nearly the same altitude as that of conductivity decrease. Several months after the eruption it appears that aerosols detected by lidar and those effective in reducing ion concentration have different profiles.     

Heterogeneous shock-induced thermal radiation in minerals

A 500 channel optical imaging intensifying and spectral digital recording system is used to record the shock-induced radiation emitted from 406 to 821 nm from transparent minerals during the time interval that a shock wave propagates through the sample. Initial results obtained for single crystals of gypsum, calcite and halite in the 30 to 40 GPa (300 to 400 kbar) pressure range demonstrate greybody emission spectra corresponding to temperatures in the 3,000 to 4,000 K range and emissivities ranging from 0.003 to 0.02. In the case of gypsum and calcite, distinctive line spectra, are superimposed on the thermal radiation. The observed color temperatures are a factor of 2 to 10 greater than the Hugoniot temperature, calculable on the basis of continuum thermodynamics and equation of state models for the shock states achieved in the three minerals. These observed high temperatures are believed to be real. We conclude that we are detecting a large number of closed spaced high temperature shear-band regions immediately behind the shock front. A shear instability model, such as proposed independently by Grady (1977, 1980), Ananin et al. (1974), and Horie (1980), in which small zones of highly deforming and melted material are produced which are the source of the detected thermal radiation and have a fractional effective area approximately numerically equal to the measured emissivity, can be used to predict an effective emissivity of 0.0065 directly behind the shock front. If shear band instability arises from viscous flow processes, Grady's model and mineral thermal properties yield for the shocked mineral viscosities values in the range of 10⁹ to 10¹⁵ P immediately behind the shock front.     

Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

New insights into the P- and S-wave velocity structure of the D" discontinuity beneath the Cocos plate

Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary.     

A new fault rupture scenario for the 2003 Mw 6.6 Bam earthquake, SE Iran: Insights from the high-resolution QuickBird imagery and field observations

The December 26, 2003 M_w 6.6 Bam earthquake is one of the most disastrous earthquakes in Iran. QuickBird panchromatic and multispectral satellite imagery with 61 cm and 2.4 m ground resolution, respectively provide new insights into the surface rupturing process associated with this earthquake. The results indicate that this earthquake produced a 2–5 km-wide surface rupture zone with a complex geometric pattern. A 10-km-long surface rupture zone developed along the pre-existing Bam fault trace. Two additional surface rupture zones, each 2–5 km long, are oblique to the pre-existing Bam fault in angles of 20–35°. An analysis of geometric and geomorphic features also shows that movement on the Bam fault is mainly right-lateral motion with some compressional component. This interpretation is consistent with field investigations, analysis of aftershocks as well as teleseismic inversion. Therefore, we suggest that the 2003 Bam earthquake occurred on the Bam fault, and that the surface ruptures oblique to the Bam fault are caused by secondary faulting such as synthetic shears (Reidel shears). Our fault model for the Bam earthquake provides a new tectonic scenario for explaining complex surface deformations associated with the Bam earthquake.     

Logic-tree Approach for Probabilistic Tsunami Hazard Analysis and its Applications to the Japanese Coasts

For Probabilistic Tsunami Hazard Analysis (PTHA), we propose a logic-tree approach to construct tsunami hazard curves (relationship between tsunami height and probability of exceedance) and present some examples for Japan for the purpose of quantitative assessments of tsunami risk for important coastal facilities. A hazard curve is obtained by integration over the aleatory uncertainties, and numerous hazard curves are obtained for different branches of logic-tree representing epistemic uncertainty. A PTHA consists of a tsunami source model and coastal tsunami height estimation. We developed the logic-tree models for local tsunami sources around Japan and for distant tsunami sources along the South American subduction zones. Logic-trees were made for tsunami source zones, size and frequency of tsunamigenic earthquakes, fault models, and standard error of estimated tsunami heights. Numerical simulation rather than empirical relation was used for estimating the median tsunami heights. Weights of discrete branches that represent alternative hypotheses and interpretations were determined by the questionnaire survey for tsunami and earthquake experts, whereas those representing the error of estimated value were determined on the basis of historical data. Examples of tsunami hazard curves were illustrated for the coastal sites, and uncertainty in the tsunami hazard was displayed by 5-, 16-, 50-, 84- and 95-percentile and mean hazard curves.     

A Simple Single-Layer Urban Canopy Model For Atmospheric Models: Comparison With Multi-Layer And Slab Models

We developed a simple, single-layer urban canopy model, and comparedit to both multi-layer and slab models. Our single-layer model has thefollowing features: (a) It is a column model of energy and momentumexchange between an urban surface and the atmosphere, (b) it includesthe influence of street canyons, which are parameterized to representthe urban geometry, (c) it includes shadowing from buildings andreflection of radiation, and (d) it estimates both the surfacetemperatures of, and heat fluxes from, three surface types: roof, wall,and road. In the simulation of the single-layer model, the roof washottest during the daytime, but coolest from midnight to early morning.This is consistent with output from the multi-layer model and fieldobservations at a residential area on a clear, summer day. The diurnalvariation of the energy budget from the single-layer model agrees wellwith that from the multi-layer model. Our single-layer model'sperformance is nearly that of a multi-layer model for studyingmesoscale heat islands. Nevertheless, it is simply parameterized,and thus easily included in larger-scale atmospheric models. The slabmodel has the largest nighttime cooling rate of the three models. Toovercome this, it needs more adjustments than for the canopy models.     

The use of sintered diamond anvils in the MA8 type high-pressure apparatus

A new multi-anvil type high-presure apparatus has been developed using sintered diamond anvils to generate pressures over 30 GPa and temperatures up to about 2000°C. A maximum sample volume of about 1 mm³ is available in this system. The pressure was confirmed by dissociation of forsterite into Mg-perovskite and periclase. The basic techniques and problems in utilizing sintered diamond in the MA8 type high-pressure apparatus are discussed with an emphasis on the future prospect of incorporating simultancous X-ray diffraction observation.     

Assessment of the Uncertainties in the NO2 and O3 Measurements by Visible Spectrometers

The column amounts of nitrogen dioxide (NO2) and ozone (O3) were measured using a visible spectrometer based on the twilight zenith-sky technique at two observatories located at similar latitudes in the northern part of Japan separated by a distance of 150 km. The measurements began in April 1991 at the Moshiri Observatory (44.4°N, 142.3°E) and in April 1994 at the Rikubetsu Observatory (43.5°N, 143.8°E). Since weather conditions and the possible influence from tropospheric pollution were not always identical at these two observatories, the overall accuracy of the measurements was studied comparing these data sets. The first year data obtained at a solar zenith angle of 90 degrees indicated that the NO2 slant column values at sunrise and sunset agreed within 0.36 and 0.54 × 1016 cm-2, respectively, corresponding to 5 % (June) and to 12 % (December) of the columns. The O3 values agreed within 0.76 × 1019 cm-2, corresponding to 4 % (March) 6 % (August) of the columns, although a part of the difference was systematic. The O3 column amounts were also compared to those obtained by the Dobson spectrometer at Sapporo (43.5°N, 143.8°E), whose latitude is similar to these observatories. When an air mass factor of 17.5 was used, the two-year Moshiri vertical column values agreed with the Dobson direct sun values to within 15 Dobson Units, or 3 6 % of the column. The difference between the two values was found to be due partly to the change in the air mass factor caused by seasonal and day-to-day changes in the shape of the O3 vertical profiles. These results confirm the reliability of the NO2 and O3 measurements by visible spectrometers at these sites for the Network for the Detection of Stratospheric Change (NDSC).