Petrology of the Chilliwack batholith,North Cascades,Washington: generation of calc-alkaline granitoids by melting of mafic lower crust with variable water fugacity

Calc-alkaline granitoid rocks of the Oligocene-Pliocene Chilliwack batholith, North Cascades, range from quartz diorites to granites (57–78% SiO₂), and are coeval with small gabbroic stocks. Modeling of major element, trace element, and isotopic data for granitoid and mafic rocks suggests that: (1) the granitoids were derived from amphibolitic lower crust having REE (rare-earth-element) and Sr-Nd isotopic characteristics of the exposed gabbros; (2) lithologic diversity among the granitoids is primarily the result of variable water fugacity during melting. The main effect of f_{H 2 O} variation is to change the relative proportions of plagioclase and amphibole in the residuum. The REE data for intermediate granitoids (quartz diorite-granodiorite; Eu/Eu^*=0.84–0.50) are modeled by melting with f_{H 2 O}<1 kbar, leaving a plagioclase + pyroxene residuum. In contrast, data for leucocratic granitoids (leuco-granodiorites and granites; Eu/Eu^* =1.0–0.54) require residual amphibole in the source and are modeled by melting with f_{H 2 O}=2–3 kbar. Consistent with this model, isotopic data for the granitoids show no systematic variation with rock type (⁸⁷Sr/⁸⁶Sr_i =0.7033–0.7043; Nd(0)=+3.3 to +5.5) and overlap significantly with data for the gabbroic rocks (⁸⁷Sr/⁸⁶Sr_i =0.7034–0.7040; Nd(0)=+3.3 to +6.9). The f_{H 2 O} variations during melting may reflect additions of H₂O to the lower crust from crystallizing basaltic magmas having a range of H₂O contents; Chillwack gabbros document the existence of such basalts. One-dimensional conductive heat transfer calculations indicate that underplating of basaltic magmas can provide the heat required for large-scale melting of amphibolitic lower crust, provided that ambient wallrock temperatures exceed 800°C. Based on lithologic and geochemical similarities, this model may be applicable to other Cordilleran batholiths.     

Ab initio valence force field calculations for quartz

We have derived valence force constants for the tetrahedral SiO₄ unit and the inter-tetrahedral SiOSi linkage from previous ab initio molecular orbital calculations on H₄SiO₄ and H₆Si₂O₇ using a split-valence polarized Gaussian basis set (6-31G^*), and used these to calculate the infrared and Raman active vibrational modes of α-quartz. The calculation gives frequencies approximately 15% greater than experiment, as expected from harmonic force constants obtained at this level of Hartree-Fock theory, but the calculation gives the correct distribution of modes within each frequency range. Calculated 28–30 Si and 16–18 O isotope shifts and pressure shifts to 6 GPa are also in reasonable agreement with experiment. We have also used our ab initio force field to calculate the vibrational spectrum for β-quartz. The results suggest either that inclusion of a torsional force constant is important for determining the stability of this high temperature polymorph, or that the β-quartz has a disordered structure with lower symmetry (P6₂) domains, as suggested by earlier diffraction studies.     

The species censorship problem: A general solution

The longevities of species constituting a statistical population have an underlying distribution whose form and parametric values reflect probabilities of origination and extinction through time. In the case that a part of the population is extant, the form of distribution and its parameters cannot be estimated directly from the longevity data without bias. Longevities of extant censored species are truncated and thus do not statistically represent the underlying distribution. The remaining uncensored species do not represent the true relative abundances of longevities. These biases can be defined from the probability densities for species longevitys and intervalr between successive originations of species. For realistic densities ofs andr, species with an intermediate longevity are preferentially censored. This simple, general result arises because the probability of censoring a species increases with its longevity, whereas the probability of censoring a given longevity varies with its relative abundance.     

Geochemical and anthropogenic enrichments of Mo in sediments from perennially oxic and seasonally anoxic lakes in Eastern Canada

We measured the vertical distributions of Mo, Fe, Mn, sulfide, sulfate, organic carbon, major ions, and pH in sediment porewater from one perennially oxic and three seasonally anoxic lacustrine basins in Eastern Canada, as well as those of Mo, acid volatile sulfide, Fe, Mn, Al, organic C, ²¹⁰Pb and ¹³⁷Cs in sediment cores from the same sites. The only input of anthropogenic Mo to these lakes comes from atmospheric deposition.The relatively monotonous distribution of Mo in the porewater of three seasonally anoxic basins suggests that Mo is not redistributed in the sedimentary column during periods of anoxia. In contrast, porewater Mo profiles obtained at three sampling dates in a perennially oxic basin display sharp Mo peaks below the sediment-water interface, indicating redistribution subsequent to deposition. Modeling of these latter porewater Mo profiles with a diagenetic reaction-transport equation coupled to comparisons among the various porewater and solid phase profiles reveal that Mo is released at 1-2 cm depth as a consequence of the reductive dissolution of Fe oxyhydroxides and scavenged both at the vicinity of the sediment-water interface, by re-adsorption onto authigenic Fe oxyhydroxides, and deeper in the sediments where dissolved sulfide concentrations are higher. The estimated rate constant for the adsorption of Mo onto Fe oxyhydroxides is 36 ± 45 cm³ mol⁻¹ s⁻¹.Diagenetic modeling indicates that authigenic Mo in sediments of the perennially oxic basin represents about one-third of the total solid phase Mo in the first cm below the sediment-water interface and only one tenth below this horizon. If we assume that no authigenic Mo is accumulated in the seasonally anoxic lake sediments we conclude that the sediment Mo concentrations, which are up to 3-16 times higher than the average lithogenic composition, depending on the lake, are mainly due to atmospheric deposition of anthropogenic Mo and not to the formation of authigenic Mo phases. Reconstructed historical records of the atmospheric Mo deposition indicate maximum values in the 1970s and 1980s and significant decreases since then. Emissions to the atmosphere associated with the smelting of non-ferrous ores and coal combustion appear to be the most important sources of anthropogenic Mo.     

Quantifying the effect of humic matter on the emission of mercury from artificial soil surfaces

A set of surface samples was created using purified laboratory grade sand treated with 0.05 μg/g Hg as the HgCl₂ salt and various concentrations of purified humic and fulvic acids. Emissions of elemental Hg from these substrates to the atmosphere were inversely correlated with the organic content of the samples (99% confidence level). The greatest differences in Hg emissions were found between samples containing the lowest concentrations of humic matter (0% versus 0.01% humic, and 0.01% versus 0.1% humic), only small differences in Hg flux were found to exist for samples with higher concentrations of humic acid (1%, 5%, and 100%). This effect was independent of the type of humic substance used, with both humic and fulvic acids showing an inhibitory effect on surface Hg emissions.     

Fluid circulation in a complex volcano-tectonic setting,inferred from self-potential and soil CO2 flux surveys: The Santa María–Cerro Quemado–Zunil volcanoes and Xela caldera (Northwestern Guatemala)

The region encompassing Santa María, Cerro Quemado, and Zunil volcanoes, close to Quetzaltenango, the second largest city of Guatemala, is volcanically and tectonically complex. In addition, the huge Xela caldera, about 20 km in diameter, crosses this area and links up to the important Zunil fault zone located between the three volcanoes. Two highly active geothermal sites, named Zunil-I and Zunil-II, are also located between these three volcanic edifices at the southeastern boundary of Xela caldera. In order to determine the permeability variations and the main structural discontinuities within this complex volcano-tectonic setting, self-potential and soil CO₂ flux measurements have been coupled, with a step of 20 m, along a 16.880 km-long profile crossing the entire area. Two shallow hydrothermal systems, with maximum lateral extensions of 1.5 km in diameter, are indicated by positive self-potential/elevation gradients below Santa María and Cerro Quemado volcanoes. Such small hydrothermal systems cannot explain the intense geothermal manifestations at Zunil-I and Zunil-II. Another minor hydrothermal system is indicated by self-potential measurements on the flank of Santa María along the edge of the Xela caldera. CO₂ flux measurements display slight variations inside the caldera and decreasing values crossing outside the caldera boundary. We hypothesize the presence of a magmatic body, inside the southeastern border of Xela caldera, to explain the deeper and more intense hydrothermal system manifested by the Zunil-I and the Zunil-II geothermal fields. This magmatic system may be independent from Santa María and Cerro Quemado volcanoes. Alternatively, the hypothesized Xela magmatic system could have a common magmatic origin with the Cerro Quemado dome complex, consistent with previous findings on regional gas emissions. Sectors bordering the Cerro Quemado dome complex also have high amplitude minima-short wavelength anomalies in self-potential, interpreted as preferential rain water infiltration along faults of major permeability, probably related with the most recent stages of Cerro Quemado dome growth.     

Potential implications of global warming and barrier island degradation on future hurricane inundation,property damages,and population impacted

Hurricane flooding is a leading natural threat to coastal communities. Recent evidence of sea level rise coupled with potential future global warming indicate that sea level rise will accelerate and hurricanes may intensify over the coming decades. In regions fronted by barrier islands, the protective capacity of these islands may diminish as they are degraded by rising sea level. Here we present a hydrodynamic and geospatial analysis of the relative role of barrier island degradation on potential future hurricane flooding. For the City of Corpus Christi, Texas, USA, hurricane flooding is projected to rise between 20% and 70% by the 2030s, resulting in an increase in property damages and impacted population. These findings indicate that adaptive management strategies should be developed and adopted for mitigating loss of natural barrier islands when these islands act as protective features for populated bayside communities. Finally, this study illustrates a method for applying models to forecast future storm protection benefits of barrier island restoration projects.     

Postglacial relative sea‐level observations from Ireland and their role in glacial rebound modelling

The British Isles have been the focus of a number of recent modelling studies owing to the existence of a high‐quality sea‐level dataset for this region and the suitability of these data for constraining shallow earth viscosity structure, local to regional ice sheet histories and the magnitude/timing of global meltwater signals. Until recently, the paucity of both glaciological and relative sea‐level (RSL) data from Ireland has meant that the majority of these glacial isostatic adjustment (GIA) modelling studies of the British Isles region have tended to concentrate on reconstructing ice cover over Britain. However, the recent development of a sea‐level database for Ireland along with emergence of new glaciological data on the spatial extent, thickness and deglacial chronology of the Irish Ice Sheet means it is now possible to revisit this region of the British Isles. Here, we employ these new data to constrain the evolution of the Irish Ice Sheet. We find that in order to reconcile differences between model predictions and RSL evidence, a thick, spatially extensive ice sheet of ～600–700 m over much of north and central Ireland is required at the LGM with very rapid deglaciation after 21 k cal. yr BP. Copyright © 2007 John Wiley & Sons, Ltd.