Quasi-static slip on the plate boundary associated with the 2003 M8.0 Tokachi-oki and 2004 M7.1 off-Kushiro earthquakes, Japan

We analyzed small repeating earthquakes recorded over a 13-year period and GPS data recorded over an 8-month period to estimate interplate quasi-static slip associated with the 2003 Tokachi-oki earthquake (M8.0) and the 2004 off-Kushiro earthquake (M7.1). The repeating-earthquake analysis revealed that the slip rate near the source region of the Tokachi-oki earthquake was relatively low (< 5 cm/year) prior to the earthquake; however, in the last 3 years leading up to the event, a minor acceleration in slip occurred upon the deeper extension of the coseismic slip area of the earthquake. Repeating-earthquake and GPS data indicate that large amounts of afterslip occurred around the rupture area following the earthquake; the afterslip mainly propagated to the east of the coseismic slip area. We also infer that the occurrence of the 2004 off-Kushiro earthquake, located about 100 km northeast of the epicenter of the Tokachi-oki earthquake, was advanced by the afterslip associated with the Tokachi-oki earthquake.     

A simple two-way coupling method of BEM and VOF model for random wave calculations

A numerical method, which combines the boundary element method (BEM) and the volume of the fluid method (VOF method), has been presented to solve wave–structure interactions; the intense wave motion at the proximity of the structure is modeled by the VOF method and the rest of the fluid region is modeled by the BEM. The combined method can considerably reduce the time-consuming VOF domain, and thus practically makes it possible to apply the VOF method for random wave calculations, in which long time computations are usually required to obtain statistically meaningful results, and therefore the use of the single-VOF model often becomes prohibitive in terms of computational time and storage memories. A VOF model CADMAS-SURF, which is based on SMAC scheme and had been constructed by a number of VOF researchers in coastal engineering in Japan, is used in the combined BEM–VOF model. The two-way coupling treatment, which enables us to deal with bidirectional wave propagations, which was originally given for the SOLA-VOF model by Yan et al. (2003a) and later improved by Kim et al. (2007), was modified for the SMAC scheme. The coupling treatments are described in detail in the paper. The validity of the combined BEM–VOF model was investigated by comparing the numerical results with the theoretical results for the propagations of Stokes 5th order waves and random waves.     

Enhancement/reduction of biological pump depends on ocean circulation in the sea-ice reduction regions of the Arctic Ocean

The biological pump is a central process in the ocean carbon cycle, and is a key factor controlling atmospheric carbon dioxide (CO₂). However, whether the Arctic biological pump is enhanced or reduced by the recent loss of sea ice is still unclear. We examined if the effect was dependent on ocean circulation. Melting of sea ice can both enhance and reduce the biological pump in the Arctic Ocean, depending on ocean circulation. The biological pump is reduced within the Beaufort Gyre in the Canada Basin because freshwater accumulation within the gyre limits nutrient supply from deep layers and shelves hence inhibits the growth of large-bodied phytoplankton. Conversely, the biological pump is enhanced outside the Beaufort Gyre in the western Arctic Ocean because of nutrient supply from shelves and greater light penetration, enhancing photosynthesis, caused by the sea ice loss. The biological pump could also be enhanced by sea ice loss in the Eurasian Basin, where uplifted isohaline surfaces associated with the Transpolar Drift supply nutrients upwards from deep layers. New data on nitrate uptake rates are consistent with the pattern of enhancement and reduction of the Arctic biological pump. Our estimates indicate that the enhanced biological pump can be as large as that in other oceans when the sea ice disappears. Contrary to a recent conclusion based on data from the Canada Basin alone, our study suggests that the biological CO₂ drawdown is important for the Arctic Ocean carbon sink under ice-free conditions.     

Energy dissipation behaviour of shear panels made of low yield steel

This paper presents the results of a pilot test conducted for evaluating the energy dissipation behaviour of shear panels made of low yield steel whose 0.2 per cent offset yield stress is 120 MPa. A total of six full-scale shear panels were tested with the loading condition, stiffener spacing, and magnitude of axial force as test variables. The shear panels tested yielded at a shear force that is approximately 1/3 of the yield shear force of equivalent shear panels made of common mild steel. Shear panels with proper stiffener arrangement exhibited stable hysteresis, thus ensuring large energy dissipation capacity. Sufficient strain hardening was observed in the shear panels tested, with their energy dissipation capacity about 1.5 times larger than that of an equivalent linear-elastic and perfect-plastic system. Plate buckling did not lead the shear panels to immediate degradation in their energy dissipation capacity. Post-buckling resistance was found to be a subject that requires further studies for quantifying the performance of shear panels made of low yield stress steel as hysteretic dampers.     

Miocene provenance change in Himalayan foreland basin and Bengal Fan sediments,with special reference to detrital garnet chemistry

Bengal Fan Miocene sediments were collected during International Ocean Discovery Program Expedition 354 and investigated using petrographic and detrital garnet chemistry analyses. The Miocene Siwalik Group, which is composed of sediments deposited in the Himalayan foreland basin, was also analyzed for comparison with the Bengal Fan data for the provenance change during the Miocene. Our petrographic analyses revealed that the Miocene sediments of the Bengal Fan and Siwalik Group consist predominantly of Higher Himalayan Crystalline (HHC)-derived detritus such as chloritoid, staurolite, sillimanite, and/or kyanite, which appear among the accessory minerals. The chemistry of the detrital garnet varies across the stratigraphy; most of the garnet is rich in almandine and poor in spessartine and pyrope. However, pyrope-rich garnet, which is considered to originate from the HHC core (granulite facies), was found in the lower to upper Miocene deposits. The deposition of HHC-derived detrital garnet began before the Middle Miocene (15 Ma) and before the Late Miocene (10–9 Ma) in the Siwalik Group. The Bengal Fan data, by contrast, indicated that pyrope-rich garnet appeared in the Early Miocene (17.3 Ma) and Late Miocene (8.5–6.5 Ma). We conclude that the Bengal Fan sediments record the erosion of the HHC zone since the Early Miocene that appears in the Siwalik sediments. Furthermore, we found that the HHC-derived inputs decreased from the late Middle Miocene (12 Ma) to the early Middle Miocene (10 Ma) in both the Nepal Himalaya foreland basin and the Bengal Fan. The disappearance of the HHC-derived detritus is probably the result of dilution by Lesser Himalayan detritus, which suggests that the Lesser Himalayan zone, which is composed of metamorphosed and unmetamorphosed sedimentary rocks, was uplifted.     

Horizontal structure of planetary-scale waves at the cloud top of Venus deduced from Galileo SSI images with an improved cloud-tracking technique

An improved cloud tracking method for deriving wind velocities from successive planetary images was developed. The new method incorporates into the traditional cross-correlation method an algorithm that corrects for erroneous cloud motion vectors by re-determining the most plausible correlation peak among all of the local maxima on the correlation surface by comparing each vector with its neighboring vectors. The newly developed method was applied to the Venusian violet images obtained by the Solid State Imaging system (SSI) onboard the Galileo spacecraft during its Venus flyby. Although the results may be biased by the choice of spatial scale of atmospheric features, the cloud tracking is the most practical mean of estimating the wind velocities with extensive spatial and temporal coverage. The two-dimensional distribution of the horizontal wind vector field over 5 days was obtained. It was found from these wind maps that the solar-fixed component in 1990 was similar to that in 1982 obtained by the Pioneer Venus orbiter. The deviation of the instantaneous zonal wind field from the solar-fixed component shows a distinct wavenumber-1 structure in the equatorial region. On the assumption that this structure is a manifestation of an equatorial Kelvin wave, the phase relationship between the zonal wind and the cloud brightness suggests a short photochemical lifetime of the violet absorber. The momentum deposition by this Kelvin wave, which is subject to radiative damping, would induce a westward mean-wind acceleration of ～0.3 m s^?1 per Earth day.     

Modelling exposure of oceanic higher trophic-level consumers to polychlorinated biphenyls: Pollution ‘hotspots’ in relation to mass mortality events of marine mammals

Marine mammals in the past mass mortality events may have been susceptible to infection because their immune systems were suppressed through the bioaccumulation of environmental pollutants such as polychlorinated biphenyls (PCBs). We compiled mortality event data sets of 33 marine mammal species, and employed a Finely-Advanced Transboundary Environmental model (FATE) to model the exposure of the global fish community to PCB congeners, in order to define critical exposure levels (CELs) of PCBs above which mass mortality events are likely to occur. Our modelling approach enabled us to describe the mass mortality events in the context of exposure of higher-trophic consumers to PCBs and to identify marine pollution ‘hotspots’ such as the Mediterranean Sea and north-western European coasts. We demonstrated that the CELs can be applied to quantify a chemical pollution Planetary Boundary, under which a safe operating space for marine mammals and humanity can exist.     

Calorimetric measurements of fusion enthalpies for Ni2SiO4 and Co2SiO4 olivines and application to olivine-liquid partitioning

Calorimetric measurements of fusion enthalpies for Ni₂SiO₄ and Co₂SiO₄ olivines were carried out using a high-temperature calorimeter, and Ni and Co partitioning between olivine and silicate liquid was analyzed using the measured heats of fusion. The fusion enthalpy of Co₂SiO₄ olivine measured by transposed-temperature drop calorimetry was 103 ± 15 kJ/mol at melting point (1688 K). The fusion enthalpy of Ni₂SiO₄ olivine was calculated based on the enthalpies of liquids in the system An₅₀Di₅₀-Ni₂SiO₄ measured by transposed-temperature drop calorimetry at 1773 K, and was 221 ± 26 kJ/mol at its metastable melting point (1923 K). The fusion enthalpy of Ni₂SiO₄ is the largest among those of olivine group, this is caused by the large crystal field stabilization energy of six-coordinated Ni²⁺ in olivine. The larger fusion enthalpy of Ni₂SiO₄ can account for the large and variable partition coefficient of Ni between olivine and silicate liquid. Based on the comparison between partition coefficients calculated from thermodynamic data and those observed in partition experiments, it is considered that the magnitude of partition coefficients is primarily dependent on the heats of fusion of the components. Furthermore, the activity coefficients for Ni-, Co- and Mn-bearing components in magmatic liquid are nearly of the same magnitude.