Visualization and ecohydrologic models: Opening the box

Earth system models synthesize the science of interactions amongst multiple biophysical and, increasingly, human processes across a wide range of scales. Ecohydrologic models are a subset of earth system models that focus particularly on the complex interactions between ecosystem processes and the storage and flux of water. Ecohydrologic models often focus at scales where direct observations occur: plots, hillslopes, streams, and watersheds, as well as where land and resource management decisions are implemented. These models complement field-based and data-driven science by combining theory, empirical relationships derived from observation and new data to create virtual laboratories. Ecohydrologic models are tools that managers can use to ask “what if” questions and domain scientists can use to explore the implications of new theory or measurements. Recent decades have seen substantial advances in ecohydrologic models, building on both new domain science and advances in software engineering and data availability. The increasing sophistication of ecohydrologic models however, presents a barrier to their widespread use and credibility. Their complexity, often encoding 100s of relationships, means that they are effectively “black boxes,” at least for most users, sometimes even to the teams of researchers that contribute to their design. This opacity complicates the interpretation of model results. For models to effectively advance our understanding of how plants and water interact, we must improve how we visualize not only model outputs, but also the underlying theories that are encoded within the models. In this paper, we outline a framework for increasing the usefulness of ecohydrologic models through better visualization. We outline four complementary approaches, ranging from simple best practices that leverage existing technologies, to ideas that would engage novel software engineering and cutting edge human–computer interface design. Our goal is to open the ecohydrologic model black box in ways that will engage multiple audiences, from novices to model developers, and support learning, new discovery, and environmental problem solving.     

Cosmogenic nuclides in core samples of the Chico L6 chondrite: Evidence for irradiation under high shielding

Abstract— Core samples were obtained from various locations of the ～ 105-kg Chico, NM, L6 chondrite in order to study the effects of large shielding on the production rates of cosmic-ray-produced nuclides. Relations between measured abundances of cosmogenic nuclides (¹⁰Be, ²⁶Al, and stable isotopes of He, Ne, and Ar) and the cosmogenic ²²Ne/²¹Ne ratio were determined and compared with recent model predictions of production rates. The measured ²²Ne/²¹Ne ratios (1.06-1.08) and significant variations observed in concentrations of cosmogenic ²¹Ne and ³He suggest an ～40-cm shielding gradient across Chico and irradiation within a large object (> 100-cm radius). Noble gas data indicate that Chico experienced greater shielding than chondrites Knyahinya or Keyes and similar to Jilin. Values of ¹⁰Be (average = 20.7 dpm/kg) and ²⁶Al (average = 71.1 dpm/kg) are nearly constant, however, and show no correlation with either ²²Ne/²¹Ne or ²¹Ne. Activities of ¹⁰Be and ²⁶Al suggest irradiation in a smaller object (～40–80 cm radius). The ²⁶Al activity and the ²⁶Al/¹⁰Be ratio (average value = 3.42) are both significantly larger than values for most other chondrites. These results could indicate a two-stage irradiation with t₁ ～ 104 Ma and t₂ ～ 4 Ma and a second-stage body the size of Knyahinya. The single stage, ¹⁰Be/²¹Ne exposure age for Chico is 65 Ma. The ²²Ne/²¹Ne ratio apparently becomes insensitive to shielding for objects the size of Chico. No substantial evidence exists for chondrites with ²²Ne/²¹Ne ratios significantly less than ～ 1.055.     

Quiescent radio emission of RS CVn binary systems

It is shown that the quiescent radio emission observed in some RS Vc binary systems cannot be ascribed to the same thermal electrons population responsible for the X-ray emission. An alternative explanation to account for the observed correlation between X-ray and radio luminosity is suggested.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

Beyond headline mitigation numbers: we need more transparent and comparable NDCs to achieve the Paris Agreement on climate change

Nationally determined contributions (NDCs) were key to reaching the Paris Agreement and will be instrumental in implementing it. Research was quick to identify the ‘headline numbers’ of NDCs: if these climate action plans were fully implemented, global mean warming by 2100 would be reduced from approximately 3.6 to 2.7°C above pre-industrial levels (Höhne et al. Climate Pol 17:1–17, 2016; Rogelj et al. Nature 534:631–639, 2016). However, beyond these headline mitigation numbers, NDCs are more difficult to analyse and compare. UN climate negotiations have so far provided limited guidance on NDC formulation, which has resulted in varying scopes and contents of NDCs, often lacking details concerning ambitions. If NDCs are to become the long-term instrument for international cooperation, negotiation, and ratcheting up of ambitions to address climate change, then they need to become more transparent and comparable, both with respect to mitigation goals, and to issues such as adaptation, finance, and the way in which NDCs are aligned with national policies. Our analysis of INDCs and NDCs (Once a party ratifies the Paris Agreement, it is invited to turn its Intended Nationally Determined Contribution (INDC) into an NDC. We refer to results from our INDC analysis rather than our NDC analysis in this commentary unless otherwise stated.) shows that they omit important mitigation sectors, do not adequately provide details on costs and financing of implementation, and are poorly designed to meet assessment and review needs.     

Two magma series and associated ore deposit types in the Permian Emeishan large igneous province, SW China

The Late Middle Permian ( 260 Ma) Emeishan large igneous province in SW China contains two magmatic series, one comprising high-Ti basalts and Fe-rich gabbroic and syenitic intrusions, the other low-Ti basalts and mafic–ultramafic intrusions. The Fe-rich gabbros are spatially and temporally associated with syenites. Each series is associated with a distinctive type of mineralization, the first with giant Fe–Ti–V oxide ore deposits such as Panzhihua and Baima, the second with Ni–Cu–(PGE) sulfide deposits such as Jinbaoshan, Limahe and Zhubu. New SHRIMP zircon U–Pb isotopic data yielded 263 ± 3 Ma for the Limahe intrusion, 261 ± 2 Ma for the Zhubu intrusion and 262 ± 2 Ma for a syenitic intrusion. These new age dates, together with previously reported SHRIMP zircon U–Pb ages, suggest that all these intrusions are contemporaneous with the Emeishan flood basalts and formed during a major igneous event at ca. 260 Ma.The oxide-bearing intrusions have higher Al₂O₃, FeO (as total iron) and total alkalis (Na₂O + K₂O) but lower MgO than the sulfide-bearing intrusions. All intrusions are variably enriched in LREE relative to HREE. The oxide-bearing intrusions display positive Nb- and Ti-anomalies and in certain cases negative Zr–Hf anomalies, whereas the sulfide-bearing intrusions have obvious negative Nb- and Ti-anomalies, a feature of crustal contamination. Individual intrusions have relatively small ranges of Nd(t) values. All the intrusions, however, have Nd(t) values ranging from − 3.9 to + 4.6, and initial ⁸⁷Sr/⁸⁶Sr ratios from 0.7039 to 0.7105. The syenites have very low MgO (< 2 wt.%) but highly variable Fe₂O₃ (2.5 to 13 wt.%) with initial ⁸⁷Sr/⁸⁶Sr ratios ranging from 0.7039 to 0.7089. Magmas from both series could have derived by melting of a heterogeneous mantle plume: the high-Ti series from a Fe-rich, more fertile source and the low-Ti series from a Fe-poor, more refractory source. In addition, the low-Ti series underwent significant crustal contamination. The two magma series evolved along different paths that led to distinct mineralization styles.     

Age and emplacement of late-Variscan granites of the western Bohemian Massif with main focus on the Hauzenberg granitoids (European Variscides, Germany)

The Variscan Hauzenberg pluton consists of granite and granodiorite that intruded late- to postkinematically into HT-metamorphic rocks of the Moldanubian unit at the southwestern margin of the Bohemian Massif (Passauer Wald). U–Pb dating of zircon single-grains and monazite fractions, separated from medium- to coarse-grained biotite-muscovite granite (Hauzenberg granite II), yielded concordant ages of 320 ± 3 and 329 ± 7 Ma, interpreted as emplacement age. Zircons extracted from the younger Hauzenberg granodiorite yielded a ²⁰⁷Pb–²⁰⁶Pb mean age of 318.6 ± 4.1 Ma. The Hauzenberg granite I has not been dated. The pressure during solidification of the Hauzenberg granite II was estimated at 4.6 ± 0.6 kbar using phengite barometry on magmatic muscovite, corresponding to an emplacement depth of 16-18 km. The new data are compatible with pre-existing cooling ages of biotite and muscovite which indicate the Hauzenberg pluton to have cooled below T = 250–400 °C in Upper Carboniferous times. A compilation of age data from magmatic and metamorphic rocks of the western margin of the Bohemian Massif suggests a west- to northwestward shift of magmatism and HT/LP metamorphism with time. Both processes started at > 325 Ma within the South Bohemian Pluton and magmatism ceased at ca. 310 Ma in the Bavarian Oberpfalz. The slight different timing of HT metamorphism in northern Austria and the Bavarian Forest is interpreted as being the result of partial delamination of mantle lithosphere or removal of the thermal boundary layer.     

Platinum-group minerals from the Jinbaoshan Pd–Pt deposit,SW China: evidence for magmatic origin and hydrothermal alteration

The Jinbaoshan Pt–Pd deposit in Yunnan, SW China, is hosted in a wehrlite body, which is a member of the Permian (∼260 Ma) Emeishan Large Igneous Province (ELIP). The deposit is reported to contain one million tonnes of Pt–Pd ore grading 0.21% Ni and 0.16% Cu with 3.0 g/t (Pd + Pt). Platinum-group minerals (PGM) mostly are ∼10 μm in diameter, and are commonly Te-, Sn- and As-bearing, including moncheite (PtTe₂), atokite (Pd₃Sn), kotulskite (PdTe), sperrylite (PtAs₂), irarsite (IrAsS), cooperite (PtS), sudburyite (PdSb), and Pt–Fe alloy. Primary rock-forming minerals are olivine and clinopyroxene, with clinopyroxene forming anhedral poikilitic crystals surrounding olivine. Primary chromite occurs either as euhedral grains enclosed within olivine or as an interstitial phase to the olivine. However, the intrusion has undergone extensive hydrothermal alteration. Most olivine grains have been altered to serpentine, and interstitial clinopyroxene is often altered to actinolite/tremolite and locally biotite. Interstitial chromite grains are either partially or totally replaced by secondary magnetite. Base-metal sulfides (BMS), such as pentlandite and chalcopyrite, are usually interstitial to the altered olivine. PGM are located with the BMS and are therefore also interstitial to the serpentinized olivine grains, occurring within altered interstitial clinopyroxene and chromite, or along the edges of these minerals, which predominantly altered to actinolite/tremolite, serpentine and magnetite. Hydrothermal fluids were responsible for the release of the platinum-group elements (PGE) from the BMS to precipitate the PGM at low temperature during pervasive alteration. A sequence of alteration of the PGM has been recognized. Initially moncheite and atokite have been corroded and recrystallized during the formation of actinolite/tremolite, and then, cooperite and moncheite were altered to Pt–Fe alloy where they are in contact with serpentine. Sudburyite occurs in veins indicating late Pd mobility. However, textural evidence shows that the PGM are still in close proximity to the BMS. They occur in PGE-rich layers located at specific igneous horizons in the intrusion, suggesting that PGE were originally magmatic concentrations that, within a PGE-rich horizon, crystallized with BMS late in the olivine/clinopyroxene crystallization sequence and have not been significantly transported during serpentinization and alteration.     

Determination of the ionospheric foF2 using a stand-alone GPS receiver

The critical frequency of ionospheric F2 layer (foF2) is a measure of the highest frequency of radio signal that may be reflected back by the F2 layer, and it is associated with ionospheric peak electron density in the F2 layer. Accurate long-term foF2 variations are usually derived from ionosonde observations. In this paper, we propose a new method to observe foF2 using a stand-alone global positioning system (GPS) receiver. The proposed method relies on the mathematical equation that relates foF2 to GPS observations. The equation is then implemented in the Kalman filter algorithm to estimate foF2 at every epoch of the observation (30-s rate). Unlike existing methods, the proposed method does not require any additional information from ionosonde observations and does not require any network of GPS receivers. It only requires as inputs the ionospheric scale height and the modeled plasmaspheric electron content, which practically can be derived from any existing ionospheric/plasmaspheric model. We applied the proposed method to estimate long-term variations of foF2 at three GPS stations located at the northern hemisphere (NICO, Cyprus), the southern hemisphere (STR1, Australia) and the south pole (SYOG, Antarctic). To assess the performance of the proposed method, we then compared the results against those derived by ionosonde observations and the International Reference Ionosphere (IRI) 2012 model. We found that, during the period of high solar activity (2011–2012), the values of absolute mean bias between foF2 derived by the proposed method and ionosonde observations are in the range of 0.2–0.5 MHz, while those during the period of low solar activity (2009–2010) are in the range of 0.05–0.15 MHz. Furthermore, the root-mean-square-error (RMSE) values during high and low solar activities are in the range of 0.8–0.9 MHz and of 0.6–0.7 MHz, respectively. We also noticed that the values of absolute mean bias and RMSE between foF2 derived by the proposed method and the IRI-2012 model are slightly larger than those between the proposed method and ionosonde observations. These results demonstrate that the proposed method can estimate foF2 with a comparable accuracy. Since the proposed method can estimate foF2 at every epoch of the observation, it therefore has promising applications for investigating various scales (from small to large) of foF2 irregularities.