边界层参数化方案对边界层热力和动力结构特征影响的比较

本文利用高分辨率中尺度WRF模式,通过改变边界层参数化方案进行多组试验,评估该模式对美国北部森林地区边界层结构的模拟能力,同时比较了五种不同边界层参数化方案模拟得出的边界层热力和动力结构.结果表明：除个别方案外,配合不同边界层方案的WRF模式都能成功模拟出白天对流边界层强湍流混合特征和夜间稳定边界层内强逆温、逆湿和低空急流等热力和动力结构.非局地YSU、ACM2方案在白天表现出强的湍流混合和卷夹,相比于局地MYJ、UW方案,模拟的对流边界层温度更高、湿度更低、混合层高度更高、感热通量更大,更接近实际观测,这表明在不稳定层结下考虑非局地大涡输送更为合理,但局地方案在风速和风向的预报上存在一定优势.TEMF方案得到的白天局地湍流混合强度为所有方案中最弱,混合层难以发展,无法体现对流边界层内气象要素垂直分布均匀的特点.对于夜间稳定边界层的模拟,不同参数化方案之间的差异较小,但是YSU方案在一定程度上高估了机械湍流,导致局地湍流混合偏强,从而影响了其对稳定边界层的模拟能力.     

壳内多震层孕震环境研究进展

壳内多震层的介质结构及其内、外环境是制约地震（尤其强震）在该层内孕育、发生的重要因素。现今多震层孕震环境研究所取得的主要进展如下：多震层位于上地壳下部至中地壳,岩石组合为角闪岩相（上部为绿片岩相）变质岩和花岗质岩石;多震层具有相对高速、高密度、高阻的介质结构,其下部对应脆 韧变形转换带和最大剪切应力带;多震层（或强震震源）之下存在壳内低速 高导体;多震层应力场具有区域性分区特点,壳内地震活动截止温度约３００～４００℃;多震层内震源断层具有区域性特点,不同地区震源断层可能具有不同的孕震机制;“坚固体孕震模式”对地震预报具有一定的指导作用     

中国不同纬度背景Na层夜间和季节变化特征的激光雷达研究

基于子午工程的北京、合肥和海南三个Na荧光激光雷达对中国不同纬度上空夜间背景Na层的长期观测,分析了我国沿东经120°N上空Na层的夜间变化和季节变化特征.对照2010年12月17日夜间三个雷达站的观测结果,发现三个地方Na层的夜间变化并不具有相关性.Na层长期变化的年加半年变化拟合结果显示,北京和合肥上空Na层柱密度具有明显的年变化特点,而海南上空Na层柱密度的半年变化特征更明显; Na层的质心高度和RMS宽度具有明显的半年变化特点,但海南地区Na层的RMS宽度的长期变化不具有半年变化特征.Na层参数的统计和对照显示,Na层柱密度的季节变化与大气温度季节变化相关,在冬季最大,在夏季最小.Na层柱密度随纬度升高而增大,同时年变化性逐渐增强;质心高度随纬度变化趋势不明显,但三个地方Na层质心高度的长期变化都具有较明显的半年变化特征;在各个月份中,北京地区Na层RMS宽度最大,合肥地区Na层RMS宽度最小,海南地区居中.     

土层对地震的随机反应分析

利用随机振动理论,本文研究了非均匀土层对地震的动力反应问题,假定土层的剪切模量随深度呈指数函数形式分布,关于基岩输入地震动加速度的功率谱密度函数,考虑了两种形式：白噪声谱和有色谱,为了比较起见,文中考虑了均匀土层的情形,数值计算表明;（１）在基岩输入地震动加速度的功率谱密度为白谱的情况下,土层的最大期望反应（相对位移,绝对加速度,剪应变）大于基岩输入地震动的功率谱为有色谱的情形;（２）在土层的上部     

土层结构对地表加速度峰值的影响

本文在研究我国数百个工程场地钻孔资料的基础上，选取和构造了若干有代表性的典型场地剖面，利用目前工程上广泛应用的场地地震反应分析的一维等效线性化波动方法，计算了不同场地在三种不同强度的地震动输入下的地表加速度峰值反应。分析研究了覆盖层厚度、软土层的埋深与厚度等对场地地表加速度峰值的影响，得到了一些有意义的结论。     

Magma ascent beneath Unzen Volcano, SW Japan, deduced from the electrical resistivity structure

The resistivity structure of Unzen Volcano has been revealed by extensive magnetotelluric surveys since the first eruption on November 17, 1990. This structure comprises a highly resistive surface layer, a low-resistive second layer at several hundred meters depth, interpreted as a water-saturated layer, a resistive third layer, and a low-resistive fourth layer at 10 km depth, possibly related to the deep magmatic activity. The structure has influenced the volcanic activity of Unzen. This activity was characterized by a series of dramatic changes in eruption type: a minor phreatic eruption on November 17, 1990; phreatic eruptions after February 12, 1991, preceded by several weeks of volcanic tremor; phreatomagmatic eruptions after April 9, and dome effusion beginning May 19, 1991. This paper presents a hypothesis in which the top of the magma column rose about 20 m/day, reached the base of the water-saturated layer at the end of January, 1991, and approached the upper boundary of this layer on April 9. Thus, the temporal change of eruption type and associated phenomena are systematically explained by an interaction between magma and groundwater contained in the saturated layer.     

Simultaneous crystallization and melting at both the roof and floor of crustal magma chambers: Experimental study using NH4Cl–H2O binary eutectic system

Thermal and compositional evolution of magmas after emplacement of basalt into continental crust has been investigated by means of fluid dynamic experiments using a cold solid mixture with eutectic composition and a hot liquid with higher salinity in the NH₄Cl–H₂O binary eutectic system. The experiments were designed to simulate cases where crystallization of a basalt magma is accompanied by melting at both the roof and floor of a crustal magma chamber. The results show that thermal and compositional convection occur simultaneously in the solution; the thermal convection is driven by cooling at the roof and the compositional convection is driven by melting and crystallization at the floor. The roof was rapidly melted by the convective heat flux, which resulted in formation of a separate eutectic melt layer (the upper liquid layer) with negligible mixing of the underlying liquid (the lower liquid layer). On the other hand, a mushy layer formed at the floor. The compositional convection at the floor carried a low heat flux, so that the heat transfer at the floor was basically explained by simple heat conduction. The thermal boundary layer in the lower liquid layer at the interface with the upper liquid layer became thicker with time and subsequently temperature decreased upward throughout the lower liquid layer. Compositional gradient with NH₄Cl content decreasing upward formed by compositional convection in the lower liquid layer. The formation of these gradients resulted in formation of double-diffusive convecting layers in the lower liquid layer. The upward heat transfer was suppressed when compared with the case where the liquid region is homogenized by vigorous convection.These experimental results imply that, when a basalt magma is emplaced in continental crust, floor melting does not always enhance the cooling of the magma, but it may even reduce the total heat loss from the magma to the crusts due to suppression of convection by formation of a stabilizing compositional gradient.     

西藏高原斜压对流边界层风、温、湿廓线特征

利用1998年夏季第2次青藏高原大气科学试验当雄观测站的边界层观测资料以及拉萨、改则和武汉等地探空资料,分析讨论西藏高原斜压对流边界层风、温、湿廓线的特征. 研究结果表明,高原地区白天对流边界层发展可高达2200m,显著超过中纬度平原地区和海面上对流边界层高度. 高原对流边界层中温度廓线具有较好的混合特征,湿度廓线有时在某一高度上出现湿度极大值. 高原对流边界层内热量和水汽收支分析表明,水平平流作用对边界层结构具有重要作用. 在对流边界层中平均风速垂直分布存在风切变现象. 水平温度梯度形成较强的斜压性是形成边界层风切变的主要原因.     

来自海底高速层径向波的理论地震图研究

本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.     

MULTIPLE REFLECTIONS AND HEAD WAVES IN THE GULF OF SUEZ*

Inferior reflection quality in the Gulf of Suez at the target depth interval is attributable in part to surficial multiple reflections. An excellent example of the latter is observed on a typical seismic line in the northern portion of the Gulf. An increase in prominence of the multiple reflections appears associated with decreasing depth to a dipping highvelocity layer. Inversion of a second-order polynomial time-distance function, fitted to the observed refraction onset time-distance values, gives the velocity-depth function for sediments between the water bottom and a high-velocity layer. Velocities thus determined increase non-linearly with depth from a value near water velocity at the water bottom. Depths to the high-velocity layer are obtained from the associated head-wave linear time-distance function and by ray tracing in the overlying sediments. As the high-velocity layer approaches the water bottom from sub-water depths exceeding 0.6 km to a depth of 56 m, intensity of the multiple reflections increases to the extent of completely dominating individual records to a time of at least 3 s. The estimated plane-wave normal-incident reflection coefficient at the top of the high-velocity layer increases with decreasing depth to this layer, approaching 0.5 at the shallowest depth. This strong reflection coefficient further substantiates the existence of multiple reflections between the high-velocity layer and water layer. However, existence of water-layer multiples cannot be ruled out. The estimated water-bottom reflection coefficient is approximately 0.3, a substantial value. Multiple reflections of considerably less intensity are apparent where the high-velocity layer is deepest, and it is likely that such are waterlayer multiple reflections. Unfortunately, water-layer multiple reflections and multiple reflections between the water surface and high-velocity layer cannot be separated by their coincidence with time-distance (normal moveout) curves, the configuration of each visibly matching the curves equally well.     

F_1层未充分发展时的电离层剖面反演

垂测电离图反演对研究电离层结构及运动、电离层波传播、空间气象应用等具有重要意义,一直以来受到十分广泛的重视.F_1层未充分发展的回波描迹是垂测电离图中观测到的最多的一种情况,此时垂测描迹上表现为F_1层到F_2层的平滑过渡,而不是F_1层充分发展时出现的尖点,然而现有电离层模型以及反演算法更多的是针对F_1层充分发展情况,即模型剖面在F_1层峰高处近似为抛物形状,此处电子浓度剖面梯度为无穷大,这不符合F_1层未充分发展时的剖面形态,即F_1层剖面未到达抛物顶点处就已经达到该层电子浓度的最大值然后到达F_2层,剖面梯度为一个有限的数值.因此,本文针对垂直探测中经常出现的F_1层未充分发展情况,引入F_1层模型设定临频这一参数,建立了F_1层未充分发展时的基于移位切比雪夫多项式的电离层模型,并从电离层剖面的光滑性考虑,在所建模型的基础上,提出了约束优化F_1层参数、F_2层参数的垂测电离图反演方法.然后,通过仿真分析,从理论上验证了所建电离层模型和反演算法的合理性,并通过反演电子浓度剖面合成的垂测、斜测描迹与实测数据的对比,对反演算法的有效性进行了进一步的试验验证.     

新安江水库(千岛湖)热力学状况及热力分层研究

利用2012年1-12月在新安江水库(千岛湖)6个点位的每月一次的水温及其他环境因子的周年观测资料,分析了水库水温逐月变化、季节变化、垂直分布及温跃层的形成与变化,探讨了温跃层特征量(温跃层深度、厚度、强度)与表层水温、水体透明度的关系.新安江水库表层和中层水温与气温存在显著的线性相关,又以表层水温线性关系最好,而下层水温与气温没有显著相关性,说明下层水温受气温的影响很小,全年处于相对恒温状态.水库表层和中层水温逐月变化明显,呈现夏季最高、春秋季次之、冬季最低的变化趋势,其中中层水温最高值出现的季节较表层水温明显后延,下层水温没有明显的逐月变化和季节变化.水温垂直分布显示,4个季节均存在不同程度的温跃层和温度分层现象,其中水深最深的大坝前水温分层最明显.小金山、三潭岛和大坝前3个典型点位从春季的4月份到冬季的2月份温跃层深度由1.61±0.47 m逐渐增加至39.37±5.35 m,而温跃层厚度和强度则在夏季最高、冬季最低,温跃层随着季节的变化呈现增强稳定减弱消失的周期变化.温跃层深度与水体透明度存在显著正相关,与表层水温存在显著负相关,并基于透明度和表层水温建立温跃层深度的多元线性回归模型.     

Near-surface layer traveltime inversion: a synthetic example1

The near-surface layer is modelled as a constant-velocity layer with varying thickness. The base of the layer is described by a B-spline curve. The optimum model is calculated by minimizing, with respect to the model parameters, the difference between traveltimes predicted by the model and those observed in the data. Once a model has been produced, corrections that are dependent on the raypath geometry through the near-surface layer can be calculated. The effect of the near-surface layer is normally considered to be consistent at each shot or geophone station for all traveltimes arriving at that location (the surface-consistent approximation). This assumption linearizes the problem, allowing timeshifts to be calculated and the traveltimes corrected to a chosen datum, representing static corrections. The single correction at each point is an averaged correction, based on an assumption that is particularly inaccurate in the presence of lateral variations of velocity or thickness of the surface layer, in the presence of large surface layer velocities or in the presence of a thick surface layer. The method presented considers the non-linear relationship between data and model explicitly, hence the correction that is dependent on the raypath. Linearization removes this dependence and reduces the problem to a surface-consistent approximation. The method is applied to synthetic data calculated from a model with surface layer variations. Comparisons are made between the corrected data resulting from the method described here and the conventional surface-consistent approach. From these results it becomes apparent that the near-surface layer inversion method presented here can reproduce accurate models and correct for near-surface layer effects in cases where conventional methods encounter difficulties. Additionally the method can be readily extended to 3D.     

A nonlinear,resistive boundary layer in rotating hydromagnetic flow

A new nonlinear boundary layer in rotating hydromagnetic flows is presented. The purpose of this layer is to provide a smooth transition between the magnetic field at a rigid, electrically insulating boundary and that far from the boundary. The boundary layer problem is solved in a cylindrical geometry, assuming that the variables obey the hydromagnetic analog of the von Karman similarity. The primary momentum balance within this boundary layer is between pressure, Coriolis and hydromagnetic forces; the viscous and nonlinear inertia forces are unimportant. Diffusion of the magnetic field is balanced by advection where the advecting fluid flow is induced by the action of the hydromagnetic body forces within the boundary layer. The structure of the layer depends upon a single parameter which measures the strength of the normal magnetic field. Steady solutions are presented and their uniqueness and temporal stability are analyzed. The relevance of this layer to the core of the Earth is discussed. It is estimated that this boundary layer would be at least as thick as one tenth of the core radius if it exists within the core.     

极区电子沉降对电离层影响的模拟研究

在极区电离层模式中,考虑了软电子沉降引起的电离,并对差分方法的应用提出了改进意见.模拟了不同特征的沉降电子对极区电离层的影响,发现平均能量低的电子束能够形成明显的电离层F层,平均能量较高的电子束能使得最大电离的高度下移,形成明显的E层,甚至其电子浓度高于F层.将电子沉降的卫星测量结果作为电离层模型的输入,所得F层临频与观测结果符合较好.通过分析中山站电离层统计结果,综合电子沉降在极隙区的分布特征和上述模拟结果,认为中山站磁中午现象主要由电子沉降所致.     

Application of a perfectly matched layer in seismic wavefield simulation with an irregular free surface

Recently, an effective and powerful approach for simulating seismic wave propagation in elastic media with an irregular free surface was proposed. However, in previous studies, researchers used the periodic condition and/or sponge boundary condition to attenuate artificial reflections at boundaries of a computational domain. As demonstrated in many literatures, either the periodic condition or sponge boundary condition is simple but much less effective than the well‐known perfectly matched layer boundary condition. In view of this, we intend to introduce a perfectly matched layer to simulate seismic wavefields in unbounded models with an irregular free surface. We first incorporate a perfectly matched layer into wave equations formulated in a frequency domain in Cartesian coordinates. We then transform them back into a time domain through inverse Fourier transformation. Afterwards, we use a boundary‐conforming grid and map a rectangular grid onto a curved one, which allows us to transform the equations and free surface boundary conditions from Cartesian coordinates to curvilinear coordinates. As numerical examples show, if free surface boundary conditions are imposed at the top border of a model, then it should also be incorporated into the perfectly matched layer imposed at the top‐left and top‐ right corners of a 2D model where the free surface boundary conditions and perfectly matched layer encounter; otherwise, reflections will occur at the intersections of the free surface and the perfectly matched layer, which is confirmed in this paper. So, by replacing normal second derivatives in wave equations in curvilinear coordinates with free surface boundary conditions, we successfully implement the free surface boundary conditions into the perfectly matched layer at the top‐left and top‐right corners of a 2D model at the surface. A number of numerical examples show that the perfectly matched layer constructed in this study is effective in simulating wave propagation in unbounded media and the algorithm for implementation of the perfectly matched layer and free surface boundary conditions is stable for long‐time wavefield simulation on models with an irregular free surface.     

Observations and model calculations of an additional layer in the topside ionosphere above Fortaleza, Brazil

Calculations using the Sheffield University plasmasphere ionosphere model have shown that under certain conditions an additional layer can form in the low latitude topside ionosphere. This layer (the F3 layer) has subsequently been observed in ionograms recorded at Fortaleza in Brazil. It has not been observed in ionograms recorded at the neighbouring station São Luis. Model calculations have shown that the F3 layer is most likely to form in summer at Fortaleza due to a combination of the neutral wind and theE × B drift acting to raise the plasma. At the location of São Luis, almost on the geomagnetic equator, the neutral wind has a smaller vertical component so the F3 layer does not form.     

Effect of the oceanic lithosphere velocity on free convection in the asthenosphere beneath mid-ocean ridges

The paper presents results obtained in experiments on a horizontal layer heated from below in its central part and cooled from above; the layer models the oceanic asthenosphere. Flow velocity and temperature profiles are measured and the flow structure under boundary layer conditions is determined (at Rayleigh numbers Ra > 5 × 10⁵). The flow in the core of a plane horizontal layer heated laterally and cooled from above develops under conditions of a constant temperature gradient averaged over the layer thickness. The flow core is modeled by a horizontal layer with a moving upper boundary and with adiabatic bounding surfaces under conditions of a constant horizontal gradient of temperature. Exact solutions of free convection equations are found for this model in the Boussinesq approximation. Model results are compared with experimental data. Temperature and flow velocity ranges are determined for the boundary layer regime. Based on the experimental flow velocity profiles, an expression is found for the flow velocity profile in a horizontal layer with a mobile upper boundary heated laterally and cooled from above. Free convection velocity profiles are obtained for the asthenosphere beneath a mid-ocean ridge (MOR) with a mobile lithosphere. An expression is obtained for the tangential stress at the top of the asthenosphere beneath an MOR and the total friction force produced by the asthenospheric flow at the asthenosphere-lithosphere boundary is determined.     

Apparent resistivity over a three layer earth with anisotropic interstratum

Summary The present work aims at investigating the effects of anisotropy on the apparent resistivity measurements over the surface of a three layer Earth. An appropriate model with the anisotropic layer embedded between two isotropic media is chosen for study. The point current source is assumed to be present at the surface of such a model. After writing the values of the potentials in each layer, expressions for apparent resistivities for Wenner and Schlumberger configurations are derived. Also, the expressions for potential for the limiting cases in which the bottom layer is infinitely resistive or infinitely conducting are derived. It is shown that these expressions can be recast into simpler forms corresponding to isotropic layer (s) on an isotropic half-space.     

Vertical structures of PM10 and PM 2.5 and their dynamical character in low atmosphere in Beijing urban areas

The vertical structures and their dynamical character of PM2.5 and PM10 over Beijing urban areas are revealed using the 1 min mean continuous mass concentration data of PM2.5 and PM10 at 8, 100, and 320 m heights of the meteorological observation tower of 325 m at Institute of Atmospheric Physics, Chinese Academy of Sciences (IAP CAS tower hereafter) on 10―26 August, 2003, as well as the daily mean mass concentration data of PM2.5 and PM10 and the continuous data of CO and NO2 at 8, 100 (low layer), 200 (middle layer), and 320 m (high layer) heights, in combination with the same period meteorological field observation data of the meteorological tower. The vertical distributions of aerosols observed on IAP CAS tower in Beijing can be roughly divided into two patterns: gradually and rapidly decreasing patterns, I.e. The vertical distribution of aerosols in calm weather or on pollution day belongs to the gradually decreasing pattern, while one on clean day or weak cold air day belongs to the rapidly decreasing pattern. The vertical distributive characters of aerosols were closely related with the dynamical/thermal structure and turbulence character of the atmosphere boundary layer. On the clean day, the low layer PM2.5 and PM10 concentrations were close to those at 8 m height, while the concentrations rapidly decreased at the high layer, and their values were only one half of those at 8 m, especially, the concentration of PM2.5 dropped even more. On the clean day, there existed stronger turbulence below 150 m, aerosols were well mixed, but blocked by the more stronger inversion layer aloft, and meanwhile, at various heights, especially in the high layer, the horizontal wind speed was larger, resulting in the rapid decrease of aerosol concentration, I.e. Resulting in the obvious vertical difference of aerosol concentrations between the low and high layers. On the pollution day, the concentrations of PM2.5 and PM10 at the low, middle, and high layers dropped successively by, on average, about 10% for each layer in comparison with those at 8 m height. On pollution days, in company with the low wind speed, there existed two shallow inversion layers in the boundary layer, but aerosols might be, to some extent, mixed below the inversion layer, therefore, on the pollution day the concentrations of PM2.5 and PM10 dropped with height slowly; and the observational results also show that the concentrations at 320 m height were obviously high under SW and SE winds, but at other heights, the concentrations were not correlated with wind directions. The computational results of footprint analysis suggest that this was due to the fact that the 320 m height was impacted by the pollutants transfer of southerly flow from the southern peripheral heavier polluted areas, such as Baoding, and Shijiazhuang of Hebei Province, Tianjin, and Shandong Province, etc., while the low layer was only affected by Beijing's local pollution source. The computational results of power spectra and periods preliminarily reveal that under the condition of calm weather, the periods of PM10 concentration at various heights of the tower were on the order of minutes, while in cases of larger wind speed, the concentrations of PM2.5 and PM10 at 320 m height not only had the short periods of minute-order, but also the longer periods of hour order. Consistent with the conclusion previously drawn by Ding et al., that air pollutants at different heights and at different sites in Beijing had the character of "in-phase" variation, was also observed for the diurnal variation and mean diurnal variation of PM2.5 and PM10 at various heights of the tower in this experiment, again confirming the "in-phase" temporal/spatial distributive character of air pollutants in the urban canopy of Beijing. The gentle double-peak character of the mean diurnal variation of PM2.5 and PM10 was closely related with the evident/similar diurnal variation of turbulent momentum fluxes, sensible heat fluxes, and turbulent kinetic energy at various heights in the urban canopy. Besides, under the condition of calm weather, the concentration of PM2.5 and PM10 declined with height slowly, it was 90% of 8 m concentration at the low layer, a little lesser than 90% at the middle layer, and 80% at the high layer, respectively. Under the condition of weak cold air weather, the concentration remarkably dropped with height, it was 70% of 8 m concentration at the low layer, and 20%―30% at the middle and high layers, especially the concentration of PM2.5 was even lower.