Assessment of hydrodynamic simulation results for eco‐hydraulic and eco‐hydrological applications: a spatial semivariance approach

Output from a three‐dimensional numerical flow model (SSIIM) is used in conjunction with high‐resolution topographic and velocity data to assess such models for eco‐hydraulic applications in river channel design and habitat appraisal. A new methodology for the comparison between field measurement and model output is detailed. This involves a comparison between conventional goodness‐of‐fit approaches applied to a spatially structured (riffle and pool) sample of model and field data, and a ‘relaxation’ method based upon the spatial semivariance of model/field departures. Conventional assessment indicates that the model predicts point‐by‐point velocity characteristics on a 0·45 m mesh to within ±0·1 m s⁻¹ over 80% of the channel area at low flow, and 50% of the area at high in‐bank flow. When a relative criterion of model fit is used, however, the model appears to perform less well: 60–70% of channel area has predicted velocities that depart from observed velocities by more than 10%. Regression analysis of observed and predicted velocities gives more cause for optimism, but all of these conventional indicators of goodness of fit neglect important spatial characteristics of model performance. Spatial semivariance is a means of supplementing model appraisal in this respect. In particular, using the relaxation approach, results are greatly improved: at a high in‐bank flow, the model results match field measurements to within 0·1 m s⁻¹ for more than 95% of the total channel area, provided that model and field comparisons are allowed within a radius of approximately 1 m from the original point of measurement. It is suggested that this revised form of model assessment is of particular relevance to eco‐hydraulic applications, where some degree of spatial and temporal dynamism (or uncertainty) is a characteristic. The approach may also be generalized to other environmental science modelling applications where the spatial attributes of model fits are of interest. Copyright © 2005 John Wiley & Sons, Ltd.     

Output and abatement effects of allocation readjustment in permit trade

In permit trading systems, free initial allocation is common practice. A recent example is the European Union Greenhouse Gas Emission Trading Scheme (EU-ETS). We investigate effects of different free allocation schemes on incentives and identify significant perverse effects on abatement and output employing a simple multi-period model. Firms have incentives for strategic action if allocation in one period depends on their actions in previous ones and thus can be influenced by them. These findings play a major role where trading schemes become increasingly popular as environmental or resource use policy instruments. This is of particular relevance in the EU-ETS where the current period is a trial-period before the first commitment period of the Kyoto protocol. Finally, this paper fills a gap in the literature by establishing a consistent terminology for initial allocation.     

Radiation-induced paramagnetic defects in natural kaolinites: Alpha dosimetry with ion beam irradiation

Three types of paramagnetic radiation-induced defects (RID), namely A, A′ (Si-O^?-centers) and B (Al-O^?-Al center), had been identified in natural kaolinites by means of electron paramagnetic resonance spectroscopy. The A-center, stable at the scale of geological periods, was thought to be of particular relevance to quantify past transits of radionuclides in the geosphere. Alpha radiation being likely the main source of RID in kaolinite, the objective of this paper is to define the role of α-particles on the creation of RID and to test the use of A-centers for an α-dosimetry. Three kaolinites with different crystalline order and containing other clays as impurities were irradiated with He⁺ ion beams. The radiation dose range (0–750 MGy) was consistent with natural radioactivity in environments from the Earth's surface. Irradiation drastically enhanced the original signals due to RID. An important increase of concentration of the unstable B-center, partly due to unrealistic dose rate provided by accelerator ion beam, was observed from the lower doses. The most stable defects remained of the Atype all along experimental irradiations. The contribution of ancillary phyllosilicates to EPR spectra was negligible. The concentration of the three types of RID was related to radiation doses up to 30 MGy. Dosimetry growth curves for the A-centers exhibited variable efficiencies and saturation levels that were related to the structural order and the chemical purity of the kaolinites: the more ordered and purer the kaolinite, the higher the efficiency and the lower the saturation plateau. Moreover, these results are of geochemical significance: dosimetry based on A-centers could be directly used to quantify past migrations of radioelements in the geosphere, by determination of the naturally-cumulated doses (paleodose) of kaolinites.     

The stress state near Spanish Peaks,colorado determined from a dike pattern

The radial pattern of syenite and syenodiorite dikes of the Spanish Peaks region is analysed using theories of elasticity and dike emplacement. The three basic components of Odé's model for the dike pattern (a pressurized, circular hole; a rigid, planar boundary; and uniform regional stresses) are adopted, but modified to free the regional stresses from the constraint of being orthogonal to the rigid boundary. Dike areal density, the White Peaks intrusion, the strike of the upturned Mesozoic strata, and the contact between these strata and the intensely folded and faulted Paleozoic rocks are used to brient the rigid boundary along a north-south line. The line of dike terminations locates the rigid boundary about 8 km west of West Peak. The location of a circular plug, Goemmer Butte, is chosen as a point of isotropic stress. A map correlating the location of isotropic stress points with regional stress parameters is derived from the theory and used to determine a regional stress orientation (N82E) and a normalized stress magnitude. The stress trajectory map constructed using these parameters mimics the dike pattern exceptionally well. The model indicates that the regional principal stress difference was less than 0.05 times the driving pressure in the West Peak intrusion. The regional stress difference probably did not exced 5 MN/m².     

Stratospheric methane measurements and predictions

A new determination of stratospheric methane from 22 km to 35 km altitude with implications on the abundance of this constituent at greater heights is presented. Previous measurements, some of which showed large discrepancies with currently admitted values, have been reinterpreted and brought into agreement. The results are in contradiction according to present theories with thein situ determined CH₄ abundances at the upper edge of the stratosphere. Implications on the vertical transport coefficient used in one dimensional models are discussed.     

Calculating the effective permeability of sandstone with multiscale lattice Boltzmann/finite element simulations

The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations.     

Morphological properties and origin of the photospheric facular granules

From time series of high resolution photographs, morphological properties of the photospheric facular granules were derived. The facular granules are cells of the common granular pattern, brighter than the normal granules when seen between cos = 0.6 and the limb. Their apparent diameter, which decreases towards the limb, is smaller than that of the normal granules: 0.65 and 1.25 respectively at cos = 0.55; their lifetime is 25 min but their bright stage lifetime is only 15 min; they are visible closer to the limb than the normal granules: 1.2 compared to 2–5; the brightening of the facular granules occurs at a faster rate than their fading. From the great similitude of both morphological properties and temperature models of facular and normal granules, it appears possible that the photospheric facular granules are convective cells modified by the presence of a magnetic field of some hundreds Gauss.     

Significance of a combined approach for replacement stones in the heritage buildings’ conservation frame

Stone substitution is a conventional operation during heritage buildings’ restoration, but becomes problematic for architects and restorers when the quarry is mined out. The compatibility of the substitution stones with the original ones has been for long mainly based on the aesthetical aspect, this resulting too often in a patchwork of original and substitution stones with different patina after several years because of differences of properties. In this study, the objective is to show how substitution stones can be selected by combining aesthetic criteria and stones properties that are relevant for analyzing their compatibility. A couple of French limestones with their potential substitution stones were selected for the study. Our results showed that potential substitution stones selected on their aesthetic criteria require to be rejected because of their differences of physical properties. On the other hand, our results showed also the possibility to select substitution stones with satisfactory aesthetic aspect and properties that enable to expect a satisfactory compatibility with the original stone.     

Size fractionation and optical properties of colloids in an organic-rich estuary (Thurso, UK)

The optical characteristics of a black water river estuary from the north coast of Scotland were examined in the filtered (0.4 µm), ultrafiltered (5 kDa) and colloid-enriched fractions of estuarine samples. The samples were collected over the full salinity range during a period when the pH was relatively constant (8.2–8.5) throughout the estuary, allowing the influence of salinity on estuarine colloidal processes to be distinguished. The properties examined in the bulk, the low molecular weight (LMW) and the colloidal fraction (HMW) were UV–visible absorption, 3-D fluorescence excitation–emission matrix (EEM) spectrum, inorganic and organic carbon, mean size (by dynamic light scattering), and size distribution by flow field-flow fractionation analysis (FlFFF). The combined results of these analyses support the view that river-borne, humic-rich colloids underwent two types of transformation upon mixing with the seawater end member. The first one resulted in an apparent increase in the abundance of LMW constituents and may be explained by coiling of the individual humic macromolecules. The second one resulted in an increase in the mean size measured in both the lower and higher colloidal size ranges, and may be explained by aggregation of colloids to form entities that were still mostly colloidal i.e., smaller than 0.4 µm. The LMW contribution to the bulk optical properties increased with increasing salinity. Very similar findings were obtained from simulated mixing experiments using a Nordic Reference NOM extract as a source of freshwater colloids. This indicates that changes in the molecular architecture and molar mass of river-borne colloids—not changes in their chemical nature—were responsible for the observed variations in the spectral characteristics of CDOM in this estuary.