PGE geochemistry of the Eagle Ni–Cu–(PGE) deposit, Upper Michigan: constraints on ore genesis in a dynamic magma conduit

The Eagle Ni–Cu–(PGE) deposit is hosted in mafic–ultramafic intrusive rocks associated with the Marquette–Baraga dike swarm in northern Michigan. Sulfide mineralization formed in association with picritic magmatism in a dynamic magma conduit during the early stages in the development of the ~1.1?Ga Midcontinent Rift System. Four main types of sulfide mineralization have been recognized within the Eagle deposit: (1) disseminated sulfides in olivine-rich rocks; (2) rocks with semi-massive sulfides located both above and below the massive sulfide zone; (3) massive sulfides; and (4) sulfide veins in sedimentary country rocks. The disseminated, massive and lower semi-massive sulfide zones are typically composed of pyrrhotite, pentlandite and chalcopyrite, whereas the upper semi-massive sulfide ore zone also contains pyrrhotite, pentlandite, and chalcopyrite, but has higher cubanite content. Three distinct types of sulfide mineralization are present in the massive sulfide zone: IPGE-rich, PPGE-rich, and PGE-unfractioned. The lower and upper semi-massive sulfide zones have different PGE compositions. Samples from the lower semi-massive sulfide zone are characterized by unfractionated PGE patterns, whereas those from the upper semi-massive sulfide zone show moderate depletion in IPGE and moderate enrichment in PPGE. The mantle-normalized PGE patterns of unfractionated massive and lower semi-massive sulfides are subparallel to those of the disseminated sulfides. The results of numerical modeling using PGE concentrations recalculated to 100% sulfide (i.e., PGE tenors) and partition coefficients of PGE between sulfide liquid and magma indicate that the disseminated and unfractionated massive sulfide mineralization formed by the accumulation of primary sulfide liquids with similar R factors between 200 and 300. In contrast, the modeled R factor for the lower semi-massive sulfide zone is <100. The fractionated sulfide zones such as those of the IPGE-rich and PPGE-rich massive sulfides and the upper semi-massive sulfide zone can be explained by fractional crystallization of monosulfide solid solution from sulfide liquids. The results of numerical modeling indicate that the sulfide minerals in the upper semi-massive sulfide zone are the products of crystallization of fractionated sulfide liquids derived from a primary sulfide liquid with an R factor similar to that for the disseminated sulfide mineralization. Interestingly, the modeled parental sulfide liquid for the IPGE-rich and PPGE-rich massive sulfide zones has a higher R factor (~400) than that for the unfractionated massive sulfide mineralization. Except one sample which has unusually high IPGE and PPGE contents, all other samples from the Cu-rich sulfide veins in the footwall of the intrusion are highly depleted in IPGE and enriched in PPGE. These signatures are generally consistent with highly fractionated sulfide liquids expelled from crystallizing sulfide liquid. Collectively, our data suggest that at least four primary sulfide liquids with different R factors (<100, 200–300, ~400) were involved in the formation of the Eagle magmatic sulfide deposit. We envision that the immiscible sulfide liquids were transported from depth by multiple pulses of magma passing through the Eagle conduit system. The sulfide liquids were deposited in the widened part of the conduit system due to the decreasing velocity of magma flow. The presence of abundant olivine in some of the sulfide ore zones indicates that the ascending magma also carried olivine crystals. Sulfide saturation was attained in the parental magma due in large part to the assimilation of country rock sulfur at depth.     

Quantitative urban climate mapping based on a geographical database: A simulation approach using Hong Kong as a case study

The urban environment has been dramatically changed by artificial constructions. How the modified urban geometry affects the urban climate and therefore human thermal comfort has become a primary concern for urban planners. The present study takes a simulation approach to analyze the influence of urban geometry on the urban climate and maps this climatic understanding from a quantitative perspective. A geographical building database is used to characterize two widely discussed aspects: urban heat island effect (UHI) and wind dynamics. The parameters of the sky view factor (SVF) and the frontal area density (FAD) are simulated using ArcGIS-embedded computer programs to link urban geometry with the UHI and wind dynamic conditions. The simulated results are synergized and classified to evaluate different urban climatic conditions based on thermal comfort consideration. A climatic map is then generated implementing the classification. The climatic map shows reasonable agreement with thermal comfort understanding, as indicated by the biometeorological index of the physiological equivalent temperature (PET) obtained in an earlier study. The proposed climate mapping approach can provide both quantitative and visual evaluation of the urban environment for urban planners with climatic concerns. The map could be used as a decision support tool in planning and policy-making processes. An urban area in Hong Kong is used as a case study.     

Toward seamless indoor-outdoor applications: Developing stakeholder-oriented location-based services

Location-Based Services (LBS),an emerging new business based on smartphone and mobile networks,are becoming more and more popular.Most of these LBSs,however,only offer non-seamless indoor/outdoor applications and simple applications without giving stakeholders the chance to play an active role.Our specific aim is to solve these issues.This paper presents concepts to solve these issues by expanding the Open Location Services Interface Standard (OpenLS) to allow seamless indoor/outdoor positioning and to extend the content of the services to include information recommended by stakeholders.     

Sulfur biogeochemical cycling and novel Fe-S mineralization pathways in a tidally re-flooded wetland

Sulfur biogeochemical cycling and associated Fe-S mineralization processes exert a major influence over acidity dynamics, electron flow and contaminant mobility in wetlands, benthic sediments and groundwater systems. While S biogeochemical cycling has been studied intensively in many environmental settings, relatively little direct information exists on S cycling in formerly drained wetlands that have been remediated via tidal re-flooding. This study focuses on a tidal wetland that was drained in the 1970s (causing severe soil and water acidification), and subsequently remediated by controlled re-flooding in 2002. We examine reduction rates and Fe-S mineralization at the tidal fringe, 7 years after the commencement of re-flooding. The initial drainage of the wetland examined here caused in-situ pyrite (FeS₂) oxidation, resulting in the drained soil layers being highly acidic and rich in -bearing Fe(III) minerals, including jarosite (KFe₃(SO₄)₂(OH)₆). Tidal re-flooding has neutralized much of the previous acidity, with the pore-water pH now mostly spanning pH 5-7. The fastest rates of in-situ reduction (up to ∼300 nmol cm⁻³ day⁻¹) occur within the inter-tidal zone in the near-surface soil layers (to ∼60 cm below ground surface). The reduction rates correlate with pore-water dissolved organic C concentrations, thereby suggesting that electron donor supply was the predominant rate determining factor. Elemental S was a major short-term product of reduction, comprising up to 69% of reduced inorganic S in the near-surface soil layers. This enrichment in elemental S can be partly attributed to interactions between biogenic H₂S and jarosite - a process that also contributed to enrichment in pore-water Fe²⁺ (up to 55 mM) and (up to 50 mM). The iron sulfide thiospinel, greigite (Fe₃S₄), was abundant in near-surface soil layers within the inter- to sub-tidal zone where tidal water level fluctuations created oscillatory redox conditions. There was evidence for relatively rapid pyrite re-formation within the re-flooded soil layers. However, the results indicate that pyrite re-formation has occurred mainly in the lower formerly drained soil layers, whereas the accumulation of elemental S and greigite has been confined towards the soil surface. The discovery that pyrite formation was spatially decoupled from that of elemental S and greigite challenges the concept that greigite is an essential precursor required for sedimentary pyrite formation. In fact, the results suggest that greigite and pyrite may represent distinct end-points of divergent Fe-S mineralization pathways. Overall, this study highlights novel aspects of Fe-S mineralization within tidal wetlands that have been drained and re-flooded, in contrast to normal, undisturbed tidal wetlands. As such, the long-term biogeochemical trajectory of drained and acidified wetlands that are remediated by tidal re-flooding cannot be predicted from the well-studied behaviour of normal tidal wetlands.     

Trace metals in porewater of surface sediments and their bioavailability in Jiaozhou Bay,Qingdao, China

Samples of porewater and the edible tissue of mottled clams, Ruditapes variegatus, were collected simultaneously at each of 12 stations from Jiaozhou Bay near Qingdao, China, in June of 2003. Chemical analysis focused on trace metals and major elements. Porewater concentrations of Cd, Cr, Pb, and Zn were distributed in a bimodal fashion among the 12 stations, with mean concentrations at six of the stations being 8–32 times mean values at the other six stations. The concentrations of the same metals in clams were remarkably similar among stations, the coefficients of variation being only 12–37%. Calculations performed with the computer program PHREEQC indicated that Pb and Cr in porewater were present only in the +2 and +3 states, respectively, and because dissolved Cd and Zn exist only in the +2 state, the bimodal distribution of these four metals likely reflects secondary effects associated with their scavenging by Fe and/or Mn under oxidizing conditions and subsequent dissolution in the reducing environment of the porewaters. Consistent with this hypothesis is the fact that the distribution of high and low metal concentrations was closely correlated with the granularity of the sediment, with lower metal concentrations associated with relatively coarse sediment. Comparison of published biological concentration factors with the ratios of metal concentrations in the clams to porewater metal concentrations indicated that the porewaters were not the primary sources of the metals in the clams. Mixing processes in the bay likely account for the rather uniform concentrations of metals in clam tissue.     

Measured and modeled wave overtopping on a natural beach

The rate of wave overtopping of a barrier beach is measured and modeled. Unique rate of wave overtopping field data are obtained from the measure of the Carmel River, California, lagoon filling during a time when the lagoon is closed-off with no river inflow. Volume changes are based on measured lagoon height changes applied to a measured hypsometric curve. Wave heights and periods are obtained from directional wave spectra data in 15 m fronting the beach. Beach morphology was measured by GPS walking surveys. Three empirical overtopping models by Van der Meer and Janssen (1995), Hedges and Reis (1998) and Pullen et al. (2007) with differing parameterizations on wave height, period and beach slope and calibrated using extensive laboratory data obtained over plane, impermeable beaches are applied in a quasi-2D manner and compared with the field observations. Three overtopping events are considered when morphology data were available less than 2 weeks prior to the event. The models are tuned to fit the data using a reduction factor to account for beach permeability, berm characteristics, non-normal wave incidence and surface roughness influence. In addition, the run-up model by Stockdon et al. (2006) based on field data is examined and found to underestimate run-up as the calculated values were too small to predict any of the observed overtopping. The three overtopping models performed similarly well with values of 0.72–0.87 for the two narrow-banded wave cases, with an average reduction factor of 0.78. The European model (Pullen et. al., 2007) performed best overall and in particular for the case of the broad-banded, double peaked wave spectrum.     

Neighborhood covenants and restrictions: an examination of fit

Covenants and restrictions (C&Rs) have increased in popularity in recent decades, being commonly applied in new residential developments. Despite the fact that developers are more commonly writing C&Rs for their developments, little research has investigated how well C&Rs actually fit the preferences of residents. Recognizing this gap in the literature, this paper examines the fit of neighborhood C&Rs with resident preferences in Porter County, Indiana, a growing suburban/exurban place in the Chicago Metropolitan Statistical Area. The researcher conducted in-depth interviews with 51 residents to ascertain how well C&Rs fit resident preferences. The majority of interviewees stated that C&Rs match resident preferences because they believed that C&Rs promoted subdivision tidiness and worked to support property values. The minority view was that C&Rs do not fit because they are too strict and residents do not obey C&Rs. Interviewees who perceived a fit were more likely to have a lower level of engagement with C&Rs and were more likely to live in developments with no C&R enforcement attempts by the homeowners’ association (HOA). The results suggest the particular content of C&Rs does not seem to fit resident preference; rather the fact that C&Rs are not enforced fits resident preference. Seemingly, residents are content being unaware of the C&Rs, then they do not know if someone is breaking the rules or not. This may imply that residents would not want to have C&Rs in the first place.     

3D near-surface soil response from H/V ambient-noise ratios

The applicability of the horizontal-to-vertical (H/V) ambient-noise spectral ratio for characterizing earthquake site effects caused by near-surface topography and velocity structures was evaluated at sites underlain by thick (i.e. >100 m) sediment deposits near the southern-end of the New Madrid seismic zone in the central United States. Three-component ambient-noise and velocity models derived from seismic (shear-wave) refraction/reflection surveys showed that a relatively horizontal, sharp shear-wave velocity interface in the soil column resulted in an H/V spectral ratio with a single well-defined peak. Observations at sites with more than one sharp shear-wave velocity contrast and horizontally arranged soil layers resulted in at least two well-defined H/V spectral ratio peaks. Furthermore, at sites where there were sharp shear-wave velocity contrasts in nonhorizontal, near-surface soil layers, the H/V spectra exhibited a broad-bandwidth, relatively low-amplitude signal instead of a single well-defined peak.