Seismic vertical array analysis for phase decomposition

We propose a vertical array analysis method that decomposes complex seismograms into body and surface wave time histories by using a velocity structure at the vertical array site. We assume that the vertical array records are the sum of vertically incident plane P and S waves, and laterally incident Love and Rayleigh waves. Each phase at the surface is related to that at a certain depth by the transfer function in the frequency domain; the transfer function is obtained by Haskell's matrix method, assuming a 1-D velocity structure. Decomposed P , S and surface waves at the surface are estimated from the vertical array records and the transfer functions by using a least-squares method in the frequency domain; their time histories are obtained by the inverse Fourier transform. We carried out numerical tests of this method based on synthetic vertical array records consisting of vertically incident plane P and S waves and laterally incident plane Love and Rayleigh waves. Perfect results of the decomposed P , S , Love and Rayleigh waves were obtained for synthetic records without noise. A test of the synthetic records in which a small amount of white noise was added yielded a reasonable result for the decomposed P , S and surface waves. We applied this method to real vertical array records from the Ashigara valley, a moderate-sized sedimentary valley. The array records from two earthquakes occurring at depths of 123 and 148 km near the array (epicentral distance of about 31 km) exhibited long-duration later phases. The analysis showed that duration of the decomposed S waves was a few seconds and that the decomposed surface waves appeared a few seconds after the direct S -wave arrival and had very long duration. This result indicated that the long-duration later phases were generated not by multireflected S waves, but by basin-induced surface waves.     

Linking the northern hemisphere sea-ice reduction trend and the quasi-decadal arctic sea-ice oscillation

The nature of the reduction trend and quasi-decadal oscillation in Northern Hemisphere sea-ice extent is investigated. The trend and oscillation that seem to be two separate phenomena have been found in data. This study examines a hypothesis that the Arctic sea-ice reduction trend in the last three decades amplified the quasi-decadal Arctic sea-ice oscillation (ASIO) due to a positive ice/ocean-albedo feedback, based on data analysis and a conceptual model proposed by Ikeda et al. The theoretical, conceptual model predicts that the quasi-decadal oscillation is amplified by the thinning sea-ice, leading to the ASIO, which is driven by the strong positive feedback between the atmosphere and ice-ocean systems. Such oscillation is predicted to be out-of-phase between the Arctic Basin and the Nordic Seas with a phase difference of 3/4, with the Nordic Seas leading the Arctic. The wavelet analysis of the sea ice data reveals that the quasi-decadal ASIO occurred actively since the 1970s following the trend starting in the 1960s (i.e., as sea-ice became thinner and thinner), as the atmosphere experienced quasi-decadal oscillations during the last century. The wavelet analysis also confirms the prediction of such out-of-phase feature between these two basins, which varied from 0.62 in 1960 to 0.25 in 1995. Furthermore, a coupled ice-ocean general circulation model (GCM) was used to simulate two scenarios, one without the greenhouse gas warming and the other having realistic atmospheric forcing along with the warming that leads to sea-ice reduction trend. The quasi-decadal ASIO is excited in the latter case compared to the no-warming case. The wavelet analyses of the simulated ice volume were also conducted to derive decadal ASIO and similar phase relationship between the Arctic Ocean and the Nordic Seas. An independent data source was used to confirm such decadal oscillation in the upper layer (or freshwater) thickness, which is consistent with the model simulation. A modified feedback loop for the sea-ice trend and ASIO was proposed based on the previous one by Mysak and Venegas and the ice/albedo and cloud/albedo feedabcks, which are responsible for the sea ice reduction trend.     

Double seismic zone beneath the middle of Hokkaido, Japan, in the southwestern side of the Kurile Arc

The vertical section of microearthquakes, determined accurately by using the Hokkaido University network, shows two dipping zones (the double seismic zone) 25–30 km apart in the depth range of 80–150 km beneath the middle of Hokkaido in the southwestern side of the Kurile arc. Hypocentral distribution of large earthquakes (m_b > 4) based on the ISC (International Seismological Centre) bulletin also shows the double seismic zone beneath the same region. The hypocentral distribution indicates that the frequency of events occurring in the lower zone is four times greater than that in the upper zone. The difference in seismic activity between the two zones beneath Hokkaido is in contrast with the region beneath northeastern Honshu in the northeastern Japan arc.Composite focal mechanisms of microearthquakes and individual mechanisms of large events mainly characterize the down-dip extension for the lower zone as is observed beneath northeastern Honshu. For the upper zone, however, the stress field is rather complex and not necessarily similar to that beneath northeastern Honshu. This may be considered to indicate the influence of slab contortion or transformation in the Hokkaido corner between the Kurile and the northeastern Japan arcs.     

Intra-oceanic island arc origin for Iratsu eclogites of the Sanbagawa belt,central Shikoku,southwest Japan

New geochemical and Sr–Nd isotopic data for the Iratsu eclogite and surrounding metamorphic rocks of the Sanbagawa belt, Japan, show that, while the protoliths of the metamorphic rocks formed in a variety of tectonic settings, the Iratsu body represents a deeply subducted and accreted island arc. The igneous protoliths of eclogites and garnet amphibolites were probably generated from a mantle source that had components of both a depleted mantle modified by slab-released fluid (as seen in a negative Nb anomaly) and an enriched mantle, similar to that of ocean island basalts (OIB). Fractional crystallization modeling indicates that the protoliths of some garnet clinopyroxenites from the Iratsu body are cumulates from a basaltic magma that crystallized under high O₂ and H₂O fugacities in the middle to lower crust. The source characteristics and crystallization conditions suggest that the protoliths of the Iratsu rocks formed in an oceanic island arc. Quartz eclogites from the marginal zone of the Iratsu body have geochemical signatures similar to turbidites from the Izu–Bonin island arc (as seen in a negative Nb anomaly and a concave REE pattern). The protoliths might be volcaniclastic turbidites that formed in a setting proximal to the oceanic island arc. Geochemical and isotopic signatures of the surrounding mafic schists are similar to normal (N-) and enriched (E-) mid-ocean-ridge basalt (MORB), and distinct from the rocks from the Iratsu body. The protoliths of the mafic schists likely formed in a plume-influenced mid-ocean ridge or back-arc basin. Pelitic schists from the surrounding rocks and pelitic gneisses from the marginal zone of the Iratsu body have evolved, continental geochemical signatures (as seen in a negative ε_Nd(t) value (~?5)), consistent with their origin as continent-derived trench-fill turbidites.     

Interannual variability of CFC-11 absorption by the ocean: an offline model study

The global ocean Chlorofluorocarbon (CFC-11) was simulated in an offline model driven by re-analysis ocean currents in order to identify the mechanisms of interannual to interdecadal variability of air?Csea CFC fluxes. The model was forced with the observed anthropogenic perturbations of atmospheric CFC-11 from the post industrial period (1938) following the OCMIP-II flux protocols along with the observed winds from 1960 to 1999 in the formulation of surface gas exchanges. The model ocean CFC-11 inventories, at the end of 1990s, accounted approximately 1% of the total atmospheric CFC-11, which is consistent with the corresponding observations. The mid-to-high latitude oceans were venue for strong (weak) oceanic sinks (sources) of CFC-11 during the winter (summer) months. The Southern Ocean (south of 40°S) and the North Atlantic (north of 35°N) provided two largest sinks of CFC-11, through which 31.4 and 14.6% of the global ocean CFC-11 entered, respectively. The eastern tropical Pacific Ocean exhibited large interannual variability of CFC-11 flux with a strong (weak) sink during La Ni?a (El Ni?o) years and represented 36% of the global CFC-11 flux variability. The North Atlantic and Southern Ocean were found as regions of large sink efficiency: a capacity to sink more CFC than outsource, although it reduced by 80 and 70%, respectively, in the last 40?years compared to 1960. The sink to source ratio of global ocean CFC-11 fluxes were reduced from 90 to 50% in the last 40?years. This indicates a saturation of CFC in the above-thermocline subsurface that makes the upper ocean less efficient in absorbing CFC in recent decades. A positive trend in CFC sink is now limited to the Southern Ocean, central tropical Pacific and western boundary current regions which possess active upwelling of old water with long time since last atmospheric contact. However, a globally averaged trend was a reduced CFC-11 sink, by emitting 30% of the total ocean CFC-11 that was absorbed during last 40?years.     

Metamorphic condition of a regional metamorphic complex in the Omuta district in northern Kyushu,southwest Japan

A high‐temperature (T) metamorphic complex occurs in the Omuta district, northern Kyushu, Japan. Three metamorphic zones are defined based on pelitic mineral assemblage, i.e. chlorite–biotite zone, muscovite–andalusite zone and sillimanite–K‐feldspar zone with ascending metamorphic grade from north to south. Two isograds trend approximately east–west, which is oblique to the boundary between the metamorphic complex and the Tamana Granodiorite located on the southeast. The metamorphic condition of two pelitic rocks that occur in the muscovite–andalusite zone and sillimanite–K‐feldspar zone are estimated as 510 ±30 °C, 300 ±60 MPa and 720 ±30 °C, 620 ±60 MPa, respectively. Thermodynamic consideration reveals that use of the same geothermobarometer enables precise determination of the difference in pressure between the samples as 320 ±10 MPa. This indicates that the pelitic samples were metamorphosed at different depth by 11–12 km that is significantly larger than the geographic distance of 6.8 km between the sample localities. This also suggests that crustal thinning took place after the high‐T metamorphism. The high‐T metamorphic complex is, therefore, not of static contact metamorphism but of dynamic regional metamorphism. The present result combined with petrological and chronological similarities implies that this complex suffered the regional Ryoke metamorphism.     

Attenuation of compressional waves in peridotite measured as a function of temperature at 200 MPa

A technique has been developed to determine attenuation in rocks at high temperature using a gas-media, high-pressure apparatus. A pulse transmission technique and a spectral ratio method are used to study compressional seismic properties of rocks. Seismic waves are transmitted to and from the sample through buffer rods of mullite. The effect of seismic wave reflections within the sample assembly are cancelled out by taking ratios of the spectra measured at different temperatures. In order to obtain good signal-to-noise ratio for resolving the attenuation at high pressure and temperature, special care is taken in the sample assembly and the ultrasonic coupling between the sample, buffer rods and transducers. A very tight connection of the sample-buffer rod-transducer is essential for obtaining high frequency signals (>300 kHz) at high temperature. A small mass is attached to each outside end of the transducer to drive low frequency signals (<250 kHz) into the sample. Before attenuation measurements, the sample and the buffer rods are tightly compacted in a platinum tube at high pressure and room temperature to ensure pressure seal of the sample assembly. The frequency range of measurement covers 50 to 450 kHz for the sample. Attenuation is very small in the buffer rod compared to the sample for the entire temperature range of the study. Because of the small attenuation, a wide frequency band of 50 kHz to 3.2 MHz can be covered for investigating the attenuation in the buffer rod. The technique has been used to measure attenuation at high confining pressure, and temperatures including sub- and hyper-solidus of upper mantle rocks. Therefore, effects of partial melting on attenuation can be studied.The method is applied to the attenuation measurement in a peridotite as a function of temperature to 1225°C at 200 MPa confining pressure. At high temperature, signal amplitude decays more rapidly at high frequency than at low frequency, from which attenuation (andQ) can be determined using a spectral ratio method. No frequency dependence ofQ is resolved for both the sample and the buffer rod over the entire temperature and frequency ranges of the measurement. The results show thatQ decreases rapidly with increasing temperature even in the temperature range below the solidus of peridotites. Such temperature sensitivity ofQ is probably more useful to probe thermal structure in the upper mantle than that of conductivity at temperatures below the solidus. The results in this study are compared with available seismic velocity, electrical conductivity and solidus data for peridotites, suggesting that there is no discontinuous change in both mechanical and electrical properties of peridotites at the solidus temperature. Even at hypersolidus temperatures, it appears that velocity drops and conductivity increases continuously (not abruptly) with increasing melt fraction. This implies that mechanical and electrical properties of the upper mantle will gradually change at the boundary where the geotherm crosses the solidus.     

Synthesis of kosmochlor and phase equilibria in the join CaMgSi2O6-NaCrSi2O6

Kosmochlor (NaCrSi₂O₆) was synthesized by the flux method from melts along the join Na₂O·2 SiO₂-Na₂O·Cr₂O₃·4 SiO₂ at 1000° C in air, and isolated by dissolving the glassy matrix with hydrofluoric and perchloric acids. The join CaMgSi₂O₆-NaCrSi₂O₆ was studied at 1 atmosphere in air by the quenching technique at temperatures between 900° and 1450° C, using mixtures of kosmochlor and diopside crystals or diopside glass as starting materials. The phases are diopside solid solution, kosmochlor, spinel (Mg-chromite), eskolaite (Cr₂O₃) and glass. The maximum solubility of kosmochlor in diopside is 24 wt percent at 1140° C, while diopside is not soluble at all in kosmochlor, resulting in the existence of a wide range of immiscibility. Petrologic significance of the results is discussed.     

Stationary Rossby-waves in a two-layer model produced by a submarine ridge

Characteristics of stationary Rossby-waves are investigated by use of equations expressing conservation of potential vorticity for a quasi-geostrophic two-layer flow of an incompressible inviscid fluid on a beta-plane. A two-dimensional problem is discussed; the basic flow is directing eastward without horizontal shear, and the stationary Rossby-waves are produced by a two-dimensional bump extending in the north-south direction which is located on a sloping bottom with a constant north-south gradient.When the slope is steep, stream-lines are scarcely displaced being influenced by the bump. When the bottom has no north-south inclination, however, there is a case where group velocity of Rossby-waves produced by the bump directs westward and the stream-lines meander sinusoidally at the west of the bump. This phenomenon is called upstream influence. These two special cases have a possibility to explain two different mean paths for the Kuroshio namely, the path along the continental slope and the path with a stationary meander at the west of the Izu-Ogasawara Ridge.     

Matrix method evaluating an internal radiation field in a plane-parallel atmosphere

Suppose that the atmosphere consists of homogeneous sublayers whose scattering and transmission matrices are known, then it is shown in this paper how to determine the intensity and polarization parameters of internal radiation field by the matrix method.