High parallelization efficiency in barotropic-mode computation of ocean models based on multi-grid boundary ghost area

In general, barotropic-mode computation requires the largest communication time in ocean models because of its iterative nature, when parallel computation is performed based on regional partitioning. In this study, reduction in the overhead time included in the communication in the parallel computation of the barotropic mode is considered to achieve a high parallelization efficiency of ocean models. We verify that the reduction in the communication frequency based on the multi-grid boundary ghost area reduces the total communication time. We find that this is because the sum of several kinds of overhead time for communication occupies a fairly large part of the total communication time. We discuss the trade-off between the decrease in communication time and the increase in computation time due to increased boundary area in such cases, leading to a determination of a suitable width of the multi-grid boundary that minimizes the total required time. We also discuss the efficiencies of the one- and two-dimensional partitioning of the model domain, when a multi-grid boundary is used. In general, two-dimensional partitioning is more efficient than one-dimensional for large partitioning numbers.     

Comparison of bromodeoxyuridine immunoassay with tritiated thymidine radioassay for measuring bacterial productivity in oceanic waters

Recently, bromodeoxyuridine (BrdU) has been successfully applied to the measurement of bacterial productivity as an alternative to tritiated thymidine (³H-TdR), which is widely used but often restricted by regulations, particularly in field settings. Here, I report improvements to existing BrdU methods to simplify procedures and increase sensitivity. The feasibility of the method was tested measuring bacterial production in low-productive waters. The method provided radioisotope-free measurements of bacterial production rates at shorter (∼1 h) on-board processing time of samples than previously reported procedures. It was applicable to the detection of rates ranging from 0.021 to 2.7 pmol BrdU l⁻¹h⁻¹. BrdU incorporation rates measured by immunoassay showed a statistically significant correlation with ³H-TdR incorporation rates measured by radioassay (r = 0.74, n = 24, p < 0.001). The linear regression obtained (BrdU = 0.80[³H-TdR] − 0.016) showed a similar relationship to previously reported regressions (BrdU = 0.65[³H-TdR] + 0.12, [³H-BrdU] = 0.69[³H-TdR] − 0.81). There were no statistically significant differences among these regression lines. These results suggest that the method described here provides a non-radioisotopic productivity measurement of bacteria in oceanic epipelagic waters, while retaining continuity of the data with other existing ³H-TdR and BrdU methods.     

A shipborne wave-recording system with digital data processing

A shipborne wave-recording system which consists of a sonic wave gauge, accelerometers, gyroscopes and a computer system is described. Signals from the measuring apparatus are fed directly into a shipborne digital computer system at a prescribed sampling rate. The time series of wave heights and the acceleration are transformed into Fourier series using an algorithm of Fast Fourier Transform. Errors contained in the observed wave heights due to ship motion are corrected in the Fourier series by using the Fourier coefficients for the vertical acceleration. Power spectra and waveforms can also be calculated in a short time with this system from Fourier coefficients. Examples of the observational results obtained in the central part of the East China Sea in 1969 are presented.     

Detection of a new interstellar molecular ion, H2COH+ (protonated formaldehyde)

A new interstellar molecular ion, H2COH+ (protonated formaldehyde), has been detected toward Sgr B2, Orion KL, W51, and possibly in NGC 7538 and DR21(OH). Six transitions were detected in Sgr B2(M). The 1(1,0)-1(0,1) transition was detected in all sources listed above. Searches were also made toward the cold, dark clouds TMC-1 and L134N, Orion (3N, 1E), and a red giant, IRC + 10216, without success. The excitation temperatures of H2COH+ are calculated to be 60-110 K, and the column densities are on the order of 10(12)-10(14) cm-2 in Sgr B2, Orion KL, and W51. The fractional abundance of H2COH+ is on the order of 10(-11) to 10-(9), and the ratio of H2COH+ to H2CO is in the range 0.001-0.5 in these objects. The values in Orion KL seem to be consistent with the "early time" values of recent model calculations by Lee, Bettens, & Herbst, but they appear to be higher than the model values in Sgr B2 and W51 even if we take the large uncertainties of column densities of H2CO into account. We suggest production routes starting from CH3OH may play an important role in the formation of H2COH+.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

Internal deformation of Sikhote Alin volcanic belt, far eastern Russia: Paleocene paleomagnetic results

We present paleomagnetic results of Paleocene welded tuffs of the 53–50 Ma Bogopol Group from the northern region (46°N, 137°E) of the Sikhote Alin volcanic belt. Characteristic paleomagnetic directions with high unblocking temperature components above 560 °C were isolated from all the sites. A tilt-corrected mean paleomagnetic direction from the northern region is D=345.8°, I=49.9°, α₉₅=14.6° (N=9). The reliability of the magnetization is ascertained through the presence of normal and reversed polarities. The mean paleomagnetic direction from the northern region of the Sikhote Alin volcanic belt reflects a counterclockwise rotation of 29° from the Paleocene mean paleomagnetic direction expected from its southern region. The counterclockwise rotation of 25° is suggested from the paleomagnetic data of the Kisin Group that underlies the Bogopol Group. These results establish that internal tectonic deformation occurred within the Sikhote Alin volcanic belt over the past 50 Ma. The northern region from 44.6° to 46.0°N in the Sikhote Alin volcanic belt was subjected to counterclockwise rotational motion through 29±17° with respect to the southern region. The tectonic rotation of the northern region is ascribable to relative motion between the Zhuravlevka terrane and the Olginsk–Taukhinsk terranes that compose the basements of the Sikhote Alin volcanic belt.     

Remote sensing parameterization of the processes of energy and water cycle over desertification areas

In order to understand the processes of land surface-atmosphere interaction over desertification area, it is indispensable to utilize of satellite remote sensing. Two scenes of Landsat TM were used to produce a set of maps of surface reflectance, MSAVI, vegetation coverage, surface temperature, net radiation, soil heat flux, sensible heat flux and latent heat flux. Statistical analysis based on these maps revealed some quantitative significant land surface characteristics. Future developments of the method are also discussed.     

Gelifluction within a solifluction lobe in the Kärkevagge valley, Swedish Lapland

Soil displacement, soil temperature, depths of thaw plane and groundwater level were continuously monitored during the period from July 1999 to June 2000 within a solifluction lobe in the Kärkevagge valley, northern Sweden. The strain–probe method was used to measure soil displacement, and we found significant soil displacements in the thawing period 2000. These displacements were the result of gelifluction. The ice content profile showed that gelifluction occurred at the same time as the thaw plane reached the layers with high ice content at shallow soil depths (0–6 and 16–25 cm deep). In contrast, gelifluction did not occur when the thaw plane reached the layers with high ice content at greater depth (46–49 cm deep). These observations indicate that thawing of ice lenses in the near–surface layer triggers gelifluction.     

Lithology and palynology of Neogene sediments on the narrow edge of the Kitakami Massif (basement rocks), northeast Japan: Significant change for depositional environments as a result of plate tectonics

Abstract A controversial stratigraphic section, the Taneichi Formation, is exposed along the Pacific Coast of northeastern Honshu, the main island of the Japanese Archipelago. Although most sediments of the formation have long been dated as late Cretaceous, the northern section of it has been assigned to (i) the Upper Cretaceous; (ii) the Paleogene; or (iii) the Neogene. In the present report, we present the data of palynological and sedimentological studies, showing that the northern section should be assigned to the Neogene. A more important point in the present study is that we invoke some basic principles of fluvial sedimentology to resolve this stratigraphic subject. The lignite layers full of Paleogene–Miocene dinoflagellate cysts and pollen assemblages drape over the boulder‐sized (>40 cm in diameter) clasts in the northern section. However, the layers totally consist of aggregates of small lignite chips, indicating that the lignites are allochthonous materials. The mega‐clasts with derived microfossils in the lignites are thought to have been deposited as Neogene fluvial (flood) sediments in the newly formed Japanese Archipelago. Prior to the Miocene, the northern Honshu was part of the Eurasian Plate, thus the boulder‐sized clasts cannot be envisaged as long river flood deposits along the continental Paleogene Pacific Coast. Instead, the mega‐clasts with the draping lignites were probably derived from nearby Miocene highlands in the newly born island arc.