Late Quaternary contourites and glaciomarine sedimentation in the Fram Strait

Two sites in the eastern Fram Strait, the Vestnesa Ridge and the Yermak Plateau, have been surveyed and sampled providing a depositional record over the last glacial‐interglacial cycle. The Fram Strait is the only deep‐water connection from the Arctic Ocean to the North Atlantic and contains a marine sediment record of both high latitude thermohaline flow and ice sheet interaction. On the Vestnesa Ridge, the western Svalbard margin, a sediment drift was identified in 1226 m of water. Gravity and multicores from the crest of the drift recovered turbidites and contourites. ¹⁴C dating indicates an age range of 8287 to 26 900 years BP (Early Holocene to Late Weichselian). The Yermak Plateau is characterized by slope sediments in 961 m of water. Gravity and multicores recovered contourites and hemipelagites. ¹⁴C ages were between 8615 and 46 437 years BP (Early Holocene to mid‐Weichselian). Downcore dinoflagellate cyst analyses from both sites provide a record of changing surface water conditions since the mid‐Weichselian, suggesting variable sea ice extent, productivity and polynyas present even during the Last Glacial Maximum. Four layers of ice‐rafted debris were also identified and correlated within the cores. These events occurred ca at 9, 24 to 25, 26 to 27 and 43 ka, asynchronous with Heinrich layers in the wider north‐east Atlantic and here interpreted as reflecting instability in the Svalbard/Barents Ice sheet and the northward advection of warm Atlantic water during the Late Weichselian. The activity of the ancestral West Spitsbergen Current is interpreted using mean sortable silt records from the cores. On the Vestnesa Ridge drift the modern mass accumulation rate, calculated using excess ²¹⁰Pb, is 0·076 g cm^?2 year^?1. On the Yermak Plateau slope the modern mass accumulation rate is 0·053 g cm^?2 year^?1.     

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高     

黄海溶解氧垂直分布最大值的成因

本文对黄海溶解氧垂直分布中最大值的成因进行了探讨．认为：浮游植物春花期光合作用产生的大量氧，是溶解氧垂直分布最大值得以在春季（5月）形成的基础；而良好的温、密跃层的存在是氧最大值得以形成的必要条件，夏季叶绿素最大值层中较强的光合作用，是氧最大值在夏季得以维持甚至增强的主要原因．提出：黄海春、夏季溶解氧垂直分布中的最大值并非主要由冬季保持而来．     

末次盛冰期以来渤、黄、东海陆架底质分布演变过程模拟研究

渤、黄、东海陆架底质的形成分布与末次盛冰期之后的海侵密切相关。末次盛冰期结束、海侵开始以来 ,潮流是渤、黄、东海陆架上的永久性主导作用应力。为从长期沉积动力演变过程的角度 ,探讨渤、黄、东海陆架底质形成分布的有关成因问题 ,利用数值模拟手段 ,再现了末次盛冰期以来 6个时期渤、黄、东海陆架潮流作用下海底的冲淤格局及底质分布。结果表明 ,扬子浅滩南侧东海外陆架的砂质沉积基本上是自 - 80 m海面以来形成的。扬子浅滩形成于 -5 2 m海面之后 ,至 - 3 0 m海面时已有一定规模 ,全新世最大海侵之后 ,逐渐形成现在规模的扬子浅滩。南黄海中部泥自 - 5 2 m海面时就已开始形成 ,- 3 0 m海面时范围很大 ,侵入北黄海 ,全新世最大海侵以来 ,逐渐调整到现在的范围。渤海中央泥、北黄海西部泥、浙闽岸外泥、辽东半岛西侧与北侧的砂质沉积、西朝鲜湾与江华湾中的砂质沉积以及苏北浅滩是自全新世最大海侵以来逐渐形成的。海州湾中砂质沉积形成的盛期在公元 8世纪之后。济洲岛西南泥、南黄海东部泥很可能分别形成于 - 3 0 m海面、- 5 2 m海面以来。全新世渤、黄、东海陆架底质分布的演变过程大致分为 2个阶段 :全新世最大海侵之前为渤、黄、东海陆架底质分布宏观格局的形成阶段 ;全新世最大海侵至今为渤     

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.     

GIS叠置图层方差分量的极大似然估计

针对GIS叠置中的同名点,以维希特分布密度为似然函数,提出了各图层方差分量的极大似然估计方法。该方法不依赖残差,不需要迭代就能估计未知参数和方差分量。     

山东海阳断裂东石兰沟段晚更新世以来地表断错特征与最大潜在地震估计

海阳断裂是胶东半岛NE向牟平 -即墨断裂带东部一条规模较大的断裂 ,尽管晚更新世以来该断裂的地表断错活动总体上已基本停息 ,但东石兰沟段在晚更新世晚期以来仍有断错地表的活动。最后一次断错地表的活动发生在距今 3 7～ 1 2万年 ,但接近 1 2万年。地表破裂长度约6 5km ,活动段长度 8km。地表断错以走滑活动为主 ,可见最大倾滑位移 0 2m ;根据断层擦痕侧伏角推测最大水平位移 1 13m。最后一次断错地表的活动若以距今 1 2万年计算 ,则最大平均倾滑速率为 0 0 17mm/a ;最大平均右旋走滑速率为 0 0 94mm/a。野外观测到该活动段的断错活动表现为突发断错 ,根据地震地表破裂参数、活动段长度与地震的关系 ,估计其最大潜在地震为 6 级