Trace element partitioning in the Fe-S-C system and its implications for planetary differentiation and the thermal history of ureilites

Element partitioning in metal-light element systems is important to our understanding of planetary differentiation processes. In this study, solid-metal/liquid-sulfide, liquid-metal/liquid-sulfide and solid-metal/troilite partition coefficients (D) were determined for 18 elements (Ag, As, Au, Co, Cr, Cu, Ge, Ir, Ni, Os, Pd, Pt, Mo, Mn, Re, Ru, Se and W) in the graphite-saturated Fe-S-C system at 1 atm. Compared at the same liquid S concentration, the solid/liquid partition coefficients are similar to those in the Fe-S system, but there are systematic differences that appear to be related to interactions with carbon dissolved in the solid metal. Elements previously shown to be “anthracophile” generally have larger solid/liquid partition coefficients in the Fe-S-C system, whereas those that are not have similar or smaller partition coefficients in the Fe-S-C system. The partitioning of trace elements between C-rich and S-rich liquids is, in most cases, broadly similar to the partitioning between solid metal and S-rich liquid. The highly siderophile elements Os, Re, Ir and W are partitioned strongly into the C-rich liquid, with D ? 100. The partition coefficients for Pt, Ge and W decrease significantly at the transition to liquid immiscibility, while the partition coefficient for Mo increases sharply. The bulk siderophile element patterns of ureilite meteorities appear to be better explained by separation of S-rich liquid from residual C-rich metallic liquid at temperatures above the silicate solidus, rather than by separation of S-rich liquid from residual solid metal at lower temperatures.     

The timing and magnitude of changes to Hortonian overland flow at the watershed scale during the post-fire recovery process

Extreme hydrologic responses following wildfires can lead to floods and debris flows with costly economic and societal impacts. Process-based hydrologic and geomorphic models used to predict the downstream impacts of wildfire must account for temporal changes in hydrologic parameters related to the generation and subsequent routing of infiltration-excess overland flow across the landscape. However, we lack quantitative relationships showing how parameters change with time-since-burning, particularly at the watershed scale. To assess variations in best-fit hydrologic parameters with time, we used the KINEROS2 hydrological model to explore temporal changes in hillslope saturated hydraulic conductivity (K_sh) and channel hydraulic roughness (n_c) following a wildfire in the upper Arroyo Seco watershed (41.5 km²), which burned during the 2009 Station fire in the San Gabriel Mountains, California, USA. This study explored runoff-producing storms between 2008 and 2014 to infer watershed hydraulic properties by calibrating the model to observations at the watershed outlet. Modelling indicates K_sh is lowest in the first year following the fire and then increases at an average rate of approximately 4.2 mm/h/year during the first 5 years of recovery. The estimated values for K_sh in the first year following the fire are similar to those obtained in previous studies on smaller watersheds (<1.5 km²) following the Station fire, suggesting hydrologic changes detected here can be applied to lower-order watersheds. Hydraulic roughness, n_c, was lowest in the first year following the fire, but increased by a factor of 2 after 1 year of recovery. Post-fire observations suggest changes in n_c are due to changes in grain roughness and vegetation in channels. These results provide quantitative constraints on the magnitude of fire-induced hydrologic changes following severe wildfires in chaparral-dominated ecosystems as well as the timing of hydrologic recovery.     

Climatological Drought in Arizona: An Analysis of Indicators for Guiding the Governor's Drought Task Force

Lacking a federal policy to address local water deficiencies within the United States, many states have developed their own methods for monitoring drought in an effort to mitigate its effects. This article provides an overview of efforts to standardize the use of drought indices in order to compare recent and historical drought both spatially and temporally for Arizona. Yearly averages of the Palmer Drought Severity Index and the Standardized Precipitation Index are placed into a frequency distribution to create standardization among the indices. The 1896–1904 drought ranked as the most severe, although the 1996–2004 drought was a close second.     

Spatial patterns in thunderstorm rainfall events and their coupling with watershed hydrological response

Weather radar systems provide detailed information on spatial rainfall patterns known to play a significant role in runoff generation processes. In the current study, we present an innovative approach to exploit spatial rainfall information of air mass thunderstorms and link it with a watershed hydrological model. Observed radar data are decomposed into sets of rain cells conceptualized as circular Gaussian elements and the associated rain cell parameters, namely, location, maximal intensity and decay factor, are input into a hydrological model. Rain cells were retrieved from radar data for several thunderstorms over southern Arizona. Spatial characteristics of the resulting rain fields were evaluated using data from a dense rain gauge network. For an extreme case study in a semi-arid watershed, rain cells were derived and fed as input into a hydrological model to compute runoff response. A major factor in this event was found to be a single intense rain cell (out of the five cells decomposed from the storm). The path of this cell near watershed tributaries and toward the outlet enhanced generation of high flow. Furthermore, sensitivity analysis to cell characteristics indicated that peak discharge could be a factor of two higher if the cell was initiated just a few kilometers aside.     

Peat extraction, conservation and the rural economy in Northern Ireland

Northern Ireland has a tradition of peat cutting; 78 per cent by area of remaining lowland bogs and 46 per cent of blanket bogs have been cut in the past. Since the early 1980s, machine peat extraction for fuel and horticulture has increased, against a background of changing economic support for agriculture, high rural unemployment and agricultural underemployment. Cutting fuel peat can reduce household fuel costs and some peat fuel is sold to gain income. Horticultural peat is cut mainly for sale. Expansion of machine extraction, with possible environmental impacts, occurred as interest in peatland conservation increased. Potential conflict arises between local economic needs and conservation. Machine cutting for fuel is widespread, but affects only 2.6 per cent of the blanket peat area; horticultural extraction is more localized and is only one-fifth of the extent of fuel extraction. Altitude, distance from roads, land ownership and rights, turbary, religious affiliation, local economic needs and the role of government departments all contribute to explaining the distribution and extent of extraction.     

The Miller Range 090340 and 090206 meteorites: Identification of new brachinite‐like achondrites with implications for the diversity and petrogenesis of the brachinite clan

Miller Range (MIL) 090340 and MIL 090206 are olivine‐rich achondrites originally classified as ureilites. We investigate their petrography, mineral compositions, olivine Cr valences, equilibration temperatures, and (for MIL 090340) oxygen isotope compositions, and compare them with ureilites and other olivine‐rich achondrites. We conclude that they are brachinite‐like achondrites that provide new insights into the petrogenesis of brachinite clan meteorites. MIL 090340,6 has a granoblastic texture and consists of ~97 modal % by area olivine (Fo = molar Mg/[Mg+Fe] = 71.3 ± 0.6). It also contains minor to trace augite, chromite, chlorapatite, orthopyroxene, metal, troilite, and terrestrial Fe‐oxides. Approximately 80% by area of MIL 090206,5 has a granoblastic texture of olivine (Fo 72.3 ± 0.1) plus minor augite and chromite, similar to MIL 090340 but also containing minor plagioclase. The rest of the section consists of a single crystal of orthopyroxene (~11 × 3 mm), poikilitically enclosing rounded grains of olivine (Fo = 76.1 ± 0.6), augite, chromite, metal, and sulfide. Equilibration temperatures for MIL 090340 and MIL 090206, calculated from olivine‐spinel, olivine‐augite, and two‐pyroxene thermometry range from ~800 to 930 °C. In both samples, symplectic intergrowths of Ca‐poor orthopyroxene + opaque phases (Fe‐oxides, sulfide, metal) occur as rims on and veins/patches within olivine. Before terrestrial weathering, the opaques were probably mostly sulfide, with minor metal. All petrologic properties of MIL 090340 and MIL 090206 are consistent with those of brachinite clan meteorites, and largely distinct from those of ureilites. Oxygen isotope compositions of olivine in MIL 090340 (δ¹⁸O = 5.08 ± 0.30‰, δ¹⁷O = 2.44 ± 0.21‰, and Δ¹⁷O = ?0.20 ± 0.12‰) are also within the range of brachinite clan meteorites, and well distinguished from ureilites. Olivine Cr valences in MIL 090340 and the granoblastic area of MIL 090206 are 2.57 ± 0.06 and 2.59 ± 0.07, respectively, similar to those of three brachinites also analyzed here (Brachina, Hughes 026, Nova 003). They are higher than those of olivine in ureilites, even those containing chromite. The valence systematics of MIL 090340, MIL 090206, and the three analyzed brachinites (lower Fo = more oxidized Cr) are consistent with previous evidence that brachinite‐like parent bodies were inherently more oxidized than the ureilite parent body. The symplectic orthopyroxene + sulfide/metal assemblages in MIL 090340, MIL 090206, and many brachinite clan meteorites have superficial similarities to characteristic “reduction rims” in ureilites. However, they differ significantly in detail. They likely formed by reaction of olivine with S‐rich fluids, with only minor reduction. MIL 090340 and the granoblastic area of MIL 090206 are similar in modal mineralogy and texture to most brachinites, but have higher Fo values typical of brachinite‐like achondrites. The poikilitic pyroxene area of MIL 090206 is more typical of brachinite‐like achondrites. The majority of their properties suggest that MIL 090340 and MIL 090206 are residues of low‐degree partial melting. The poikilitic area of MIL 090206 could be a result of limited melt migration, with trapping and recrystallization of a small volume of melt in the residual matrix. These two samples are so similar in mineral compositions, Cr valence, and cosmic ray exposure ages that they could be derived from the same lithologic unit on a common parent body.     

Fractional melting and smelting on the ureilite parent body

We investigate petrologic and physical aspects of melt extraction on the parent asteroid of the ureilite meteorites (UPB). We first develop a petrologic model for simultaneous melting and smelting (reduction of FeO by C) at various depths. For a model starting composition, determined from petrologic constraints to have been CV-like except for elevated Ca/Al (2.5 × CI), we determine (1) degree of melting, (2) the evolution of mg, (3) production of CO + CO₂ gas and (4) the evolution of mineralogy in the residue as a function of temperature and pressure. We then use these relationships to examine implications of fractional vs. batch melt extraction.In the shallowest source regions (∼30 bars), melting and smelting begin simultaneously at ∼1050 °C, so that mg and the abundance of low-Ca pyroxene (initially pigeonite, ultimately pigeonite + orthopyroxene) begin to increase immediately. However, in the deepest source regions (∼100 bars), smelting does not begin until ∼1200 °C, so that mg begins to increase and low-Ca pyroxene (pigeonite) appears only after ∼21% melting. The final residues in these two cases, obtained just after the demise of augite, match the end-members of the ureilite mg range (∼94-76) in pyroxene abundance and type. In all source regions, production of CO + CO₂ by smelting varies over the course of melting. The onset of smelting results in a burst of gas production and very high incremental gas/melt ratios (up to ∼2.5 by mass); after a few % (s)melting, however, these values drastically decline (to <0.05 in the final increments).Physical modelling based on these relationships indicates that melts would begin to migrate upwards after only ∼1-2% melting, and thereafter would migrate continuously (fractionally) and rapidly (reaching the surface in < a year) in a network of veins/dikes. All melts produced during the smelting stage in each source region have gas contents sufficient to cause them to erupt explosively and be lost. However, since in all but the shallowest source regions part of the melting sequence occurs without smelting, fractional melting implies that a significant fraction of UPB melts may have erupted more placidly to form a thin crust (∼3.3 km thick for a 100 km radius body).Our calculations suggest that melt extraction was so rapid that equilibrium trace element partitioning may not have been attained. We present a model for disequilibrium fractional melting (in which REE partitioning is limited by diffusion) on the UPB, and demonstrate that it produces a good match to the ureilite data. The disequilibrium model may also apply to trace siderophile elements, and might help explain the “overabundance” of these elements in ureilites relative to predictions from the smelting model.Our results suggest that melt extraction on the UPB was a rapid, fractional process, which can explain the preservation of a primitive oxygen isotopic signature on the UPB.     

Interannual variations in the local spatial autocorrelation of tropospheric temperatures

Many dimensions of the global temperature pattern have been explored intensely over the past few decades. In this investigation, we explore the underlying spatial autoregressiveness of annual tropospheric temperature anomalies measured by polar-orbiting satellites. We found that the percent of the Earth covered by areas of significantly high local autocorrelation ranged from 11.92% to 25.90% over the 1979 to 2008 study period. We gathered 13 different teleconnection indices that have been linked to regional to global temperatures and found two (Polar/Eurasia pattern and North Atlantic Oscillation) that were positively correlated to the percent area with high local spatial autocorrelation. The magnitude but not the sign of the El Ni?o Southern Oscillation is also an important factor in the variation of spatial autocorrelation. The physical mechanism driving high local spatial autocorrelation is thought to be zonal and hemispheric flow governed by the subtropical and polar jet streams.