Electron microprobe study of a ‘mysterite’-bearing inclusion from the Krymka LL-chondrite

The black inclusion from the Krymka LL3 chondrite previously found to contain ‘mysterite’ by Lewiset al. (1979) belongs to a hitherto unknown class of carbonaceous chondrites. Its olivine and pyroxene compositions. Fo 97–99 and En 96, respectively, are characteristic of carbonaceous chondrites and its plagioclase composition. An 100, is characteristic of C3's. It contains a peculiar group of Co-, Cr-rich metal grains whose compositions are similar, but not identical, to those in C2 chondrites and which also bear some similarities to those in Renazzo. Its weight ratios of total $\text{Fe}\text{SiO}{}_2$ and $\text{solSiO}{}_2\text{MgO}$ are 0.74 and 1.43, respectively, and its atomic ratio of $\text{Si}\text{Al}$ is 10.7, exactly the same as in carbonaceous chondrites. Its bulk chemical composition is very close to that of the Murchison C2 chondrite. The association of mysterite with a special type of carbonaceous chondrite material suggests that mysterite formed by low-temperature condensation in a different region of the nebula from other carbonaceous chondrites.     

pH influence on sorption characteristics of heavy metal in the vadose zone

Sorption is an important process in the modelling and prediction of the movement of heavy metals in unsaturated clay barriers. This experimental study investigates the effect of pH changes in the acidic range on the sorption characteristics of heavy metals such as: lead, copper and zinc in an unsaturated soil. A series of one-dimensional coupled solute and moisture leaching column tests, using different heavy metal solutions, were conducted on an unsaturated illitic soil at varying pH values. Variations of volumetric water content (VWC) with distance were measured for different time durations, and concentrations of heavy metals in the liquid and solid phases were analysed. Partitioning coefficient profiles of contaminants along the soil column were determined for each individual layer in the soil.

Results from column leaching tests showed that the sorption characteristics of heavy metals are controlled by many factors which should be taken into consideration, i.e. the VWC, time of wetting, soil pH, and the influent heavy metal concentrations. Simplification of K_d as a constant and of the VWC as a linear function cannot be considered a good assumption and may lead to an improper evaluation of the sorption phenomena and also to serious errors in predicting contaminant transport through unsaturated soils.     

艾比湖盐尘对周边地区土壤盐分及景观变化的影响

通过对准噶尔盆地西部艾比湖地区盐尘的扩散和堆积进行观测和取样分析,鉴别出粉尘中可溶盐的化学成分、pH和盐尘堆积强度,分析研究了盐尘的形成、分布特征和灾害强度以及盐尘暴与温度、降水、风等气候因子之间的相关关系.结果表明:近50 a来艾比湖的沉积环境总体是比较稳定的,但由于气候波动与人类活动引起艾比湖水位、水域面积曾发生较明显的变化,并产生了盐尘暴的发生、湖泊和土壤盐度以及湖周景观的变化.根据艾比湖周边盐漠区盐尘的时空分布规律、活动特点和危害强度分析,盐尘对准噶尔盆地西部土壤盐分及景观变化的影响强度范围和影响程度划分为3级区.     

The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Long-Term Trends of Nutrients and Phytoplankton in Chesapeake Bay

Climate effects on hydrology impart high variability to water-quality properties, including nutrient loadings, concentrations, and phytoplankton biomass as chlorophyll-a (chl-a), in estuarine and coastal ecosystems. Resolving long-term trends of these properties requires that we distinguish climate effects from secular changes reflecting anthropogenic eutrophication. Here, we test the hypothesis that strong climatic contrasts leading to irregular dry and wet periods contribute significantly to interannual variability of mean annual values of water-quality properties using in situ data for Chesapeake Bay. Climate effects are quantified using annual freshwater discharge from the Susquehanna River together with a synoptic climatology for the Chesapeake Bay region based on predominant sea-level pressure patterns. Time series of water-quality properties are analyzed using historical (1945–1983) and recent (1984–2012) data for the bay adjusted for climate effects on hydrology. Contemporary monitoring by the Chesapeake Bay Program (CBP) provides data for a period since mid-1984 that is significantly impacted by anthropogenic eutrophication, while historical data back to 1945 serve as historical context for a period prior to severe impairments. The generalized additive model (GAM) and the generalized additive mixed model (GAMM) are developed for nutrient loadings and concentrations (total nitrogen—TN, nitrate?+?nitrate—NO₂?+?NO₃) at the Susquehanna River and water-quality properties in the bay proper, including dissolved nutrients (NO₂?+?NO₃, orthophosphate—PO₄), chl-a, diffuse light attenuation coefficient (K _D (PAR)), and chl-a/TN. Each statistical model consists of a sum of nonlinear functions to generate flow-adjusted time series and compute long-term trends accounting for climate effects on hydrology. We present results identifying successive periods of (1) eutrophication ca. 1945–1980 characterized by approximately doubled TN and NO₂?+?NO₃ loadings, leading to increased chl-a and associated ecosystem impairments, and (2) modest decreases of TN and NO₂?+?NO₃ loadings from 1981 to 2012, signaling a partial reversal of nutrient over-enrichment. Comparison of our findings with long-term trends of water-quality properties for a variety of estuarine and coastal ecosystems around the world reveals that trends for Chesapeake Bay are weaker than for other systems subject to strenuous management efforts, suggesting that more aggressive actions than those undertaken to date will be required to counter anthropogenic eutrophication of this valuable resource.     

H2S fluxes from Mt. Etna, Stromboli, and Vulcano (Italy) and implications for the sulfur budget at volcanoes

We present here new measurements of sulfur dioxide and hydrogen sulfide emissions from Vulcano, Etna, and Stromboli (Italy), made by direct sampling at vents and by filter pack and ultraviolet spectroscopy in downwind plumes. Measurements at the F0 and FA fumaroles on Vulcano yielded SO₂/H₂S molar ratios of ≈0.38 and ≈1.4, respectively, from which we estimate an H₂S flux of 6 to 9 t · d⁻¹ for the summit crater. For Mt. Etna and Stromboli, we found SO₂/H₂S molar ratios of ≈20 and ≈15, respectively, which combined with SO₂ flux measurements, suggest H₂S emission rates of 50 to 113 t · d⁻¹ and 4 to 8 t · d⁻¹, respectively. We observe that “source” and plume SO₂/H₂S ratios at Vulcano are similar, suggesting that hydrogen sulfide is essentially inert on timescales of seconds to minutes. This finding has important implications for estimates of volcanic total sulfur budget at volcanoes since most existing measurements do not account for H₂S emission.     

Sedimentology of an ancient erg margin: the Lower Jurassic Aztec Sandstone, southern Nevada and southern California

The Lower Jurassic Aztec Sandstone is an aeolian-deposited quartzose sandstone that represents the western margin of the southerly-migrating Navajo-Nugget sand sea (or erg). Vertical and lateral facies relations suggest that the erg margin encroached upon volcanic highlands, alluvial fan, wadi and sabkha environments. In southern Nevada, 700 m thick facies successions record the arrival of the Aztec sand sea. Initial erg sedimentation in the Valley of Fire consists of lenticular or tongue-shaped aeolian sand bodies interstratified with fluvially-deposited coarse sandstone and mudstone. Above, evaporite-rich fine sandstone and mudstone are overlain by thick, cross-stratified aeolian sandstone that shows an upsection increase in set thickness. The lithofacies succession represents aeolian sand sheets and small dunes that migrated over a siliciclastic sabkha traversed by ephemeral wadis. These deposits were ultimately buried by large dunes and draas of the erg. In the Spring Mountains, a similar facies succession also contains thin, lenticular volcaniclastic conglomerate and sandstone. These sediments represent the distal margin of an alluvial fan complex sourced from the west. Thin aeolian sequences are interbedded with volcanic flow rocks, ash-flow tuffs, debris flows, and fluvial deposits in the Mojave Desert of southern California. These aeolian strata represent erg migration up the eastern flanks of a magmatic arc. The westward diminution of aeolian-deposited units may reflect incomplete erg migration, thin accumulation of aeolian sediment succeptible to erosion, and stratigraphic dilution by arc-derived sediment. A two-part division of the Aztec erg is suggested by lithofacies associations, the size and geometry of aeolian cross-strata, and sediment dispersal data. The leading or downwind margin of the erg, here termed the fore-erg, is represented by a 10–100 m thick succession of isolated pods, lenses, and tongues of aeolian-deposited sediment encased in fluvial and sabkha deposits. Continued sand-sea migration brought large dunes and draas of the erg interior into the study area; these 150–500 m thick central-erg sediments buried the fore-erg deposits. The trailing, upwind margin of the erg is represented by back-erg deposits in northern Utah and Wyoming.