The effects of the spinel—garnet phase transition on the formation of rifted sedimentary basins

We have examined the effects of the spinel-garnet phase transition on subsidence of extensional sedimentary basins. For a constant positive Clapeyron slope ( dP/dT ), the phase boundary moves downwards in the syn-rift and upwards in the post-rift phase. For a non-linear Clapeyron curve ( dP/dT > 0 above 900°C and dP/dT < 0 below 900°C), theory predicts for the reaction of the spinel-garnet phase transition, the direction of phase boundary movement is dependent on the stretching factor, the position of the Clapeyron curve and the lithospheric thickness. A smaller syn-rift and larger post-rift subsidence are predicted for a deeper phase boundary and a thicker lithosphere. The model with a non-linear Clapeyron curve is applied to the subsidence histories of a young extensional basin (Gulf of Lion) and an old continental margin (eastern Canada). The observed syn-rift uplift and the larger post-rift subsidence can be reasonably explained by this model, where the optimum depth of the phase boundary for eastern Canada (˜90 km) is consistent with the estimate from seismic observations and is larger than that for the Gulf of Lion (˜ 50 km). The depth of the spinel-garnet phase boundary is sensitive to the composition of mantle rocks and increases with the extraction of basaltic components from the lithosphere, compatible with our result that the phase boundary is deeper for an older and thicker lithosphere. Thus the surface movement associated with the rifting for these areas may reflect the chemical evolution of the continental lithosphere.     

K–Ar geochronology of the temporal change of eruptive style in the eastern Izu peninsula, central Japan

Abstract   A recent K–Ar study elucidated that eruptive style in the eastern Izu peninsula changed from polygenetic to monogenetic volcano at 0.3–0.2 Ma. To narrow down the time of change, we determined 10 K–Ar ages on Togasayama Andesite of Amagi volcano, the youngest polygenetic volcano in the area, and Togasayama Monogenetic Volcano, one of the oldest monogenetic volcanoes in the area, which overlies a part of the Togasayama Andesite. Dating results showed that the Togasayama Andesite effused at least from 0.34 to 0.20 Ma, whereas the Togasayama Monogenetic Volcano erupted at 0.26–0.29 Ma, suggesting that the northern part of the Togasayama Andesite effused after the eruption of the Togasayama Monogenetic Volcano. Considering previous data, it is therefore inferred that change of eruptive style in the eastern Izu area occurred during the period 0.29–0.20 Ma, with considerable overlap of both polygenetic and monogenetic volcanism.     

Measures for groundwater security during and after the Hanshin-Awaji earthquake (1995) and the Great East Japan earthquake (2011), Japan

Many big earthquakes have occurred in the tectonic regions of the world, especially in Japan. Earthquakes often cause damage to crucial life services such as water, gas and electricity supply systems and even the sewage system in urban and rural areas. The most severe problem for people affected by earthquakes is access to water for their drinking/cooking and toilet flushing. Securing safe water for daily life in an earthquake emergency requires the establishment of countermeasures, especially in a mega city like Tokyo. This paper described some examples of groundwater use in earthquake emergencies, with reference to reports, books and newspapers published in Japan. The consensus is that groundwater, as a source of water, plays a major role in earthquake emergencies, especially where the accessibility of wells coincides with the emergency need. It is also important to introduce a registration system for citizen-owned and company wells that can form the basis of a cooperative during a disaster; such a registration system was implemented by many Japanese local governments after the Hanshin-Awaji Earthquake in 1995 and the Great East Japan Earthquake in 2011, and is one of the most effective countermeasures for groundwater use in an earthquake emergency. Emphasis is also placed the importance of establishing of a continuous monitoring system of groundwater conditions for both quantity and quality during non-emergency periods.     

Pressure and temperature dependence of cation distribution in Mg-Mn olivine

Synthetic (Mg_0.51, Mn_0.49)₂SiO₄ olivine samples are heat-treated at three different pressures; 0, 8 and 12 GPa, all at the same temperature (～500° C). X-ray structure analyses on these single crystals are made in order to see the pressure effect on cation distribution. The intersite distribution coefficient of Mg and Mn in M1 and M2 sites, K _D = (Mn/Mg) _M1/(Mn/Mg) _M2, of these samples are 0.192 (0 GPa), 0.246 (8 GPa) and 0.281 (12 GPa), indicating cationic disordering with pressure. The small differences of cell dimensions between these samples are determined by powder X-ray diffraction. Cell dimensions b and c decrease, whereas a increases with pressure of equilibration. Cell volume decreases with pressure as a result of a large contraction of the b cell dimension. The effect of pressure on the free energy of the cation exchange reaction is evaluated by the observed relation between the cell volume and the site occupancy numbers. The magnitude of the pressure effect on cation distribution is only a fifth of that predicted from the observed change in volume combined with thermodynamic theory. This phenomenon is attributed to nonideality in this solid solution, and nonideal parameters are required to describe cation distribution determined in the present and previous experiments. We use a five-parameter equation to specify the cationic equilibrium on the basic of thermodynamic theory. It includes one energy parameter of ideal mixing, two parameters for nonideal effects, one volume parameter, and one thermal parameter originated from the lattice vibrational energy. The present data combined with some of the existing data are used to determine the five parameters, and the cation distribution in Mg-Mn olivine is described as a function of temperature, pressure, and composition. The basic framework of describing the cationic behavior in olivine-type mineral is worked out, although the result is preliminary: each of the determined parameters is not accurate enough to enable us to make a reliable prediction.     

Pressure effect on the divalent cation distribution in nonideal solid solution of forsterite and fayalite

The divalent cation distribution in olivine (Mg, Fe)₂SiO₄ under high pressure and temperature was studied to clarify the detailed state of olivine in the mantle. Single crystal samples were heated for a sufficient period of time for the cations to migrate and quenched fast enough to preserve the equilibrated state under high pressures, and the crystal structure was determined with X-ray method. The pressure effect on the distribution coefficient K _D[= (Fe/Mg) _M1/(Fe/Mg) _M2] was determined for the first time; dK _D/dP?0.02 GPa^?1. A set of five thermodynamic parameters required to describe the regular solution model was determined from data concerning the pressure dependence and the known temperature and compositional effects. As a result we have shown how K _D depends on pressure, temperature, and composition. The notable feature clarified is the very large contribution of nonideality in the olivine solid solution. The K _D of olivine crystals in the mantle is predicted; K _D increases to ～ 2.2 at the depth of 400 km, in contrast to 0.9 ～ 1.2 of natural samples available at the surface of the Earth.     

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.     

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.     

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.     

Tectonic implications of damage zone-related fault-fracture networks revealed in drill core through the Nojima fault, Japan

1800 m of drill core through the Nojima fault zone, Japan, reveals subsidiary fault and fracture networks that developed in the fault zone that triggered the 1995 Ms 7.2 Kobe earthquake. The subsidiary fault zones contain a fault gouge of < 1 cm bounded by thin zones of foliated cataclasite or breccia. Fractures are filled with calcite veins, calcite-cemented breccias, clay, and iron-oxide and carbonate alternation of the granitic host rock. These features are typical of extensional fractures that form the conduit network for fluid flux close to a major fault zone. The zone of distributed deformation surrounding the main fault is 50 m in width, and the dip of the Nojima fault at > 1 km depth is 75°. The fault-fracture networks associated with the Nojima fault zone are coseismic and were filled with carbonate and fine-grained material during repeated seismic-related infiltration of the fault zone by carbonate-bearing subsurface water. This study shows that fault-related fracture networks plays an important role as fluid flow conduits within seismically active faults, and can change in character from zones of high permeability to low permeability due to cementation and/or pore collapse.