Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.     

A hypothesis of the Senoumi submarine megaslide in Suruga Bay in Japan—based on the undrained dynamic-loading ring shear tests and computer simulation

The distinctive bathymetric feature exists in the Suruga Bay, Japan. It has been called as Senoumi (Stone flower sea) from old times. Senoumi is a 30?km wide and 20?km long concave feature. Its origin has not been explained yet; however, the feature might be a combined consequence of intensive tectonic activity in the plate border, landslides, and a submarine flow coming from the Oi River. If the Senoumi was caused by a landslide, the latter would be larger than any on-land landslide in Japan. The downshelf “exit” from this feature is much narrower than its central part. This is not usual shape of landslides, but it is similar to the liquefied landslides such as those in quick clays which mobilize great strength reduction after failure. To study Senoumi as a landslide, the shear behaviors of the following three soil samples were investigated by the cyclic and seismic undrained stress control ring shear tests. One sample is volcanic ash taken from the base of landslide deposits (mass transport deposits), from 130 to 190?m deep layer below the submarine floor which was drilled and cored by the Integrated Ocean Drilling Program Expedition 333. Another two samples are the Neogene silty–sand and silt taken from the Omaezaki hill adjacent to the Senoumi, because the shear zone might have been formed in Neogene layers extending from on-land to the continental shelf. The largest strength reduction from peak to steady-state shear resistance in the undrained cyclic loading test was found in volcanic ash. The strength reduction in Neogene silty–sand was smaller than volcanic ash, while the Neogene silt mobilized the least post-failure strength reduction. An integrated model simulating the initiation and motion of earthquake-induced rapid landslides (landslide simulation (LS)-RAPID, Sassa et al. Landslides 7–3:219–236, 2010) was applied to this study. The steady-state shear resistance and other geotechnical parameters measured by the undrained ring shear tests and the greatest strong motion record in the 2011 off-the-Pacific Coast of Tohoku earthquake (M _w 9.0), also known as “2011 Tohoku Earthquake” at the observation point MYG004 (2,933?gal) were input to this model. As the result, it was found that landslides would be triggered by 0.30–1.0 times of MYG004 in volcanic ash, 0.4–1.0 times of MYG004 in Neogene silty–sand and Neogene silt, though the depth and area of triggered landslides were different in soils and intensity of shaking. Feature, created by LS-RAPID using the parameters of volcanic ash, was most similar to the Senoumi in depth and extent. The result obtained from this study includes a hypothesis to be proved, but presents the strong need to investigate the risk of the large-scale submarine landslides which could enhance tsunami wave and possibly enlarge the submarine landslide retrogressively into the adjacent coastal plain by the upcoming mega earthquake in the Nankai Trough.     

Geochemical Characteristics,Origin, and Evolution of Ore‐Forming Fluids of the Khut Copper Skarn Deposit,West of Yazd in Central Iran

The Khut copper skarn deposit is located at about 50 km northwest of Taft City in Yazd province in the middle part of the Urumieh‐Dokhtar magmatic arc. Intrusion of granitoid of Oligocene–Miocene age into carbonate rocks of the Triassic Nayband Formation led to the formation of marble and a calcic skarn. The marble contains high grade Cu mineralization that occurs mainly as open space filling and replacement. Cu‐rich sulfide samples from the mineralized marble are also anomalous in Au, Zn, and Pb. In contrast, the calcic skarn is only weakly anomalous in Cu and W. The calcic skarn is divided into garnet skarn and garnet–pyroxene skarn zones. Paragenetic relationships and microthermometric data from fluid inclusions in garnet and calcite indicate that the compositional evolution of skarn minerals occurred in three main stages as follows. (i) The early prograde stage, which is characterized by Mg‐rich hedenbergite (Hd_53.7Di_42.3–Hd_86.1Di_9.5) with Al‐bearing andradite (69.8–99.5 mol% andradite). The temperature in the early prograde skarn varies from 400 to 500°C at 500 bar. (ii) The late prograde stage is manifested by almost pure andradite (96.2–98.4 mol% andradite). Based on the fluid inclusion data from garnet, fluid temperature and salinity in this stage is estimated to vary from 267 to 361°C and from 10.1 to 21.1 wt% NaCl equivalent, respectively. Pyrrhotite precipitation started during this stage. (iii) The retrograde stage occurs in an exoskarn, which consists of an assemblage of ferro‐actinolite, quartz, calcite, epidote, chlorite, sphalerite, pyrite, and chalcopyrite that partially replaces earlier mineral assemblages under hydrostatic conditions during fracturing of the early skarn. Fluids in calcite yielded lower temperatures (T < 260°C) and fluid salinity declined to ～8 wt% NaCl equivalent. The last stage mineralization in the deposit is supergene weathering/alteration represented by the formation of iron hydroxide, Cu‐carbonate, clay minerals, and calcite. Sulfur isotope data of chalcopyrite (δ³⁴S of +1.4 to +5.2‰) show an igneous sulfur source. Mineralogy and mineral compositions of the prograde assemblage of the Khut skarn are consistent with deposition under intermediately oxidized and slightly lower fS₂ conditions at shallow crustal levels compared with those of other typical Fe‐bearing Cu–Au skarn systems.     

Mineralogy and geochemistry of granitoids from Kinnaur region,Himachal Higher Himalaya,India: Implication on the nature of felsic magmatism in the collision tectonics

Felsic magmatism in the southern part of Himachal Higher Himalaya is constituted by Neoproterozoic granite gneiss (GGn), Early Palaeozoic granitoids (EPG) and Tertiary tourmaline-bearing leucogranite (TLg). Magnetic susceptibility values (<3 ×10^?3 SI), molar Al₂ O ₃/(CaO + Na₂ O + K ₂O) (≥1.1), mineral assemblage (bt–ms–pl–kf–qtz ± tur ± ap), and the presence of normative corundum relate these granitoids to peraluminous S-type, ilmenite series (reduced type) granites formed in a syncollisional tectonic setting. Plagioclase from GGn (An₁₀–An₃₁) and EPG (An₁₅–An₃₃) represents oligoclase to andesine and TLg (An₂–An₁₅) represents albite to oligoclase, whereas compositional ranges of K-feldspar are more-or-less similar (Or₈₈ to Or₉₅ in GGn, Or₈₆ to Or₉₇ in EPG and Or₈₇ to Or₉₄ in TLg). Biotites in GGn (Mg/Mg + Fe^t= 0.34–0.45), EPG (Mg/Mg + Fe^t= 0.27–0.47), and TLg (Mg/Mg + Fe^t= 0.25–0.30) are ferribiotites enriched in siderophyllite, which stabilised between FMQ and HM buffers and are characterised by dominant 3Fe\(\rightleftharpoons \)2Al, 3Mg\(\rightleftharpoons \)2Al substitutions typical of peraluminous (S-type), reducing felsic melts. Muscovite in GGn (Mg/Mg + Fe^t=0.58–0.66), EPG (Mg/Mg + Fe^t=0.31?0.59), and TLg (Mg/Mg + Fe^t=0.29–0.42) represent celadonite and paragonite solid solutions, and the tourmaline from EPG and TLg belongs to the schorl-elbaite series, which are characteristics of peraluminous, Li-poor, biotite-tourmaline granites. Geochemical features reveal that the GGn and EPG precursor melts were most likely derived from melting of biotite-rich metapelite and metagraywacke sources, whereas TLg melt appears to have formed from biotite-muscovite rich metapelite and metagraywacke sources. Major and trace elements modelling suggest that the GGn, EPG and TLg parental melts have experienced low degrees (～13, ～17 and ～13%, respectively) of kf–pl–bt fractionation, respectively, subsequent to partial melting. The GGn and EPG melts are the results of a pre-Himalayan, syn-collisional Pan-African felsic magmatic event, whereas the TLg is a magmatic product of Himalayan collision tectonics.     

Dynamic evolution of nutrient discharge under stormflow and baseflow conditions in a coastal agricultural watershed in Ishigaki Island,Okinawa, Japan

Excessive terrestrial nutrient loadings adversely impact coral reefs by primarily enhancing growth of macroalgae, potentially leading to a phase‐shift phenomenon. Hydrological processes and other spatial and temporal factors affecting nutrient discharge must be examined to be able to formulate effective measures for reducing nutrient export to adjacent reefs. During storm events and baseflow periods, water samples were obtained from the tropical Todoroki River, which drains an intensively agricultural watershed into Shiraho coral reef. In situ nutrient analyzers were deployed for 6 months to hourly measure dissolved nutrient (NO₃⁻‐N and PO₄³⁻‐P) concentrations. Total phosphorus (TP) and suspended solid concentration (TSS) were increased by higher rainfall intensity (r = 0·94, p < 0·01) and river discharge Q (r = 0·88, p < 0·01). In contrast, NO₃⁻‐N concentration tends to decrease drastically (e.g. from 3 to 1 mg l⁻¹) during flood events. When base flow starts to dominate afterwards, NO₃⁻‐N manifested an increasing trend, but decreases when baseflow discharge becomes low. This counter‐clockwise hysteresis for NO₃⁻‐N highlights the significant influence of groundwater discharge. N delivery can therefore be considered a persistent process compared to sediment and P discharge, which are highly episodic in nature. Based on GIS analysis, nutrient concentration along the Todoroki River was largely affected by the percentage of sugarcane/bare areas and bedrock type. The spatial distribution of N concentration in the river reflects the considerable influence of subsurface geology—higher N levels in limestone‐dominated areas. P concentrations were directly related to the total length of artificial drainage, which enhances sediment transport. The use of high‐resolution monitoring data coupled with GIS‐based spatial analysis therefore enabled the clarification of control factors and the difference in the spatio‐temporal discharge characteristics between N and P. Thus, although erosion‐reduction schemes would reduce P discharge, other approaches (e.g. minimize fertilizer) are needed to reduce N discharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Natural and anthropogenic aquatic environmental changes reconstructed by paleolimnological analyses in Lake Kitaura,central Japan

In order to assess the recent anthropogenic environmental changes in Lake Kitaura, central Japan, changes during the past few centuries were reconstructed from results of radiometric and tephrochlonological age determination, magnetic susceptibility measurements, total organic carbon analyses, total nitrogen analyses and fossil diatom analyses on a sediment core from the lake. A total of six major and sub-zones are recognized according to the diatom fossil assemblages, and we discuss aquatic environmental change in Lake Kitaura mainly based on these diatom assemblage change. Zone Ia and Zone Ib (older than AD 1707) are marine to brackish. In Zone IIa (AD␣1707–AD 1836), most of the brackish diatoms disappeared, and were replaced by freshwater species indicating a decrease in salinity. We interpret the salinity decrease in Zone I–IIa as a sea-level fall during the Little Ice Age. The salinity of the lake decreased to near freshwater conditions in Zone IIb (AD 1836–AD 1970), which could arise from alteration in River Tone or development of a sandspit in the mouth of River Tone in addition to sea-level change. In Zone IIIa (AD 1970–AD 1987), the diatom assemblage indicates a freshwater environment, and sedimentation rates increase rapidly. These changes reflect sedimentary environment change and an ecosystem transition due to the construction of the tide gate. In Zone IIIb (AD 1987–AD 2002), the diatom flux (valves cm⁻² y⁻¹) increased and species composition changed. The changes in Zone IIIb show a good agreement with limnological monitoring data gathered from the lake. These paleolimnological data suggest that the recent human-induced changes of the aquatic environment of the lake after the 1970s exceed rates during the period concerned in this study.     

Simple formulations of ground impedance functions for rigid surface foundations

A differential equation is formulated for the dynamic response of ground medium by using a simplified ground model. Applying Galerkin's procedure for weighted residual, this equation leads to a governing equation only at the ground surface. The equation indicates that the ground surface behavior can be computed even further by a simplified model. By solving the governing equation for the boundary conditions along the surface, expressions in simple closed forms are developed for the dynamic response analysis of a massless rigid foundation that rests on the ground surface. Despite their significant simplicity, the developed expressions compute the values very close to those computed by far more complex rigorous solutions. They are found to be capable of capturing the important characteristics of the dynamic ground behavior well.     

Assessment of Debris-Flow Hazards of Alluvial Fans

A strategy is presented for the assessment ofdebris-flow hazards on alluvial fans on the basis ofa case study carried out on the southern foot of astratovolcano named Mt. Yatsugatake. Transformation ofcommercial forests into a golf course was planned ata corner of the Kikkakezawa fan. The case studyinvolves an assessment of hazards due to debris flowsof different magnitudes and their recurrenceintervals. The environmental conditions for therecurrence of these debris flows are discussed as wellas the extent of the areas affected by them. In orderto generalize the case study, concepts of hazardpotential, hazard and risk for debris flows onalluvial fans are established. Accordingly, the hazardpotential is the possible hazards at any location onan alluvial fan in an indefinitely long time period,which can be assessed on the basis of hydrological andother geoscientific investigations. Hazards associatedwith a particular land use can be evaluated on thebasis of the hazard potential considering the locationand the time period associated with the land use. Riskcan be further assessed considering the life styles inthis land and the social conditions.     

Coronal Mass Ejections Observed at the Total Solar Eclipse on 13 November 2012

White-light observations of the total solar eclipse on 13 November 2012 were made at two sites, where the totality occurred 35 min apart. The structure of the corona from the solar limb to a couple of solar radii was observed with a wide dynamic range and a high signal-to-noise ratio. An ongoing coronal mass ejection (CME) and a pre-CME loop structure just before the eruption were observed in the height range between 1?–?2 R_⊙. The source region of CMEs was revealed to be in this height range, where the material and the magnetic field of CMEs were located before the eruption. This height range includes the gap between the extreme ultraviolet observations of the low corona and the spaceborne white-light observations of the high corona, but the eclipse observation shows that this height range is essential for the study of CME initiation. The eclipse observation is basically just a snapshot of CMEs, but it indicates the importance of a continuous coverage of CME observations in this height range in the future.