Influence of calibration methodology on ground water flow predictions

We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).     

The Pioneer anomaly and new physics

The Pioneer anomaly is one of the most important problems in modern physics. The observed blueshift of the Doppler signals coming back from the space probes Pioneer 10 and 11 is interpreted as being due to an anomalous acceleration a_p = (8.74 ± 1.33) × 10^–8 cm s^–2 towards the Sun. In this paper the blueshift is explained by the frequency shifts of the receivers. These frequency shifts result from an increase in elementary particle masses in time, the rate of increase being tied up with the present‐day Hubble parameter H₀. The result is that the seeming acceleration a_p is the product of H₀ and the velocity of light. Taking new physics into consideration, this paper presents a new explanation of the Pioneer anomaly based on the assumption that the Universe is eternal and infinite without expansion or contraction (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Late Mesozoic magmatism and Cenozoic tectonic deformations of the Barents Sea continental margin: Effect on hydrocarbon potential distribution

The paper is focused on the two tectonic-geodynamic factors that made the most appreciable contribution to the transformation of the lithospheric and hydrocarbon potential distribution at the Barents Sea continental margin: Jurassic-Cretaceous basaltic magmatism and the Cenozoic tectonic deformations. The manifestations of Jurassic-Cretaceous basaltic magmatism in the sedimentary cover of the Barents Sea continental margin have been recorded using geological and geophysical techniques. Anomalous seismic units related to basaltic sills hosted in terrigenous sequences are traced in plan view as a tongue from Franz Josef Land Archipelago far to the south along the East Barents Trough System close to its depocentral zone with the transformed thinned Earth’s crust. The Barents Sea igneous province has been contoured. The results of seismic stratigraphy analysis and timing of basaltic rock occurrences indicate with a high probability that the local structures of the hydrocarbon (HC) fields and the Stockman-Lunin Saddle proper were formed and grew almost synchronously with intrusive magmatic activity. The second, no less significant multitectonic stress factor is largely related to the Cenozoic stage of evolution, when the development of oceanic basins was inseparably linked with the Barents Sea margin. The petrophysical properties of rocks from the insular and continental peripheries of the Barents Sea shelf are substantially distinct as evidence for intensification of tectonic processes in the northwestern margin segment. These distinctions are directly reflected in HC potential distribution.     

Radium uptake during barite recrystallization at 23 ± 2 °C as a function of solution composition: An experimental Ba and Ra tracer study

High-purity synthetic barite powder was added to pure water or aqueous solutions of soluble salts (BaCl₂, Na₂SO₄, NaCl and NaHCO₃) at 23 ± 2 °C and atmospheric pressure. After a short pre-equilibration time (4 h) the suspensions were spiked either with ¹³³Ba or ²²⁶Ra and reacted under constant agitation during 120-406 days. The pH values ranged from 4 to 8 and solid to liquid (S/L) ratios varied from 0.01 to 5 g/l. The uptake of the radiotracers by barite was monitored through repeated sampling of the aqueous solutions and radiometric analysis. For both ¹³³Ba and ²²⁶Ra, our data consistently showed a continuous, slow decrease of radioactivity in the aqueous phase.Mass balance calculations indicated that the removal of ¹³³Ba activity from aqueous solution cannot be explained by surface adsorption only, as it largely exceeded the 100% monolayer coverage limit. This result was a strong argument in favor of recrystallization (driven by a dissolution-precipitation mechanism) as the main uptake mechanism. Because complete isotopic equilibration between aqueous solution and barite was approached or even reached in some experiments, we concluded that during the reaction all or substantial fractions of the initial solid had been replaced by newly formed barite.The ¹³³Ba data could be successfully fitted assuming constant recrystallization rates and homogeneous distribution of the tracer into the newly formed barite. An alternative model based on partial equilibrium of ¹³³Ba with the mineral surface (without internal isotopic equilibration of the solid) could not reproduce the measured activity data, unless multistage recrystallization kinetics was assumed. Calculated recrystallization rates in the salt solutions ranged from 2.8 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ mol m⁻² s⁻¹ (2.4-16 μmol m⁻² d⁻¹), with no specific trend related to solution composition. For the suspensions prepared in pure water, significantly higher rates (∼5.7 × 10⁻¹⁰ mol m⁻² s⁻¹ or ∼49 μmol m⁻² d⁻¹) were determined.Radium uptake by barite was determined by monitoring the decrease of ²²⁶Ra activity in the aqueous solution with alpha spectrometry, after filtration of the suspensions and sintering. The evaluation of the Ra uptake experiments, in conjunction with the recrystallization data, consistently indicated formation of non-ideal solid solutions, with moderately high Margules parameters (W_AB = 3720-6200 J/mol, a₀ = 1.5-2.5). These parameters are significantly larger than an estimated value from the literature (W_AB = 1240 J/mol, a₀ = 0.5).In conclusion, our results confirm that radium forms solid solutions with barite at fast kinetic rates and in complete thermodynamic equilibrium with the aqueous solutions. Moreover, this study provides quantitative thermodynamic data that can be used for the calculation of radium concentration limits in environmentally relevant systems, such as radioactive waste repositories and uranium mill tailings.     

Some cosmochemical consequences of a turbulent protoplanetary cloud

Physical processes affecting dust and vapor in a turbulent disk, which is hot in a small central region (T > 1500°K) and cold at its outer edge (T ～ 50°K) are discussed. Such disks are obtained from hydrodynamic collapse calculations of rotating gas clouds. While these calculations give the result that most of the disk, including the planetary region, remains cold (?400°K) throughout the evolution of the system, it can be shown, nevertheless, that the information of the hot interior is preserved in solids. Therefore it is suggested that the basic cosmochemical evidence from meteorite analysis, in favor of a hot protoplanetary cloud, can be understood. Detailed calculations involving estimates of sizes, both for inclusions and chondrules, as well as rim structures are given.     

Observation of chromospheric evaporation during the solar maximum mission

A sample of flares detected in 1980 with the Bent Crystal Spectrometer and the Hard X-Ray Burst Spectrometer on the Solar Maximum Mission satellite has been analysed to study the upward motions of part of the soft X-ray emitting plasma. These motions are inferred from the presence of secondary blue-shifted lines in the Ca XIX and Fe XXV spectral regions during the impulsive phase of disk flares. Limb flares do not show such blue-shifted lines indicating that the direction of the plasma motion is mainly radial and outward. The temporal association of these upward motions with the rise of the thermal phase and with the impulsive hard X-ray burst, as well as considerations of the plasma energetics, favour the interpretation of this phenomenon in terms of chromospheric evaporation. The two measureable parameters of the evaporating plasma, emission measure and velocity, depend on parameters related to the energy deposition and to the thermal phase. The evaporation velocity is found to be correlated with the spectral index of the hard X-ray flux and with the rise time of the thermal emission measure of the coronal plasma. The emission measure of the rising plasma is found to be correlated with the total energy deposited by the fast electrons in the chromosphere by collisions during the impulsive phase and with the maximum emission measure of the coronal plasma.     

Variations of cosmic-ray energy in interplanetary space

A consistent theory of energy exchange between high-energy charged cosmic-ray particles and the random inhomogeneities of a magnetic field frozen in the moving solar wind plasma is developed. It is shown that the mode of the particle energy variations at a given law of plasma velocity variation in space is determined by the specific form of the particle distribution function. The equation for the density of cosmic-ray energy is obtained. Consideration is given to the generation of a charged particle energy spectrum in the course of multiple scatterings by the random inhomogeneities of the magnetic field.     

Sediment budget of the Napo River,Amazon basin,Ecuador and Peru

The upstream‐downstream sediment budget along the Napo River (100 520 km², 6300 m³ s^?1) was studied in the Andean foothills of Ecuador, at the west of the Amazon basin. A comparative study was made during four hydrological cycles (2001–2005) for three hydrological stations located upstream, and during one hydrological cycle (2004–2005) for the fourth one located near the mouth of the Napo River (region of Iquitos in Peru). This analysis showed an unusual increase in the concentration of suspended sediment recorded for the western part of the Amazon plain. Like the runoff (81 l s^?1 km²), which is a world's maximum, the erosion rate (1160 t km^?2 year^?1, i.e. 47% of total suspended solid (TSS) export at the exit of Ecuador), one of the highest for a floodplain basin is the result of a stepper slope than in the rest of the Andean foothills, where typically sedimentation phenomena are predominant, and can be explained in part by a greater tectonic activity. Similar phenomenes were evidenced in small mountainous rivers in New Guinea (Milliman and Syvitski, 1992; Milliman, 1995). On the headwaters of the Napo River drainage basin, the tectonic uplift causes the Pastaza Megafan's existence. This progressively diverts the course of Napo River towards north and also provokes the remobilization of fine fluvial deposits. Moreover, this geodynamic trend is completed by the impact of volcanic eruption, earthquakes and landslides. The combination of these phenomena, so common in the region, has provided a large sediment transfer, not only at present but also in the past, as can be confirmed by the presence of incised terraces, mainly formed by volcanic materials. Then, these results were compared with a similar study carried out further south in the Madeira basin at the Bolivian foothills. These studies show the spatio‐temporal variability of the relation between sediment transfer and geodynamic processes at the Andean Piedmont. Copyright © 2009 John Wiley & Sons, Ltd.     

The linkage between velocity patterns and sediment entrainment in a forced‐pool and riffle unit

A field‐based project was initiated in order to characterize velocities and sediment entrainment in a forced‐pool and riffle sequence. Three‐dimensional velocities and turbulence intensities were measured with an acoustic Doppler velocimeter at 222 different points at three similar flows that averaged approximately 4·35 m³ s⁻¹ within a large pool–riffle unit on North Saint Vrain Creek, Colorado. Sediment‐sorting patterns were observed with the introduction of 500 tracer particles painted according to initial seeding location. Tracer particles moved sporadically during a 113 day period in response to the annual snowmelt peak flow, which reached a maximum level of 14·8 m³ s⁻¹. Velocity data indicate high instantaneous velocities and turbulence levels in the centre of pools. Patterns of sediment deposition support the notion that stream competence is higher in the pool than the downstream riffle. Flow convergence around a large channel constriction appears to play a major role in multiple processes that include helical flow development and sediment routing, and backwater development with low velocities and turbulence levels above the constriction that may locally limit sediment supply. Jet flow, flow separation, vortex scour and turbulence generation enhance scour in the centre of pools. Ultimately, multiple processes appear to play some role in maintenance of this forced pool and the associated riffle. Copyright © 2008 John Wiley & Sons, Ltd.     

Radiative transfer in an infinite cylinder

The problem of radiation transfer in a cylinder with diffuse reflectivity and containing an energy source is connected with the source-free radiation transfer problem with isotropic boundary condition. Equation for the radiation heat flux is obtained for a polynomial source. In the special case of isotropic scattering, the radiation heat flux is given in terms of the albedo of the second problem. An expression is also given for the net radiation heat flux.