Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

Empirical correlations of soil parameters based on Cone Penetration Tests (CPT) for Greek soils

This study presents the results of Cone Penetration Tests (CPT) and estimated correlations with the results of other in-situ tests (standard penetration test, cross-hole, seismocone) and laboratory tests (classification and undrained shear strength tests). These tests cover nearly all types of penetrable soils in Greece and, together with the formulation of empirical correlations, complement and improve on previous studies due to the significant number of CPT measurements.     

Thermodynamic properties of the Sb(III) hydroxide complex Sb(OH)3(aq) at hydrothermal conditions

The standard thermodynamic properties and Helgeson-Kirkham-Flowers (HKF) parameters for Sb(OH)_3(aq) have been estimated. For this purpose, the available solubility data for senarmontite, valentinite, stibnite, and native Sb in a wide range of temperatures (15 to 450°C) and pressures (1 to 1000 bar), and thermodynamic properties of Sb oxides (senarmontite and valentinite) have been critically analyzed. Published data were complimented by results from new experiments performed by solubility and solid-state galvanic cell methods. Both experimental data and thermodynamic calculations show that the hydroxide complex Sb(OH)_3(aq) is primarily responsible for hydrothermal transport of antimony, especially at temperatures above 250°C.     

Isolation and water-level fluctuations of Lake Kachishayoot, Northern Québec, Canada

This study combines different methods, including grain size, macrofossil, and pollen analyses, to reconstruct paleogeographical and paleoclimatological conditions for periods before, during, and after the isolation of a small lake (Lake Kachishayoot) in northern Québec. After the retreat of the Laurentide Ice Sheet around 8000 ¹⁴C yr B.P., the area was submerged by the Tyrrell Sea. The transition from marine to lacustrine environment occurred about 5400 yr B.P. Two major periods of water-level fluctuations were inferred from organic and mineral sediments: a high water level that occurred after 3200 yr B.P. and a low water level that started before 2200 yr B.P. Our chronological data for the first period are consistent with those from nearby Lac des Pluviers and from other lakes in east central Canada and in the northeastern United States. During the low-water-level period, however, there is no evidence for minor fluctuations, whereas other lakes in northern Québec and east-central Canada underwent several brief lowerings. Long-term changes in atmospheric circulation caused by changing global boundary conditions likely explained long-term water-level fluctuations of Lake Kachishayoot.     

Estuary-specific variation in the air-water gas exchange coefficient for oxygen

Oxygen air-water gas exchange was measured using floating chambers in two shallow tidal estuaries of differing bathymetry and local terrain, near Waquoit Bay, Massachusetts (United States). The specific chamber design permitted measurements of gas flux in 15 min, allowing analysis of the relationship with wind speed and tidal stage. Exchange coefficients ranged from 0.5 to 2.5 g O₂·m^?2 h^?1 atm^?1 (equivalent to piston velocities of 1.5 to 7 cm h^?1) for wind speeds of 0.3 to 9 m s^?1 at 10 m elevation. While the relationships for each estuary appear linear (significant linear regressions with wind speed were shown for each estuary, and the slopes were different at the 99.5% confidence level), the range of speeds differed at the two sites and an exponential function of wind speed was consistent with the combined data from both estuaries. A power function of wind speed was not an acceptable model. The exchange coefficients for our estuaries are from 57% to as low as 9% of that predicted by previously published generic equations. Because the atmospheric correction can be significant in shallow, metabolically active coastal waters, we suggest that empirically determined relationships for gas exchange versus wind for a specific estuary are preferable to the predictions of the general equations. While the floating chamber method should be used cautiously, at low winds speeds (below 8 m s^?1) and in slowly flowing waters, it provides a convenient approach for quantifying these site-specific differences. The differences, especially those between shallow sheltered systems and the open waters best fit by some published relationships, are ecologically important and do not appear yet to be measurable by other methods.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal     

Regional effects of large-scale extreme wind events over orographically structured terrain

Summary ?Above orographically structured terrain considerable differences of the regional wind field may be identified during large-scale extreme wind events. So far, these regional differences could not be resolved by climate models. To determine the relationships between large-scale atmospheric conditions, the influence of orography, and the regional wind field, data measured in the upper Rhine valley within the framework of the REKLIP Regional Climate Project were analyzed and calculations were made using the KAMM mesoscale model. In the area of the upper Rhine valley, ratios of the wind velocity in the Rhine valley at 10 m above ground level, ν_val, and the large-scale flow velocity, ν_lar, are between ν_val/ν_lar ≈ 0.1 and ν_val/ν_lar ≈ 1. The ν_val/ν_lar ratio exhibits a strong dependence on thermal stratification, δ, and decreases from ν_val/ν_lar ≈ 1 at δ = 0 K m⁻¹ to ν_val/ν_lar ≈ 0.2 at δ = 0.0075 K m⁻¹. In areas, where the lateral mountainous border of the Rhine valley is interrupted, the ν_val/ν_lar ratio increases again with increasing stability or decreasing Froude number. This is obviously due to flow around the Black Forest under stable stratification. It is demonstrated by model calculations that a complex wind field develops in the Rhine valley at small Froude numbers (Fr < 1) irrespective of the direction of large-scale flow. The ν_val/ν_lar ratio is characterized by small values in the direct lee side (ν_val/ν_lar ≈ 0.2) and high values on the windward side of the lateral mountainous border of the Rhine valley (ν_val/ν_lar ≈ 0.8). Received October 22, 2001; revised June 18, 2002; accepted June 23, 2002     

A coupled model study on ENSO, MJO and Indian summer monsoon rainfall relationships

Summary In this paper, we have tried to understand the ENSO, MJO and Indian summer monsoon rainfall relationships from observation as well as from coupled model results. It was the general feeling that El-Niño years are the deficient in Indian monsoon rainfall and converse being the case for the La-Niña years. Recent papers by several authors noted the failure of this relationship. We find that the model output does confirm a breakdown of this relationship. In this study we have seen that a statistically defined modified Indian summer monsoon rainfall (MISMR) index, a linearly regressed ISMR index and dynamical Webster index (WBSI), shows an inverse relationship with ENSO index during the entire period of integration (1987 to 1999). It is also seen from this study that the amplification of the MJO signals were large and the ENSO signals were less pronounced during the years of above normal ISMR. The MJO signal amplitudes were small and ENSO signals were strong during the years of deficient ISMR. It has been noted that here is a time lag between the MJO and ENSO signal in terms of their modulation aspect. If time lag is added with the ENSO signal then both signals maintain the amplitude modulation theory. A hypothesis is being proposed here to define a relationship between MJO and ENSO signals for the entire period between 1987 and 1999.Received September 18, 2002; revised November 22, 2002; accepted December 20, 2002 Published online: May 8, 2003     

Description and evaluation of the bergen climate model: ARPEGE coupled with MICOM

A new coupled atmosphere–ocean–sea ice model has been developed, named the Bergen Climate Model (BCM). It consists of the atmospheric model ARPEGE/IFS, together with a global version of the ocean model MICOM including a dynamic–thermodynamic sea ice model. The coupling between the two models uses the OASIS software package. The new model concept is described, and results from a 300-year control integration is evaluated against observational data. In BCM, both the atmosphere and the ocean components use grids which can be irregular and have non-matching coastlines. Much effort has been put into the development of optimal interpolation schemes between the models, in particular the non-trivial problem of flux conservation in the coastal areas. A flux adjustment technique has been applied to the heat and fresh-water fluxes. There is, however, a weak drift in global mean sea-surface temperature (SST) and sea-surface salinity (SSS) of respectively 0.1 °C and 0.02 psu per century. The model gives a realistic simulation of the radiation balance at the top-of-the-atmosphere, and the net surface fluxes of longwave, shortwave, and turbulent heat fluxes are within observed values. Both global and total zonal means of cloud cover and precipitation are fairly close to observations, and errors are mainly related to the strength and positioning of the Hadley cell. The mean sea-level pressure (SLP) is well simulated, and both the mean state and the interannual standard deviation show realistic features. The SST field is several degrees too cold in the equatorial upwelling area in the Pacific, and about 1 °C too warm along the eastern margins of the oceans, and in the polar regions. The deviation from Levitus salinity is typically 0.1 psu – 0.4 psu, with a tendency for positive anomalies in the Northern Hemisphere, and negative in the Southern Hemisphere. The sea-ice distribution is realistic, but with too thin ice in the Arctic Ocean and too small ice coverage in the Southern Ocean. These model deficiencies have a strong influence on the surface air temperatures in these regions. Horizontal oceanic mass transports are in the lower range of those observed. The strength of the meridional overturning in the Atlantic is 18 Sv. An analysis of the large-scale variability in the model climate reveals realistic El Niño – Southern Oscillation (ENSO) and North Atlantic–Arctic Oscillation (NAO/AO) characteristics in the SLP and surface temperatures, including spatial patterns, frequencies, and strength. While the NAO/AO spectrum is white in SLP and red in temperature, the ENSO spectrum shows an energy maximum near 3 years.     

Comparison of the Langrangian Footprint

We compare flux and concentration footprint estimates of athree-dimensional Lagrangian stochastic dispersion modelapplying backward trajectories with the results of ananalytical footprint model by Kormann and Meixner.The comparison is performed for varying stability regimesof the surface layer as well as for different measurementheights. In general, excellent correspondence is found.