Microtremor array method using spatial autocorrelation analysis of Rayleigh-wave data

Journal of Seismology - Microtremor array measurements, and passive surface wave methods in general, have been increasingly used to non-invasively estimate shear-wave velocity structures for...     

Transitions in thermal convection with strongly temperature-dependent viscosity in a wide box

Numerical models are systematically presented for time-dependent thermal convection of Newtonian fluid with strongly temperature-dependent viscosity in a two-dimensional rectangular box of aspect ratio 3 at various values of the Rayleigh number Ra_b defined with viscosity at the bottom boundary up to 1.6×10⁸ and the viscosity contrast across the box r_η up to 10⁸. We found that there are two different series of bifurcations that take place as r_η increases. One series of bifurcations causes changes in the behavior of the thermal boundary layer along the surface boundary from small-viscosity-contrast (SVC) mode, through transitional (TR) mode, to stagnant-lid (ST) mode, or from SVC mode directly to ST mode, depending on Ra_b. Another series of bifurcations causes changes in the aspect ratio of convection cells; convection with an elongated cell can take place at moderate r_η (10³–10^5.5 at Ra_b=6×10⁶), while only convection of aspect ratio close to 1 takes place at small r_η and large r_η. The parameter range of r_η and Ra_b for elongated-cell convection overlaps the parameter range for SVC and ST modes and include the entire parameter range for TR mode. In the elongated-ST regime, the lid of highly viscous fluid along the top boundary is not literally ‘stagnant’ but can horizontally move at a velocity high enough to induce a convection cell with aspect ratio much larger than 1.     

Consequences of implementing a reservoir operation algorithm in a global hydrological model under multiple meteorological forcing

We investigated dam behaviours during high-flow events and their robustness against perturbations in meteorological conditions using the H08 global hydrological model. Differences in these behaviours were examined by comparing simulation runs, with and without dams and using multiple meteorological datasets, at a case-study site, Fort Peck Dam on the Missouri River, USA. The results demonstrated that dam-regulated river flow reduced temporal variability over large time periods and also dampened inter-forcing discrepancies in river discharge (smoothing effects). However, during wet years, differences in peak flow were accentuated downstream of the dam, resulting in divergence in simulated peak flow across the meteorological forcing (pulsing effect). The pulsing effect was detected at other major dams in global simulations. Depending upon the meteorological forcing, the dams act as a selective filter against high-flow events. Synergy between a generic dam scheme and differences in meteorological forcing data might introduce additional uncertainties in global hydrological simulations.     

Vector Magnetic Field Synoptic Charts from the <Emphasis Type="Italic">Helioseismic and Magnetic Imager</Emphasis> (HMI)

Vector magnetic field synoptic charts from the Helioseismic and Magnetic Imager (HMI) are now available for each Carrington Rotation (CR) starting from CR 2097 in May 2010. Synoptic charts are produced using 720-second cadence full-disk vector magnetograms remapped to Carrington coordinates. The vector field is derived from the Stokes parameters (\(I, Q, U, V\)) using a Milne–Eddington-based inversion model. The \(180^{\circ}\) azimuth ambiguity is resolved using the minimum energy algorithm for pixels in active regions and for strong-field pixels (the field is greater than about 150 G) in quiet-Sun regions. Three other methods are used for the rest of the pixels: the potential-field method, the radial acute-angle method, and the random method. The vector field synoptic charts computed using these three disambiguation methods are evaluated. The noise in the three components of the vector magnetic field is generally much higher in the potential-field method charts. The component noise levels are significantly different in the radial-acute charts. However, the noise levels in the random-method charts are lower and comparable. The assumptions used in the potential-field and radial-acute methods to disambiguate the weak transverse field introduce bias that propagates differently into the three vector-field components, leading to unreasonable pattern and artifacts, whereas the random method appears not to introduce any systematic bias. The current sheet on the source surface, computed using the potential-field source-surface model applied to random-method charts, agrees with the best solution (the result computed from the synoptic charts with the minimum energy algorithm applied to each and every pixel in the vector magnetograms) much better than the other two. Differences in the synoptic charts determined with the best method and the random method are much smaller than those from the best method and the other two. This comparison indicates that the random method is better for vector field synoptic maps computed from near-central meridian data. The vector field synoptic charts provided by the Joint Science Operations Center (JSOC) are therefore produced with the random method.     

Relation Between Coronal Hole Areas and Solar Wind Speeds Derived from Interplanetary Scintillation Measurements

We investigate the relation between coronal hole (CH) areas and solar wind speeds during 1995?–?2011 using the potential field (PF) model analysis of magnetograph observations and interplanetary scintillation (IPS) observations by the Institute for Space-Earth Environmental Research (formerly Solar-Terrestrial Environment Laboratory) of Nagoya University. We obtained a significant positive correlation between the CH areas (\(A\)) derived from the PF model calculations and solar wind speeds (\(V\)) derived from the IPS observations. The correlation coefficients between them are usually high, but they drop significantly in solar maxima. The slopes of the \(A\)?–?\(V\) relation are roughly constant except for the period around solar maximum, when flatter or steeper slopes are observed. The excursion of the correlation coefficients and slopes at solar maxima is ascribed partly to the effect of rapid structural changes in the coronal magnetic field and solar wind, and partly to the predominance of small CHs. It is also demonstrated that \(V\) is inversely related to the flux expansion factor (\(f\)) and that \(f\) is closely related to \(A^{-1/2}\); hence, \(V \propto A^{1/2}\). A better correlation coefficient is obtained from the \(A^{1/2}\)?–?\(V\) relation, and this fact is useful for improving space weather predictions. We compare the CH areas derived from the PF model calculations with He i 1083 nm observations and show that the PF model calculations provide reliable estimates of the CH area, particularly for large \(A\).     

Survival of eclogite xenolith in a Cretaceous granite intruding the Central Dabieshan migmatite gneiss dome (Eastern China) and its tectonic implications

Investigation of an eclogite xenolith, discovered in a Cretaceous granite from the Central Domain of the Dabieshan massif in eastern China, yields new petrological insights into the high to ultrahigh-pressure metamorphism, experienced by the Qinling-Dabie orogen. Prior to inclusion as a xenolith in the granite during the Early Cretaceous, this eclogite xenolith had recorded a complex metamorphic evolution that complies with subduction and exhumation processes experienced by the continental crust of the South China Block. Well-preserved mineral parageneses substantiate the prograde and retrograde stages revealed by inclusions in porphyroblastic garnet and zoned minerals such as garnet, omphacite and amphibole in the matrix. The relatively low P/T re-equilibration during a late metamorphic stage was textually inferred by the presence of aluminous and calcic-subcalcic amphiboles such as katophorite, edenite, taramite and pargasite as main matrix phases. According to our U/Pb, Rb/Sr and new ⁴⁰Ar/³⁹Ar geochronological results, namely109 ± 1 and 112 ± 2 Ma plateau ages for muscovite and amphiboles, respectively, two successive but distinct cooling stages account for the thermal history of the granite–migmatite gneiss dome that forms the Central Dabieshan Domain. We argue that prior to the Cretaceous doming, the Central Dabieshan Domain experienced a tectono-metamorphic evolution similar to that observed in the high-pressure to ultra high-pressure units developed in the Southern Dabieshan Domain and Hong’an massif.     

Trace-element characteristics of east–west mantle geochemical hemispheres

Recent statistical analyses on the isotopic compositions of oceanic, arc, and continental basalts have revealed that the Earth's mantle is broadly divided into eastern and western hemispheres. The present study aimed to characterize the isotopically defined east–west geochemical hemispheres using trace-element concentrations. Basalt data with Rb, Sr, Nd, Sm, Pb, Th, and U in addition to the isotopic ratios ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb were selected mostly from the GEOROC and PetDB databases. A total of 4787 samples were used to investigate the global geochemical variations. The results show that the wide trace-element variations are broadly explained by the melting of melt-metasomatized and fluid-metasomatized mantle sources. The larger amount of the fluid component derived from subducted plates in the eastern hemisphere than that in the western hemisphere is inferred from the basalts. These characteristics support the hypothesis that focused subduction towards the supercontinent created the mantle geochemical hemispheres.     

Measurement of persistent organic pollutants (POPs) in plastic resin pellets from remote islands: toward establishment of background concentrations for International Pellet Watch

Plastic resin pellets collected from remote islands in the Pacific, Atlantic, and Indian Oceans and the Caribbean Sea were analyzed for polychlorinated biphenyls (PCBs), dichloro-diphenyltrichloroethane and its degradation products (DDTs), and hexachlorocyclohexanes (HCHs). Concentrations of PCBs (sum of 13 congeners) in the pellets were 0.1-9.9 ng/g-pellet. These were 1-3 orders of magnitude smaller than those observed in pellets from industrialized coastal shores. Concentrations of DDTs in the pellets were 0.8-4.1 ng/g-pellet. HCH concentrations were 0.6-1.7 ng/g-pellet, except for 19.3 ng/g-pellet on St. Helena, where current use of lindane is likely influence. This study provides background levels of POPs (PCBs<10 ng/g-pellet, DDTs <4 ng/g-pellet, HCHs <2 ng/g-pellet) for International Pellet Watch. Sporadic large concentrations of POPs were found in some pellet samples from remote islands and should be considered in future assessments of pollutants on plastic debris.     

Subduction of mantle wedge peridotites: Evidence from the Higashi-akaishi ultramafic body in the Sanbagawa metamorphic belt

The Higashi-akaishi garnet-bearing ultramafic body in the Sanbagawa metamorphic belt, Southwest Japan, represents a rare example of oceanic-type ultrahigh-pressure metamorphism. The body of 2 km × 5 km is composed mostly of anhydrous dunite with volumetrically minor lenses of clinopyroxene-rich rocks. Dunite samples contain high Ir-type platinum group elements (PGE) and Cr in bulk rocks, high Mg and Ni in olivine, and high Cr in spinel. On the other hand, clinopyroxene-rich rocks contain low concentrations of Ir-type PGE and Cr, high concentrations of fluid-mobile elements in bulk rocks, and low Ni and Mg in olivine. Clinopyroxene is diopsidic with low Al₂O₃. The compositions of bulk rocks and mineral chemistry of spinel, olivine, and clinopyroxene suggest that the olivine-dominated rocks are residual mantle peridotites after high degrees of influx partial melting, and that the clinopyroxene-rich rocks are cumulates of subduction-related melts. Thus, the Higashi-akaishi ultramafic body originated from the interior of the mantle wedge, most likely the forearc upper mantle. It was then incorporated into the Sanbagawa subduction channel by a mantle flow, and underwent high pressure metamorphism to a depth greater than 100 km. Such a strong active flow in the mantle wedge is likely facilitated by the lack of serpentinites along the interface between the slab and the overlying mantle, as it was too hot for serpentine. These unusually hot conditions and strong active mantle flow may reflect conditions in the earliest stage of development of subduction, and may have been maintained by massive upwelling and subsequent eastward flow of asthenospheric mantle in the northeastern Asian continent in Cretaceous time when the Sanbagawa belt began to form.     

Stratigraphic records of tsunamis along the Japan Sea,southwest Hokkaido,northern Japan

The stratigraphy of tsunami deposits along the Japan Sea, southwest Hokkaido, northern Japan, reveals tsunami recurrences in this particular area. Sandy tsunami deposits are preserved in small valley plains, whereas gravelly deposits of possible tsunami origin are identified in surficial soils covering a Holocene marine terrace and a slope talus. At least five horizons of tsunami events can be defined in the Okushiri Island, the youngest of which immediately overlies the Ko‐d tephra layer (1640 AD) and was likely formed by the historical Oshima‐Ohshima tsunami in 1741 AD. The four older tsunami deposits, dated using accelerator mass spectrometry ¹⁴C, were formed at around the 12th century, 1.5–1.6, 2.4–2.6, and 2.8–3.1 ka, respectively. Tsunami sand beds of the 1741 AD and circa 12th century events are recognized in the Hiyama District of Hokkaido Island, but the older tsunami deposits are missing. The deposits of these two tsunamis are found together at the same sites and distributed in regions where wave heights of the 1993 tsunami (Hokkaido Nansei‐oki earthquake, Mw = 7.7) were less than 3 m. Thus, the 12th century tsunami waves were possibly generated near the south of Okushiri Island, whereas the 1993 tsunami was generated towards the north of the island. The estimated recurrence intervals of paleotsunamis, 200–1100 years with an average of 500 years, likely represents the recurrence interval of large earthquakes which would have occurred along several active faults offshore of southwest Hokkaido.