Pull-down basin in the central part of Japan due to subduction-induced mantle flow

The driving force for the basin subsiding against isostatic balance in and around Lake Biwa in the Kinki district, Japan is discussed. The lake region is characterized by strong negative Bouguer anomalies, especially by a steep horizontal gradient zone of gravity anomaly running along the western margin of the lake. The large negative anomaly (>50 mgal) cannot be explained by low-density sediments beneath it. A down-warping structure extending to the Moho depth should be taken into account. This conjecture has been strongly supported by a short-period receiver function imaging, which shows a clear offset of about 8 km for the Moho discontinuity under the steep gravity gradient zone.A question arises as to what is the driving force to create such a large down-warping structure. We consider that the subduction of the shallow-dipping slab under the region (Philippine Sea Slab) may cause crustal deformation by dragging the viscous mantle downward. In order to verify this model, we simulated the induced mantle flow due to the subduction of the Philippine Sea Slab and the pressure distribution on the crust–mantle boundary. This numerical experiment showed that the induced flow makes a strong negative pressure zone under the lake region if the slab has a vertical offset along the direction of subduction. This offset of the slab is consistent with plate models deduced from hypocentral distributions and Sp phases of the deep-focus earthquakes.     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.     

Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.     

一次区域性大暴雨环境特征分析

通过对冀东南、鲁西北一次区域性大暴雨的分析表明：除有利的大尺度环流背景以外，超低空急流是大暴雨不稳定能量和水汽的主要输送者，高能区不稳定能量释放，产生大暴雨所需的强烈上升运动。暴雨落区与超低空急流前部最大风速切变区的位置密切相关。     

新形势下县局天气预报思路浅谈

在9210工程业务环境下，结合县局的预报工作实际，建立了以数值预报指导产品为基础，以人机交互系统为平台，以提高县局预报业务效率为目的，制作的县级预报思路，并投入实际业务运行中，收到良好的效果。     

上地幔变黏度小尺度对流的数值研究

基于二维模型，利用有限元方法，研究上地幔-岩石圈系统的变黏度小尺度对流. 考虑该系统的黏度随温度以指数形式变化，数值结果表明，当黏度随温度变化较剧烈时，由于低温高黏度，系统的最上部物质不参与对流，系统发育形成一个类似于岩石圈的静止盖层. 计算表面热流、地形起伏及重力异常与对流格局有较好的相关性，高热流、上升地形对应于对流的上升区，反之低热流、下降地形与对流的下降区对应.     

多相流动电磁波成像测井测量敏感场计算

在多相流动电磁波成像测井研究中，测量敏感场是设计测量仪器和研究图像重建算法的前提工作. 本文用数值方法计算多相流动电磁波成像测井测量敏感场. 根据测量探头电极阵列的对称性，将三维敏感场问题划分为轴向和横向两个平面问题，然后应用有限元算法求解两个平面电势分布问题，对于不同介质模型，计算各种测量条件下测量区域内电势分布，进而根据电势分布计算测量敏感场. 结果表明测量敏感区域为发射电极到测量电极的弧形区域，在发射电极到测量电极的电势梯度线上，测量敏感函数近似于一个指数函数.     

河北省黄壁庄水库坝基渗漏单井示踪试验研究

以河北省黄壁庄水库副坝坝基的渗漏塌坑为例，利用“九五”国家重点科技成果推广项目“智能化地下水动态参数测量仪”在天然流场下，测量出水库自然渗漏量地下水的渗透流速流向、垂直向流速流向，每米含水层的渗透系数等工程必须的水文地质参数。确定性地诊断和检测出黄壁水库坝基渗漏隐患出现的具体位置和高程，从而为大坝的渗漏处理和安全处理提供详细的现场测量数据。     

利用数字图像处理技术提高地震剖面图像信噪比

提出了利用数字图像处理技术提高地震剖面信噪比的新方法，首先根据数字图像处理要求的格式，对地震剖面数据进行转换，得到地震剖面图像，分析了地震数据特点和初步地震图像的实验结果后，设计了新的预处理方法——“二维沿层滤波”，在此基础上，利用可以计算帧间运动速度及其变化都较大的改进的光流分析技术，计算出多幅地震剖面对应点的偏移量，然后应用图像积累技术对这多幅地震剖面进行积累，实现对三维地震数据体提高信噪比的处理，该方法充分利用了三维地震信息，不但可以提高整个数据体的信噪比，而且可以减少信号能量的损失，并保持原来的信号能量关系，使地震剖面的质量得到明显提高，为地震解释奠定良好的基础。     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.