Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Historic Dead Sea Level Fluctuations Calibrated with Geological and Archaeological Evidence

The Dead Sea, the Holocene terminal lake of the Jordan River catchment, has fluctuated during its history in response to climatic change. Biblical records, calibrated by radiocarbon-dated geological and archaeological evidence, reinforce and add detail to the chronology of the lake-level fluctuations. There are three historically documented phases of the Dead Sea in the Biblical record: low lake levels ca. 2000–1500 B.C.E. (before common era); high lake levels ca. 1500–1200 B.C.E.; and low lake levels between ca. 1000 and 700 B.C.E. The Biblical evidence indicates that during the dry periods the southern basin of the Dead Sea was completely dry, a fact that was not clear from the geological and archaeological data alone.     

Most likely locations of large earthquakes in the Kanto and Tokai areas, Japan, based on the local recurrence times

We examined the hypothesis that minima in local recurrence time, T_L, or equivalently maxima in local probability, P_L, may map asperities in the Kanto and Tokai areas of Japan, where the earthquake catalog of the National Research Institute for Earth Science and Disaster Prevention (NIED) is complete at the M=1.5 (M1.5) level. We mapped T_L (P_L) based on the a- and b-values of the nearest earthquakes within 20 km of every node of a grid spaced 0.01° for M7 target events. Only earthquakes within the top 33 km were used. The b-values increase strongly with depth, in several areas. Therefore, some of the T_L (P_L) anomalies are not revealed if data from the entire crustal seismogenic zone are mixed. Thus, we mapped T_L (P_L) separately for the top 15 km and the rest of the depth range, as well as for the entire seismogenic crust. The resulting T_L- and P_L-maps show that approximately 12% of the total area shows anomalously short recurrence times. Out of six shallow target events with M≥6.5 and which occurred since 1890, five are located within the anomalous areas with T_L <450 years. We interpret this to mean that areas with anomalously short T_L map asperities, which are more likely than other areas to generate future target events. The probability that this result is due to chance is vanishingly small. The great Kanto rupture of 1923 appears to have initiated in the most significant asperity we mapped in the study area. One anomaly is located in the northeastern part of the area of the proposed future rupture of the Tokai earthquake, and another one at its southwestern corner. The absolute values of T_L calculated are uncertain because they depend on the size of the volume used for the calculation.     

Modelling of meteorological conditions at an urban scale for the PUMA winter campaign

For the winter 2000 campaign of the Pollution of Urban Midlands Atmosphere project, observation and numerical modelling of meteorological conditions over the West Midlands conurbation, UK, was undertaken. Modelling was performed using the regional atmospheric meteorological system (RAMS). This paper presents a comparison of modelled and observed wind and temperature for 25 and 26 January 2000. The RAMS model uses two nested grids with a mesh size of 2 km for the inner grid which is embedded in the outer grid with a mesh size of 8 km. Statistical evaluation of the model results against the observational data of wind speed, direction and temperature at 10 m was conducted. In general, the modelling results are in a reasonable agreement with observation. The statistical evaluation suggests that model performance is poorer for the inner grid than the outer grid as the model uncertainties (mainly mean bias) transfer from the outer to inner one. The low indices of agreement of temperature and wind are mainly associated with the systematic root-mean-square-difference values. For temperature, the systematic bias may also be affected by representation of cloud amount by the model. For wind, the model tends to have a poor performance for calm conditions, as under a stable anti-cyclonic situation local wind patterns associated with topography may develop, although the topography of the region is relatively flat. The results for the inner grid reveal some subtle spatial patterns at a scale smaller than 10 km near hills and valleys with differences in elevation of a few hundred metres.     

Development and testing of numerical methods for two-way nested air pollution modelling

A new 3-D model REGINA (REGIonal high resolutioN Air pollution model) is under development at the National Environmental Research Institute (NERI). The model is based on pieces from several models developed over the last decades at NERI. The aim of the work is to develop a nested model which can operate with very high resolution in both space and time. To fulfill this aim the choice and implementation of accurate numerical methods is crucial. The model will be applied for studying air pollution phenomena (both monitoring, forecasting and scenarios) over Denmark. The present paper is focussed on the modification, implementation and testing of a numerical method for treating the horizontal advection in the model as well as the implementation of two-way nesting techniques. The horizontal transport in the model is solved using an accurate space derivative algorithm. This method is traditionally implemented with periodic boundary conditions, however, this is not an option for nested modelling. A new method for calculating non-periodic boundary conditions has been developed in order to overcome this problem. Extensive testing of the numerical solution of the advection and the coupling of the solution of advection and chemistry in the model using Molenkamp–Crowley rotation tests have been carried out. The results show that the model with the current implementation of numerical methods is suitable for calculating air pollution levels with high resolution.     

Nonlinear earthquake response analysis and energy calculation for seismic slit shear wall structures

Based on the concept of structural passive control, a new type of slit shear wall, with improved seismic performance when compared to an ordinary solid shear wall, was proposed by the authors in 1996. The idea has been verified by a series of pseudo-static and dynamic tests. In this paper a macro numerical model is developed for the wall element and the energy dissipation device. Then, nonlinear time history analysis is carried out for a 10-story slit shear wall model tested on a shaking table. Furthermore, the seismic input energy and the individual energy dissipated by the components are calculated by a method based on Newmark-β assumptions for this shear wall model, and the advantages of this shear wall are further demonstrated by the calculation results from the viewpoint of energy. Finally, according to the seismic damage criterion on the basis of plastic accumulative energy and maximum response, the optimal analysis is carried out to select design parameters for the energy dissipation device.     

DAS24-3B数据采集系统采集控制模块的工作原理

DAS24-3B是中国科学院地质与地球物理研究所研制的宽频带地震数据采集系统。采集控制模块是该系统的重要组成部分。本文阐述了控制模块的工作原理，介绍了DAS24-3B系统的基本组成和采集器工作过程，并给出该系统如何与笔记本端应用程序界面通讯和联络的操作建议，使该系统能够在工作中发挥更大的作用。     

连续小波变换在新疆地震活动周期分析中的应用

利用连续小波变换分析和研究了 1 970年 1月— 2 0 0 1年 1 2月新疆地区 30年来的地震活动情况。结果表明 ,新疆地震活动存在比较稳定的 1 0年左右活动周期 ,并存在 5年和 1 7年左右的不稳定的准活动周期 ,也就是说新疆地区地震活动既存在比较稳定的优势周期 ,也存在一定的时变性。根据 MS≥ 4.7地震时间序列的连续小波变换结果可推测 ,2 0 0 2— 2 0 0 5年新疆地区地震活动相对偏弱 ,中、强地震的发生次数偏少 ,有可能发生 5级或 6级地震 ,而到 2 0 0 6— 2 0 0 7年新疆有可能再次发生 7级左右的大地震。