Recharge estimation for transient ground water modeling

Reliable ground water models require both an accurate physical representation of the system and appropriate boundary conditions. While physical attributes are generally considered static, boundary conditions, such as ground water recharge rates, can be highly variable in both space and time. A practical methodology incorporating the hydrologic model HELP3 in conjunction with a geographic information system was developed to generate a physically based and highly detailed recharge boundary condition for ground water modeling. The approach uses daily precipitation and temperature records in addition to land use/land cover and soils data. The importance of the method in transient ground water modeling is demonstrated by applying it to a MODFLOW modeling study in New Jersey. In addition to improved model calibration, the results from the study clearly indicate the importance of using a physically based and highly detailed recharge boundary condition in ground water quality modeling, where the detailed knowledge of the evolution of the ground water flowpaths is imperative. The simulated water table is within 0.5 m of the observed values using the method, while the water levels can differ by as much as 2 m using uniform recharge conditions. The results also show that the combination of temperature and precipitation plays an important role in the amount and timing of recharge in cooler climates. A sensitivity analysis further reveals that increasing the leaf area index, the evaporative zone depth, or the curve number in the model will result in decreased recharge rates over time, with the curve number having the greatest impact.     

Trace metals in the North Sea

The distribution of zinc, manganese, nickel, copper and cadmium in water and biological material has been measured in the southern North Sea and Straits of Dover, and off the north-east coast of England. These investigations provide baseline values against which future pollution changes can be measured.     

Preparation,collisional decay and transfer of alignment in C2H2 21

Stimulated Raman scattering is used to prepare ensembles of aligned acetylene molecules in thev ₂=1 excited vibrational level. Decay and transfer of the alignment is followed by laser induced fluorescence. Measurements of the initially excited alignment agree well with theoretical calculations. The rate constants for decay decrease as the rotational quantum numberJ of the excited level increases.     

A search for faint giant radio galaxies in the NVSS survey

Giant radio galaxies are the most extended radio sources in the universe having the sizes of groups and clusters of galaxies (about 1 Mpc). Their total number is comparable with the amount of clusters, revealing the Zeldovich-Sunyaev effect, which may lead to a biased estimation of the angular power spectrum for the ZS effect. To assess a possible contribution to the power spectrum, we need to collect the data on all the observed giant radio galaxies. We propose an algorithm for the selection of large sources (sized over 4′) from the NVSS survey catalog, using the parameters of possible components of radio sources, cataloged as compact objects. As a result of the first stage of our work, we have selected 61 new giant radio galaxy candidates, of which 35 have a pronounced morphological type FRII, and 26 galaxies—the FRI type.     

Phase equilibrium constraints on a deep crustal metamorphic field gradient: metapelitic rocks from the Ivrea Zone (NW Italy)

The metamorphic rocks of the Ivrea Zone in NW Italy preserve a deep crustal metamorphic field gradient. Application of quantitative phase equilibria methods to metapelitic rocks provides new constraints on the P–T conditions recorded in Val Strona di Omegna, Val Sesia and Val Strona di Postua. In Val Strona di Omegna, the metapelitic rocks show a structural and mineralogical change from mica‐schists with the common assemblage bi–mu–sill–pl–q–ilm ± liq at the lowest grades, through metatexitic migmatites (g–sill–bi–ksp–pl–q–ilm–liq) at intermediate grades, to complex diatexitic migmatites (g–sill–ru–bi–ksp–pl–q–ilm–liq) at the highest grades. Partial melting in the metapelitic rocks is consistent with melting via the breakdown of first muscovite then biotite. The metamorphic field gradient in Val Strona di Omegna is constrained to range from conditions of ～3.5–6.5 kbar at ≈650 °C to ～10–12 kbar at >900 °C. The peak P–T estimates, particularly for granulite facies conditions, are significantly higher than those of most earlier studies. In Val Sesia and Val Strona di Postua, cordierite‐bearing rocks record the effects of contact metamorphism associated with the intrusion of a large mafic body (the Mafic Complex). The contact metamorphism occurred at lower pressures than the regional metamorphic peak and overprints the regional metamorphic assemblages. These relationships are consistent with the intrusion of the Mafic Complex having post dated the regional metamorphism and are inconsistent with a model of magmatic underplating as the cause of granulite facies metamorphism in the region.     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.