Simulating present climate of the global ocean–ice system using the Meteorological Research Institute Community Ocean Model (MRI.COM): simulation characteristics and variability in the Pacific sector

A long-term spin-up and a subsequent interannual simulation are conducted for the ocean–ice component of the climate model intercomparison project (CMIP)-class earth system model of the Japan Meteorological Agency/Meteorological Research Institute. This experiment has three purposes: first is to assess the ability of our model with the Coordinated Ocean–ice Reference Experiments (COREs) forcing in reproducing the present ocean-climate; second is to understand the ocean-climate variability for the past 60 years; third is to present an example of evaluating an ocean–ice interannual variability simulation. The Pacific Ocean is focused on for the last two purposes. After integrating for about 1500 years with repeated use of a detrended CORE interannual forcing, the model reaches a quasi-steady state where the present climate is reproduced satisfactorily. Then, the interannual variability simulation is conducted with the retrieved forcing trend and the result is analyzed. The simulation is successful at reproducing the long-term variability in the Pacific and surrounding oceans. Brief analyses of the tropical and mid-latitude upper layer, deep circulation, and the Arctic sea ice are presented. A caveat in treating other parts of the globe is due to the recent intense convection in the Southern Ocean caused by a remarkably increasing trend of the Southern Hemisphere westerly. Overall, the current simulation with our CMIP-class ocean–ice model is shown to be useful for studying the present ocean-climate variability, specifically in the Pacific sector. It could also be used as a benchmark control experiment that facilitates further research, model development, and intercomparison.     

Exploration of S-wave velocity profiles at strong motion stations in Eskisehir,Turkey, using microtremor phase velocity and S-wave amplification

We have explored 1D S-wave velocity profiles of shallow and deep soil layers over a basement at strong motion stations in Eskisehir Province, Turkey. Microtremor array explorations were conducted at eight strong motion stations in the area to know shallow 1D S-wave velocity models. Rayleigh wave phase velocity at a frequency range from 3 to 30 Hz was estimated with the spatial autocorrelation analysis of array records of vertical microtremors at each station. Individual phase velocity was inverted to a shallow S-wave velocity profile. Low-velocity layers were identified at the stations in the basin. Site amplification factors from S-wave parts of earthquake records that had been estimated at the strong motion stations by Yamanaka et al. (2017) were inverted to the S-wave velocities and Q-values of the sedimentary layers. The depths to the basement with an S-wave velocity of 2.2 km/s are about 1 km in the central part of the basin, while the basement becomes shallow as 0.3 km in the marginal part of the basin. We finally discussed the effects of the shallow and deep sedimentary layers on the 1D S-wave amplification characteristics using the revealed profiles. It is found that the shallow soil layers have no significant effects in the amplification at a frequency range lower than 3 Hz in the area.     

The three-dimensional shear wave structure in the mantle by overtone waveform inversion - I. Radial seismogram inversion

Summary. The three-dimensional (3-D) shear wave structure of the mantle, down to the depth of about 900 km, is obtained by inverting waveforms of radial component seismograms. Radial component seismograms contain large amplitude overtone signals which circle the Earth as wave packets and are sometimes called X1, X2, X3, … We use data which contain R1, X1 and X2 and filtered between 2 and 10mHz. It is shown that, unless each seismogram is weighted, all seismograms are not fitted uniformly. Only data from large earthquakes are fitted and the final velocity anomalies are biased by the small number of large earthquake data. Resolution is good at shallow depths, becomes worse in the intermediate depth range between about 400 and 500 km and then becomes better at greater depth ranges (600–900km). Even though we use only spheroidal mode data, velocity anomalies in the shallow structure show excellent correlation with the age of the surface rocks of the Earth. In the deeper regions, between about 600 and 900km, South America shows a fast velocity anomaly which may indicate the slab penetration beyond 700 km there. Another region which shows a fast velocity anomaly is the Mariana trench, but other subduction regions do not show such features.     

Effects of large-scale wind on the Kuroshio path south of Japan in a 60-year historical OGCM simulation

The effects of large-scale wind forcing on the bimodality of the Kuroshio path south of Japan, the large meander (LM) and non-large meander (NLM), were studied by using a historical simulation (1948–2007) with a high-resolution Ocean general circulation models (OGCM). The Kuroshio in this simulation spent much time in the NLM state, and reproduced several aspects of its long-term path variability for the first time in historical OGCM simulation, presumably because the eddy kinetic energy was kept at a moderate level. By using the simulated fields, the relationships between wind forcing (or Kuroshio transport) and path variation proposed by past studies were tested, and specific roles of eddies in those variations were investigated. The long-term variation of the simulated net Kuroshio transport south of Japan was largely explained by the linear baroclinic Rossby wave adjustment to wind forcing. In the simulated LM events, a triggering meander originated from the interaction of a wind-induced positive sea surface height (SSH) anomaly with the upstream Kuroshio and was enlarged by cyclonic eddies from the recirculation gyre. The cyclonic eddy of the trigger meander was followed by a sizable anticyclonic eddy on the upstream side. Subsequently, a weak (strong) Kuroshio favored the LM (NLM). The LM tended to be maintained when the Kuroshio transport off southern Japan was small, and increasing Kuroshio transport promoted decay of an existing LM. The supply of disturbances from upstream, which is related to the wind-induced SSH variability at low latitudes, contributed to the maintenance of an existing LM.     

ESR analysis of the Nojima fault gouge, Japan, from the DPRI 500 m borehole

Abstract Electron spin resonance (ESR) analyses of quartz grains in fault gouge were performed for a core sample taken from the Nojima Fault that moved during the 1995 Kobe earthquake (Hyogo-ken Nanbu earthquake). Distribution of radiation-induced defects in the gouge at a depth of 389.4 m was measured by extracting quartz grains from seven discrete positions within 30 mm of the fault plane on the granite side. The decrease in E'₁ and Al centers was observed within 2 mm of the fault plane, suggesting partial annealing due to faulting. Partial annealing even at that depth suggested that conventional ESR dating, which is based on the hypothesis of complete annealing during faulting, was not applicable. Theoretical calculations of the temperature rise and of the thermal annealing of defects have been made by assuming a simple annealing model of heat generation on the fault plane. Thermal energy was calculated to have been approximately 8 MJ/m² to explain the profile of the heat-affected region. Thermal energy was much larger than the one estimated from hydrothermal solution, and corresponded to the frictional heat calculated for a normal stress of 20 MPa, a displacement of 2 m, and a frictional coefficient of 0.2.     

Sequential high-pressure transformations of FeGeO3 high-P clinopyroxene (C 2 /c) at temperatures up to 365 ¯C

 In order to elucidate high-pressure transformations of high-P clinopyroxene (C2/c) at kinetically low temperature where atoms are not thermally activated, the transformation processes of FeGeO₃ clinopyroxene (C2/c) have been investigated at pressures up to 20 GPa and 365 °C by powder X-ray diffraction using a synchrotron radiation source and TEM observation. With increasing pressure up to 20 GPa at room temperature, FeGeO₃ high-P clinopyroxene (C2/c) reversibly transforms into a new high-pressure phase, FeGeO₃(II). On increasing the temperature up to 365 °C, this phase rapidly transforms into FeGeO₃ ilmenite within about 2 h. Intensity analysis of the X-ray diffraction pattern reveals that the high-pressure phase of FeGeO₃(II) has an intermediate structure between clinopyroxene and ilmenite: the cation arrangement is similar to that of clinopyroxene and the oxygen arrangement is similar to that of ilmenite. The comparison of the crystal structures of these polymorphs suggests that clinopyroxene to FeGeO₃(II) and FeGeO₃(II) to ilmenite transformations are performed by the slight deformation of the oxygen packing and the short-range movement of cations, respectively. It is shown that this high-P clinopyroxene transforms into ilmenite through a low-activation energy path under the low-temperature condition. Received: 30 August 2000 / Accepted: 10 February 2001