Roles of Biogeochemical Processes in the Oceanic Carbon Cycle Described with a Simple Coupled Physical-Biogeochemical Model

To evaluate the contribution of biogeochemical processes to the oceanic carbon cycle and to calculate the ratio of calcium carbonate to organic carbon downward export, we have incorporated biological and alkalinity pumps in the yoked high-latitude exchange/interior diffusion-advection (YOLDA) model. The biogeochemical processes are represented by four parameters. The values of the parameters are tuned so that the model can reproduce the observed phosphate and alkalinity distributions in each oceanic region. The sensitivity of the model to the biogeochemical parameters shows that biological production rates in the euphotic zone and decomposition depths of particulate matters significantly influence horizontal and vertical distributions of biogeochemical substances. The modeled vertical fluxes of particulate organic phosphorus and calcium carbonate are converted to vertical carbon fluxes by the biological pump and the alkalinity pump, respectively. The downward carbon flux from the surface layer to the deep layer in the entire region is estimated to be 3.36 PgC/yr, which consists of 2.93 PgC/yr from the biological pump and 0.43 PgC/yr from the alkalinity pump, which is consistent with previous studies. The modeled rain ratio is higher with depth and higher in the Pacific and Indian Oceans than in the Atlantic Ocean. The global rain ratio at the surface layer is calculated to be 0.14 to 0.15. This value lies between the lower and higher ends of the previous estimates, which range widely from 0.05 to 0.25. This study indicates that the rain ratio is unlikely to be higher than 0.15, at least in the surface waters.     

Processes Causing the Temporal Changes in Si/N Ratios of Nutrient Consumptions and Export Flux during the Spring Diatom Bloom

Two processes are generally explained as causes of temporal changes in the stoichiometric silicon/nitrogen (Si/N) ratios of sinking particles and of nutrient consumption in the surface water during the spring diatom bloom: (1) physiological changes of diatom under the stress of photosynthesis of diatom and (2) differences of regeneration between silicon and nitrogen. We investigated which process plays an important role in these changes using a one-dimensional ecosystem model that explicitly represents diatom and the other non-silicious phytoplankton. The model was applied to station A7 (41°30′ N, 145°30′ E) in the western North Pacific, where diatom regularly blooms in spring. Model simulations show that the Si/N ratios of the flux exported by the sinking particles at 100 m depth and of nutrient consumptions in the upper 100 m surface water have their maxima at the end of the spring diatom bloom, the values and timings of which are significantly different from each other. Analyses of the model results show that the differences of regeneration between silicon and nitrogen mainly cause the temporal changes of the Si/N ratios. On the other hand, the physiological changes of diatoms under stress can hardly cause these temporal changes, because the effect of the change in the diatom's uptake ratio of silicon to nitrogen is cancelled by that in its sinking rate.     

An Ecosystem Model Coupled with Nitrogen-Silicon-Carbon Cycles Applied to Station A7 in the Northwestern Pacific

A model based on that of Kishi et al. (2001) has been extended to 15 compartments including silicon and carbon cycles. This model was applied to Station A7 off Hokkaido, Japan, in the Northwestern Pacific. The model successfully simulated the observations of: 1. a spring bloom of diatoms; 2. large seasonal variations of nitrate and silicate concentrations in the surface water; and 3. large inter-annual variations in chlorophyll-a. It also reproduced the observed features of the seasonal variations of carbon dioxide partial pressure (pCO₂)—a peak in pCO₂ in winter resulting from deep winter convection, a rapid decrease in pCO₂ as a result of the spring bloom, and an almost constant pCO₂ from summer through fall (when the effect of increasing temperature cancels the effect of biological production). A comparison of cases with and without silicate limitation shows that including silicate limitation in the model results in: 1. decreased production by diatoms during summer; and 2. a transition in the dominant phytoplankton species, from diatoms to other species that do not take up silicate. Both of these phenomena are observed at Station A7, and our results support the hypothesis that they are caused by silicate limitation of diatom growth. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Comparison of Simulated Particle Fluxes Using NEMURO and Other Ecosystem Models in the Western North Pacific

JGOFS has revealed the importance of marine biological activity to the global carbon cycle. Ecological models are valuable tools for improving our understanding of biogeochemical cycles. Through a series of workshops, the North Pacific Marine Science Organization (PICES) developed NEMURO (North Pacific Ecosystem Model Understanding Regional Oceanography) a model, specifically designed to simulate the lower trophic ecosystem in the North Pacific Ocean. Its ability to simulate vertical fluxes generated by biological activities has not yet been validated. Here compare NEMURO with several other lower trophic level models of the northern North Pacific. The different ecosystem models are each embedded in a common three-dimensional physical model, and the simulated vertical flux of POM and the biomass of phytoplankton are compared. The models compared are: (1) NEMURO, (2) the Kishi and Nakata Model (Kishi et al., 1981), (3) KKYS (Kawamiya et al., 1995, 2000a, 2000b), and (4) the Denman model (Denman and Peña, 2002). With simple NPZD models, it is difficult to describe the production of POM (Particulate Organic Matter) and hence the simulations of vertical flux are poor. However, if the parameters are properly defined, the primary production can be well reproduced, even though none of models we used here includes iron limitation effects. On the whole, NEMURO gave a satisfactory simulation of the vertical flux of POM in the northern North Pacific.     

Effects of presence of a circumpolar region on buoyancy-driven circulation

Effects of the presence of a circumpolar region on buoyancy-driven circulation are investigated by using an idealized numerical ocean model. Comparison of circulation and meridional density (heat) transport is made between a closed ocean and an ocean with a cyclic gap near its southern boundary. The presence of the circumpolar region leads to disconnection of the meridional overturning across the circumpolar region. And the circumpolar eastward flow reaches the bottom of the ocean. It is essential for this that the pycnocline is deeper than the bottom of the gap. Since the amount of the mass transported northward must return southward at the levels deeper than the bottom of the cyclic gap, the weak stratification, hence weak vertical geostrophic shear, at the deeper levels leads to inactive communication across the circumpolar region. Meridional heat transport across the circumpolar region is made mainly by horizontal diffusion for the ocean with the cyclic gap, while the contribution of the advection is dominant for the closed ocean. Sensitivity of meridional heat transport to change in horizontal diffusivity is studied. The meridional heat transport for the ocean with the cyclic gap is more sensitive than for the closed ocean. The change in heat transport occurs not only in the circumpolar region but also in the rest of the ocean. It is suggested that subgrid scale phenomena, especially mesoscale eddies, in the circumpolar region controls the whole ocean to a great extent.     

Fluid–Sediment Interaction in a Marine Shallow-Water Hydrothermal System in the Wakamiko Submarine Crater, South Kyushu, Japan

The Wakamiko submarine crater is a small depression located in Kagoshima Bay, southwest Japan. Marine shallow‐water hydrothermal activity associated with fumarolic gas emissions at the crater sea floor (water depth 200 m) is considered to be related with magmatic activity of the Aira Caldera. During the NT05‐13 dive expedition conducted in August 2005 using remotely operated vehicle Hyper‐Dolphine (Japan Agency for Marine‐Earth Science and Technology), an active shimmering site was discovered (tentatively named the North site) at approximately 1 km from the previously known site (tentatively named the South site). Surface sediment (up to 30 cm) was cored from six localities including these active sites, and the alteration minerals and pore fluid chemistry were studied. The pore fluids of these sites showed a drastic change in chemical profile from that of seawater, even at 30 cm below the surface, which is attributed to mixing of the ascending hydrothermal component and seawater. The hydrothermal component of the North site is estimated to be derived from a hydrothermal aquifer at 230°C based on the hydrothermal end‐member composition. Occurrence of illite/smectite interstratified minerals in the North site sediment is attributed to in situ fluid–sediment interaction at a temperature around 150°C, which is in accordance with the pore fluid chemistry. In contrast, montmorillonite was identified as the dominant alteration mineral in the South site sediment. Together with the significant low potassium concentration of the hydrothermal end‐member, the abundant occurrence of low‐temperature alteration mineral suggests that the hydrothermal aquifer in the South site is not as high as 200°C. Moreover, the montmorillonite is likely to be unstable with the present pore fluid chemistry at the measured temperature (117°C). This disagreement implies unstable hydrothermal activity at the South site, in contrast to the equilibrium between the pore fluid and alteration minerals in the North site sediment. This difference may reflect the thermal and/or hydrological structure of the Wakamiko Crater hydrothermal system.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

Performance of Mortar and Chemical Grout Injection into Surrounding Soil When Slurry Pipe-jacking Method is Used

Small-diameter shallow tunnels are often being built by using the slurry pipe-jacking method. This system involves the pushing or thrusting of a drivage machine and concrete pipes into the ground. Chemical grout injection into the surrounding soil around the tunnel is carried out after the drivage and pushing processes are finished. The purpose of the chemical grout injection is to maintain permanent stability of the surrounding soil. However, the behavior of the chemical grouting material in the surrounding soil around the tunnel and the amount of optimum injection is not clearly understood. From these points of view, this paper discusses the performance of the chemical grouting material, when it is injected into the surrounding soil around the tunnel, by means of 2-D Eulerian–Lagrangian seepage analysis. Moreover, the effectiveness of the chemical grout injection was evaluated by using the non-linear finite element method. This investigation show when the range of the grouted zone is designed; it is necessary that the relationship between Young’s modulus of the soil/grouted zone and the confining stress be taken into consideration in order to establish effective, economical and safe chemical grout injection system. Understanding the performance of the seepage/dispersion behavior of the chemical grout and the characteristics of soil/ grouted zone is also important.