Spatial variation of overbank aggradation rate and its influence on avulsion frequency

Abstract River avulsions are commonly considered to be driven by the aggradation and growth of alluvial ridges, and the associated increase in cross‐valley slope relative to either the down‐channel slope or the down‐valley slope (the latter is termed the slope ratio in the present paper). Therefore, spatial patterns of overbank aggradation rate over stratigraphically relevant time scales are critical in avulsion‐dominated models of alluvial architecture. Detailed evidence on centennial‐ to millennial‐scale floodplain deposition has, to date, been largely unavailable. New data on such long‐term overbank aggradation rates from the Rhine–Meuse and Mississippi deltas demonstrate that the rate of decrease of overbank deposition away from the channel belt is much larger than has been supposed hitherto, and can be similar to observations for single overbank floods. This leads to more rapid growth of alluvial ridges and more rapid increase in slope ratios, potentially resulting in increased avulsion frequencies. A revised input parameter for overbank aggradation rate was used in a three‐dimensional model of alluvial architecture to study its effect on avulsion frequency. Realistic patterns of avulsion and interavulsion periods (≈1000 years) were simulated with input data from the Holocene Rhine River, with avulsions occurring when the slope ratio is in the range 3–5. However, caution should be practised with respect to uncritical use of these numbers in different settings. Evidence from the two study areas suggests that the avulsion threshold cannot be represented by one single value, irrespective of whether critical slope ratios are used, as in the present study, or superelevation as has been proposed by other investigators.     

A model of oscillatory zoning in solid solutions grown from aqueous solutions: Applications to the (Ba,Sr)SO4 system

Barite-celestite crystals can be synthesized from aqueous solutions during counter-diffusion in a gel column connecting two reservoirs. It is known that such crystals may exhibit oscillatory zoning, whereby the barium composition in the crystal fluctuates more or less regularly from the core of the crystal to its rim. We present here a simple model of oscillatory zoning in such binary solid solutions A₁A₂ grown from aqueous solutions. The model combines diffusive transport of the relevant ions with an autocatalytic growth process. The latter is formulated as a continuous growth in which the probability of finding a kink site on the growing surface depends on the chemical composition of that surface. Thus, an A₁-rich surface favors the growth of A₁ over A₂, as long as A₁ is present in the vicinity of the surface. Precipitation results in a local depletion of A₁ in the aqueous solution, and the system may switch to a A₂ growth mode, until diffusion replenishes the amount of A₁, and so on. We use a dynamical equation for the molar fraction of component A₁ in the crystal, which results from mass conservation across the rough crystal-solution interface. Linear stability analysis and direct numerical solutions show that the system exhibits oscillatory behavior. Using the barite-celestite system as a framework, the scaling is consistent with the experimental observations. We discuss the variety of zoning patterns and textures numerically obtained as the concentrations of reactants in the reservoirs vary. This model might help in understanding the formation of oscillatory zoning in hydrothermal environments.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Erratum: Comparison of Seismic Dispersion and Attenuation Models: Stud. Geophys. Geod., 46(2002), 293-320

Studia Geophysica et Geodaetica -     

Estimates of lichen growth–rate in northern Sweden

Two lichenometric techniques were compared in a study of lichen growth–rate in northern Sweden. The first technique, based on the maximum lichen diameter on glacier moraines, was identical to the technique used in the 1970s, whereas the other utilized the lichen diameter measured on 100 randomly selected boulders. The results indicate that it does not matter which technique is chosen, as long as the technique is used consistently on both the calibration surfaces and the surfaces to be dated. The use of data from both the 1970s and the 2000s increased the number of calibration surfaces available. The new calibration curve indicates that the age of Little Ice Age moraines was underestimated by up to about 30 years in the study conducted in the 1970s.     

Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.     

Reactive Nitrogen Compounds at Spitsbergen in the Norwegian Arctic

Simultaneousindependent measurements of NO_y and NO_x(NO_x= NO + NO₂) by high-sensitivitychemiluminescence systems and of PAN (peroxyacetylnitrate) and PPN (peroxypropionyl nitrate) by GC-ECDwere made at Spitsbergen in the Norwegian Arcticduring the first half year of 1994. The average mixingratio of the sum of PAN and PPN (denoted PANs)increased from around 150 pptv in early winter to amaximum of around 500 pptv in late March, whereasepisodic peak values reached 800 pptv. This occurredsimultaneously with a maximum in ozone which increasedto 45–50 ppbv in March–April. The average NO_xmixing ratio was 27 pptv and did not show any cyclethrough the period. The NO_y mixing ratio showeda maximum in late March, while the difference betweenNO_y and PAN decreased during spring. This is anindication of the dominance of PAN in the NO_ybudget in the Arctic, but possible changes in theefficiency of the NO_y converter could alsocontribute to this. Although most PAN in theArctic is believed to be due to long range transport,the observations indicate local loss and formationrates of up to 1–2 pptv h^-1 in April–May.Measurements of carbonyl compounds suggest thatacetaldehyde was the dominant, local precursor ofPAN.Now at 1.     

On the variability of the flow along the Meso-American Barrier Reef system: a numerical model study of the influence of the Caribbean current and eddies

A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize, Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly (SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies, and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef, the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation. The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region.