Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

The Preboreal climate reversal and a subsequent solar‐forced climate shift

Accurate chronologies are essential for linking palaeoclimate archives. Carbon‐14 wiggle‐match dating was used to produce an accurate chronology for part of an early Holocene peat sequence from the Borchert (The Netherlands). Following the Younger Dryas–Preboreal transition, two climatic shifts could be inferred. Around 11 400 cal. yr BP the expansion of birch (Betula) forest was interrupted by a dry continental phase with dominantly open grassland vegetation, coeval with the PBO (Preboreal Oscillation), as observed in the GRIP ice core. At 11 250 cal. yr BP a sudden shift to a humid climate occurred. This second change appears to be contemporaneous with: (i) a sharp increase of atmospheric ¹⁴C; (ii) a temporary decline of atmospheric CO₂; and (iii) an increase in the GRIP ¹⁰Be flux. The close correspondence with excursions of cosmogenic nuclides points to a decline in solar activity, which may have forced the changes in climate and vegetation at around 11 250 cal. yr BP. Copyright © 2004 John Wiley & Sons, Ltd.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

Wind-driven interannual variability over the northeast Pacific Ocean

Interannual variability of the sea surface height (SSH) over the northeast Pacific Ocean is hindcast with a reduced-gravity, quasi-geostrophic model that includes linear damping. The model is forced with monthly Ekman pumping fields derived from the NCEP reanalysis wind stresses. The numerical solution is compared with SSH observations derived from satellite altimeter data and gridded at a lateral resolution of 1 degree. Provided that the reduced gravity parameter is chosen appropriately, the results demonstrate that the model has significant hindcast skill over interior regions of the basin, away from continental boundaries. A damping time scale of 2 to 3 years is close to optimal, although the hindcast skill is not strongly dependent on this parameter.A simplification of the quasi-geostrophic model is considered in which Rossby waves are eliminated, yielding a Markov model driven by local Ekman pumping. The results approximately reproduce the hindcast skill of the more complete quasi-geostrophic model and indicate that the interannual SSH variability is dominated by the local response to wind forcing. There is a close correspondence the two leading empirical orthogonal modes of the local model and those of the observed SSH anomalies. The latter account for over half of the variance of the interannual signal over the region.     

潮汐对瞬态气候响应模拟的影响研究

In this study, the impact of oceanic processes on the sensitivity of transient climate change is investigated using two sets of coupled experiments with and without tidal forcing, which are termed Exp_Tide and Exp_Control,respectively. After introducing tidal forcing, the transient climate response(TCR) decreases from 2.32 K to 1.90 K,and the surface air temperature warming at high latitudes decreases by 29%. Large ocean heat uptake efficiency and heat storage can explain the low TCR in Exp_Tide. Approximately 21% more heat is stored in the ocean in Exp_Tide(1.10×10~(24) J) than in Exp_Control(0.91×10~(24) J). Most of the large ocean warming occurs in the upper 1 000 m between 60°S and 60°N, primarily in the Atlantic and Southern Oceans. This ocean warming is closely related to the Atlantic Meridional Overturning Circulation(AMOC). The initial transport at mid-and high latitudes and the decline in the AMOC observed in Exp_Tide are both larger than those observed in Exp_Control. The spatial structures of AMOC are also different with and without tidal forcing in present experiments. The AMOC in Exp_Tide has a large northward extension. We also investigated the relationship between AMOC and TCR suggested by previous studies using the present experiments.     

大洋碳循环与气候演变的热带驱动

20世纪气候演变研究的最大突破,在于地球轨道变化驱动冰期旋回的米兰柯维奇理论。然而近年来学术界对热带过程和大气CO2浓度变化的研究进展,暴露了传统的轨道驱动理论存在着对低纬区和碳循环在全球气候系统中作用估计不足的严重缺陷。国家重点基础研究发展计划项目"大洋碳循环与气候演变的热带驱动"拟以南海与西太平洋暖池的深海记录为依据,进行全球性对比和跨越地球圈层的探索,通过观测分析结果与数值模拟的结合、地质记录与现代过程的结合,检验和论证大洋碳储库长周期变化机制的假说,对于不同时间尺度上低纬过程如何通过碳循环在全球气候环境演变中的作用,实现理论上的突破。同时简要介绍了该项目的目的、科学意义、关键科学问题及预期目标等。     

Income distribution dynamics and cross-region convergence in Europe

This paper presents a continuous version of the model of distribution dynamics to analyse the transition dynamics and implied long-run behaviour of the EU-27 NUTS-2 regions over the period 1995–2003. It departs from previous research in two respects: first, by introducing kernel estimation and three-dimensional stacked conditional density plots as well as highest density regions plots for the visualisation of the transition function, based on Hyndman et al. (J Comput Graph Stat 5(4):315–336, 1996), and second, by combining Getis’ spatial filtering view with kernel estimation to explicitly account for the spatial dimension of the growth process. The results of the analysis indicate a very slow catching-up of the poorest regions with the richer ones, a process of shifting away of a small group of very rich regions, and highlight the importance of geography in understanding regional income distribution dynamics.

Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.     

Late Holocene linkages between decade–century scale climate variability and productivity at Lake Tanganyika, Africa

Microlaminated sediment cores from the Kalya slope region of Lake Tanganyika provide a near-annually resolved paleoclimate record between ∼∼2,840 and 1,420 cal. yr B.P. demonstrating strong linkages between climate variability and lacustrine productivity. Laminae couplets comprise dark, terrigenous-dominated half couplets, interpreted as low density underflows deposited from riverine sources during the rainy season, alternating with light, planktonic diatomaceous ooze, with little terrigenous component, interpreted as windy/dry season deposits. Laminated portions of the studied cores consist of conspicuous dark and light colored bundles of laminae couplets. Light and dark bundles alternate at decadal time scales. Within dark bundles, both light and dark half couplets are significantly thinner than within light bundles, implying slower sediment accumulation rates during both seasons over those intervals.Time series analyses of laminae thickness patterns demonstrate significant periodicities at interannual–centennial time scales. Longer time scale periodicities (multidecadal to centennial scale) of light and dark half couplet thicknesses are coherent and in some cases are similar to solar cycle periods on these time scales. Although laminae thickness cycles do not strongly covary with the actual Δ¹⁴C record for this same time period, two large Δ¹⁴C anomalies are associated with substantial decreases in both light and dark laminae thickness. In contrast to the multidecadal– centennial time scale, significant annual to decadal periodicities, which are broadly consistent with ENSO/PDO forcing and their impact on East African climate, are not coherent between light and dark half couplets. The coherency of light–dark couplets at decadal–centennial time scales, but not at shorter time scales, is consistent with a model of a long-term relationship between precipitation (recorded in wet season dark laminae thickness) and productivity (light laminae thickness), which is not manifest at shorter time scales. We hypothesize that this coupling results from long-term recharging of internal nutrient loading during wet periods (higher erosion of soil P) and reduced loading during drought intervals. The relationship is not expressed on short time scales during which the dominant control on productivity is wind-driven, dry season upwelling, which is uncorrelated with wet-season precipitation. Our record greatly extends the temporal record of this quasi-periodic behavior throughout the late Holocene and provides the first evidence linking decade- to century-scale episodes of enhanced productivity to enhanced precipitation levels and nutrient recharge in a productive tropical lake.