A new form of the cokriging equations

Myers developed a matrix form of the cokriging equations, but one that entails the solution of a large system of linear equations. Large systems are troublesome because of memory requirements and a general increase in the matrix condition number. We transform Myers’s system into a set of smaller systems, whose solution gives the classical kriging results, and provides simultaneously a nested set of lower dimensional cokriging results. In the course of developing the new formulation we make an interesting link to the Cauchy-Schwarz condition for the invertibility of a system, and another to a simple situation of coregionalization. In addition, we proceed from these new equations to a linear approximation to the cokriging results in the event that the crossvariograms are small, allowing one to take advantage of a recent results of Xie and others which proceeds by diagonalizing the variogram matrix function over the lag classes.     

Observations and simulations of an unsteady island wake in the Firth of Forth, Scotland

Beamer Rock, a 50-m-wide island in the Firth of Forth, produces a distinctive von Kármán vortex street wake, the characteristics of which depend on the speed and direction of the tidal flow. ADCP (acoustic Doppler current profiler), CTD (conductivity–temperature–depth) and aerial photograph data were collected from the region during flood and ebb tidal flow under neap and spring conditions. Good agreement was found between the observations and the results from a fixed-grid depth-averaged numerical tidal model. The island wake parameter correctly predicted the unsteady nature of the wake, and the Strouhal number (defined in terms of flow past a circular cylinder) was found to give excellent predictions of the wake wavelength when scaled on the island width. Contrary to published results, the study shows that it is possible to accurately simulate an unsteady island wake using a relatively coarse fixed-grid numerical model.Responsible Editor: Jens Kappenberg     

第二松花江流域水电厂综合中长期水文预报的分析研究

基于对第二松花江流域上游小山、松山、两江水电站中长期水文预报研究成果，论述了综合中长期水文预报的研究思路和定性预报、定量预报的分析研究方法。该研究通过2002年实践检验，具有较高的预报精度。     

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.     

Mechanism of arrested charnockite formation at Nemmara, Palghat region, southern India

The formation of arrested charnockite is an excellent example of structurally controlled channellised fluid flow along specific sites accompanied by selective elemental mobility and mineralogical changes. The present paper recognises and focuses study on three types of arrested charnockite formation from Palghat region, namely, shear-controlled, foliation parallel and boudin-neck types, and address their spatio-temporal relations to regional-scale charnoenderbite. The shear-controlled and foliation parallel types post-date deformation and migmatisation. The boudin-neck type, on the contrary, is coeval with partial melting and followed the path of cooling and decreasing water activity in the gneiss. K-feldspar veining around plagioclase and quartz, symplectitic intergrowth of biotite+quartz after orthopyroxene and K-feldspar, and fluid inclusion data suggests the presence of alkalic supercritical brine and low-density CO₂-rich fluid during charnockite formation. Charnockite domains developed following the breakdown of hornblende, biotite and quartz are characterised by a more or less pronounced depletion of Fe, Ca, Mg and Ti and trace elements Y and Zr, compared to their counterpart gneiss. REE spectra indicate a subtle depletion in the HREE near the centre of the charnockite domain. Although close-pair samples of gneiss–charnockite are isochemical, on a scale of a few millimetres, bi-directional element movement, related to the formation of new mineral was noted. It is postulated that arrested charnockite formation developed in situ on local scale within the granitic domains of the hornblende-biotite gneiss, in the presence of CO₂-rich fluids and alkalic supercritical saline brine. This process post-dated the time of regional granulite (charnoenderbite) and large regional scale retrogression and migmatisation.     

Channel Bed Mobility Downstream from the Elwha Dams,Washington*

Dams are a major source of fragmentation and degradation of rivers. Although substantial research has been conducted on the environmental impacts of large structures in the United States, smaller dams have received less attention. This study evaluated the impact of two dams of moderate size, the Elwha Dams, on the downstream channel system using field data collection at river cross‐sections. The relationship of average boundary shear stress (τ_o) to critical shear stress (τ_cr) served as the basis for determining channel bed material mobility under the two‐year and ten‐year flood events. The channel had the greatest channel bed mobility at the natural cross‐section upstream from the dams, low bed mobility between the structures, and an increase in channel bed mobility in the low gradient river segment near the mouth of the river. Low bed mobility tended to be associated with a lack of channel system complexity, including reduction or loss of bars and low alluvial terraces and their associated young riparian communities. Although these run‐of‐the‐river dams do not modify streamflow greatly, the loss of sediment from the channel system has had a substantial impact on bed mobility and geomorphic and biotic complexity of the Elwha River.     

The use of heavy metal concentrations and dendrochronology in the reconstruction of sediment accumulation, Ma a Panew River Valley, southern Poland

Heavy metal concentrations were investigated in overbank sediments of the Mala Panew River, southern Poland. Samples were collected from seven vertical profiles located within channel infills of a 20th century floodplain at three sites, each up to 50 m wide. In each profile, 15–24 samples were collected and analysed for Ba, Cd, Cu, Pb, and Zn. Sequential extraction of these elements was carried out in the 0.063-mm fraction of selected samples. Additionally, the age of the oldest trees growing close to the profiles has been used to estimate the initiation of sediment accumulation there. Ba, Cu, and Pb, which occur mostly in less mobile, moderately reducible, and residual fractions, were used for sediment dating. Zn and Cd, which in 50–75% occur in the mobile exchangeable fraction, were not suitable for dating. Correlation of Ba, Cu, and Pb concentrations in vertical profiles with changes in the load of effluents discharged to the river showed abrupt changes in the thickness of the strongly polluted sediments across the floodplains. A comparison of the relative changes between heavy metal peaks in sediments of similar age in the different profiles suggests a variable rate of downward metal migration. In general, none of the heavy metals investigated seems to have been mobilised within the stratigraphic layers above the water table. In layers located at stratigraphically lower levels, the Zn and Cd peaks seem to migrate several centimetres to several decimetres down in the profile. In profiles inundated for several weeks every year, Zn and Cd, as well as the relatively less mobile Ba, Cu, and Pb, have migrated downward by several decimetres. The investigation shows that frequent fluctuations of the water table have blurred the original depositional metal patterns of metal concentrations within a period of less than 40 years.     

Turbulent flow over a dune: Green River,Colorado

Detailed echo‐sounder and acoustic Doppler velocimeter measurements are used to assess the temporal and spatial structure of turbulent flow over a mobile dune in a wide, low‐gradient, alluvial reach of the Green River. Based on the geometric position of the sensor over the bedforms, measurements were taken in the wake, in transitional flow at the bedform crest, and in the internal boundary layer. Spatial distributions of Reynolds shear stress, turbulent kinetic energy, turbulence intensity, and correlation coefficient are qualitatively consistent with those over fixed, two‐dimensional bedforms in laboratory flows. Spectral and cospectral analysis demonstrates that energy levels in the lee of the crest (i.e. wake) are two to four times greater than over the crest itself, with minima over the stoss slope (within the developing internal boundary layer). The frequency structure in the wake is sharply defined with single, dominant peaks. Peak and total spectral and cross‐spectral energies vary over the bedform in a manner consistent with wave‐like perturbations that ‘break’ or ‘roll up’ into vortices that amalgamate, grow in size, and eventually diffuse as they are advected downstream. Fluid oscillations in the lee of the dune demonstrate Strouhal similarity between laboratory and field environments, and correspondence between the peak frequencies of these oscillations and the periodicity of surface boils was observed in the field. Copyright © 2005 John Wiley & Sons, Ltd.     

The effect of bed mobility on resistance to overland flow

When sediment grains are transported as bed load in overland flow, there is a net transfer of momentum from the flow to the grains. When these grains collide with other grains, whether on the bed or in the flow, streamwise flow velocity decreases and resistance to flow increases. Resistance to flow generated in this manner is termed bed‐load transport resistance. Resistance to flow f over a plane bed may be partitioned into grain resistance f_g and bed‐load transport resistance f_bt. We use the symbols f_btf and f_btm to denote f_bt for flows over fixed beds and over mobile beds, respectively, and we compute the effect of bed mobility on flow resistance f_mob by subtracting f_btf from f_btm. The data for this study come from 54 flume experiments with fixed beds and 38 with mobile beds. On average f_mob is approximately equal to half of f_btm, which is about one‐quarter of f. Hence, f_mob is about one‐tenth of f. Predictive equations are developed for f_btf, f_btm and f_mob using dimensional analysis to identify the relevant independent variables and regression analysis to evaluate the coefficients associated with these variables. Values of f_mob are always positive which implies that mobile beds offer greater resistance to flow than do fixed beds. Evidently bed‐load grains colliding with mobile beds lose more momentum to the bed than do grains colliding with fixed beds. In other words, grain collisions with mobile beds are less elastic than those with fixed beds. Copyright © 2005 John Wiley & Sons, Ltd.     

The effective number of cycles of earthquake ground motion

The seismic response of any system that accumulates damage under cyclic loading is dependent not only on the maximum amplitude of the motion but also its duration. This is explicitly recognized in methods for estimating the liquefaction potential of soil deposits. Many researchers have proposed that the effective number of cycles of the ground motion is a more robust indicator of the destructive capacity of the shaking than the duration. However, as is the case with strong‐motion duration, there is no universally accepted approach to determining the effective number of cycles of motion, and the different methods that have been proposed can give widely varying results for a particular accelerogram. Definitions of the effective number of cycles of motion are reviewed, classified and compared. Measurements are found to differ particularly for accelerograms with broad‐banded frequency content, which contain a significant number of non‐zero crossing peaks. The key seismological parameters influencing the number of cycles of motion and associated equations for predicting this quantity for future earthquakes are identified. Correlations between cycle counts and different duration measures are explored and found to be rather poor in the absence of additional parameters. Copyright © 2004 John Wiley & Sons, Ltd.