Determination of the elastic modulus set of foliated rocks from ultrasonic velocity measurements

Ultrasonic measurements of compressional and shear wave velocities under hydrostatic pressure up to 70 MPa were carried out on cylindrical specimens cored across and along the foliation planes. Our measurements revealed that the foliation of the metamorphic rocks induces a clear velocity anisotropy between two orthogonal directions; faster along the foliation plane and slower across the plane in most rock types. All velocity components monotonically increase with the confining pressure, probably due to the closure of microcracks distributed in rock specimens. We determined the complete set of dynamic moduli of foliated metamorphic rocks with two assumptions; transverse isotropy due to the foliation and ellipsoidal seismic energy propagation from a point source. The calculated elastic moduli referring to different directions could be valuable for the design of various engineering structures in planar textured rock mass.     

Changes of precipitation extremes over South Korea projected by the 5 RCMs under RCP scenarios

The change of extreme precipitation is assessed with the HadGEM2-AO - 5 Regional Climate Models (RCMs) chain, which is a national downscaling project undertaken cooperatively by several South Korean institutes aimed at producing regional climate change projection with fine resolution (12.5 km) around the Korean Peninsula. The downscaling domain, resolution and lateral boundary conditions are held the same among the 5 RCMs to minimize the uncertainties from model configuration. Climatological changes reveal a statistically significant increase in the mid-21st century (2046- 2070; Fut1) and the late-21st century (2076-2100; Fut2) precipitation properties related to extreme precipitation, such as precipitation intensity and average of upper 5 percentile daily precipitation, with respect to the reference period (1981-2005). Changes depending on the intensity categories also present a clear trend of decreasing light rain and increasing heavy rain. In accordance with these results, the change of 1-in-50 year maximum precipitation intensity over South Korea is estimated by the GEV method. The result suggests that the 50-year return value (RV50) will change from -32.69% to 72.7% and from -31.6% to 96.32% in Fut1 and from -31.97% to 86.25% and from -19.45% to 134.88% in Fut2 under representative concentration pathway (RCP) 4.5 and 8.5 scenarios, respectively, at the 90% confidence level. This study suggests that multi-RCMs can be used to reduce uncertainties and assess the future change of extreme precipitation more reliably. Moreover, future projection of the regional climate change contains uncertainties evoked from not only driving GCM but also RCM. Therefore, multi-GCM and multi-RCM studies are expected to provide more robust projection.     

Some statistical characteristics of large deepwater waves around the Korean Peninsula

The relationship between significant wave height and period, the variability of significant wave period, the spectral peak enhancement factor, and the directional spreading parameter of large deepwater waves around the Korean Peninsula have been investigated using various sources of wave measurement and hindcasting data. For very large waves comparable to design waves, it is recommended to use the average value of the empirical formulas proposed by Shore Protection Manual in 1977 and by Goda in 2003 for the relationship between significant wave height and period. The standard deviation of significant wave periods non-dimensionalized with respect to the mean value for a certain significant wave height varies between 0.04 and 0.21 with a typical value of 0.1 depending upon different regions and different ranges of significant wave heights. The probability density function of the peak enhancement factor is expressed as a lognormal distribution, with its mean value of 2.14, which is somewhat smaller than the value in the North Sea. For relatively large waves, the probability density function of the directional spreading parameter at peak frequency is also expressed as a lognormal distribution.     

Generation of random waves in time-dependent extended mild-slope equations using a source function method

We develop techniques of numerical wave generation in the time-dependent extended mild-slope equations of Suh et al. [1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Engineering 32, 91–117] and Lee et al. [2003. Extended mild-slope equation for random waves. Coastal Engineering 48, 277–287] for random waves using a source function method. Numerical results for both regular and irregular waves in one and two horizontal dimensions show that the wave heights and the frequency spectra are properly reproduced. The waves that pass through the wave generation region do not cause any numerical disturbances, showing usefulness of the source function method in avoiding re-reflection problems at the offshore boundary.     

Particle shape effect on thermal conductivity and shear wave velocity in sands

This study presents the correlations between quantified shape parameters and geotechnical properties for nine sand specimens. Four shape parameters, sphericity, convexity, elongation and slenderness, were quantified with two-dimensional microscopic images with the aid of image processing techniques. An instrumented oedometer cell is used to measure compressibility, thermal conductivity and shear wave velocity during loading, unloading and reloading stages. As the particle shape inherently determines the initial loose packing condition, initial void ratio and shape parameters are well correlated with compressibility. The applied stress in soils increases the interparticle contact area and contact quality; round particles tend to achieve higher thermal conductivity and shear wave velocity during stress-induced volume change. Multiple linear regression is implemented to capture the relative contributions of involved variables, revealing that the thermal evolution is governed by the initial packing density and particle shape. The experimental observations underscore the predominant effect that particle shape has on the geomechanical and physical properties upon stress-induced soil behavior.     

关于黄海深部构造的地球物理认识

利用黄海海区重磁观测数据,以最新的地震层析成像和浅层反射地震探测结果作为约束,并利用小波分解、纹理特征图像处理等手段,对研究区进行了地球物理场特征分析、岩石物性总结、断裂信息提取和磁性基底埋深反演计算,同时对南黄海海域三条剖面进行了广义逆重磁数据拟合反演和地震P波速度成像.根据地球物理数据的各种处理结果,提出在南黄海西部存在一条串珠状地球物理线性构造带,并将这条NNW向断续延伸的构造带称之为南黄海西缘断裂带.该断裂带延伸长,断裂两侧前新生代地层差异较大,有可能是造成陆区和南黄海南部盆地区前新生代油气远景差异的原因之一.文中得到南黄海磁性基底埋藏分布具有“中间浅四周深”的分布特征,与地震层析成像结果相吻合.磁性基底的局部起伏和局部构造的边界断层共同控制了前新生代残留盆地的格架和残余厚度分布,反映出“区域控制局部,深层约束浅层”的规律.     

利用卫星长波长地球物理位场对东亚地区深部结构的研究

有研究表明从中国北部经黄海、朝鲜半岛、日本海至日本本州岛存在一个大型断裂带,被称为“北纬40°断裂带”.为了验证该断裂带的存在性,本文通过GRACE和CHAMP卫星的位场数据中的大尺度区域重力和磁力异常来研究该区的大型构造特征.区域性研究范围为24°N~56°N,90°E~150°E,详细研究范围为32°N~42°N,122°E~132°E.2008年Taylor等利用CHAMP卫星2005年6月至12月收集的轨道数据,从测量数据中消除地核和外部产生磁场效应,经化磁极之后得到磁异常图,从中来追踪构造特征,并根据地质和地球物理资料建立了数学模型来进行解释.本文在此基础上,增加了GRACE卫星2003年10月份测量的重力异常场数据.实测差分位场数据利用高斯-勒让德求积法转换为垂向重力异常值,通过减去最新地球重力场模型（EGM96）的场值来去除波长超过1100 km的长波长异常,所得的重力异常数据和磁场数据进行比较.为了能从中得到研究区的共同构造特征,建立一个球面棱柱模型进行正演模拟、并进行重、磁异常场的波数相关性分析.研究发现在重、磁异常中均发现一个呈负相关的纬向构造特征,东亚地区“北纬40°断裂带”横穿了整个朝鲜半岛.     

利用远震接收函数研究黄海周边地区的地壳上地幔结构

利用宽频带流动台站(YSBSN)记录的远震波形数据和远震接收函数方法,反演了黄海东、西两侧地壳上地幔的S波速度结构.结果表明,莫霍面深度在30~38 km之间变化,位于中方一侧的JNN台下方地壳厚度最大,可以归因于华北板块和扬子板块的碰撞;韩方一侧的地壳厚度自北向南逐渐变厚,但仍然难以厘定朝鲜半岛南部潜在碰撞带的位置,这些问题的解决需要更大范围的流动台站观测.由于部分台站位于巨厚的沉积层和多孔的火山岩之上,与浅部构造的相关性使得接收函数表现出较大振幅的混响,从而影响了来自深部结构的转换震相.     

格列菲列岛的地质地球物理研究

本文是中、韩国际合作项目“黄海地质地球物理场特征研究”的一部分 .具体研究了格列菲列岛的基础地质问题及地球物理场特征 .格列菲列岛是由三个小岛组成的位于黄海中部的一组群岛 ,位于朝鲜半岛布格重力异常正值区 ,莫霍面埋深在 2 9km左右 .三个岛屿主要由安山集块岩构成 .安山集块岩是由大小不同的火成岩碎块和熔岩混合而成 .火成岩的碎块成分主要为石英岩、片岩、片麻岩、玄武岩、安山岩和流纹岩 .主要节理方向为NE10°～NE2 0°,倾角近乎垂直 ;另一组节理为NW70°～NW80° ,倾角也近垂直 .本研究为我们了解这个地区的基础地质情况 ,岩石物性和随后进行的黄海海域典型剖面的地球物理物性反演工作以及地质模型的建立 ,打下重要基础 .     

Three-dimensional topography corrections of magnetotelluric data

Topographic effects due to irregular surface terrain may prevent accurate interpretation of magnetotelluric (MT) data. Three-dimensional (3-D) topographic effects have been investigated for a trapezoidal hill model using an edge finite-element method. The 3-D topography generates significant MT anomalies, and has both galvanic and inductive effects in any polarization. This paper presents two different correction algorithms, which are applied to the impedance tensor and to both electric and magnetic fields, respectively, to reduce topographic effects on MT data. The correction procedures using a homogeneous background resistivity derived from a simple averaging method effectively decrease distortions caused by surface topography, and improve the quality of subsurface interpretation. Nonlinear least-squares inversion of topography-corrected data successfully recovers most of structures including a conductive or resistive dyke.