Numerical simulations of very large impacts on the Earth

Vertical impacts on the Earth of asteroids 500-3000 km in diameter at 15 km/s have been numerically modelled using the hydrodynamic SOVA code. This code has been modified for the spherical system of coordinates well suited for simulations of very large impacts when the entire Earth is involved in motion. The simulations include cratering process, upward motion of deep mantle layers, fall of ejecta on the Earth, escape of matter to space, and formation of rock vapour atmospheres. The calculations were made for the period preceding disappearance of rock vapour atmospheres caused by radiation several years after the largest impacts. For very large vertical impacts at 15 km/s, escaping masses proved to be negligibly small. Quantities of kinetic, internal, potential, and radiated away energies are obtained as functions of time and space. After the impacts, a global layer of condensed ejecta covers the whole of the Earth's surface and the ejecta energy is sufficient to vaporise an ocean 3 km deep. The mass of rock vapour atmosphere is 10-23% of the impactor mass. This atmosphere has a greater mass than the water atmosphere if impactor is 2000 km in diameter or larger.     

Biodegradation and transport of benzene,toluene, and xylenes in a simulated aquifer: comparison of modelled and experimental results

Both laboratory experiments and numerical modelling were conducted to study the biodegradation and transport of benzene–toluene–xylenes (BTX) in a simulated semi‐confined aquifer. The factors incorporated into the numerical model include advection, hydrodynamic dispersion, adsorption, and biodegradation. The various physico‐chemical parameters required by the numerical model were measured experimentally. In the experimental portion of the study, BTX compounds were introduced into the aquifer sand. After the contaminants had been transported through the system, BTX concentrations were measured at 12 equally spaced wells. Subsequently, microorganisms obtained from the activated sludge of a sewage treatment plant and cultured in BTX mixtures were introduced into the aquifer through the 12 sampling wells. The distribution data for BTX adsorption by the aquifer sand form a nonlinear isotherm. The degree of adsorption by the sand varies, depending on the composition of the solute. The degradation time, measured from the time since the bacteria were added to the aquifer until a specific contaminant was no longer detectable, was 35–42 h for BTX. The dissolved oxygen, after degradation by BTX compounds and bacteria, was consumed by about 40–60% in the entire simulated aquifer; thus the aerobic conditions were maintained. This study provides insights for the biodegradation and transport of BTX in aquifers by numerical modelling and laboratory experiments. Experimental and numerical comparisons indicate that the results by Monod degradation kinetics are more accurate than those by the first‐order degradation kinetics. Copyright © 2002 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

形态参数下地表采动裂缝三维建模及可视化方法

实现对地表采动裂缝的三维建模及可视化能促进裂缝的形状、成因研究,并为开采沉陷监测和生态环境恢复提供科学支持。目前三维建模方法和软件主要面向建筑物、管道、地质体等地物地貌,若直接应用于地表裂缝的建模与表达会存在野外测绘工作量大、形态规律表达不准确、呈现不逼真等问题。本文结合三维地形建模、分形、空间插值等理论,提出了基于几何形态参数的地表裂缝三维建模方法,并利用ArcEngine平台实现了对裂缝的三维可视化;选取甘肃省东峡矿区为试验区,根据裂缝分布与形态预计参数,实现了对区域内裂缝的模拟,并通过与实测数据和高分辨率影像进行对比,验证了该方法能真实逼真呈现裂缝的形态、延伸与细节信息。     

交叉裂隙岩体裂纹扩展试验及混合有限-离散元数值模拟研究

对预制交叉裂隙岩石试样进行裂纹扩展试验,研究了裂纹萌生、扩展、聚合过程,分析了主裂隙和轴向载荷夹角及主次裂隙夹角对裂纹起裂应力和聚合应力的影响,并用混合有限元-离散元程序,即图形处理器并行化的3D Y-HFDEM代码对试验进行了仿真计算,实现了岩石破坏从连续介质向非连续介质的过渡,对裂纹类及损伤破坏模式进行了识别,捕捉到了试验中难以发现的现象。研究表明：随主裂隙与轴向载荷夹角增加,裂纹聚合区的拉伸裂纹数量增加;裂纹起裂和聚合应力与主裂隙与轴向载荷夹角成正比;主次裂隙夹角增加,岩石的破坏模式由拉伸破坏转为剪切破坏,交叉裂隙加剧岩石破碎程度;主裂隙尖端萌生扩展的拉伸-剪切混合裂缝引起的破坏在岩石破坏中占主导地位,是导致岩体失去承载能力的主控裂纹;混合有限元-离散元仿真软件GPGPU并行化的3D Y-HFDEM IDE在岩石裂纹扩展研究中具有优势,可以捕捉实验室难以发现的损伤断裂类型,可以作为岩石裂纹扩展研究的有力工具。     

Maximized erosion at the last glacial maximum revealed by thermoluminescence thermochronometry

Increased erosion rates during the last 2 Myr of Earth history are commonly associated with widespread glaciation and global cooling. However, whether erosion rates actually increased during that period of time remains debated. This is in part because the respective role of fluvial and glacial erosion is difficult to assess, and because many existing techniques used to estimate erosion rates may be affected by averaging biases that may produce increased erosion rates towards the present. Here, we apply thermoluminescence thermochronometry to granitic bedrock samples of the southern-central Chilean Andes, a region which experienced a full transition from glacial-humid to fluvial-dry conditions during the Quaternary. Inverse modelling of the thermochronological data enables us to extract unbiased time series of erosion rates during the last ∼100 kyr. We find that erosion rates were systematically higher during around the last glacial maximum than at present, with differences exceeding an order of magnitude. Our findings indicate that glacial erosion is probably more efficient than fluvial erosion and imply that erosion has remained transient during the Quaternary.     

Seismic performance of two-story subway station structures with different isolating systems

When the sliding bearing is fixed only at the top of the middle column of the underground structure, the cracks at the side end of the middle plate should be aggravated while the seismic damage of the mid-column should be alleviated. To enhance the seismic performance of the mid-plate, a new isolation design method has been mentioned while the elastic sliding bearings are set at the top of the mid-columns and between the side end of the mid-plate and the side wall at the same time. By establishing a nonlinear finite element analysis model for the static-dynamic coupling interaction system, the seismic response characteristics of the cast-in-place station structure without a sliding bearing have been analyzed and compared with those of the station structure with the sliding bearing fixed only at the top of the middle columns, and those of the station structure with sliding bearing be fixed between the mid-plate and the sidewall at the same time. The results show that the new isolation station structures suffer fewer earthquake damages at the mid-plate and mid-columns at the same time, which can improve the overall seismic performance of the subway station structure.     

Tracing exhumation record in high-pressure micaschists: A new tectonometamorphic model of the evolution of the eastern part of the Fore Sudetic Block,Kamieniec Metamorphic Belt,NE Bohemian Massif,SW Poland

The Kamieniec Metamorphic Belt comprises a volcano-sedimentary succession exposed within a collision zone between the Saxothuringian and Brunovistulian crustal domains of the European Variscides. The studied rocks recorded two metamorphic episodes. The first episode, M₁, occurred at conditions of c. 485 ± 25 °C and 18 ± 1.8 kbar related to burial within a subduction zone. The subsequent episode, M₂, was linked to the final phases of exhumation to mid-crustal level, associated with pressure and temperature (P–T) conditions ranging from c. 520 ± 26 °C and 6 ± 0.6 kbar through 555 ± 28 °C and 7 kbar ± 0.7 to ~590 ± 30 °C and 3–4 ± 0.4 kbar. The documented deformation record is ascribed to three events, D₁ to D₃, interpreted as related to the burial and subsequent exhumation of the Kamieniec Metamorphic Belt. The D₁ event must have witnessed the subduction of the Kamieniec Metamorphic Belt rock succession whereas the D₂ event was associated with the exhumation and folding of the Kamieniec Metamorphic Belt in an E-W-directed shortening regime. A subsequent folding related to the D₂ event was initiated at HP conditions, however, the planar fabric produced during a late stage of the D₂ event, defined by a low-pressure mineral assemblage M₂, indicates that the D₂ final stage was synchronous with the onset of the M₂ episode. Consequently, the entire D₂ event seems to have been associated with the exhumation of the Kamieniec Metamorphic Belt to mid crustal level. The third deformation event D₃, synchronous with the M₂ episode, marked the last stage of the exhumation, and was linked to emplacement of granitoid veins and lenses. The latter resulted in heating and rheological weakening of the entire rock succession and in the formation of non-coaxial shear zones.     

Grampian migmatites in the Buchan Block,NE Scotland

Rocks exposed along the Scottish coast between Fraserburgh and Inzie Head contain information critical to understanding the evolution of the Buchan Block, the type locality for low‐P, high‐T regional metamorphism, and its relationship with the rest of the Grampian terrane, one of the major tectonostratigraphic components of the Scottish Caledonides. The ~8 km long section traverses a regional network of shear zones and, at the highest grades around Inzie Head, passes into the core of the Buchan Anticline, a large‐scale open fold that is commonly regarded as a late structure, post‐dating metamorphism. The metasedimentary rocks increase in grade from upper amphibolite to granulite facies and preserve unequivocal evidence for partial melting. The diatexite migmatites around Inzie Head, along with other gneissose units within the Buchan Block, have been regarded as allochthonous Precambrian basement rocks that were thrust into their current position during the Grampian orogenesis. However, field observations show that the onset of in situ partial melting in metapelitic rocks, which was associated with the formation of garnet‐bearing aplites and associated pegmatites, occurred around Fraserburgh, where shear fabrics are absent. Thus, the rocks preserve a continuous metamorphic field gradient that straddles the shear zone network. This observation supports an alternative interpretation that anatexis was the result of mid‐Ordovician (Grampian) metamorphism, rather than an older tectonothermal event, and that the Inzie Head gneisses are autochthonous. Using an average mid‐Dalradian pelite as a plausible representative protolith, phase equilibria modelling satisfactorily reproduces the observed appearance and disappearance of key minerals providing that peritectic garnet produced with the first formed melts (represented by the garnet‐bearing aplites) depleted the source rocks in Mn. The modelled metamorphic field gradient records a temperature increase of at least 150 °C (from ~650 °C near Fraserburgh to in excess of 800 °C at Inzie Head) but is isobaric at pressures of 2.7–2.8 kbar, suggesting the Buchan Anticline developed synchronous with partial melting. The Buchan Anticline is likely an expression of crustal thinning and asthenospheric upwelling, which produced voluminous gabbroic intrusions that supplied the heat for Buchan metamorphism.