Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

Assessment of the hydraulic slope flow approach using a mesoscale model

The simplified hydraulic two-layer model for a katabatic flow is analysed using the outputs from a high-resolution mesoscale simulation. A stably stratified night is simulated for the Duero basin, a complex terrain area located in the northern Spanish plateau, with large vertical and horizontal spatial resolution. Well-defined katabatic flows on the basin slopes are generated by the simulation, that are relatively stationary and quasi-bidimensional for some areas in the central part of the night. The bulk quantities used in the two-layer approach as well as the different terms in the equations are computed from the three-dimensional information provided by the mesoscale simulation. This method allows to inspect how well the simplified approach represents the katabatic flow generated by the mesoscale model. The study shows that the hydraulic model allows for a comprehensive analysis of the basic mechanisms of the slope flows but is not able to close the budget equations, since the residuals are large.     

Ar-Ar and I-Xe ages and thermal histories of three unusual metal-rich meteorites

Portales Valley, Sombrerete, and Northwest Africa (NWA) 176 are three unrelated meteorites, which consist of silicate mixed with substantial amounts of metal and which likely formed at elevated temperatures as a consequence of early impacts on their parent bodies. Measured ³⁹Ar-⁴⁰Ar ages of these meteorites are 4477 ± 11 Ma and 4458 ± 16 Ma (two samples of Portales Valley), 4541 ± 12 Ma, and 4524 ± 13 Ma, respectively (Ma = million years; all one-sigma errors). The Ar-Ar data for Portales Valley show no evidence of later open system behavior suggested by some other chronometers. Measured ¹²⁹I-¹²⁹Xe ages of these three meteorites are 4559.9 ± 0.5 Ma, 4561.9 ± 1.0 Ma, and ∼4544 Ma, respectively (relative to Shallowater = 4562.3 ± 0.4 Ma). From stepwise temperature release data, we determined the diffusion characteristics for Ar and Xe in our samples and calculated approximate closure temperatures for the K-Ar and I-Xe chronometers. Adopting results and interpretations about these meteorites from some previous workers, we evaluated all these data against various thermal cooling models. We conclude that Portales Valley formed 4560 Ma ago, cooled quickly to below the I-Xe closure temperature, then cooled deep within the parent body at a rate of ∼4 °C/Ma through K-Ar closure. We conclude that Sombrerete formed 4562 Ma ago and cooled relatively quickly. NWA 176 likely formed and cooled quickly ∼4544 Ma ago, or later than formation times of most meteorite parent bodies. For all three meteorites, the Ar-Ar ages are in better agreement with I-Xe ages and preferred thermal models if we increase these Ar-Ar ages by ∼20 Ma. Such age corrections would be consistent with probable errors in ⁴⁰K decay parameters in current use, as suggested by others. The role of impact heating and possible disruption and partial reassembly of meteorite parent bodies to form some meteorites likely was an important process in the early solar system.     

Opportunities accruing from the application of public private partnerships in down market urban housing in Kenya

GeoJournal - The purpose of this study was to investigate the opportunities which accrue from application of PPPs in development of down market urban housing in Kenya, hence attracting private...     

Renewed uplift at the Yellowstone Caldera measured by leveling surveys and satellite radar interferometry

 A first-order leveling survey across the northeast part of the Yellowstone caldera in September 1998 showed that the central caldera floor near Le Hardy Rapids rose 24±5 mm relative to the caldera rim at Lake Butte since the previous survey in September 1995. Annual surveys along the same traverse from 1985 to 1995 tracked progressive subsidence near Le Hardy Rapids at an average rate of –19±1 mm/year. Earlier, less frequent surveys measured net uplift in the same area during 1923–1976 (14±1 mm/year) and 1976–1984 (22±1 mm/year). The resumption of uplift following a decade of subsidence was first detected by satellite synthetic aperture radar interferometry, which revealed approximately 15 mm of uplift in the vicinity of Le Hardy Rapids from July 1995 to June 1997. Radar interferograms show that the center of subsidence shifted from the Sour Creek resurgent dome in the northeast part of the caldera during August 1992 to June 1993 to the Mallard Lake resurgent dome in the southwest part during June 1993 to August 1995. Uplift began at the Sour Creek dome during August 1995 to September 1996 and spread to the Mallard Lake dome by June 1997. The rapidity of these changes and the spatial pattern of surface deformation suggest that ground movements are caused at least in part by accumulation and migration of fluids in two sill-like bodies at 5–10 km depth, near the interface between Yellowstone's magmatic and deep hydrothermal systems. Received: 30 November 1998 / Accepted: 16 April 1999     

A meta-analysis of global crop water productivity of three leading world crops (wheat,corn, and rice) in the irrigated areas over three decades

ABSTRACT

The overarching goal of this study was to perform a comprehensive meta-analysis of irrigated agricultural Crop Water Productivity (CWP) of the world’s three leading crops: wheat, corn, and rice based on three decades of remote sensing and non-remote sensing-based studies. Overall, CWP data from 148 crop growing study sites (60 wheat, 43 corn, and 45 rice) spread across the world were gathered from published articles spanning 31 different countries. There was overwhelming evidence of a significant increase in CWP with an increase in latitude for predominately northern hemisphere datasets. For example, corn grown in latitude 40–50° had much higher mean CWP (2.45?kg/m³) compared to mean CWP of corn grown in other latitudes such as 30–40° (1.67?kg/m³) or 20–30° (0.94?kg/m³). The same trend existed for wheat and rice as well. For soils, none of the CWP values, for any of the three crops, were statistically different. However, mean CWP in higher latitudes for the same soil was significantly higher than the mean CWP for the same soil in lower latitudes. This applied for all three crops studied. For wheat, the global CWP categories were low (≤0.75?kg/m³), medium (>0.75 to <1.10?kg/m³), and high CWP (≥1.10?kg/m³). For corn the global CWP categories were low (≤1.25?kg/m³), medium (>1.25 to ≤1.75?kg/m³), and high (>1.75?kg/m³). For rice the global CWP categories were low (≤0.70?kg/m³), medium (>0.70 to ≤1.25?kg/m³), and high (>1.25?kg/m³). USA and China are the only two countries that have consistently high CWP for wheat, corn, and rice. Australia and India have medium CWP for wheat and rice. India’s corn, however, has low CWP. Egypt, Turkey, Netherlands, Mexico, and Israel have high CWP for wheat. Romania, Argentina, and Hungary have high CWP for corn, and Philippines has high CWP for rice. All other countries have either low or medium CWP for all three crops. Based on data in this study, the highest consumers of water for crop production also have the most potential for water savings. These countries are USA, India, and China for wheat; USA, China, and Brazil for corn; India, China, and Pakistan for rice. For example, even just a 10% increase in CWP of wheat grown in India can save 6974 billion liters of water. This is equivalent to creating 6974 lakes each of 100?m³ in volume that leads to many benefits such as acting as ‘water banks’ for lean season, recreation, and numerous ecological services. This study establishes the volume of water that can be saved for each crop in each country when there is an increase in CWP by 10%, 20%, and 30%.     

Electrodynamic Disaggregation: Does it Affect Apatite Fission-Track and (U-Th)/He Analyses?

Apatite fission-track and (U-Th)/He analyses require the liberation of intact idiomorphic apatite grains from rock samples. While routinely being carried out by mechanical methods, electrodynamic disaggregation (ED) offers an alternative approach. The high-voltage discharges produced during the ED process create localised temperature peaks (10000 K) along a narrow plasma channel. In apatite, such high temperatures could potentially reduce the length of fission tracks, which start to anneal at temperatures > 60 °C, and could also enhance He diffusion, which becomes significant at 30–40 °C over geological time scales. A comparison of fission-track analyses and (U-Th)/He ages of apatites prepared both by mechanical (jaw crusher, disk mill) and ED processing provides a way of determining whether heating during the latter method has any significant effect. Apatites from three samples of different geological settings (an orthogneiss from Madagascar, the Fish Canyon Tuff, and a muscovite-gneiss from Greece) yielded statistically identical track length distributions compared to samples prepared mechanically. Additionally, (U-Th)/He ages of apatites from a leucogranite from Morocco prepared by both methods were indistinguishable. These first results indicated that during electrodynamic disaggregation apatite crystals were not heated enough to partially anneal the fission tracks or induce significant diffusive loss of He.     

Dynamics of lava flow fronts, Arenal Volcano, Costa Rica

Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ₀) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time t_x can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.