Efficiency of a two-step upscaling method for permeability evaluation at Darcy and pore scales

This work presents a new subdivision method to upscale absolute permeability fields. This process, called two-step method, consists in (i) solving micro-scale equations on subdomains obtained from the full domain regular decomposition and (ii) solve a second upscaling with Darcy’s law on the permeability fields obtained in the first step. The micro-scale equations used depend on the case studied. The two-step upscaling process is validated on randomly generated Darcy-scale permeability fields by measuring the numerical error induced by upscaling. The method is then applied to real domains obtained from sandstone micro-tomographic images. The method specificities due to pore-space structure are discussed. The main advantage of the two-step upscaling method resides in the drastic reduction of computational costs (CPU time and memory usage) while maintaining a numerical error similar to that of other upscaling procedures. This new upscaling method may improve permeability predictions by the use of finer meshes or larger sample volumes.     

Scheduling the EChO survey with known exoplanets

The Exoplanet Characterization Observatory (EChO) is a concept of a dedicated space telescope optimized for low-resolution transit and occultation spectroscopy to study the exoplanet diversity through the composition of their atmospheres. The scope of this paper is to answer the following question: Can we schedule a nominal EChO mission, with targets known today (in mid 2013), given the science requirements, realistic performances and operational constraints? We examine this issue from the point of view of duration of the mission and the scheduling restrictions with a sample of exoplanet systems known nowadays. We choose different scheduling algorithms taking into account the science and operational constraints and we verified that it is fairly straightforward to schedule a mission scenario over the lifetime of EChO compliant with the science requirements. We identified agility as a critical constraint that reduces significantly the efficiency of the survey. We conclude that even with known targets today the EChO science objectives can be reached in the 4.5 years duration of the mission. We also show that it is possible to use gaps between exoplanet observations, to fit the required calibration observations, data downlinks and station keeping operations or even to observe more exoplanet targets to be discovered in the coming years.     

Estimating the error in vertically interpolated forecast‐observation differences: U‐ and v‐wind components

Abstract

Data assimilation in numerical weather forecasting corrects current forecast values by subtracting a portion of interpolated forecast‐minus‐observation differences at the points of a three‐dimensional grid. Deviations used in updating a forecast data field are forecast errors obtained or derived from observations available at update time. When observations are missing at mandatory levels, construction of full vertical soundings by interpolation introduces extraneous errors. The present paper is concerned with determination of the error in vertical extrapolations of surface winds, and of aircraft and satellite cloud‐tracked winds. In addition it examines the effect on accuracy of using location‐specific statistics compared to averaged statistics as the basis for the interpolation weighting scheme and compares errors of one‐ and two‐variable interpolations.

Interpolation accuracy tests demonstrate the influence of the interpolation scheme on the quality of interpolated information used in forecast updating. The results show that the level of accuracy exceeds the benchmark provided by monthly mean forecast error values only with bivariate interpolation of wind components from off‐level data sources.     

Canadian middle atmosphere model: Preliminary results from the chemical transport module

Abstract

An important objective of middle atmosphere global climate modelling is the development of the capability of predicting the response of the middle atmosphere to natural or anthropogenic perturbations. To achieve this, a comprehensive chemistry package interactively coupled with radiative and dynamical modules is required. This paper presents preliminary results obtained with a photochemistry module which has been incorporated in the Canadian Middle Atmosphere Model (CMAM). The module contains 42 species including necessary oxygen, hydrogen, nitrogen, chlorine, bromine and methane oxidation cycle species. Photochemical balance equations are solved on‐line throughout the middle atmosphere at every dynamical time step. A full diurnal cycle is simulated with photolysis rates provided by a look‐up table. The chemistry solver is a mass conserving, fully implicit, backward difference scheme which currently uses less than 10% of the GCM run time. We present the results obtained from short integrations and compare them with UARS measurements. The model ozone distribution appears in quantitative agreement with observations showing peak values near 10 ppmv and confined to the 35‐km region. The abundance of nitrogen, chlorine, bromine oxides and their respective contributions to the overall ozone budget is realistic. The study illustrates the capability of the model to simulate middle atmosphere photochemistry for the disparate conditions occurring throughout the region.     

Climate and climate change in western canada as simulated by the Canadian regional climate model

Abstract

A þrst climate simulation performed with the novel Canadian Regional Climate Model (CRCM) is presented. The CRCM is based on fully elastic non‐hydrostatic þeld equations, which are solved with an efþcient semi‐implicit semi‐Lagrangian (SISL) marching algorithm, and on the parametrization package of subgrid‐scale physical effects of the second‐generation Canadian Global Climate Model (GCMII). Two 5‐year integrations of the CRCM nested with GCMII simulated data as lateral boundary conditions are made for conditions corresponding to current and doubled CO₂ scenarios. For these simulations the CRCM used a grid size of 45 km on a polar‐stereographic projection, 20 scaled‐height levels and a time step of 15 min; the nesting GCMII has a spectral truncation of T32, 10 hybrid‐pressure levels and a time step of 20 min. These simulations serve to document: (1) the suitability of the SISL numerical scheme for regional climate modelling, (2) the use of GCMII physics at much higher resolution than in the nesting model, (3) the ability of the CRCM to add realistic regional‐scale climate information to global model simulations, and (4) the climate of the CRCM compared to that of GCMII under two greenhouse gases (GHG) scenarios.