Acid-base status of soils in groundwater discharge zones — relation to surface water acidification

Critical load calculations have suggested that groundwater at depth of 2 m in Sweden is very sensitive to acid load. As environmental isotope studies have shown that most of the runoff in streams has passed through the soil, there is a risk in the near future of accelerated acidification of surface waters.

To assess the importance of the last soil horizon of contact before discharge, the upper 0–0.2m of soils in seven discharge zones were analysed for pools of base cations, acidity and base saturation. The sites were about 3–4 m² in size and selected from two catchments exposed to different levels of acid deposition.

The soils in the seven sites had high concentrations of exchangeable base cations and consequently high base saturation. The high correlation (r² = 0.74) between base saturation in the soils of the discharge zones and mean pH of the runoff waters suggested that the discharge zone is important for surface water acidification. The high pool of exchangeable base cations will buffer initially against the acid load. As the cation exchange capacity (meq dm⁻³) and base saturation were lower in the sites from the catchment receiving lower deposition, these streams may be more vulnerable to acidification in the near future. The high concentration of base cations in non-exchangeable fractions may also buffer against acidification as it is likely that some of these pools will become exchangeable with time.     

The diurnal and semidiurnal barometric oscillations global distribution and annual variation

Summary The global distributions of the annual and seasonal means of the diurnal (S ₁) and semidiurnal (S ₂) surface pressure oscillations are investigated by spherical harmonic analysis. The main waves are,S ₁ ¹ (with wave number 1) forS ₁ andS ₂ ² forS ₂.S ₁ ¹ is much less predominant among the waves ofS ₁ thanS ₂ ² among those ofS ₂. As in the case of the lunar semidiurnal barometric tideL ₂ the pressure maxima occur earlier in the Southern than in the Northern Hemisphere. In the case ofS ₂ the standing waveS ₂ ⁰ and the waveS ₂ ³ are also of interest besidesS ₂ ². Although the present analysis extends only from 60°N to 60°S, whileS ₂ ⁰ is largest at polar latitudes, its results show thatS ₂ ⁰ should be smaller at high southerly than at high northerly latitudes, as has been observed. Thus this observed asymmetrical distribution ofS ₂ ⁰ may be due to causes outside the polar regions rather than to their geographical differences. The best approximation to the observed distribution ofS ₂ ⁰ is obtained by including a mode representing an oscillation independent of longitude and latitude indicating a small semidiurnal variation of the mean global surface presure, which is an unlikely result on physical grounds.The seasonal variation ofS ₁ ¹ expressed in percent of the annual mean is smaller than that ofS ₂ ², and both are less than the unexplained seasonal variation ofL ₂ ².The main wavesS ₁ ¹ andS ₂ ² are expressed not only by associated Legendre functions, but also by Hough functions.National Center for Atmospheric Research, Boulder, Colorado, U.S.A., sponsored by the National Science Foundation.     

Studying Regions by Studying Firms*

This paper presents a method for inferring regional economic structure and prospects from key informant interviews. It describes assumptions which must be made and steps to be taken in moving from individuals to establishment, firm, and industry aggregations. A technique is offered for mapping core relationships among firms and other regional and external actors, and for aggregating firm-centered maps into industrial and regional maps. The mapping technique can be expanded to incorporate historicity and dynamic change.     

The Uppsala urban meteorology project

An extensive urban meteorology project in Uppsala, Sweden, is described. The city itself is considered suitable as a kind of model city, being almost circular, having sharp urban/rural boundaries and being situated in a relatively flat area. The project incorporates numerous measurements. A main station with a 100-m mast is situated at the NE urban/rural border. Its instrumentation consists of a slow-response (profile) system and a turbulence system. The profile instrumentation is described in detail. In addition to the main station, the project comprises: measurements from two 14-m masts, one on the top of a centrally located building, the other mobile; pilot balloon ascents; tethered balloon soundings, car-borne temperature traverses etc. Some preliminary results are presented: analyses of wind profiles from the mast and three-dimensional heat island studies.     

Upscaling Permeability Measurements Within Complex Heterolithic Tidal Sandstones

We investigate numerically the effect of sample volume on the effective single-phase permeability of heterolithic tidal sandstones, using three-dimensional models reconstructed directly from large rock specimens measuring 45 × 30 × 15 cm. We find that both individual and averaged effective permeability values vary as a function of sample volume, which indicates that permeability data obtained from core-plugs will not be representative at the scale of a reservoir model grid-block regardless of the number of measurements taken. However, the error introduced by averaged data may be minimized using the appropriate averaging scheme for a given facies type and flow direction.     

Effects of seawater carbonate ion concentration and temperature on shell U, Mg, and Sr in cultured planktonic foraminifera

We investigate the sensitivity of U/Ca, Mg/Ca, and Sr/Ca to changes in seawater [CO₃²⁻] and temperature in calcite produced by the two planktonic foraminifera species, Orbulina universa and Globigerina bulloides, in laboratory culture experiments. Our results demonstrate that at constant temperature, U/Ca in O. universa decreases by 25 ± 7% per 100 μmol [CO₃²⁻] kg⁻¹, as seawater [CO₃²⁻] increases from 110 to 470 μmol kg⁻¹. Results from G. bulloides suggest a similar relationship, but U/Ca is consistently offset by ∼+40% at the same environmental [CO₃²⁻]. In O. universa, U/Ca is insensitive to temperature between 15°C and 25°C. Applying the O. universa relationship to three U/Ca records from a related species, Globigerinoides sacculifer, we estimate that Caribbean and tropical Atlantic [CO₃²⁻] was 110 ± 70 μmol kg⁻¹ and 80 ± 40 μmol kg⁻¹ higher, respectively, during the last glacial period relative to the Holocene. This result is consistent with estimates of the glacial-interglacial change in surface water [CO₃²⁻] based on both modeling and on boron isotope pH estimates. In settings where the addition of U by diagenetic processes is not a factor, down-core records of foraminiferal U/Ca have potential to provide information about changes in the ocean’s carbonate concentration.Below ambient pH (pH < 8.2), Mg/Ca decreased by 7 ± 5% (O. universa) to 16 ± 6% (G. bulloides) per 0.1 unit increase in pH. Above ambient pH, the change in Mg/Ca was not significant for either species. This result suggests that Mg/Ca-based paleotemperature estimates for the Quaternary, during which surface-ocean pH has been at or above modern levels, have not been biased by variations in surface-water pH. Sr/Ca increased linearly by 1.6 ± 0.4% per 0.1 unit increase in pH. Shell Mg/Ca increased exponentially with temperature in O. universa, where Mg/Ca = 0.85 exp (0.096*T), whereas the change in Sr/Ca with temperature was within the reproducibility of replicate measurements.     

Magnetic anisotropy of calcite at room-temperature

The intrinsic room temperature magnetic properties of pure calcite were determined from a series of natural crystals, and they were found to be highly dependent on the chemical composition. In general, dia-, para-, and ferromagnetic components contribute to the magnetic susceptibility and the anisotropy of magnetic susceptibility (AMS). With a combination of magnetic measurements and chemical analysis these three contributions were determined and related to their mineralogical sources. The intrinsic diamagnetic susceptibility of pure calcite is − 4.46 ± 0.16 × 10^− 9 m³/kg (− 12.09 ± 0.5 × 10^− 6 SI) and the susceptibility difference is 4.06 ± 0.03 × 10^− 10 m³/kg (1.10 ± 0.01 × 10^− 6 SI). These diamagnetic properties are easily dominated by other components. The paramagnetic contribution is due to paramagnetic ions in the crystal lattice that substitute for calcium; these are mainly iron and manganese. The measured paramagnetic susceptibility agrees with the values calculated from the known concentration of paramagnetic ions in the crystals according to the Curie law of paramagnetic susceptibility. Substituted iron leads to an increase in the AMS. The paramagnetic susceptibility difference was found to correlate linearly with the iron content for concentrations between 500 and 10,000 ppm. An empirical relation was determined: (k₁ − k₃)^para (kg/m³) = Fe-content (ppm) × (1 ± 0.1) × 10^− 12 (kg/m³/ppm). The maximum susceptibility difference (Δk = k₁ − k₃) was found to be unaffected by iron contents below 100 ppm. Ferromagnetic contributions due to inclusions of ferromagnetic minerals can dominate the susceptibility. They were detected by acquisition of isothermal remanent magnetization (IRM) and their contribution to the AMS was separated by high-field measurements.     

Optimization strategies for sediment reduction practices on roads in steep,forested terrain

Many forested steeplands in the western United States display a legacy of disturbances due to timber harvest, mining or wildfires, for example. Such disturbances have caused accelerated hillslope erosion, leading to increased sedimentation in fish‐bearing streams. Several restoration techniques have been implemented to address these problems in mountain catchments, many of which involve the removal of abandoned roads and re‐establishing drainage networks across road prisms. With limited restoration funds to be applied across large catchments, land managers are faced with deciding which areas and problems should be treated first, and by which technique, in order to design the most effective and cost‐effective sediment reduction strategy. Currently most restoration is conducted on a site‐specific scale according to uniform treatment policies. To create catchment‐scale policies for restoration, we developed two optimization models – dynamic programming and genetic algorithms – to determine the most cost‐effective treatment level for roads and stream crossings in a pilot study basin with approximately 700 road segments and crossings. These models considered the trade‐offs between the cost and effectiveness of different restoration strategies to minimize the predicted erosion from all forest roads within a catchment, while meeting a specified budget constraint. The optimal sediment reduction strategies developed by these models performed much better than two strategies of uniform erosion control which are commonly applied to road erosion problems by land managers, with sediment savings increased by an additional 48 to 80 per cent. These optimization models can be used to formulate the most cost‐effective restoration policy for sediment reduction on a catchment scale. Thus, cost savings can be applied to further restoration work within the catchment. Nevertheless, the models are based on erosion rates measured on past restoration sites, and need to be updated as additional monitoring studies evaluate long‐term basin response to erosion control treatments. Copyright © 2006 John Wiley & Sons, Ltd.     

Dispersion and cost analysis of some finite difference schemes in one-parameter acoustic wave modeling

Given a precision threshold to be imposed on the group velocity error and a finite difference scheme for the acoustic wave equation, it is possible to determine time-step and grid-spacing in an optimal manner, i.e., so as to minimize the computational cost. Using this optimal cost as a criterion, it becomes easy to compare schemes for efficiency in homogeneous media. Heterogeneous media with constant density can be accommodated to a certain extent by minimizing the cost over a range of Courant numbers. Such analysis shows that, amongst the second-order Taylor series schemes in time, higher-order schemes are generally more efficient than lower-order schemes. However, this result does not extend to very high order schemes.