Two-liquid partition coefficients: Experimental data and geochemical implications

Partition coefficients for Cs, Ba, Sr, Ca, Mg, La, Sm, Lu, Mn, Ti, Cr, Ta, Zr, and P between immiscible basic and acidic liquids in the system K₂O-Al₂O₃-FeO-SiO₂ were experimentally determined at 1,180 °C and 1 atm. Phosphorus is most strongly enriched in the basic melt (by a factor of 10), followed by rare earth elements, Ta, Ca, Cr, Ti, Mn, Zr, Mg, Sr, and Ba (enriched by a factor of 1.5). Of the elements studied, only Cs is enriched in the acidic melt. The two-liquid partition coefficients of Zr, Ta, Sm, and Mn are constant for concentrations ranging from <0.1% to as high as 1 wt.-%, suggesting that Henry's law is applicable in silicate melts (at least for these elements) to concentrations well above typical trace element levels in rocks. The strong relative preference of many elements for the basic melt implies that the structural characteristics of basic melts more readily permit stable coordination of cations by oxygen. Partitioning of elements between crystal and liquid in a magma must therefore be influenced by the composition (and consequent structure) of the liquid.Application of the two-liquid partition coefficients to possible occurrences of liquid immiscibility in magmas reveals that typical basalt-rhyolite associations are probably not generated by two-liquid phase separation. However, liquid immiscibility cannot be discounted as a possible origin for lamprophyric rocks containing felsic segregations.     

Crystal size distribution (CSD) in rocks and the kinetics and dynamics of crystallization

Crystal-size in crystalline rocks is a fundamental measure of growth rate and age. And if nucleation spawns crystals over a span of time, a broad range of crystal sizes is possible during crystallization. A population balance based on the number density of crystals of each size generally predicts a log-linear distribution with increasing size. The negative slope of such a distribution is a measure of the product of overall population growth rate and mean age and the zero size intercept is nucleation density. Crystal size distributions (CSDs) observed for many lavas are smooth and regular, if not actually linear, when so plotted and can be interpreted using the theory of CSDs developed in chemical engineering by Randolph and Larson (1971). Nucleation density, nucleation and growth rates, and orders of kinetic reactions can be estimated from such data, and physical processes affecting the CSD (e.g. crystal fractionation and accumulation, mixing of populations, annealing in metamorphic and plutonic rocks, and nuclei destruction) can be gauged through analytical modeling. CSD theory provides a formalism for the macroscopic study of kinetic and physical processes affecting crystallization, within which the explicit affect of chemical and physical processes on the CSD can be analytically tested. It is a means by which petrographic information can be quantitatively linked to the kinetics of crystallization, and on these grounds CSDs furnish essential information supplemental to laboratory kinetic studies. In this three part series of papers, Part I provides the general CSD theory in a geological context, while applications to igneous and metamorphic rocks are given, respectively, in Parts II and III.     

Stable isotope and petrologic evidence for open-system degassing during the climactic and pre-climactic eruptions of Mt. Mazama, Crater Lake, Oregon

Evaluation of the extent of volatile element recycling in convergent margin volcanism requires delineating likely source(s) of magmatic volatiles through stable isotopic characterization of sulfur, hydrogen and oxygen in erupted tephra with appropriate assessment of modification by degassing. The climactic eruption of Mt. Mazama ejected approximately 50 km³ of rhyodacitic magma into the atmosphere and resulted in formation of a 10-km diameter caldera now occupied by Crater Lake, Oregon (lat. 43°N, long. 122°W). Isotopic compositions of whole-rocks, matrix glasses and minerals from Mt. Mazama climactic, pre-climactic and postcaldera tephra were determined to identify the likely source(s) of H₂O and S. Integration of stable isotopic data with petrologic data from melt inclusions has allowed for estimation of pre-eruptive dissolved volatile concentrations and placed constraints on the extent, conditions and style of degassing.Sulfur isotope analyses of climactic rhyodacitic whole rocks yield δ³⁴S values of 2.8-14.8‰ with corresponding matrix glass values of 2.4-13.2‰. δ³⁴S tends to increase with stratigraphic height through climactic eruptive units, consistent with open-system degassing. Dissolved sulfur concentrations in melt inclusions (MIs) from pre-climactic and climactic rhyodacitic pumices varies from 80 to 330 ppm, with highest concentrations in inclusions with 4.8-5.2 wt% H₂O (by FTIR). Up to 50% of the initial S may have been lost through pre-eruptive degassing at depths of 4-5 km. Ion microprobe analyses of pyrrhotite in climactic rhyodacitic tephra and andesitic scoria indicate a range in δ³⁴S from −0.4‰ to 5.8‰ and from −0.1‰ to 3.5‰, respectively. Initial δ³⁴S values of rhyodacitic and andesitic magmas were likely near the mantle value of 0‰. Hydrogen isotope (δD) and total H₂O analyses of rhyodacitic obsidian (and vitrophyre) from the climactic fall deposit yielded values οf −103 to −53‰ and 0.23-1.74 wt%, respectively. Values of δD and wt% H₂O of obsidian decrease towards the top of the fall deposit. Samples with depleted δD, and mantle δ¹⁸O values, have elevated δ³⁴S values consistent with open-system degassing. These results imply that more mantle-derived sulfur is degassed to the Earth’s atmosphere/hydrosphere through convergent margin volcanism than previously attributed. Magmatic degassing can modify initial isotopic compositions of sulfur by >14‰ (to δ³⁴S values of 14‰ or more here) and hydrogen isotopic compositions by 90‰ (to δD values of −127‰ in this case).     

A-type granites: geochemical characteristics,discrimination and petrogenesis

New analyses of 131 samples of A-type (alkaline or anorogenic) granites substantiate previously recognized chemical features, namely high SiO₂, Na₂O+K₂O, Fe/Mg, Ga/Al, Zr, Nb, Ga, Y and Ce, and low CaO and Sr. Good discrimination can be obtained between A-type granites and most orogenic granites (M-, I and S-types) on plots employing Ga/Al, various major element ratios and Y, Ce, Nb and Zr. These discrimination diagrams are thought to be relatively insensitive to moderate degrees of alteration. A-type granites generally do not exhibit evidence of being strongly differentiated, and within individual suites can show a transition from strongly alkaline varieties toward subalkaline compositions. Highly fractionated, felsic I- and S-type granites can have Ga/Al ratios and some major and trace element values which overlap those of typical A-type granites.A-type granites probably result mainly from partial melting of F and/or Cl enriched dry, granulitic residue remaining in the lower crust after extraction of an orogenic granite. Such melts are only moderately and locally modified by metasomatism or crystal fractionation. A-type melts occurred world-wide throughout geological time in a variety of tectonic settings and do not necessarily indicate an anorogenic or rifting environment.Geological Survey of Canada contribution no. 18886     

Resolving the contributing factors to Mississippi Delta subsidence: Past and Present

To date, quantification of individual components that contribute to shallow and deep‐seated subsidence in passive margin deltas worldwide has proven problematic. A new, regional gridded chronostratigraphic dataset for the Lower Mississippi Delta region, derived from 80,928 well reports across the northern Gulf of Mexico (GOM), has bridged the disparity between geodetic mean rates measuring total land surface subsidence across annual‐to‐decadal timescales and the deep‐seated stratigraphic subsidence rates that record isostatic response over timescales of >10⁴ years. Through a quantitative assessment of gridded chronostratigraphic surfaces, sections, and subsidence rates extending from the Middle Pleistocene (0.58 Ma) to the Late Pliocene (3.85 Ma), we identify both temporal and spatial variability in deep‐seated subsidence across the northern GOM. Targeted deep‐seated subsidence data extracted across prior GOM Holocene sea‐level sample locations have revealed more than an order of magnitude greater rates of isostatic compensation in the Mississippi depocentre versus similar GOM sea‐level control sites in Florida and Alabama, casting doubt on efforts towards a representative Holocene sea‐level curve. Spatial variability in subsidence was also assessed locally in both the strike and dip directions to assess the contributions of growth faults. Fault throw displacement magnitude was discovered to decrease with depth, accounting for less than half of the total deep‐seated subsidence record of the Middle Pleistocene. Temporal subsidence complexities were also revealed including a direct, inverse logarithmic relationship between subsidence rate and time indicating variable subsidence component controls across different timescales. Despite the spatial and temporal complexities, this dataset serves as the first regional baseline for deep‐seated subsidence rates across the northern GOM.     

Using hyperspectral imagery to predict post-wildfire soil water repellency

A principal task of evaluating large wildfires is to assess fire's effect on the soil in order to predict the potential watershed response. Two types of soil water repellency tests, the water drop penetration time (WDPT) test and the mini-disk infiltrometer (MDI) test, were performed after the Hayman Fire in Colorado, in the summer of 2002 to assess the infiltration potential of the soil. Remotely sensed hyperspectral imagery was also collected to map post-wildfire ground cover and soil condition. Detailed ground cover measurements were collected to validate the remotely sensed imagery and to examine the relationship between ground cover and soil water repellency. Percent ash cover measured on the ground was significantly correlated to WDPT (r = 0.42; p-value < 0.0001), and the MDI test (r = − 0.37; p-value < 0.0001). A Mixture Tuned Matched Filter (MTMF) spectral unmixing algorithm was applied to the hyperspectral imagery, which produced fractional cover maps of ash, soil, and scorched and green vegetation. The remotely sensed ash image had significant correlations to the water repellency tests, WDPT (r = 0.24; p-value = 0.001), and the MDI test (r = − 0.21; p-value = 0.005). An iterative threshold analysis was also applied to the ash and water repellency data to evaluate the relationship at increasingly higher levels of ash cover. Regression analysis between the means of grouped data: MDI time vs. ash cover data (R² =0.75) and vs. Ash MTMF scores (R² = 0.63) yielded significantly stronger relationships. From these results we found on-the-ground ash cover greater than 49% and remotely sensed ash cover greater than 33% to be indicative of strongly water repellent soils. Combining these results with geostatistical analyses indicated a spatial autocorrelation range of 15 to 40 m. Image pixels with high ash cover (> 33%), including pixels within 15 m of these pixel patches, were used to create a likelihood map of soil water repellency. This map is a good indicator of areas where soil experienced severe fire effects—areas that likely have strong water repellent soil conditions and higher potential for post-fire erosion.     

Vesiculation of high fountaining Hawaiian eruptions: episodes 15 and 16 of 1959 Kīlauea Iki

The 1959 summit eruption of Kīlauea volcano produced the highest recorded Hawaiian fountain in Hawai‘i. Quantitative analysis of closely spaced samples from the final two high-fountaining episodes of the eruption result in a fine-scale textural study of pyroclasts and provide a record of postfragmentation processes. As clast vesicularity increases, the vesicle number density decreases and vesicle morphology shifts from small and round to larger and more irregular. The shift in microtexture corresponds to greater degrees of postfragmentation expansion of clasts with higher vesicularity. We suggest the range of clast morphologies in the deposit is related to thermal zonation within a Hawaiian fountain where the highest vesicularity clasts traveled in the center and lowest traveled along the margins. Vesicle number densities are greatest in the highest fountaining episode and therefore scale with intensity of activity. Major element chemical analyses and fasciculate crystal textures indicate microlite-rich zones within individual clasts are portions of recycled lava lake material that were incorporated into newly vesiculating primary melt.     

Historical variability and feedbacks among land cover,stream power,and channel geometry along the lower Canadian River floodplain in Oklahoma

In 1820, the lower Canadian River meandered through a densely forested floodplain. By 1898, most of the floodplain had been cleared for agriculture and changes in channel geometry and specific stream power followed, particularly channel widening and straightening with a lower potential specific stream power. In 1964, a large upstream hydropower dam was constructed, which changed the flow regime in the lower Canadian River and consequently the channel geometry. Without destructive overbank floods, the channel narrowed rapidly and considerably due to encroachment by floodplain vegetation. The lower Canadian River, which was once a highly dynamic floodplain‐river system, has now been transformed into a relatively static river channel. These changes over the past 200 years have not been linear or independent. In this article, we use a variety of data sources to assess these historical changes along the lower Canadian River floodplain and identify feedbacks among floodplain cultivation, dam construction, specific stream power, and channel width, slope, and sinuosity. Finally, we combine the results of our study with others in the region to present a biogeomorphic response model for large Great Plains rivers that characterizes channel width changes in response to climate variability and anthropogenic disturbances. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of historical landslide time series in the Emilia‐Romagna region,northern Italy

A catalogue of historical landslides, 1951–2002, for three provinces in the Emilia‐Romagna region of northern Italy is presented and its statistical properties studied. The catalogue consists of 2255 reported landslides and is based on historical archives and chronicles. We use two measures for the intensity of landsliding over time: (i) the number of reported landslides in a day (D_L) and (ii) the number of reported landslides in an event (S_event), where an event is one or more consecutive days with landsliding. From 1951–2002 in our study area there were 1057 days with 1 ≤ D_L ≤?45 landslides per day, and 596 events with 1 ≤ S_event ≤ 129 landslides per event. In the first set of analyses, we find that the probability density of landslide intensities in the time series are power‐law distributed over at least two‐orders of magnitude, with exponent of about ?2·0. Although our data is a proxy for landsliding built from newspaper reports, it is the first tentative evidence that the frequency‐size of triggered landslide events over time (not just the landslides in a given triggered event), like earthquakes, scale as a power‐law or other heavy‐tailed distributions. If confirmed, this could have important implications for risk assessment and erosion modelling in a given area. In our second set of analyses, we find that for short antecedent rainfall periods, the minimum amount of rainfall necessary to trigger landslides varies considerably with the intensity of the landsliding (D_L and S_event); whereas for long antecedent periods the magnitude is largely independent of the cumulative amount of rainfall, and the largest values of landslide intensity are always preceded by abundant rainfall. Further, the analysis of the rainfall trend suggests that the trigger of landslides in the study area is related to seasonal rainfall. Copyright © 2010 John Wiley & Sons, Ltd.     

龙门山褶皱逆冲带中前山逆冲带之新期构造及断层运动学(英文)

The disastrous Wenchuan, Sichuan, earthquake (MS=8.0) on 12 May 2009 ruptured several major thrust faults of the Longmenshan fold-and-thrust belt, along the western margin of the Sichuan Basin.