南海东北部末次冰期以来的沉积环境演变

利用南海东北部上陆坡处所取的D孔柱样进行硅藻分析,结合粒度与碎屑矿物分析结果,探讨该海域末次冰期以来沉积环境演变过程.结果表明：粒度与硅藻分析结果较为吻合,共同显示D孔柱样可分为上下两层.上表层（0～2 cm）的砂质沉积层为冰后期海进过程改造的晚更新世残留沉积层,即变余沉积;下层（2～130 cm）主要是末次冰期形成的浅海沉积.整个柱样可反映出海洋同位素MIS 4期海退至MIS 3期海进与MIS 2期海退至冰后期海进的沉积变化韵律.     

Surface Elevation Distribution of Sea Waves Based on the Maximum Entropy Principle

A probability density function of surface elevation is obtained through improvement of the method introduced by Cieslikiewicz who employed the maximum entropy principle to investigate the surface elevation distribution. The density function can be easily extended to higher order according to demand and is non-negative everywhere, satisfying the basic behavior of the probability. Moreover because the distribution is derived without any assumption about sea waves, it is found from comparison with several accepted distributions that the new form of distribution can be applied in a wider range of wave conditions. In addition, the density function can be used to fit some observed distributions of surface vertical acceleration although something remains unsolved.     

长江口枯季悬沙粒度与浓度之间的关系

2003和2004年枯季在长江口采集水样并作水文观测，对所获水样进行过滤和粒度分析，以计算悬沙浓度和悬沙粒度分布。结果表明，2003年11月小潮期间，悬沙中值粒径与悬沙浓度存在着显著的指数关系，在大潮期间没有显著关系；在2004年2月小湖期间，两者之间没有显著关系，但在大潮期间存在着显著的指数关系。枯季水体悬沙以粉砂组份为主，并且随着向口外的推移，细颗粒组份逐渐增加，但在拦门沙最大浑浊带附近，由于絮凝作用，沉积物粒度变幅较大，可产生粒径粗化的现象。小潮期间，砂含量较低，但与悬沙浓度之间有显著相关关系；大潮期间，悬沙粒径粗化，但砂含量与悬沙浓度之间的关系不显著。上述分布趋势与沉积物来源、当地的水动力条件和絮凝作用等因素有密切关系。     

电磁波信号在对流层中路径弯曲改正的级数展开式

大气折射延迟是空间大地测量中的一个主要误差源。本文严格地给出了中性大气路径弯曲改正的原理性公式。通过对各种因素量级的分析,在亚毫米量级的精度上给出了一个以视天顶距正切和正割乘积为参数的级数展开模型。新模型在球对称大气模型下,级数展开模型的精度达到了毫米量级,理论上可以展开到任意阶次,且可允许观测高度角小到5°。     

随机Dirichlet问题解的调和展开

针对重力学随机Dirichlet问题,通过适当地对边界检验函数的分解,并在随机边界样本空间中提取确定性部分的对偶基,本文将随机Dirichlet问题的一般解展开为一随机系数的调和级数形式。     

Seismic hazard of Egypt

Earthquake hazard parameters such as maximum expected magnitude,M _max, annual activity rate,, andb value of the Gutenberg-Richter relation have been evaluated for two regions of Egypt. The applied maximum likelihood method permits the combination of both historical and instrumental data. The catalogue used covers earthquakes with magnitude 3 from the time interval 320–1987. The uncertainties in magnitude estimates and threshold of completeness were taken into account. The hazard parameter determination is performed for two study areas. The first area, Gulf of Suez, has higher seismicity level than the second, all other active zones in Egypt.b-values of 1.2 ± 0.1 and 1.0 ± 0.1 are obtained for the two areas, respectively. The number of annually expected earthquakes with magnitude 3 is much larger in the Gulf of Suez, 39 ± 2 than in the other areas, 6.1 ± 0.5. The maximum expected magnitude is calculated to be 6.5 ± 0.4 for a time span of 209 years for the Gulf of Suez and 6.1 ± 0.3 for a time span of 1667 years for the remaining active areas in Egypt. Respective periods of 10 and 20 years were reported for earthquakes of magnitude 5.0 for the two subareas.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.     

Lake-Level Fluctuations and Small-Scale Vegetation Patterns During the Late Glacial in The Netherlands

Late glacial changes in the vegetation were studied in and around a former lake on the southeastern side of a coversand ridge near Milheeze (southern Netherlands). Analyses of microfossils and macroremains and AMS ¹⁴C dating were performed on four sediment cores along a transect from sand ridge to the lake centre. Small-scale vegetation patterns and lake-level fluctuations were reconstructed in detail based on the information provided by the transect. For the first time in The Netherlands, cores along a transect within one lake were used to reconstruct the amplitude of late glacial lake-level fluctuations. Near Milheeze, a small and shallow lake was formed during the Bølling. The large increase in the water level during the Bølling and early Allerød, and the transition to more eutrophic conditions at the start of the Allerød, were probably related to the disappearance of permafrost. During the Allerød, open birch and pine woodlands developed in the area. In the lake, organic deposits accumulated, and the lake size and depth fluctuated. At the start of the Younger Dryas, higher lake water levels were recorded and woodlands became more open as a result of both a drop in the temperature and an increase in the effective precipitation. During the late Younger Dryas the lake water level dropped as the climate became drier and temperatures slightly increased. Accumulation of organic deposits in the lake ceased at the end of the Younger Dryas, which was caused by a drop in the water level in combination with the hydroseral succession process within the lake itself. The climatic signal reflected in the late glacial flora and lake-level fluctuations agree well with other published data from The Netherlands.