Use of precipitation data for diabatic initialization to improve the tropical analysis of divergence and moisture

Summary First, we review the present status of diabatic initialization used for numerical weather prediction and conclude that the deficiency of diabatic initialization mostly stems from the shortcoming in evaluating diabatic heating rates accurately, particularly the release of latent heat by cumulus convection. This indicates the need to adjust the initial conditions for physical processes, and Krishnamurti and his colleagues introduced in 1984 the concept of physical initialization. Since cumulus convection is most sensitive to input data among many physical processes, the adjustment of atmospheric input data to a prediction model to produce desired initial precipitation rates is referred to as cumulus initialization.In this article we describe a general approach to diabatic initialization with a special emphasis on cumulus initialization. We present the results of forecasting experiments with a version of the NCAR Community Climate Model (CCM) to demonstrate the efficacy of a cumulus initialization procedure to ameliorate the spinup problern of precipitation. Finally, we discuss application of the present methodology of cumulus initialization for a stability-dependent mass-flux cumulus parameterization of CCM2 to pave the way to complete the diabatic normal mode initialization package for CCM2. Note that the present cumulus initialization scheme can be used to assimilate into the atmospheric analysis of the tropics the precipitation rates estimated by satellite radiometric imagery data.With 8 FiguresThe National Center for Atmospheric Research is sponsored by the National Science Foundation     

IR active orientation of OH bending mode in topaz

Infrared (IR) and nearinfrared (NIR) absorption spectra of hydrous and F-rich topazes were measured to assign an OH bending mode of topaz. Three absorption peaks at 1165, 3650, and 4803 cm⁻¹ are assigned to OH related absorption peaks. Since a peak at 4803 cm⁻¹ can be assigned to a combination mode of 1165 and 3650 cm⁻¹, the 1165 cm⁻¹ peak is harmonic with the 3650 cm⁻¹ peak. Polarized IR absorption spectra of (100), (010), and (001) planes of the hydrous topaz were measured to examine IR active orientation of the 1165 cm⁻¹ OH related mode. Three pleochroic distributions of the absorption peak at 1165 cm⁻¹ on (100), (010), and (001) planes indicate an active orientation of the 1165 cm⁻¹ OH related mode. The IR active orientation of the 1165 cm⁻¹ OH related mode in topaz is normal to the OH dipole. The orthogonality and harmonic combination mode indicate that the 1165 cm⁻¹ peak is OH bending mode. The active orientation of OH bending mode is polarized in the plane normal to the OH dipole. The polarization suggests that anisotropic thermal vibration of protons on the hydroxyl is maximum along the IR active orientation. Received: August 16, 1996 / Revised, accepted: April 20, 1997     

Candidate precursors: pulse-like geoelectric signals possibly related to recent seismic activity in Japan

Telemetric network observations of pulse-like geoelectric charge signals using a vertical dipole buried under the ground were undertaken at various observation sites over a wide area in Japan from April 1996. From continuous records of the signals during the six months following that, we attempted to select anomalous signals, possibly related to seismic electric activity. Special attention was paid to the elimination of noise resulting from industrial and meteorological electric activity, comparison with other electromagnetic signals in the VLF band and the relation between the precursor period and the distance from the eventual main shock that occurred in Japan. Four candidate precursor electric signals, which were not contaminated by industrial and meteorological electric activity, were then selected, of which the second appeared before the Akita-ken Nairiku-nanbu earthquake swarm, for which the maximum M = 5.9 occurred on 1996 August 11, and the third and fourth before the Chiba-ken Toho-oki earthquake, M = 6.6, on 1996 September 11. A tentative qualitative model explaining why the candidate precursory signal is related to stress building up before an earthquake is presented in terms of the electrification of gases released from fracturing rocks immediately prior to the main shock.     

A study of microearthquake seismicity and focal mechanisms within the Sea of Marmara (NW Turkey) using ocean bottom seismometers (OBSs)

We have carried out seismological observations within the Sea of Marmara (NW Turkey) in order to investigate the seismicity induced after Gölcük–İzmit (Kocaeli) earthquake (M_w 7.4) of August 17, 1999, using ocean bottom seismometers (OBSs). High-resolution hypocenters and focal mechanisms of microearthquakes have been investigated during this Marmara Sea OBS project involving deployment of 10 OBSs within the Çınarcık (eastern Marmara Sea) and Central-Tekirdağ (western Marmara Sea) basins during April–July 2000. Little was known about microearthquake activity and their source mechanisms in the Marmara Sea. We have detected numerous microearthquakes within the main basins of the Sea of Marmara along the imaged strands of the North Anatolian Fault (NAF). We obtained more than 350 well-constrained hypocenters and nine composite focal mechanisms during 70 days of observation. Microseismicity mainly occurred along the Main Marmara Fault (MMF) in the Marmara Sea. There are a few events along the Southern Shelf. Seismic activity along the Main Marmara Fault is quite high, and focal depth distribution was shallower than 20 km along the western part of this fault, and shallower than 15 km along its eastern part. From high-resolution relative relocation studies of some of the microearthquake clusters, we suggest that the western Main Marmara Fault is subvertical and the eastern Main Marmara Fault dips to south at 45°. Composite focal mechanisms show a strike-slip regime on the western Main Marmara Fault and complex faulting (strike-slip and normal faulting) on the eastern Main Marmara Fault.     

U-Pb and Pb-Pb ages of zircons from basaltic eucrites: Implications for early basaltic volcanism on the eucrite parent body

We have undertaken petrologic and SHRIMP U-Th-Pb isotopic studies on zircons from basaltic eucrites (Yamato [Y]-75011, Y-792510, Asuka [A]-881388, A-881467 and Padvarninkai) with different thermal and shock histories. Eucritic zircons are associated with ilmenite in most cases and have subhedral shapes in unmetamorphosed and metamorphosed eucrites. Some zircons in highly metamorphosed eucrites with granulitic texture occur alone in pyroxene, and typically have rounded to subrounded shapes due to recrystallization. Superchondritic Zr/Hf ratios of eucritic zircons indicate that they crystallized from incompatible element-rich melts after crystallization of ilmenite. Concentrations of uranium and thorium in zircons in the unmetamorphosed eucrite Y-75011 are higher than those in metamorphosed eucrites.The U-Pb systems of eucritic zircons are almost concordant but some zircon grains show reverse discordance. Radiogenic lead-loss up to 48% from zircons is observed in the shock-melted eucrite Padvarninkai. The ²⁰⁷Pb-²⁰⁶Pb ages of zircon in Y-75011 (4550 ± 9 Ma, n = 5) are nearly identical, within analytical uncertainty, to the ages of zircons from the metamorphosed eucrite Y-792510 (4545 ± 15 Ma, n = 13), the highly metamorphosed eucrites A-881388 (4555 ± 54 Ma, n = 5) and A-881467 (4558 ± 13 Ma, n = 8), and the shock-melted eucrite Padvarninkai (4555 ± 13 Ma, n = 18). The averaged ²⁰⁷Pb-²⁰⁶Pb age of zircon from five eucrites analyzed in this study is 4554 ± 7 Ma (95% confidence limits, n = 49), indistinguishable from the averaged U-Pb age (4552 ± 9 Ma) of the same samples. Because of the high closure temperature of lead in zircon (T_closure = ∼1050°C with a cooling rate of 0.2°C/yr), the ²⁰⁷Pb-²⁰⁶Pb ages of eucritic zircon do not represent metamorphic ages but crystallization ages of extrusive lavas.This fact strongly suggests that volcanism of the eucrite parent body occurred at a very early stage of the Solar System history, 7-20 Ma after CAI formation (4567.2 ± 0.6 Ma), thus basaltic eucrites crystallized from parental magmas within a short interval following the differentiation of their parent body. The U-Pb ages of eucritic zircons are older than the U-Pb, Sm-Nd and Rb-Sr ages of some basaltic eucrites, which is consistent with differences in closure temperatures of each isotopic system, and suggests that thermal and shock metamorphism affected the isotopic systems of pyroxene, plagioclase and phosphates.     

Deducing the spatial variability of exchange within a longitudinal channel water balance

Developing an appropriate data collection scheme to infer stream–subsurface interactions is not trivial due to the spatial and temporal variability of exchange flowpaths. Within the context of a case study, this paper presents the results from a number of common data collection techniques ranging from point to reach scales used in combination to better understand the spatial complexity of subsurface exchanges, infer the hydrologic conditions where individual influences of hyporheic and groundwater exchange components on stream water can be characterized, and determine where gaps in information arise. We start with a tracer‐based, longitudinal channel water balance to quantify hydrologic gains and losses at a sub‐reach scale nested within two consecutive reaches. Next, we look at groundwater and stream water surface levels, shallow streambed vertical head gradients, streambed and aquifer hydraulic conductivities, water chemistry, and vertical flux rates estimated from streambed temperatures to provide more spatially explicit information. As a result, a clearer spatial understanding of gains and losses was provided, but some limitations in interpreting results were identified even when combining information collected over various scales. Due to spatial variability of exchanges and areas of mixing, each technique frequently captured a combination of groundwater and hyporheic exchange components. Ultimately, this study provides information regarding technique selection, emphasizes that care must be taken when interpreting results, and identifies the need to apply or develop more advanced methods for understanding subsurface exchanges. Copyright © 2013 John Wiley & Sons, Ltd.     

Anisotropic phase transition of rutile under shock compression

Shock recovery experiments for single crystal and powdered specimens of TiO₂ with the rutile structure were performed in the pressure range up to 72 GPa. Single crystal specimens were shocked parallel to [100], [110] and [001] directions. X-ray powder diffraction analysis showed that the amount of -PbO₂ type TiO₂ produced by shock-loading depended strongly on the shock propagation direction. The maximum yield (about 70%) was observed for shock loading to 36 GPa parallel to the [100] direction. In the [001] shock direction, the yield is much smaller than that of the [100] direction. This anisotropic yield was consistent with the observed anisotropy of the phase transition pressure in shock compression measurements. However, transformation to the -PbO₂ type cannot explain the large volume change observed above about 20 GPa. On the basis of the high pressure behavior of MnF₂, we assumed that the high pressure phase was either fluorite or distorted fluorite type and that the phase conversion to the -PbO₂ type was induced spontaneously in the pressure reduction process.We present a displacive mechanism of phase transition under shock compression from the rutile structure to the fluorite structure, in which the rutile [100] is shown to correspond to the fluorite [001] or [110] and the rutile [001] to the fluorite [110]. Direct evidence is obtained by examining the [100] shocked specimen by high resolution electron microscopy.     

Reactions to define the biotite isograd in the Ryoke metamorphic belt,Kii Peninsula,Japan

Two types of biotite isograd are defined in the low-grade metamorphism of the Wazuka area, a Ryoke metamorphic terrain in the Kii Peninsula, Japan. The first, BI₁, is defined by the reaction of chlorite+K-feldspar= biotite+muscovite+quartz+H₂O that took place in psammitic rocks, and the second, BI₂, by the continuous reaction between muscovite, chlorite, biotite and quartz in pelitic rocks. The Fe/Mg ratios of the host rocks do not significantly affect the reactions. From the paragenesis of pelitic and psammitic metamorphic rocks, the following mineral zones were established for this low-pressure regional metamorphic terrain: chlorite, transitional, chlorite-biotite, biotite, and sillimanite. The celadonite content of muscovite solid solution in pelitic rocks decreases systematically with the grade of metamorphism from 38% in the chlorite zone to 11% in the biotite zone. Low pressure does not prohibit muscovite from showing the progressive change of composition, if only rocks with appropriate paragenesis are chosen. A qualitative phase diagram of the AKF system relevant to biotite formation suggests that the higher the pressure of metamorphism, the higher the celadonite content of muscovite at BI₁, which is confirmed by comparing the muscovites from the Barrovian and Ryoke metamorphism.     

The influence of weathering on the geotechnical properties and slope angles of mudstone in the Mineoka earth-slide area,Japan

Three active earth-slide slopes of Tertiary mudstone were investigated: Slope 1 has an angle of 17 ?4°, Slope 2 of 12.9°, and Slope 3 of 11.6°. Infinite slope analysis indicated that the instability of these three slopes can be well explained by using the residual strength parameters of earth-slide soils near the sliding surface in conjunction with the highest ground water table. The residual angle of shearing resistance, ø'_r, plays an important role in the determination of slope angle because it differs greatly among slopes. Mineralogical studies and X-ray diffraction analyses were performed for the clay minerals included in the slope material. The results showed that illite and chlorite were found in Slope 1, and also in Slopes 2 and 3 together with interstratified illite/montmorillonite and montmorillonite. The degree of weathering is progressive in order of Slopes 1, 2 and 3. The alteration of clay minerals by weathering causes the reduction in ø'_r -values, i.e., 19.4° in Slope 1 (steep, less weathered) and 12.1-9.2° in Slopes 2 and 3 (gentle, much more weathered). This result indicates that the degree of weathering has a great influence on the value of ø'_r, which in turn determines the slope angle.     

Radiation induced lattice defects in natural zircon (ZrSiO4) observed at atomic resolution

Fission tracks and point defects in natural zircon are directly observed by a 1 MV electron microscope at atomic resolution for three types of samples adjusted to the 100 orientation. Lattice planes intersecting the fission tracks at high angles are distorted in a rather irregular manner over a wide region up to more than 100 Å wide. Diameter of the tracks, ranging from 25 Å to 40Å, is much narrower than those so far reported for the U-doped synthetic zircon (100–200 Å), UO₂ thin film (100 Å), mica (66 Å, 240 Å) or fluorophlogopite (150 Å). The fact that fairly long tracks thousands of angstroms in length are observed in thin 100-oriented sample hundreds of angstroms in thickness and that some of them are nearly parallel to a low index lattice plane suggest a possible occurrence of channelling in the process of track formation. Parallel tracks often observed in chemically etched specimens support the idea of channelling. Slightly bent tracks are sometimes observed. It is concluded from computer simulation that many contrast anomalies of bright and dark spots in the lattice image are due to point defects of vacancies and interstitial atoms, mainly produced by the direct atomic collision with α-particles or by passage of ionizing nuclear particles. Optimum conditions of the observation of point defects with highest contrast are studied. One interstitial Zr atom or one Zr ion vacancy will give very low contrast and will be not detectable unless the crystal is less than two unit cells thick. A pair of Zr ion vacancies, however, yields extended detectable limit of thickness. Some of the observed defects are in good accordance with those simulated.