Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model

We apply a coupled thermomechanical ice sheet—self-gravitating viscoelastic solid Earth model (SGVEM), allowing for the dynamic exchange of ice thickness and bedrock deformation, in order to investigate the effect of viscoelastic deformation on ice dynamics and vice versa. In a synthetic glaciation scenario, we investigate the interaction between the ice sheet and the solid Earth deformation, the glacial-isostatic adjustment (GIA), accounting for an atmospheric forcing depending on the ice sheet surface altitude. We compare the results from the coupled model to runs with the common elastic lithosphere/relaxing asthenosphere (ELRA) model, where the lithosphere is represented by a thin plate and the mantle relaxes with one characteristic relaxation time, as well as to a rigid Earth without any deformation. We find that the deformational behaviour of the SGVEM on ice dynamics (i.e. stored ice volume, ice thickness and velocity field) is comparable to the ELRA for an optimal choice of the parameters in steady state, but exhibits differences in the transient behaviour. Beyond the ice sheet, in the region of peripheral forebulge, the differences in the transient surface deformation between ELRA and SGVEM are substantial, demonstrating the inadequacy of the ELRA model for interpreting constraints on GIA in the periphery of the ice sheet, such as sea-level indicators and GPS uplift rates.     

中国日降雨量变化特征分析

将旋转主成分分析和复变量莫莱特小波分析应用于14a中国境内实测与ECMWF再分析日降雨量资料,研究其年际、季际与季节内变化的时空分布特征和再分析日降雨量资料的有效性。结果表明:从实测资料中分析的降水特征,都可以由再分析资料再现出来。除周内振荡大有出入外,其余频率的振荡二者大致相符或基本吻合。因此,周内振荡还不能被再分析数据较好地刻划,但是与旱涝灾害和季风爆发联系密切的周际以上振荡,在再分析日降雨量中基本上能被正确地表现出来。     

Indication for precession of comet Halley's nucleus from dust tail analysis

The distinct structures ta comet Halley's dust tail around the beginning of March 1986 are analyzed by means of a computer simulation based on nucleus data obtained by the Giotto mission. It is shown that the assumption of a considerable free precession is required to understand the ground based dust tail observations from that time supposing a rotational period of some 50 hr. But a precession-free rotation with a period of about 7 days does not contradict an analysis of the dust tail structure. In both cases, an asymmetric distribution of the relevant emission sources is required.     

Projection of Future Precipitation Extremes Change (2001-2050) in the Yangtze River Basin

 Daily maximum rainfall (R1D) was higher in the Jialing River basin, the Taihu Lake area and the mid-lower main stream section of the Yangtze River basin in the 1990s, and there was a good relationship between ECHAM5/MPI-OM model simulation and the observed data about extreme precipitation (R1D). Under the IPCC SRES A2, A1B, and B1 scenarios, R1Ds are all projected to be in increasing trends in the upper Yangtze River basin during 2001-2050, and R1D shows a more significant increasing tendency under the A2 scenario when compared with the A1B scenario before 2020. With respect to the middle and lower Yangtze River basin, an increasing tendency is projected before 2025, and since then the increasing tendency will become insignificant. There might be more floods to the south of the Yangtze River and more droughts to the north in the next decades.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.     

Empirical assessment of rockfall and debris flow risk along the Karakoram Highway,Pakistan

Natural Hazards - The Karakoram Highway links north Pakistan with southwest China. It passes through unique geomorphological, geological and tectonic setting. This study focused 200-km-long section...     

Sequence stratigraphy and architecture of a late Early–Middle Aptian carbonate platform succession from the western Maestrat Basin (Iberian Chain, Spain)

The attributes of a ‘four-systems-tract’ sequence are at times difficult to identify in outcrop-scale carbonate successions. Poor exposure conditions, variable rates of sediment production, erosion and/or superposition of surfaces that are intrinsic to the nature of carbonate systems frequently conceal or remove its physical features. The late Early–Middle Aptian platform carbonates of the western Maestrat Basin (Iberian Chain, Spain) display facies heterogeneity enabling platform, platform-margin and slope geometries to be identified, and provide a case study that shows all the characteristics of a quintessential four systems tract-based sequence. Five differentiated systems tracts belonging to two distinct depositional sequences can be recognized: the Highstand Systems Tract (HST) and Forced Regressive Wedge Systems Tract (FRWST) of Depositional Sequence A; and the Lowstand Prograding Wedge Systems Tract (LPWST), Transgressive Systems Tract (TST) and subsequent return to a highstand stage of sea-level (HST) of Depositional Sequence B. An extensive carbonate platform of rudists and corals stacked in a prograding pattern marks the first HST. The FRWST is constituted by a detached, slightly cross-bedded calcarenite situated at the toe of the slope in a basinal position. The LPWST is characterized by a small carbonate platform of rudists and corals downlapping over the FRWST and onlapping landwards. The TST exhibits platform backstepping and marly sedimentation. Resumed carbonate production in shelf and slope settings characterizes the second HST. A basal surface of forced regression, a subaerial unconformity, a correlative conformity, a transgressive surface and a maximum flooding surface bound these systems tracts, and are well documented and widely mappable across the platform-to-basin transition area analyzed. Moreover, the sedimentary succession studied is made up of four types of parasequence that constitute stratigraphic units deposited within a higher-frequency sea-level cyclicity. Ten lithofacies associations form these basic accretional units. Each facies assemblage can be ascribed to an inferred depositional environment in terms of bathymetry, hydrodynamic conditions and trophic level. The architecture of the carbonate platform systems reflects a flat-topped non-rimmed depositional profile. Furthermore, these carbonate shelves are interpreted as having been formed in low hydrodynamic conditions. The long-term relative fall in sea-level occurred during the uppermost Early Aptian, which subaerially exposed the carbonate platform established during the first HST and resulted in the deposition of the FRWST, is interpreted as one of global significance. Moreover, a possible relationship between this widespread sea-level drop and glacio-eustasy seems plausible, and could be linked to the cooling event proposed in the literature for the late Early Aptian. Because of the important implications in sequence stratigraphy of this study, the sedimentary succession analyzed herein could serve as an analogue for the application of the four-systems-tract sequence stratigraphic methodology to carbonate systems.