Distribution of local Cs anomalies on the seafloor near the Fukushima Dai-ichi Nuclear Power Plant

An estimated 3.5 ± 0.7 × 10¹⁵ Bq of ¹³⁷Cs is thought to have been discharged into the ocean following the melt down at Fukushima Dai-ichi Nuclear Power Plant (F1NPP). While efforts have been made to monitor seafloor radiation levels, the sampling techniques used cannot capture the continuous distribution of radionuclides. In this work, we apply in situ measurement techniques using a towed gamma ray spectrometer to map the continuous distribution of ¹³⁷Cs on the seafloor within 20 km of the F1NPP. The results reveal the existence of local ¹³⁷Cs anomalies, with levels of ¹³⁷Cs an order of magnitude higher than the surrounding seafloors. The sizes of the anomalies mapped in this work range from a few meters to a few hundreds of meters in length, and it is demonstrated that the distribution of these anomalies is strongly influenced by meter scale features of the terrain.     

Successful application of multiscale methods in a real reservoir simulator environment

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.     

Simple determination of dimethylsulfide (DMS) and dimethylsulfoniopropionate (DMSP) using solid-phase microextraction and gas chromatography-mass spectrometry

Solid-phase microextraction (SPME) is a simple, sensitive and less destructive method for the determination of dimethylsulfide (DMS) in seawater. Combined with detection by gas chromatography-mass spectrometry (GC-MS), the method had sufficient sensitivity (minimum detectable concentration of DMS was 0.05 nM), and practical levels of reproducibility (relative standard deviation ≤7%) and linearity (r ² > 0.995) over a wide concentration range (0.5 to 910 nM). The protocol developed was applied to a Sagami Bay water sample to determine concentrations of DMS and DMSP, and in situ DMSP-lyase activity.     

The effect of a hydrous phase on P‐wave velocity anisotropy within a detachment shear zone in the slow‐spreading oceanic crust: A case study from the Godzilla Megamullion,Philippine Sea

We studied the contributions of plagioclase, clinopyroxene, and amphibole to the P‐wave velocity properties of gabbroic mylonites of the Godzilla Megamullion (site KH07‐02‐D18) in the Parece Vela Rift of the central Parece Vela Basin, Philippine Sea, based on their crystal‐preferred orientations (CPOs), mineral modes, and elastic constants and densities of single crystals. The gabbroic mylonites have been classified into three types based on their microstructures and temperature conditions: HT1, HT2 and medium‐temperature (MT) mylonites. The P‐wave velocity properties of the HT1 mylonite are dominantly influenced by plagioclase CPOs. Secondary amphibole occurred after deformation in the HT1 mylonite, so that its effect on P‐wave velocity anisotropy is minimal due to weak CPOs. Although the HT2 mylonite developed deformation microstructures in the three minerals, the P‐wave velocity properties of the HT2 mylonite are essentially isotropic, resulting from the destructive interference of different P‐wave velocity anisotropy patterns produced by the distinct CPOs of the three constituent minerals (i.e., plagioclase, clinopyroxene, and amphibole). The P‐wave velocity properties of the MT mylonite are influenced mainly by amphibole CPOs, whereas the effect of plagioclase CPOs on P‐wave velocity anisotropy becomes very small with a decrease in the intensity of plagioclase CPOs. As a result, the gabbroic mylonites tend to have weak P‐wave velocity anisotropy in seismic velocity, although their constituent minerals show distinct CPOs. Such weakness in the whole‐rock P‐wave velocity anisotropy could result from the destructive contributions of the different mineral CPOs with respect to the structural framework (foliation and lineation). These results show that amphibole has a high potential for P‐wave velocity anisotropy by aligning both crystallographically and dimensionally during deformation in the hydrous oceanic crust. The results also suggest that the effect of a hydrous phase on P‐wave velocity anisotropy within the detachment shear zone in a slow‐spreading oceanic crust varies depending on the degree of deformation and on the timing of hydrothermal activity.     

Long-term changes in the fish community structure from the Tsushima warm current region of the Japan/East Sea with an emphasis on the impacts of fishing and climate regime shift over the last four decades

Japanese fisheries production in the Japan/East Sea between 1958 and 2003 increased to their peak (1.76 million tons) in the late 1980s and decreased abruptly with the collapse of Japanese sardine. Catch results for 58 fisheries and various environmental time-series data sets and community indices, including mean trophic level (MTL) and Simpson’s diversity index (DI), were used to investigate the impacts of fishing and climate changes on the structure of the fish community in the Tsushima warm current (TWC) region of the Japan/East Sea. The long-term trend in fisheries production was largely dependent on the Japanese sardine that, as a single species, contributed up to 60% of the total production in the Japanese waters of the Japan/East Sea during the late 1980s. Excluding Japanese sardine, production of the small pelagic species was higher during 1960s and 1990s but lower during 1970s and 1980s. This variation pattern generally corresponds with the trend in water temperature, warmer before early 1960s and after 1990s but colder during 1970s and 1980s. The warm-water, large predatory fishes and cold water demersal species show opposite responses to the water temperature in the TWC region, indicating the significant impact of oceanic conditions on fisheries production of the Japan/East Sea. Declines in demersal fishes and invertebrates during 1970s and 1980s suggested some impact of fishing. MTL and DI show a similar variation pattern: higher during 1960s and 1990s but lower during 1970s and 1980s. In particular, the sharp decline during the 1980s resulted from the abundant sardine catches, suggesting that dominant species have a large effect on the structure of the fish community in the Japan/East Sea. Principal component analysis for 58 time-series data sets of fisheries catches suggested that the fish community varied on inter-annual to inter-decadal scales; the abrupt changes that occurred in the mid-1970s and late 1980s seemed to correspond closely with the climatic regime shifts in the North Pacific. These results strongly suggest that the structure of the fish community in the Japan/East Sea was largely affected by climatic and oceanic regime shifts rather than by fishing. There is no evidence showing “fishing down food webs” in the Japan/East Sea. However, in addition to the impacts of abrupt shifts that occurred in the late 1980s, the large predatory and demersal fishes seem to be facing stronger fishing pressure with the collapse of the Japanese sardine.     

Deuteron whistler and trans-equatorial propagation of the ion cyclotron whistler

New ion cyclotron whistlers which have the asymptotic frequency of one half the local proton gyrofrequency, $\text{G}{}_{\mathrm{p}}\text{2}$, and the minimum (or equatorial) proton gyrofrequency, G_pm, along the geomagnetic field line passing through the satellite have been found in the low-latitude topside ionosphere from the spectrum analysis of ISIS VLF electric field data received at Kashima, Japan. Ion cyclotron whistlers with asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are observed only in the region of $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{2}$ or rarely $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{4}$, where B is the local magnetic field and B_m is the mini magnetic field along the geomagnetic field line passing through the satellite.The particles with one half the proton gyrofrequency may be the deuteron or alpha particle. Theoretical spectrograms of the electron whistlers (R-mode) and the ion cyclotron whistlers (L-mode) propagating along the geomagnetic field lines are computed for the appropriate distributions of the electron density and the ionic composition, and compared with the observed spectrograms.The result shows that the ion cyclotron whistler with the asymptotic frequency of $\text{G}{}_{\mathrm{p}}\text{2}$ is the deuteron whistler, and that the ion cyclotron whistlers with the asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are caused by the trans-equatorial propagation of the proton or deuteron whistler from the other hemisphere.     

A Preliminary Study of Oceanic Bomb Radiocarbon Inventory in the North Pacific during the Last Two Decades

Radiocarbon and total carbonate data were obtained near the 1973 GEOSECS stations in the North Pacific along 30°N and along 175°E between 1993 and 1994. In these stations, we estimated radiocarbon originating from atomic bomb tests using tritium, trichlorofluoromethane and silicate contents. The average penetration depth of bomb radiocarbon during the two decades has deepened from 900 m to 1300 m. Bomb radiocarbon inventories above the average value for the whole North Pacific were found widely in the western subtropical region around 30°N both in the 1970s and 1990s, and its area in the 1990s was broader than that in the 1970s. In most of the North Pacific, while the bomb radiocarbon has decreased above 25.4, the bomb radiocarbon flux below 25.4 was over 1 × 10¹² atom m^-2yr^-1 in the subtropical region around 30°N. In the tropical area south of 20°N, the bomb radiocarbon inventory below 25.4 increased from zero to over 10 × 10¹² atom m^-2 during the last three decades. These distributions suggest that the bomb radiocarbon removed from the surface is currently accumulated with bomb ¹⁴C flux of over 1 × 10¹² atom m^-2yr^-1 below 25.4 in the subtropical region, mainly by advection from the higher latitude, and that part of the accumulated bomb ¹⁴C gradually spread southward with about 30 years.