A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

大连滨海湿地景观格局变化及其驱动机制

以2000年和2006年TM卫星影像为主要数据源,配合其它非遥感数据,在遥感与地理信息系统技术支持下,运用景观生态学原理,选取反映景观空间结构和景观异质性的指数,对大连地区湿地的整体景观格局和类型景观格局及其动态变化进行定量分析.结果表明:6年间,大连湿地面积减少了97.62 km2;整体景观多样性指数和均匀度指数降低,优势度指数增加;各类景观格局时间序列上也存在明显差异性变化.湿地景观格局指数的变化, 反映了移山填海工业园区的扩大及养殖业的大力发展等人为活动对景观格局的深刻影响.人为活动已成为大连市湿地景观格局变化的主要驱动因子.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

热液成矿作用的地质地球化学动力学——以胶东牟平-乳山金矿带为例

通过对胶东牟平—乳山热液脉状金,铅锌矿床成矿动力学控制规律的研究,说明了建立热液成矿反应体系的方法,并得出以下认识:（1）区域断裂构造活动,通过影响成矿反应体系的热力学性质和条件,控制元素的富集和分散。（2）热液成矿过程中容矿断裂活动可划分为两种作用方式:脆性破裂和韧—脆性张开,构成热液成矿的两种构造动力学环境。（3）热液成矿反应体系是一种开放的动态的热力学体系。断裂的脆性破裂阶段,使体系处于强烈的过饱和状态,矿质在远离平衡状态下快速结晶,加剧物质分异;在脆—韧性张开阶段,矿质在接近平衡态体系中缓慢晶出。（4）热液演化晚期,金在残余溶液中富集;当断裂再次发生构造脉动震颤时,早期形成的块状黄铁矿矿石碎裂,富金溶液充填其中,形成含裂隙金和晶隙金的富矿石。成矿体系热力学演化与构造动力学条件有利的匹配控制着富矿石的形成。     

鲁23井“动水位”观测的分析

本文介绍了鲁23井“动水位”观测情况。提出了一个分析有泄流的“动水位”的新方法,即把流量换算成水位,以等效静水位来表示“动水位”与流量。由于把“动水位”与流量用一个量来计算,从而为“动水位”观测定量的、统一的分析研究提供了新的方法。这对提高“动水位”井的监测预报地震能力起一定的作用。 本文还以鲁23井“动水位”观测资料为例,利用新方法进行处理的观测分析结果与理论研究结果是一致的。     

Room temperature solid-state detectors for Gamma ray burst spectroscopy

We have evaluated several solid state detectors which offer excellent energy resolution at room temperature for soft X-rays. For soft X-rays (< 1 keV to 20 keV), silicon P-intrinsic-N (PIN) and avalanche-mode photodiodes (APD's) have been studied. Using commercially available charge sensitive pre-amplifiers, these photodiodes provide 1 keV resolution without cooling. Their detection efficiencies are limited to about 20 keV and 15 keV, respectively. To overcome this constraint, we have studied thick (1.5 mm) PIN detectors made by Micron Semiconductor Ltd., U.K., as well as compound semiconducting materials with high Z constituents such as CZT and PbI₂. PbI₂ allows high detection efficiencies of photons up to 100 keV with detectors 100–300 microns thick. These new detectors offer the capability to study the low-energy spectral evolution of Gamma ray bursts (GRBs). A matrix of these detectors could be used as an image plane detector with moderate spatial resolution for imaging.     

Effects of dynamic shear strength and time-histories stress analysis on newmark sliding block analyses

A simplified procedure for evaluating aseismic stability of slope subjected to earthquake shaking, in which the effect of dynamic shear strength and time-history stress on the yielding angular acceleration of sliding block is taken into account, is presented. The fundamental feature of this procedure is the dynamic shear strength. The numerical computations are performed by using the proposed method. It is shown that the computed sliding displacement for a given core dam, with either method of dynamic shear strength properly considered, is more rational compared with the conventional computational results based on static shear strength.     

林草复合经营模式养分动态关系

林草复合经营模式在目前的退耕还林工程中起着重要的作用，研究林草复合模式的林与草之间的营养元素关系，通过对不同林草养分的研究得出林木、草、土壤之间的养分动态关系，对矿质养分相互关系的研究为合理搭配林草提供重要依据，为退化山地生态系统的恢复和重建提供重要的模式。模式Ⅰ在退耕还林初期，土壤养分供给氮素虽未亏缺，但从三者问的消长关系看出：氮磷素供应并不富足，特别是水解性氮。土壤交换性钙供应不足，引起杉木与黑麦草对钙质竞争。而土壤全钙含量很高，钙质向交换性钙转化的速度缓慢。模式Ⅱ与模式Ⅰ得到相似结果，所以在这两种模式中应加强人为管理，在退耕初期应适当增施氮磷钙肥，辅以土壤结构改良，加速土壤养分向植物可利用态转化。模式Ⅲ牛鞭草的生物量极高，磷、钙索相对缺乏。由于牛鞭草被不断的采割用于牲畜的饲料，对该模式中必须进行一定养分的输入，以保持杂交竹与牛鞭草这一林草生态系统的养分平衡。     

1993～2003年厦门市湿地动态变化及其驱动因素分析

通过对厦门市1993年TM遥感影像和2003年ETM 遥感影像的解译分析和数据统计,从湿地面积动态变化和湿地消长情况等方面,用单一土地利用动态度等指标对1993~2003年期间厦门市的湿地变化进行了研究。结果表明,1993~2003年期间厦门市天然湿地面积减少了2103hm2,其中有954hm2(滩涂864hm2)转变为非湿地,有1149hm2转变为人工湿地;而人工湿地面积增加了2775hm2,其中养殖场面积从1993年的5778hm2迅猛增加到2003年的9768hm2,年均增长率为6.91%,水库与水塘和盐田则分别减少了961hm2和254hm2。有27.74%、11.41%、9.33%和3.63%的养殖场面积分别是由非湿地、滩涂、盐田和水库与水塘转变而来的;有54.37%的盐田转变成了养殖场,有33.24%、30.24%的红树林转变为非湿地、养殖场,有10.35%、8.64%的滩涂转变为养殖场、非湿地。对社会经济驱动因素的分析表明,经济发展、人口增长和农业生产比较利益因素极大影响着厦门市湿地的动态变化。     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.