利用MODIS fAPAR傅立叶时间序列分析研究植被光合作用活动对净辐射和降雨的响应：青藏高原个例研究

气候变化对植被动力学有非常大的影响。为了定量描述气候变化对植被的影响，文章利用MODIS fAPAR 数据和NCEP 的净辐射和降雨再分析数据对青藏高原地区气候变化对植被的影响进行了时间序列分析。研究所用的数据时间跨度为2000年至2005年。首先利用NCEP 再分析数据建立了干旱度因子的时间序列，为了与MODIS fAPAR 具有相同的时间采样间隔，由NCEP的日净辐射和日降雨量得到每8天的平均净辐射和8日降雨的和。根据一定时间间隔的净辐射与降雨量的比可以用来衡量相对于可利用水分的剩余能量，因此该比值也是干旱灾害的度量。其次，对MODIS fAPAR 的傅立叶时间序列分析提供了两个植被光合作用对干旱相应的因子，即fAPAR的年平均值及其年振幅值。在时间和空间尺度上对植被光合作用活动与干旱指数之间的关系进行了定量分析。对湿年和干年之间的响应差异进行了比较。研究表明较干地区对气候变化的响应最为显著。分析应该扩展到更长的时间跨度以便更加有效地在时间和空间尺度上评估气候变化对植被动力学的影响。     

Earthquake source model using strong motion displacement as response of finite elastic media

The strong motion displacement records available during an earthquake can be treated as the response of the earth as the a structural system to unknown forces acting at unknown locations. Thus, if the part of the earth participating in ground motion is modelled as a known finite elastic medium, one can attempt to model the source location and forces generated during an earthquake as an inverse problem in structural dynamics. Based on this analogy, a simple model for the basic earthquake source is proposed. The unknown source is assumed to be a sequence of impulses acting at locations yet to be found. These unknown impulses and their locations are found using the normal mode expansion along with a minimization of mean square error. The medium is assumed to be finite, elastic, homogeneous, layered and horizontal with a specific set of boundary conditions. Detailed results are obtained for Uttarkashi earthquake. The impulse locations exhibit a linear structure closely associated with the causative fault. The results obtained are shown to be in good agreement with reported values. The proposed engineering model is then used to simulate the acceleration time histories at a few recording stations. The earthquake source in terms of a sequence of impulses acting at different locations is applied on a 2D finite elastic medium and acceleration time histories are found using finite element methods. The synthetic accelerations obtained are in close match with the recorded accelerations.     

Measurement of vertical environmental-head profiles in unconfined sand aquifers using a multi-channel manometer board

The vertical hydraulic gradient in an unconfined sand aquifer frequently exceeds the horizontal gradient when measured close to a surface discharge boundary, around an abstraction bore or close to a saline interface. The vertical gradient can be efficiently investigated using a multi-channel manometer board connected to a bundled piezometer installed using hollow-stem-auger-drilling techniques. The design and construction of a manometer board is described. Manometer board data from a site in a 30-m-deep fresh-water aquifer adjacent to a pumping bore are described. The bore is screened in the lower part of the aquifer and the impacts of partial penetration are evident on the manometer board results. Data are then described from the centre of a coastal sand spit where fresh water overlies saline water. Fluid electrical conductivity measurements from the individual tubes of the bundled piezometer are used to derive fluid density at each of the bundled piezometer ports. The average fluid density of water above the mini-piezometer port at each depth is then calculated using Gaussian quadrature. The environmental-water head profiles are compared to point-water and fresh-water head profiles and the implications of the incorrect use of point-water or fresh-water head vertical profiles are clearly demonstrated.     

北京市永定河流域地下水^14C年龄的初步分析

应用同位素方法初步分析北京郊区永定河流域地下水的演化特点。沿永定河冲洪积扇地下水流动方向布置取样剖面，共有取样点14个，对采集的水样进行^14C和氚含量分析，并确定地下水同位素年龄。浅层孔隙水的^14C年龄的变化范围为730～4900a，深层孔隙水为13420-22480a；^14C年龄在垂直方向上由浅部至深部逐渐增大，最大变化幅度为从3010a增至22480a；浅层孔隙水的氚含量都在14．99～30．56TU之间，深层孔隙水大部分在0．51～4．71TU之间。运用地下水^14C和氚年龄在垂向和水平方向变化的结果，验证了地下水的流向并计算了地下水的流速变化范围为5．02～62．63m／a，从山前至平原浅层地下水径流速度逐渐变小，反映了地下水水平径流强度逐渐减弱，地下水交替逐渐变差，浅层孔隙水以垂向交替为主，深层孔隙水以水平径流为主。     

Mechanics of a sandy subsoil subjected to cyclic loadings

A simple model for compaction of a sand, that may be useful in various geotechnical applications is presented. The model has been formulated in terms of the cyclic stress and strain amplitudes. The compaction properties of a dry sand are characterized by a common compaction curve described by two coefficients. The second aim of this paper is to show some applications of the theory proposed. The attention has been restricted to the two extreme cases of practical importance, namely to the behaviour of a dry sand (or a saturated sand but in free draining conditions), and to the behaviour of saturated sand in undrained conditions. Some numerical algorithms showing how to deal with those problems are presented and illustrated on the following examples: settlement and pore pressure generation in a soil stratum subjected to an earthquake, settlement of a foundation, pore pressure generation and liquefaction in a sea-bed. The results obtained suggest that the model proposed, as well as a method of dealing with boundary value problems can serve as a useful tool for the analysis of a sandy subsoil subjected to cyclic loadings.     

The mechanics of a magma chamber-fault system in trans-tension with application to Coso

The mechanical interaction between an elliptically shaped magma chamber and a fault subject to transtension is investigated with particular reference to the Coso geothermal field. The geologic setting of the Coso field is interpreted as a releasing bend step-over structure formed by the Airport Lake and Owens Valley dextral strike-slip fault system. The role of the Coso volcano-magmatic center in the development of the “over-step” structure is examined by treating the magma chamber as a liquid inclusion in a viscoelastic crust containing a fault (Airport Lake). The problem is numerically solved using a 2D viscoelastic finite element model with thermally activated viscosity to account for thermal weakening of the rock. The temperature distribution around the magma body is calculated based on a 3D steady-state approach and using the mesh-less numerical method. The fault is modeled as a frictionless contact. The simulated distributions of stress and strain around the inclusion display a rotation caused by the shearing component of the applied transtension. The results indicate that the fault tends to overstep the chamber in a geometric pattern similar to a step-over. There is good correspondence between the computed distributions of the maximum shear stress in the vicinity of the magma chamber and the map of earthquake epicenters at a depth of 7–10 km in Coso.     

自适应滤波器

自适应滤波器是一个空间域的数字滤波器。它对不应作滤波处理的点不作滤波,使滤波过程对区域成分的畸变降到最小。它采用一元三点内插法作插值,使滤波公司简单,计算速度快而又能保证精度。判别局部成分的标志及滤波计算都只用局部成分前二个点及后三个点,因而边缘损失小。     

Burial and infilling of a karst in Papua New Guinea by road erosion sediments

The anthropogenic impact on karst in Papua New Guinea is briefly introduced and a specific case is presented detailing the effect of road erosion sediments on a small karst. The karst is in the perennially humid tropics and covered with primary rain forest. The road was placed high above the karst on steep friable rock and traverses several of its catchments. The changes to and the rate of burial of parts of the karst and the infilling of the caves are described. The karst drainage has altered, and there is increased water storage. The sediment build-up ceased in less than a year due to vegetation and stabilization of the road embankments. It is concluded that any construction within a catchment leading to a karst should be assessed as to its impact on the karst.     

Modelling fluid flow through a single fracture with different contacts using cellular automata

The representation of the surface of a rock fracture and a numerical method to simulate fluid flow in single fractures are the keys to understanding the hydraulic behaviour of rock fractures. In this paper, a cellular automaton (CA) approach is used to generate the single fracture structure, which is assumed to be composed of contacts and voids. We develop a CA evolution rule to produce a contact area, and randomly model a single rock fracture with different contact ratios to reflect natural fracture properties such as dead voids, islands and tortuous flow path. Then, based on the localisation theory of a CA, a numerical method to simulate fluid flow in single fractures with contacts is developed. In this method, the fracture is discretised into a system composed of cell elements. Different apertures, i.e., zero for contacts and non-zero for voids, are assigned to each cell element. Therefore, the contribution of the cell elements in a contact on a cell’s transmissivity can be ignored completely. The local transmissivity is assumed to conform to the cubic law. The fluid flow in a fracture with different contact situations is then modelled using the method established in this paper. The fluid flow path, flow velocity and fluid head distributions as well as the channel flow in the fracture are well-modelled. The flow behaviour of the fracture strongly depends on the effective fluid flow path.     

Design and simulation of a vertical broadband seismometer

In this paper, a new vertical broadband seismometer is designed, simulated, and evaluated. The proposed seismometer has a bandwidth from 40 s to 50 Hz and sensitivity to ground velocity of 1500 V’/m/s. The mechanical suspension used in the proposed design is a mass-leaf spring rotational pendulum. The displacement of the seismic mass is sensed by a laser interferometer and then forwarded to an electronic force feedback system which feeds a linear actuator to keep the mass almost fixed in its position. The current needed to fix the mass in its position is an indication of the ground acceleration. The performance of the simulated seismometer is evaluated by applying a data record of two actual earthquakes to the simulated model and comparing its response to the earthquakes with that of a real seismometer (Trillium-40) having the same bandwidth and sensitivity. The proposed seismometer has performed within 0.4 % residual variance for the local event and 0.03 % residual variance for the regional one compared to Trillium-40.     

Discrete fracture model for coupled flow and geomechanics

We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well.     

Clay fine fissuring monitoring using miniature geo-electrical resistivity arrays

This article describes a miniaturised electrical imaging (resistivity tomography) technique to map the cracking pattern of a clay model. The clay used was taken from a scaled flood embankment built to study the fine fissuring due to desiccation and breaching process in flooding conditions. The potential of using a miniature array of electrodes to follow the evolution of the vertical cracks and number them during the drying process was explored. The imaging technique generated two-dimensional contoured plots of the resistivity distribution within the model before and at different stages of the desiccation process. The change in resistivity associated with the widening of the cracks were monitored as a function of time. Experiments were also carried out using a selected conductive gel to slow down the transport process into the cracks to improve the scanning capabilities of the equipment. The main vertical clay fissuring network was obtained after inversion of the experimental resistivity measurements and validated by direct observations.     

Storm-influenced siliciclastic and carbonate ramp deposits, the Lower Ordovician Dumugol Formation, South Korea

The Dumugol Formation (Lower Ordovician) in the southern part of the Baegunsan syncline, South Korea, contains mixed siliciclastic and carbonate ramp deposits. The ramp sediments were frequently influenced by storm events resulting in tempestites of sandstone-mudstone couplets, bioclastic grainstones to packstones, flat-pebble conglomerates, a skeletal lag layer and laminated calcisiltites. All tempestites are characterized by an erosive to sharp base, poor grading and a transitional upper boundary. The difference in lithology of tempestites appears to have been controlled by the nature of substrates and by proximality. For example, laminated calcisiltites have developed on the shallow carbonate ramp, flat-pebble conglomerates are closely associated with nodular limestones on shallow and deep ramps, and thin skeletal lag layers from fossiliferous argillaceous sediments formed in a basinal setting. The stratigraphic succession of the Dumugol Formation represents an initial transgression followed by a regression. The vertical facies change records the transition from a shallow siliciclastic ramp to a deep carbonate ramp, to a basin, shallowing to a deep carbonate ramp, and to a shallow carbonate ramp. Storm effects are mostly well preserved in shallow to deep ramp deposits.     

Transportation and the recent development of the Queensland coal industry

The role of transportation in the development of the Queensland coal industry since 1955 is reviewed in this paper. It is found that transportation has played a concomitant rather than a prerequisite role in the industry's development. The commitment of the Japanese to Queensland coal has been the most important stimulus to the development. The paper ends with a view to future developments, concluding that it would be unwise for Queensland to develop a transportation infrastructure as a stimulus to expanding the coal industry.     

Evidence of multiple sources involved in the genesis of the neoproterozoic itapetim granitic complex,NE Brazil,based on geochemical and isotopic data

The Itapetim Complex is a multiple facies intrusion of porphyritic monzogranite hosting isolated, or swarms of, dioritic enclaves and cut by late dikes of biotite granodiorite. It is a syn-tectonic intrusion in relation to the D3 regional deformation phase of the Brasiliano Orogeny. The complex has contacts with metagreywackes and gneissic granites associated with the Mesoproterozoic Cariris Velhos event and, to the northwest, with a belt of gneisses of Paleoproterozoic ages. At least two different sources seem to be involved in the evolution of the Itapetim Complex. The diorite source appears to be a depleted mantle that underwent some small degree of mixture with a Paleoproterozoic enriched mafic crust. The monzogranite source is probably metagreywacke, a mixture of enriched crust of Paleoproterozoic age and volcanic rocks related to the Cariris Velhos event. The granodiorites were generated by partial melting of a source compositionally similar to those suggested for the monzogranite. U–Pb in zircon gave a crystallization age of 638±4.9 Ma for the porphyritic monzogranites and associated mafic rocks. The granodiorite dikes represent a late magmatic event, under intermediate fO₂ conditions. The dikes, intruded into brittle fractures, are probably related to 570 Ma intrusions within the Pajeú Paraı́ba Terrane.     

Limit analysis of 2-D and 3-D structures based on an ellipsoid yield criterion

In this paper, a nonlinear numerical technique is developed to calculate the limit load and failure mode of structures obeying an ellipsoid yield criterion by means of the kinematic limit theorem, nonlinear programming theory and displacement-based finite element method. Using an associated flow rule, a general yield criterion expressed by an ellipsoid equation can be directly introduced into the kinematic theorem of limit analysis. The yield surface is not linearized and instead a nonlinear purely kinematic formulation is obtained. The nonlinear formulation has a smaller number of constraints and requires less computational effort than a linear formulation. By applying the finite element method, the kinematic limit analysis with an ellipsoid yield criterion is formulated as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the dissipation power which is to be minimized and an upper bound to the plastic limit load of a structure can then be calculated by solving the minimum optimization problem. An effective, direct iterative algorithm has been developed to solve the resulting nonlinear programming formulation. The calculation is based purely on kinematically admissible velocities. The stress field does not need to be calculated and the failure mode of structures can be obtained. The proposed method can be used to calculate the bearing capacity of clay soils in a direct way. Some examples are given to illustrate the validity and effectiveness of the proposed method.     

变参数DCD算法在岩土工程位移反分析中的应用

常规DCD(dynamic canonical descent)算法具有全局优化能力且无需考虑目标函数的可微性,只要预先定义优化空间即可,但是该算法的收敛速度有限,为此提出了变参数DCD算法,并建立了其相应的算法迭代格式。在数值试验和工程应用中将该变参数DCD算法与常规的DCD算法进行比较,其结果均表明：变参数DCD算法在全局优化能力和收敛速度上找到了一个均衡点,该算法不仅具有DCD算法的全局优化能力,而且收敛时所需的目标函数评估次数少,在优化过程中该算法展示出了稳定性强且优化结果可靠度高的一面。     

An example of high-pressure low-temperature metamorphic rocks from an island-arc: the Paikon Series (Innermost Hellenides, Greece)

Abstract The Paikon Series is considered to be a volcanic arc sequence with a mainly neritic sedimentary sequence and bimodal tholeiitic volcanism of early Mesozoic age. The metamorphic assemblages are syn- to post-kinematic with respect to a pre-Tithonian tectonic phase and range from the lawsonite-chlorite-albite facies through transitional Na-amphibole-greenschist facies to the chlorite sub-zone of the greenschist facies. The metamorphic imprint of the Paikon Series corresponds to a temperature range from less than 330° C to ± 450° C under a total pressure from 3 kbar to 6–7 kbar. The overprinting of these facies on an earlier blueschist assemblage, related either to a subduction zone or to a tectonic overpressure caused by thrusting, is suspected.     

Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.     

并行组合模拟退火算法在边坡稳定分析中的应用

采用自适应有限元法分析边坡的应力场,然后应用并行组合模拟退火算法建立基于应力场的边坡任意形状临界滑动面的全局搜索方法。该方法适用于任意几何形状,不同岩土材料分层以及多种荷载作用的复杂边坡。并行组合模拟退火算法是将模拟退火算法和遗传算法相结合的优化算法,它可以从多个初始点开始并行寻优,能以较快的速度找到全局最优解。通过工程实例分析,证明这种基于自适应有限元分析和并行组合模拟退火算法搜索边坡临界滑动面的方法可行、高效。