Comparison of alternative upscaled model formulations for simulating DNAPL source dissolution and biodecay

This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution.     

Source process of very-long-period events accompanying long-period signals at Cotopaxi Volcano,Ecuador

Renewed seismic activity of Cotopaxi, Ecuador, began in January 2001 with the increased number of long-period (LP) events, followed by a swarm of volcano-tectonic (VT) earthquakes in November 2001. In late June 2002, the activity of very-long-period (VLP) (2 s) events accompanying LP (0.5–1 s) signals began beneath the volcano. The VLP waveform was characterized by an impulsive signature, which was accompanied by the LP signal showing non-harmonic oscillations. We observed temporal changes of both the VLP and LP signals from the beginning until September 2003: The VLP signal gradually disappeared and the LP signal characterized by decaying harmonic oscillations became dominant. Assuming possible source geometries, we applied a waveform inversion method to the observed waveforms of the largest VLP event. Our inversion and particle motion analyses point to volumetric changes of a sub-vertical crack as the VLP source, which is located at a depth of 2–3 km beneath the northeastern flank. The spectral analysis of the decaying harmonic oscillations of LP events shows frequencies between 2.0 and 3.5 Hz, with quality factors significantly above 100. The increased VT activity and deformation data suggest an intrusion of magma beneath the volcano. A release of gases with small magma particles may have repetitively occurred due to the pressurization, which was caused by sustained bubble growth at the magma ceiling. The released particle-laden gases opened a crack above the magma system and triggered the resonance of the crack. We interpret the VLP and LP events as the gas-release process and the resonance of the crack, respectively.     

A stochastic ground motion model with separable temporal and spectral nonstationarities

A fully nonstationary stochastic model for strong earthquake ground motion is developed. The model employs filtering of a discretized white‐noise process. Nonstationarity is achieved by modulating the intensity and varying the filter properties in time. The formulation has the important advantage of separating the temporal and spectral nonstationary characteristics of the process, thereby allowing flexibility and ease in modeling and parameter estimation. The model is fitted to target ground motions by matching a set of statistical characteristics, including the mean‐square intensity, the cumulative mean number of zero‐level up‐crossings and a measure of the bandwidth, all expressed as functions of time. Post‐processing by a second filter assures zero residual velocity and displacement, and improves the match to response spectral ordinates for long periods. Copyright © 2008 John Wiley & Sons, Ltd.     

The behavior of specific sediment yield in different grain size fractions in the tributaries of the middle Yellow River as influenced by eolian and fluvial processes

Based on data from 35 stations on the tributaries of the Yellow River, annual specific sediment yield (Y_s) in eight grain size fractions has been related to basin‐averaged annual sand–dust storm days (D_ss) and annual precipitation (P_m) to reveal the influence of eolian and fluvial processes on specific sediment yield in different grain size fractions. The results show that Y_s in fine grain size fractions has the highest values in the areas dominated by the coupled wind–water process. From these areas to those dominated by the eolian process or to those dominated by the fluvial process, Y_s tends to decrease. For relatively coarse grain size fractions, Y_s has monotonic variation, i.e. with the increase in D_ss or the decrease in P_m, Y_s increases. This indicates that the sediment producing behavior for fine sediments is different from that for relatively coarse sediments. The results all show that Y_s for relatively coarse sediments depends on the eolian process more than on the fluvial process, and the coarser the sediment fractions the stronger the dependence of the Y_s on the eolian process. The Y_s–D_ss and Y_s–P_m curves for fine grain size fractions show some peaks and the fitted straight lines for Y_s–D_ss and Y_s–P_m relationships for relatively coarse grain size fractions show some breaks. Almost all these break points may be regarded as thresholds. These thresholds are all located in the areas dominated by the coupled wind–water process, indicating that these areas are sensitive for erosion and sediment production, to which more attention should be given for the purpose of erosion and sediment control. A number of regression equations were established, based which the effect of rainfall, sand–dust storms and surface material grain size on specific sediment yield can be assessed. Copyright © 2007 John Wiley & Sons, Ltd.     

2004年苏门答腊大地震激发的T-波

T-波是由海底地震或者海陆边界俯冲带附近地震激发,并在海洋低速层中传播的声波.2004年12月26日,在印度洋东部印尼苏门答腊岛附近发生MW=9.3级大地震,其产生的能量在印度洋中激发了巨大的海啸,造成了严重的人员伤亡和财产损失,受到了世界科学家们极大的关注.本文从台站(PALK)及台站(DGAR)记录到的地震的信号中,提取出了清晰的高频T-波,并在频率域内分析,最终得到了T-波的频谱已及频率随时间变化图像.另外,通过对大地震时间域和频率域内T-波信号的分析,了解到此次大地震断层破裂过程持续的时间大致为500 s,其间伴随有两次明显的能量释放过程.分析数据表明两次能量释放过程的间隔大致为80~100 s.T-波分析将为推断海洋地震以及海陆边界俯冲带附近地震的特征,提供一种独立的研究手段和方法.     

Observation and study on the whole process of cloud-to-ground lightning using narrowband radio interferometer

A narrowband radio interferometer has been developed and used to locate the entire sources of VHF radiations from a negative cloud-to-ground (CG) lightning discharge which contains 19 strokes. This system uses five antennas to form an array consisting of short- and long-baselines along two or- thogonal directions. The system error which comes from frequency conversion is reduced by phase detection through direct high frequency amplifying. An interactive graphic analysis procedure is used to remove the fringe ambiguities which exist inherently in interferometry and to determine the direction of lightning radiation sources in two dimensions (azimuth and elevation) as a function of time at a time resolution of microsecond orders. With the developed system, the whole progression process in time and space of a lightning flash can be reconstructed. In this paper, combining the synchronous data of electric filed change and VHF radiation, the whole processes of an example negative CG flash have been studied in detail. It is found that the preliminary breakdown event of the CG flash started from negative charge region and exhibited firstly a downward pregression and then an upward propagation. There were very intense and continuous radiations during stepped leaders which became much stronger when the first return stroke began. In contrast, there were less and only discrete radiations during dart leaders. Stepped leader and dart leader may transform to each other depending on the state of the ionization of the path. The progression speed of initial stepped leaders was about 105 ms?1, while that was about 4.1×106 and 6.0×106 ms?1 for dart leaders and dart-stepped leaders, respectively. M events produced hook-shaped field changes accompanied by active burst of radiations at their begin- nings. Followed these active radiation processes, M events appeared to contact finally into conducting main discharge channels. The mean progression speed of M events was about 7×107 ms?1, greater than that of the dart leaders and dart-step leaders. K events and attempted leaders were essentially the same as dart leaders except that they could not reach the ground and initiate return strokes.     

CRG桩在花岗岩持力层中的试验研究

CRG桩(预应力高强混凝土管注浆桩)取得了国家发明专利。介绍了CRG桩施工方法及试验桩施工情况,讨论了在花岗岩地层施工中遇到的孔内沉渣清理、管桩上浮及灌浆问题。对CRG桩的灌浆充填体进行了钻芯检测,发现充填体胶结质量不高,充填体强度低,影响桩身强度的提高。但静载试验表明,CRG桩的侧阻力和端阻力都得到充分发挥,单桩承载力比同直径的钻孔灌注桩有较大的提高。为在花岗岩持力层中CRG桩的试验积累了经验。     

改革开放以来中国城市等级体系演变的基本特征

选取372个设市城市作为样本,在构建1987年、2005年城市规模数据集的基础上,对比分析了城市规模位序变动的特征,以此刻画近30年来中国城市等级体系演变的基本轮廓.研究发现:城市规模序列呈现大规模、大幅度的变动;城市位序变动随城市规模等级呈现分异特征,城市规模等级越小,城市位序变动越明显;城市规模增长还表现为明显的区域差异性.对变动幅度较大的138个城市进行典型观察,结果显示:成长明显的城市具有地域集中性,在交通沿线和沿海地区尤为突出,城市规模位序的提升与地区经济发展密切相关,部分旅游城市上升明显,而资源型城市则明显下降,城市的成长受行政性调整力量影响较大.     

The complex nonlinear systems with fractal as well as chaotic dynamics of annual runoff processes in the three headwaters of the Tarim River

This paper attempted to identify fractal and chaotic characteristics of the annual runoff processes in headwaters of the Tarim River. Methods of fractal analyses were used to explore several aspects of the temporal changes from 1957 to 2002. The main findings are as follows: (1) The annual runoff processes of the three headwaters of the Tarim River are complex nonlinear systems with fractal as well as chaotic dynamics. (2) The correlation dimensions of attractor derived from the time series of the annual runoff for the Hotan, Yarkand and Aksu rivers are all greater than 3.0 and non-integral, implying that all three rivers are chaotic dynamical systems that are sensitive to initial conditions, and the dynamic modeling of their annual runoff process requires at least four independent variables. (3) The time series of annual runoff in each river presents a long-term correlation characteristic. The Hurst exponent for the period of 1989 to 2002 suggests that we may expect to see an increasing trend in the annual runoff of the Aksu and Yarkand rivers in the years after 2002, but a decreasing tendency for the Hotan River in the same period.     

干旱区土壤植被系统恢复的生态水文学研究进展

干旱区约占地球陆表面积的40%以上,所支撑的复杂多样的生态系统是全球陆地生态系统的重要组成,水是干旱区诸多生态系统过程的驱动力和关键的非生物限制因子。生态水文学作为新兴的交叉学科,为客观、全面地诠释干旱区植被与土壤系统相互作用与反馈机理提供了新的理念和途径,同时,也为干旱区以植被建设为主的生态恢复实践提供了理论基础。以沙坡头沙漠研究试验站在腾格里沙漠东南缘荒漠化草原和草原化荒漠过渡带的长期定位监测研究为例,分析、综述了我国在干旱区生态水文学研究领域所取得的主要成绩和相关进展,重点讨论了沙地植被从组成、结构和功能群演变等多层次上对土壤生态水文过程的响应,以及植被对土壤环境和水循环的长期反馈作用和调控机理;将生物土壤结皮首次作为影响植被-土壤系统水文过程的重要环节来考虑,对其发挥的生态水文功能的相关研究进行了总结;归纳了基于水量平衡的干旱区生态恢复理论模式及其在实践中的应用。