Upscaling of hydraulic conductivity and telescopic mesh refinement   总被引：1，自引：0，他引：1  

Walker DD  Gylling B  Selroos JO 《Ground water》2005,43(1):40-51

Performance assessments of repositories for the underground disposal of nuclear fuel and waste include models of ground water flow and transport in the host rocks. Estimates of hydraulic conductivity, K, based on field measurements may require adjustment (upscaling) for use in numerical models, but the choice of upscaling approach can be complicated by the use of nested modeling, large-scale fracture zones, and a high degree of heterogeneity. Four approaches to upscaling K are examined using a reference case based on exhaustive site data and an application of nested modeling to evaluate performance assessment of a waste repository. The upscaling approaches are evaluated for their effects on the flow balance between nested modeling domains and on simple measures of repository performance. Of the upscaling approaches examined in this study, the greatest consistency of boundary flows was achieved using the observed scale dependence for the rock domains, measured values from the large-scale interference test for the conductor domain, and a semivariogram regularization based on the Moye model for packer test interpretation. Making the assumption that large fracture zones are two-dimensional media results in the greatest changes to the median of travel time and improves the flow balance between the nested models. The uncertainty of upscaling methods apparently has a small impact on median performance measures, but a significant impact on the variances and earliest arrival times.  相似文献   

Representing Grounding Line Dynamics in Numerical Ice Sheet Models: Recent Advances and Outlook     

David Docquier  Laura Perichon  Frank Pattyn 《Surveys in Geophysics》2011,32(4-5):417-435

Recent satellite observations of the Antarctic and Greenland ice sheets show accelerated ice flow and associated ice sheet thinning along coastal outlet glaciers in contact with the ocean. Both processes are the result of grounding line retreat due to melting at the grounding line (the grounding line is the contact of the ice sheet with the ocean, where it starts to float and forms an ice shelf or ice tongue). Such rapid ice loss is not yet included in large-scale ice sheet models used for IPCC projections, as most of the complex processes are poorly understood. Here we report on the state-of-the art of grounding line migration in marine ice sheets and address different ways in which grounding line migration can be attributed and represented in ice sheet models. Using one-dimensional ice flow models of the ice sheet/ice shelf system we carried out a number of sensitivity experiments with different spatial resolutions and stress approximations. These are verified with semi-analytical steady state solutions. Results show that, in large-scale finite-difference models, grounding line migration is dependent on the numerical treatment (e.g. staggered/non-staggered grid) and the level of physics involved (e.g. shallow-ice/shallow-shelf approximation).  相似文献   

Thermal effects of thrust faulting     

Jon Brewer 《Earth and Planetary Science Letters》1981

Calculations based on simple models of overthrust sheets in crystalline basement rocks show that significant thermal effects may result from their movements. If rates are sufficiently high (e.g. plate tectonic rates), the thrust sheets sufficiently thick (5, 10 and 15 km are modelled here), the distances moved sufficiently large, and for reasonable values of the coefficient of friction along the thrust plane overthrusting can cause metamorphic mineral zonations and heat flow anomalies observable in the field. Regions where large-scale overthrusting has occurred should be characterized by a decrease with depth of grade of metamorphic mineral assemblages and anomalously low heat flow. The theoretical effects are presented as a series of maximum temperature vs. depth and heat flow vs. time plots.  相似文献   

Death valley regional ground-water flow model calibration using optimal parameter estimation methods and geoscientific information systems     

Frank A. D'Agnese  Claudia C. Faunt  Mary C. Hill  A. Keith Turner 《Advances in water resources》1999,22(8):2555

A regional-scale, steady-state, saturated-zone ground-water flow model was constructed to evaluate potential regional ground-water flow in the vicinity of Yucca Mountain, Nevada. The model was limited to three layers in an effort to evaluate the characteristics governing large-scale subsurface flow. Geoscientific information systems (GSIS) were used to characterize the complex surface and subsurface hydrogeologic conditions of the area, and this characterization was used to construct likely conceptual models of the flow system. Subsurface properties in this system vary dramatically, producing high contrasts and abrupt contacts. This characteristic, combined with the large scale of the model, make zonation the logical choice for representing the hydraulic-conductivity distribution. Different conceptual models were evaluated using sensitivity analysis and were tested by using nonlinear regression to determine parameter values that are optimal, in that they provide the best match between the measured and simulated heads and flows. The different conceptual models were judged based both on the fit achieved to measured heads and spring flows, and the plausibility of the optimal parameter values. One of the conceptual models considered appears to represent the system most realistically. Any apparent model error is probably caused by the coarse vertical and horizontal discretization.  相似文献   

Variability of the topographic core-mantle torque calculated from core surface flow models     

S. Asari  H. Shimizu  H. Utada 《Physics of the Earth and Planetary Interiors》2006,154(1):85-111

With the prospect of studying the relevance of the topographic core-mantle coupling to the variations of the Earth’s rotation and also its applicability to constraining the core surface flow, we investigate the variability of the topographic torque estimated by using core surface flow models accompanied by (a) uncertainty due to the non-uniqueness problem in the flow inversion, and (b) variance originating in that of geomagnetic secular variation models employed in the inversion. Various flow models and their variances are estimated by inverting prescribed geomagnetic models at the epoch 1980. The subsequent topographic torque is then calculated by using a core-mantle boundary topography model obtained by seismic tomography. The calculated axial and equatorial torques are found subject to the variability of order 10¹⁹ and 10²⁰  Nm, respectively, on which (b) is more effective than (a). The variability of the torque is attributed even to (a) and (b) of the large-scale flows (degrees 2 and 3). Yet, it still seems unlikely for the decadal polar motion with the observed amplitude to be excited exclusively by the equatorial topographic torque associated with any of reasonable core surface flow models. It is also confirmed that, with the topography model adopted here, the axial topographic torque on a rigid annulus in the core (coaxial with the Earth’s rotation axis) associated with any of reasonable flow models is larger by two orders of magnitude than the plausible inertial torque on such cylinders. This implies that any core surface flow model consistent with the topographic coupling does not exist, unless the topography model is appropriately modified. Nevertheless, the topographic coupling might provide not only a weak constraint for explaining the decadal LOD variations, but also the possibility to probe the core surface flow and the core dynamics.  相似文献   

Transport in chemically and mechanically heterogeneous porous media. I: Theoretical development of region-averaged equations for slightly compressible single-phase flow     

《Advances in water resources》1996,19(1):29-47

Transport processes in heterogeneous porous media are often treated in terms of one-equation models. Such treatment assumes that the velocity, pressure, temperature, and concentration can be represented in terms of a single large-scale averaged quantity in regions having significantly different mechanical, thermal, and chemical properties. In this paper we explore the process of single-phase flow in a two-region model of heterogeneous porous media. The region-averaged equations are developed for the case of a slightly compressible flow which is an accurate representation for a certain class of liquid-phase flows. The analysis leads to a pair of transport equations for the region averaged pressures that are coupled through a classic exchange term, in addition to being coupled by a diffusive cross effect. The domain of validity of the theory has been identified in terms of a series of length and timescale constraints.In Part II the theory is tested, in the absence of adjustable parameters, by comparison with numerical experiments for transient, slightly compressible flow in both stratified and nodular models of heterogeneous porous media. Good agreement between theory and experiment is obtained for nodular and stratified systems, and effective transport coefficients for a wide range of conditions are presented on the basis of solutions of the three closure problems that appear in the theory. Part III of this paper deals with the principle of large-scale mechanical equilibrium and the region-averaged form of Darcy's law. This form is necessary for the development and solution of the region-averaged solute transport equations that are presented in Part IV. Finally, in Part V we present results for the dispersion tensors and the exchange coefficient associated with the two-region model of solute transport with adsorption.  相似文献   

Mantle-driven geodynamo features – Effects of compositional and narrow D″ anomalies     

Hagay Amit  Gaël Choblet 《Physics of the Earth and Planetary Interiors》2012

Lower mantle heterogeneity could cause deviations from axial symmetry in geodynamo properties. Global tomography models are commonly used to infer the pattern of core–mantle boundary heat flux via a linear relation that corresponds to a purely thermal interpretation of lower mantle seismic anomalies, ignoring both non-thermal origins and non-resolved small scales. Here we study the possible impact on the geodynamo of narrow thermal anomalies in the base of the mantle, originating from either compositional heterogeneity or sharp margins of large-scale features. A heat flux boundary condition composed of a large-scale pattern and narrow ridges separating the large-scale positive and negative features is imposed on numerical dynamos. We find that hot ridges located to the west of a positive large-scale core–mantle boundary heat flux anomaly produce a time-average narrow elongated upwelling, a flow barrier at the top of the core and intensified low-latitudes magnetic flux patches. When the ridge is located to the east of a positive core–mantle boundary heat flux anomaly, the associated upwelling is weaker and the homogeneous dynamo westward drift leaks, precluding persistent intense low-latitudes magnetic flux patches. These signatures of the core–mantle boundary heat flux ridge are evident in the north–south component of the thermal wind balance. Based on the pattern of lower mantle seismic tomography (Masters et al., 2000), we hypothesize that hot narrow thermal ridges below central Asia and the Indian Ocean and below the American Pacific coast produce time-average fluid upwelling and a barrier for azimuthal flow at the top of the core. East of these ridges, below east Asia and Oceania and below the Americas, time-average intense geomagnetic flux patches are expected.  相似文献   

Calculation of primary and secondary flow and boundary shear stresses in a meandering channel     

Thorsten Stoesser  Nils Ruether  Nils Reidar Boe Olsen 《Advances in water resources》2010

Turbulent flow in a meandering channel is computed with two Computational Fluid Dynamics (CFD) codes solving the Navier–Stokes equations by employing different turbulence closure approaches. The first CFD code solves the steady Reynolds-Averaged Navier–Stokes equations (RANS) using an isotropic turbulence closure. The second code is based on the concept of Large Eddy Simulation (LES). LES resolves the large-scale turbulence structures in the flow and is known to outperform RANS models in flows in which large-scale structures dominate the statistics. The results obtained from the two codes are compared with experimental data from a physical model study. Both, LES and RANS simulation, predict the primary helical flow pattern in the meander as well as the occurrence of an outer-bank secondary cell. Computed primary as well as secondary flow velocities are in reasonably good agreement with experimental data. Evidence is given that the outer-bank secondary cell in a meander bend is the residual of the main secondary cell of the previous bend. However, the RANS code, regardless of the turbulence model employed, overpredicts the size and strength of the outer-bank secondary cell. Furthermore, only LES is able to uphold the outer-bank second secondary cell beyond the bend apex until the exit of the bend as turbulence anisotropy contributes to its persistence. The presence of multiple secondary cells has important consequences for the distribution of shear stresses along the wetted perimeter of the channel, and thereby the sediment transport in meandering channels. Consequently, even though LES is expected to compute the bed-shear stresses along the wetted perimeter of the channel with a higher degree of accuracy than the RANS model, comparisons between LES and RANS computed wall shear stresses agree well. These findings are useful for practitioners who need to rely on RANS model predictions of the flow in meandering channels at field scale.  相似文献   

Seasonal Prediction for the Indian Monsoon Region with FSU Ocean-atmosphere Coupled Model: Model Mean and 2002 Anomalous Drought     

A. K. Mitra  L. Stefanova  T. S. V. Vijaya Kumar  T. N. Krishnamurti 《Pure and Applied Geophysics》2005,162(8-9):1431-1454

Seasonal climate prediction for the Indian summer monsoon season is critical for strategic planning of the region. The mean features of the Indian summer monsoon and its variability, produced by versions of the ‘Florida State University Coupled Ocean-Atmosphere General Circulation Model’ (FSUCGCM) hindcasts, are investigated for the period 1987 to 2002. The coupled system has full global ocean and atmospheric models with coupled assimilation. Four member models were created by choosing different combinations of parameterizations of the physical processes in the atmospheric model component. Lower level wind flow patterns and rainfall associated with the summer monsoon season are examined from this fully coupled model seasonal integrations. By comparing with observations, the mean monsoon condition simulated by this coupled model for the June, July and August periods is seen to be reasonably realistic. The overall spatial low-level wind flow patterns and the precipitation distributions over the Indian continent and adjoining oceanic regions are comparable with the respective analyses. The anomalous below normal large-scale precipitation and the associated anomalous low-level wind circulation pattern for the summer monsoon season of 2002 was predicted by the model three months in advance. For the Indian summer monsoon, the ensemble mean is able to reproduce the mean features better compared to individual member models.  相似文献   

Use of upscaled elevation and surface roughness data in two-dimensional surface water models   总被引：1，自引：0，他引：1  

J.D. Hughes  J.D. DeckerC.D. Langevin 《Advances in water resources》2011,34(9):1151-1164

In this paper, we present an approach that uses a combination of cell-block- and cell-face-averaging of high-resolution cell elevation and roughness data to upscale hydraulic parameters and accurately simulate surface water flow in relatively low-resolution numerical models. The method developed allows channelized features that preferentially connect large-scale grid cells at cell interfaces to be represented in models where these features are significantly smaller than the selected grid size. The developed upscaling approach has been implemented in a two-dimensional finite difference model that solves a diffusive wave approximation of the depth-integrated shallow surface water equations using preconditioned Newton-Krylov methods. Computational results are presented to show the effectiveness of the mixed cell-block and cell-face averaging upscaling approach in maintaining model accuracy, reducing model run-times, and how decreased grid resolution affects errors. Application examples demonstrate that sub-grid roughness coefficient variations have a larger effect on simulated error than sub-grid elevation variations.  相似文献   

Comparison of direct statistical and indirect statistical-deterministic frameworks in downscaling river low-flow indices     

Deepti Joshi  Andre St-Hilaire  Taha B. M. J. Ouarda  Anik Daigle  Nathalie Thiemonge 《水文科学杂志》2013,58(11):1996-2010

ABSTRACT

This work explores the ability of two methodologies in downscaling hydrological indices characterizing the low flow regime of three salmon rivers in Eastern Canada: Moisie, Romaine and Ouelle. The selected indices describe four aspects of the low flow regime of these rivers: amplitude, frequency, variability and timing. The first methodology (direct downscaling) ascertains a direct link between large-scale atmospheric variables (the predictors) and low flow indices (the predictands). The second (indirect downscaling) involves downscaling precipitation and air temperature (local climate variables) that are introduced into a hydrological model to simulate flows. Synthetic flow time series are subsequently used to calculate the low flow indices. The statistical models used for downscaling low flow hydrological indices and local climate variables are: Sparse Bayesian Learning and Multiple Linear Regression. The results showed that direct downscaling using Sparse Bayesian Learning surpassed the other approaches with respect to goodness of fit and generalization ability.