Phytoplankton dynamics in Lake Biwa during the 20th century: complex responses to climate variation and changes in nutrient status

We examined algal remains and fossil pigments in ²¹⁰Pb-dated sediment cores from Lake Biwa to explore historical changes in the phytoplankton community of the lake over the past 100 years and to identify environmental factors that caused those changes. Fluxes of fossil pigments and algal remains were very low before the 1960s, but increased through the 1960s and 1970s, indicating that the lake had eutrophied in the 20 years since 1960. After 1980, however, fluxes of all fossil pigments and algal remains decreased or stabilized. Redundancy analysis with meteorological and limnological variables explained more than 70% of the variation of these fluxes and showed that the decrease in fluxes of most algal taxa that occurred in the 1980s was related to changes in meteorological variables such as wind velocity, rather than changes in the lake’s trophic state. Sedimentary records of algal remains also revealed that Aulacoseira nipponica, an endemic diatom species that grows in winter, decreased dramatically after 1980, while Fragilaria crotonensis, a cosmopolitan spring diatom species, became dominant. Replacement of one dominant diatom species by another could not be explained simply by changes in the lake trophic state, but was reasonably strongly related with an increase in winter water temperature. These results suggest that the phytoplankton community in Lake Biwa was influenced by changes in local environmental conditions (nutrient loading) through the 1960s and 1970s, but more so by regional (meteorological) and global (climate warming) factors since 1980.     

Effects of plant roots on the soil erosion rate under simulated rainfall with high kinetic energy

ABSTRACT

This study examined the effects of herbaceous plant roots on interrill erosion using two herbaceous species: clover (Trifolium repens) and oats (Avena sativa). We developed a simple rainfall simulator with relatively high normalized kinetic energy (KE; 23.2 J m^?2 mm^?1). Under simulated rainfall, we measured eroded soil for 42 boxes with various amounts of aboveground and belowground biomass. Aboveground vegetation had a significant effect on the soil erosion rate (SER). We found a clear negative relationship between the percent vegetation cover (c) and the SER. In contrast, plant roots showed no effects on the SER. The SER was not significantly different between the boxes with and without plant roots under similar c conditions. Thus, plant roots could have less of an effect on the SER under higher KE conditions.

Editor M.C. Acreman Associate editor N. Verhoest     

A precise bathymetric map of the world’s deepest seafloor,Challenger Deep in the Mariana Trench

Data from three bathymetric surveys by R/V Kairei using a 12-kHz multibeam echosounder and differential GPS were used to create an improved topographic model of the Challenger Deep in the southwestern part of the Mariana Trench, which is known as the deepest seafloor in the world. The strike of most of the elongated structures related to plate bending accompanied by subduction of the Pacific plate is N70°E and is not parallel to the trench axis. The bending-related structures were formed by reactivation of seafloor spreading fabric. Challenger Deep consists of three en echelon depressions along the trench axis, each of which is 6–10 km long, about 2 km wide, and deeper than 10,850 m. The eastern depression is the deepest, with a depth of 10,920 ± 5 m.     

3-D internal architecture of methane hydrate-bearing turbidite channels in the eastern Nankai Trough, Japan

The reservoir architecture of methane hydrate (MH) bearing turbidite channels in the eastern Nankai Trough, offshore Japan is evaluated using a combination of 3-D seismic and well data. On the 3-D seismic section, the MH-bearing turbidite channels correspond to complex patterns of strong seismic reflectors, which show the 3-D internal architecture of the channel complex. A seismic-sequence stratigraphic analysis reveals that the channel complex can be roughly classified into three different stages of depositional sequence (upper, middle, and lower). Each depositional sequence results in a different depositional system that primarily controls the reservoir architecture of the turbidite channels. To construct a 3-D facies model, the stacking patterns of the turbidite channels are interpreted, and the reservoir heterogeneities of MH-bearing sediments are discussed. The identified channels at the upper sequence around the β1 well exhibit low-sinuosity channels consisting of various channel widths that range from tens to several hundreds of meters. Paleo-current flow directions of the turbidite channels are typically oriented along the north-northeast-to-south-southwest direction. High-amplitude patterns were identified above the channels along the north-to-south and north-northeast-to-south-southeast directions. These roughly coincide with the paleo-current flow of the turbidite channels. An interval velocity using high-density velocity analysis shows that velocity anomalies (>2000 m/s) are found on the northeastern side of the turbidite channels. The depositional stage of the northeastern side of the turbidite channels exhibits slightly older sediment stages than the depositional stages of the remaining channels. Hence, the velocity anomalies of the northeastern side of the channels are related to the different stages of sediment supply, and this may lead to the different reservoir architectures of the turbidite channels.     

Formation process of pyrite polyframboid based on the heavy-metal analysis by micro-PIXE

Pyrite framboid in sedimentary rocks could be concerned with arsenic contamination in groundwater of acidic environment and has been studied for the formation process of its unique morphology. However, not much has been discussed about the formation process based on the heavy-metal distribution in pyrite framboids and their aggregates. To reveal the formation process of pyrite framboids and their aggregates, mudstone from the Late Cretaceous Hakobuchi Group, central Hokkaido, Japan, are investigated for mode, petrographical, mineralogical, and micro-PIXE (particle-induced X-ray emission) analysis in this study. Spherical and sub-spherical pyrite framboids observed in polyframboids can be divided into two types based on the diameter of framboids (D) and the average diameter of microcrystals (d) within framboids: type 1, ranging from 2 to 9 μm and from 0.4 to 0.9 μm, respectively, and type 2, ranging from 8 to 50 μm and from 0.5 to 1.8 μm, respectively. Based on the quantitative traverse analyses and 2-D elemental mapping results by micro-PIXE, heterogeneities in the concentrations of heavy metals within the analyzed aggregations of pyrite are exhibited. On the basis of the As and Pb zoning patterns, the time range and chronological stages of pyrite-aggregation growth are revealed.     

A stability analysis of a conduit flow model for lava dome eruptions

Periodic variations in magma discharge rate and ground deformation have been commonly observed during lava dome eruptions. We performed a stability analysis of a conduit flow model by Barmin et al. [Barmin, A., Melnik, O., Sparks, R.S.J., 2002. Periodic behavior in lava dome eruptions. Earth and Planetary Science Letters 199 (1-2), 173–184], in which the periodic variations in magma flow rate and chamber pressure are reproduced as a result of the temporal and spatial changes of the magma viscosity controlled by the kinetics of crystallization. The model is reduced to a dynamical system where the time derivatives of the magma flow rate (dQ/dt) and the chamber pressure (dP/dt) are functions of Q and P evaluated at a shifted time t − t_?. Here, the time delay t_? represents the time for the viscosity of fluid particle to increase in a conduit. The dynamical system with time delay is approximated by a simple two-dimensional dynamical system of Q and P where t_? is given as a parameter. The results of our linear stability analyses for these dynamical systems indicate that the transition from steady to periodic flow depends on nonlinearities in the steady state relation between Q and P. The steady state relation shows a sigmoidal curve in Q − P phase plane; its slope has negative values at intermediate flow rates. The steady state solutions become unstable, and hence P and Q oscillate periodically, when the negative slope of the steady state relation ([dP/dQ]_S) exceeds a critical value; that is [dP/dQ]_S < − t_?γ/(2V_ch), where V_ch is the chamber volume and γ is an elastic constant which is related to the rigidity of chamber wall. We also found that the period and the pattern of oscillation of the conduit flow primarily depend on a quantity defined by LV_ch/r⁴, where L is the conduit length and r is the conduit radius.     

Gravity and S-wave modelling across the Jan Mayen Ridge,North Atlantic; implications for crustal lithology

The horizontal components from fourteen Ocean Bottom Seismometers deployed along four profiles focused along the western margin of the Jan Mayen microcontinent, North Atlantic, have been modelled with regard to S-waves, based on P-wave models obtained earlier. The seismic models have furthermore been constrained by 2D gravity modelling. High V _p/V _s-ratios (2.3–7.9) within the Cenozoic sedimentary section are attributed to significant porosities, whereas V _p/V _s-ratios in the order of 1.9–2.2 for the Mesozoic and Paleozoic sedimentary rocks indicate shale-dominated lithology throughout the area. The eastern side of the Jan Mayen Ridge is interpreted as a passive, volcanic margin, based on relatively high crustal V _p/V _s-ratios (1.9), whereas lower V _p/V _s-ratios (1.75–1.8) suggest the presence of intermediate composition crust and non-volcanic margin on the western side of the ridge. In the westernmost part of the Jan Mayen Basin, slightly increased upper mantle V _p/V _s-ratios may indicate some degree of serpentization of upper mantle peridotites.     

Application of a 3-D chemical fate prediction model (FATE3D) to predict dioxin concentrations in the Tokyo Bay

A 3-D chemical fate prediction model (FATE3D) was applied to predict the dioxin concentrations in the seawater of Tokyo Bay, Japan. The simulations were carried out for a period of one year (from September 2002 to August 2003). Parameters such as meteorological data, flow field conditions, concentrations and sinking rates of organic particulate matter, initial and boundary conditions, and loading fluxes and physico-chemical properties of dioxins were used as the model inputs.The simulation results compared favorably with the field measurements of dioxin concentrations in the bay for both the particulate and dissolved phases, indicating the validity and predictive capability of the model. Furthermore, the differences in the seasonal cycles and distributions between the particulate- and dissolved-phase dioxins in the bay were estimated from the simulation results.However, the particulate-phase dioxin concentrations in the bottom layers (+1 m from the bottom) were underestimated, probably because the resuspension process was not taken into account in the model. The improvement of the model's predictive capability, including the resuspension process, shall be the focus of our next study.