U–Pb ages of granitoids around the Kofu basin: Implications for the Neogene geotectonic evolution of the South Fossa Magna region,central Japan

The Japanese archipelago underwent two arc–arc collisions during the Neogene. Southwest Honshu arc collided with the Izu‐Bonin‐Mariana arc and the northeast Honshu arc collided with the Chishima arc. The complicated geological structure of the South Fossa Magna region has been attributed to the collision between the Izu‐Bonin‐Mariana arc and the southwest Honshu arc. Understanding the geotectonic evolution of this tectonically active region is crucial for delineating the Neogene tectonics of the Japanese archipelago. Many intrusive granitoids occur around the Kofu basin, in the South Fossa Magna region. Although the igneous ages of these granitoids have been mainly estimated through biotite and hornblende K–Ar dating, here, we perform U–Pb dating of zircon to determine the igneous ages more precisely. In most cases, the secondary post‐magmatic overprint on the zircon U–Pb system was minor. Based on our results, we identify four groups of U–Pb ages: ca 15.5 Ma, ca 13 Ma, ca 10.5 Ma, and ca 4 Ma. The Tsuburai pluton belongs to the first group, and its age suggests that the granite formation within the Izu‐Bonin‐Mariana arc dates back to at least 15.5 Ma. The granitoids of the second group intruded into the boundary between the Honshu arc and the ancient Izu‐Bonin‐Mariana arc, suggesting that the arc–arc collision started by ca 13 Ma. As in the case of the Kaikomagatake pluton, the Chino pluton likely corresponds to a granodiorite formed in a rear‐arc setting in parallel with the other granodiorites of the third group. The U–Pb age of the Kogarasu pluton, which belongs to the fourth group, is the same as those of the Tanzawa tonalitic plutons. This might support a syncollisional rapid granitic magma formation in the South Fossa Magna region.     

A Historical Note on the Study of Ocean Surface Waves

The modern study of ocean surface waves started with a pioneer study by Sverdrup and Munk (1947). More than half a century has passed since then and the study of ocean surface waves has greatly advanced. The current numerical wave models, supported by many fundamental studies, enable us to compute ocean surface waves on a global scale with sufficient accuracy for practical purposes. However, physical process controlling the energy balance of ocean surface waves is still not completely understood. The present note is a rough sketch of the historical development of the study of ocean surface waves in the latter half of the twentieth century when the Oceanographic Society of Japan was founded and grew. This revised version was published online in August 2006 with corrections to the Cover Date.     

Ba, Si, U, Al, Sc, La, Th, C and C/C in a sediment core in the western subarctic Pacific as proxies of past biological production

A sediment core covering the last 145 kyrs was collected in the western subarctic Pacific (WSAP), and analyzed for Ba, U, Al, Sc, La, Yb, Th, biogenic opal (Opal) and organic carbon (C_org) as well as its isotopic ratio (δ¹³C). This study examined the change of past biological production in WSAP with multiple proxies, together with understanding the relation between Loess from the Asian continent and the biological production. The Loess content was estimated from the metal components, Al, Sc, La, Yb and Th. In this high latitude core (50°N), the Loess content was generally high during the glacial periods, but it was also high even in some interglacial periods. The excess amount of Ba relative to the detrital material composition, Ba_ex, showed the best correlation with the Vostok δD (r = 0.72, p < 0.001), indicating that the biological production was lower in the glacial periods than in the interglacials. This corroborates the pervasive correlation between Ba_ex in the polar region, WSAP and the Antarctic Sea, and Antarctic temperature, combined with previous research. This correlation might be explained by the stratification caused by cooling. In addition, the time variations of Ba_ex in WSAP were similar to those of Ba_ex in the Okhotsk Sea and of other proxies (C_org and Opal) in both the Okhotsk and the Bering Sea, indicating the spatial homogeneity of Ba_ex in WSAP including proximal marginal seas. The Opal content was more weakly correlated with the Vostok δD (r = 0.46, p < 0.001) than Ba_ex, reflecting that Opal in WSAP including proximal marginal seas was spatially heterogeneous compared to Ba_ex. While both the C_org content and U_ex, the excess amount of U relative to the detritus composition, were not positively correlated with the Vostok δD, they behaved similarly in the sediments. The positive correlation between δ¹³C and the Vostok δD (r = 0.42, p < 0.001), between δ¹³C and Ba_ex (r = 0.60, p < 0.001) and between δ¹³C and Opal (r = 0.36, p < 0.01) indicates that δ¹³C in WSAP may give some information on the phytoplankton growth rate. There was not a significant correlation between the spatially homogeneous Ba_ex in WSAP and Loess (r = − 0.16, p > 0.01), suggesting that the increase of biological production with the increase of Loess supply during the glacial periods did not occur.     

Application of a two-way nested model to the seamount problem

A two-way nested model has been constructed and applied to the idealized ocean where a uniform mean flow impinges on the isolated Gaussian-shaped seamount and produces two eddies (cold and warm) in the depths. The performance of the nested model has been evaluated subjectively and objectively. Both subjective and objective analyses confirm the traditional view that the nested model can well capture the performance of isolated eddies. Objective analysis, however, reveals some quantitatively important features of a two-way nested model. One is penetration of improved features into the coarse domain and another is the deterioration of mean flow field inside the nested area, neither of which is clear from subjective analysis. With successful application of two-way nested model to the seamount problem, we expect that such a nested model will also be applicable to other oceanic phenomena, particularly to some coastal problems whose time scale is short and where the topographic effects are dominant.     

On the relation between the growth rate of water waves and wind speed

Various wind velocitiesu _*,U _/2,U andU ₁₀ are correlated to the measured growth rate of water waves , whereu _* is the friction velocity of the wind, andU _/2,U andU ₁₀ are the wind speeds respectively at the heights /2, and 10m above sea surface (: wave length). It is shown that within a range of the dimensionless wind speed, 0.1<u _* /C<0.6, there are no appreciable differences in the correlations, whereC is the phase velocity of water waves. The present relation between andU shows qualitatively similar properties as the one obtained by Al'Zanaidi and Hui (1984); the growth rate for waves with rough surface is larger than that with smooth surface. However, our present relations give, for the both waves with different surface roughness, larger values by factors 1.71.8 than those given by Al'Zanaidi and Hui's relation.     

Seafloor electromagnetic induction studies in the Bay of Bengal

Seafloor magnetometer array experiments were conducted in the Bay of Bengal to delineate the subsurface conductivity structure in the close vicinity of the 85°E Ridge and Ninety East Ridge (NER), and also to study the upper mantle conductivity structure of the Bay of Bengal. The seafloor experiments were conducted in three phases. Array 1991 consisted of five seafloor stations across the 85°E Ridge along 14°N latitude with a land reference station at Selam (SLM). Array 1992 also consisted of five seafloor stations across 85°E Ridge along 12°N latitude. Here we used the data from Annamalainagar Magnetic Obervatory (ANN) as land reference data. Array 1995 consisted of four seafloor stations across the NER along 9°N latitude with land reference station at Tirunelveli (TIR). OBM-S4 magnetometers were used for seafloor measurements. The geomagnetic Depth Sounding (GDS) method was used to investigate the subsurface lateral conductivity contrasts. The vertical gradient sounding (VGS) method was used to deliniate the depth-resistivity structure of the oceanic crust and upper mantle. 1-D inversion of the VGS responses were conducted and obtained a 3-layer depth-resistivity model. The top layer has a resistivity of 150–500 m and a thickness of about 15–50 km. The second layer is highly resistive (2000–9000 m) followed by a very low resistive (0.1–50 m) layer at a depth of about 250–450 km. The 3-component magnetic field variations and the observed induction arrows indicated that the electromagnetic induction process in the Bay of Bengal is complex. We made an attempt to solve this problem numerically and followed two approaches, namely (1) thin-sheet modelling and (2) 3-D forward modelling. These model calculations jointly show that the observed induction arrows could be explained in terms of shallow subsurface features such as deep-sea fans of Bay of Bengal, the resistive 85°E Ridge and the sea water column above the seafloor stations. VGS and 3-D forward model responses agree fairly well and provided depth-resistivity profile as a resistive oceanic crust and upper mantle underlained by a very low resistive zone at a depth of about 250–400 km. This depth-range to the low resistive zone coincide with the seismic low velocity zone of the northeastern Indian Ocean derived from the seismic tomography. Thus we propose an electrical conductivity structure for the oceanic crust and upper mantle of the Bay of Bengal.     

船舶螺旋桨尾流场的数值分析

利用基于速度势的低阶面元法计算船舶螺旋桨的尾流场。采用计算较为简捷的关于扰动速度势的基本积分微分方程，并采用双曲面形状的面凶以消除面元间的缝隙。Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件，使桨叶面上表面的压力在随边有良好的一致性。在计算面元的影响系数时，应用了Morino导出的解析计算公式，加快了数值计算的速度。从解面元法的基本积分方程得到的偶极强度和源汇强度，直接求得尾流场的速度分布。     

Late Cenozoic volcanic activity in the Chugoku area, southwest Japan arc during back-arc basin opening and reinitiation of subduction

Abstract Temporal–spatial variations in Late Cenozoic volcanic activity in the Chugoku area, southwest Japan, have been examined based on 108 newly obtained K–Ar ages. Lava samples were collected from eight Quaternary volcanic provinces (Daisen, Hiruzen, Yokota, Daikonjima, Sambe, Ooe–Takayama, Abu and Oki) and a Tertiary volcanic cluster (Kibi Province) to cover almost all geological units in the province. Including published age data, a total of 442 Cenozoic radiometric ages are now available. Across‐arc volcanic activity in an area approximately 500 km long and 150 km wide can be examined over 26 million years. The period corresponds to syn‐ and post‐back‐arc basin opening stages of the island arc. Volcanic activity began in the central part of the rear‐arc ca 26 Ma. This was followed by arc‐wide expansion at 20 Ma by eruption at two rear‐arc centers located at the eastern and western ends. Expansion to the fore‐arc occurred between 20 and 12 Ma. This Tertiary volcanic arc was maintained until 4 Ma with predominant alkali basalt centers. The foremost‐arc zone activity ceased at 4 Ma, followed by quiescence over the whole arc between 4 and 3 Ma. Volcanic activity resumed at 3 Ma, covering the entire rear‐arc area, and continued until the present to form a Quaternary volcanic arc. Adakitic dacite first occurred at 1.7 Ma in the middle of the arc, and spread out in the center part of the Quaternary volcanic arc. Alkali basalt activities ceased in the area where adakite volcanism occurred. Fore‐arc expansion of the volcanic arc could be related to the upwelling and expansion of the asthenosphere, which caused opening of the Japan Sea. Narrowing of the volcanic zone could have been caused by progressive Philippine Sea Plate subduction. Deeper penetration could have caused melting of the slab and resulted in adakites. Volcanic history in the Late Cenozoic was probably controlled by the history of evolution of the upper mantle structure, coinciding with back‐arc basin opening and subsequent reinitiation of subduction.     

网式大地电磁阻抗张量及其影响因素分析

网式大地电磁（Network-MT,N-MT）法采用长数公里至数十公里的电话线为电极线测量电场,很难形成两条笔直且相互垂直的电极线,因此阻抗张量的计算不如大地电磁法中直接.本文依据阻抗张量的旋转规则提出了一种计算N-MT阻抗张量的简便算法.依据该算法计算了中国东北地区5个N-MT测站的阻抗张量,获得了基于阻抗张量的视电阻率、相位曲线和最佳主轴方位角分布图像,为我国东北N-MT资料的进一步处理和解释提供了基础性数据.此外,本文对比分析了朝阳测站中6条N-MT观测电极线上阻抗张量旋转值与观测值之间的差异,重点讨论了产生这种系统性偏差的各种因素,提出其主要因素可能来自“电场等效各向异性”效应,即测站附近的地壳内部存在与观测电极线尺度相比拟的横向非均匀构造,而测站各电极线沿不同方向跨越不均匀构造,此时各电极线上的电场分量不遵循同一电场矢量的分解准则,导致地表观测三角形内阻抗分量不满足统一的阻抗张量旋转规则.