Earth Observation Scientific Workflows in a Distributed Computing Environment

Geospatially Enabled Scientific Workflows offer a promising toolset to help researchers in the earth observation domain with many aspects of the scientific process. One such aspect is that of access to distributed earth observation data and computing resources. Earth observation research often utilizes large datasets requiring extensive CPU and memory resources in their processing. These resource intensive processes can be chained; the sequence of processes (and their provenance) makes up a scientific workflow. Despite the exponential growth in capacity of desktop computers, their resources are often insufficient for the scientific workflow processing tasks at hand. By integrating distributed computing capabilities into a geospatially enabled scientific workflow environment, it is possible to provide researchers with a mechanism to overcome the limitations of the desktop computer. Most of the effort on extending scientific workflows with distributed computing capabilities has focused on the web services approach, as exemplified by the OGC's Web Processing Service and by GRID computing. The approach to leveraging distributed computing resources described in this article uses instead remote objects via RPyC and the dynamic properties of the Python programming language. The Vistrails environment has been extended to allow for geospatial processing through the EO4Vistrails package ( http://code.google.com/p/eo4vistrails/ ). In order to allow these geospatial processes to be seamlessly executed on distributed resources such as cloud computing nodes, the Vistrails environment has been extended with both multi‐tasking capabilities and distributed processing capabilities. The multi‐tasking capabilities are required in order to allow Vistrails to run side‐by‐side processes, a capability it does not currently have. The distributed processing capabilities are achieved through the use of remote objects and mobile code through RPyC.     

The influence of negative emission technologies and technology policies on the optimal climate mitigation portfolio

Combining policies to remove carbon dioxide (CO₂) from the atmosphere with policies to reduce emissions could decrease CO₂ concentrations faster than possible via natural processes. We model the optimal selection of a dynamic portfolio of abatement, research and development (R&D), and negative emission policies under an exogenous CO₂ constraint and with stochastic technological change. We find that near-term abatement is not sensitive to the availability of R&D policies, but the anticipated availability of negative emission strategies can reduce the near-term abatement optimally undertaken to meet 2°C temperature limits. Further, planning to deploy negative emission technologies shifts optimal R&D funding from ??carbon-free?? technologies into ??emission intensity?? technologies. Making negative emission strategies available enables an 80% reduction in the cost of keeping year 2100 CO₂ concentrations near their current level. However, negative emission strategies are less important if the possibility of tipping points rules out using late-century net negative emissions to temporarily overshoot the CO₂ constraint earlier in the century.     

Preparing for climate change in Washington State

Climate change is expected to bring potentially significant changes to Washington State’s natural, institutional, cultural, and economic landscape. Addressing climate change impacts will require a sustained commitment to integrating climate information into the day-to-day governance and management of infrastructure, programs, and services that may be affected by climate change. This paper discusses fundamental concepts for planning for climate change and identifies options for adapting to the climate impacts evaluated in the Washington Climate Change Impacts Assessment. Additionally, the paper highlights potential avenues for increasing flexibility in the policies and regulations used to govern human and natural systems in Washington.     

Adaptive co-management and the paradox of learning

Much emphasis has been placed on the importance of learning to support collaborative environmental management and achieve sustainability under conditions of social–ecological change. Yet, on-going struggles to learn from experience and respond to complex social–ecological conditions reflect an emerging paradox. Despite widespread support of learning as a normative goal and process, core concepts, assumptions and approaches to learning have been applied in vague and sometimes uncritical ways. Greater specificity with respect to learning goals, approaches and outcomes is required. In response to this gap, we examine five dimensions of the learning paradox in the context of adaptive co-management, where the learning and linking functions of governance are stressed: (i) definitions of learning; (ii) learning goals and expectations; (iii) mechanisms by which learning takes place; (iv) questions regarding who is involved in the process of learning; and (v) the risks and ethical ambiguities faced by different actors expected to willingly participate in a learning process, whether formal or informal. Lessons from experience with a series of cases from the global North and South illustrate the implications of these dimensions. Resolving the dimensions of this learning paradox will require greater attention to capacity-building, recognition of the role of risk, and consideration of how incentives could be used to encourage learning. Further consideration of the role of power and marginality among groups participating in the learning process is also needed, as is more systematic evaluation to monitor and measure learning outcomes.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

Mass-independent fractionation of oxygen isotopes in the mesostasis of a chondrule from the Semarkona LL3.0 ordinary chondrite

A large chondrule from Semarkona, the most primitive ordinary chondrite known, has been discovered to contain a record of mass transport during its formation. In most respects, it is a normal Type I, group A1, low-FeO chondrule that was produced by reduction and mass-loss during the unidentified flash-heating event that produced the chondrules, the most abundant structural component in primitive meteorites. We have previously measured elemental abundances and abundance profiles in this chondrule. We here report oxygen isotope ratio abundances and ratio abundance profiles. We have found that the mesostasis is zoned in oxygen isotope ratio, with the center of the chondrule containing isotopically heavier oxygen than the outer regions, the outer regions being volatile rich from the diffusion of volatiles into the chondrule during cooling. The δ¹⁷O values range from −2.0‰ to 9.9‰, while δ¹⁸O range from −1.9‰ to 9.6‰. More importantly, a plot of δ¹⁷O against δ¹⁸O has a slope of 1.1 ± 0.2 (1σ) and 0.88 ± 0.10 (1σ) when measured by two independent methods. Co-variation of δ¹⁷O with δ¹⁸O that does not follow mass-dependent fractionation has often been seen in primitive solar system materials and is usually ascribed to the mixing of different oxygen reservoirs. We argue that petrographic and compositional data indicate that this chondrule was completely melted at the time of its formation so that relic grains could not have survived. Furthermore, there is petrographic and compositional evidence that there was no aqueous alteration of this chondrule subsequent to its formation. Although it is possible to formulate a series of exchanges between the chondrule and external ¹⁶O-rich and ¹⁶O-poor reservoirs that may explain the detailed oxygen isotope systematics of this chondrule, such a sequence of events looks very contrived. We therefore hypothesize that reduction, devolatilization, and crystallization of the chondrule melt may have produced ¹⁶O-rich olivines and ¹⁶O-poor mesostasis plotting on a slope-one line as part of the chondrule-forming process in an analogous fashion to known chemical mass-independent isotopic fractionation mechanisms. During cooling, volatiles and oxygen near the terrestrial line in oxygen isotope composition produced the outer zone of volatile rich and ¹⁶O-rich mesostasis. The chondrule therefore not only retains a record of considerable mass transport accompanying formation, but also may indicate that the isotopes of oxygen underwent mass-independent fractionation during the process.     

Investigation of acidic dissolution of mixed clays between pH 1.0 and −3.0 using Si and Al X-ray absorption near edge structure

Although widely investigated in relation to acid mine drainage systems at pH > 1.0, we know little about the impact of sulfuric acid (H₂SO₄) on the geochemistry and mineralogy of clays at pH < 1.0 (including negative pH values). Thus, laboratory batch experiments were conducted on three mixed clay samples with different mass ratios of phyllosilicates (smectite, illite, and kaolinite) to investigate the impact of H₂SO₄ from pH 1.0 to −3.0 for exposure periods of 14, 90, 180, and 365 days. Si and Al K- and L_2,3-edge X-ray absorption near edge structure (XANES) spectroscopy were employed on these samples to determine the chemical and structural changes that occur during acidic dissolution of phyllosilicates that cannot be distinguished using X-ray diffraction analyses. A series of silicate, phyllosilicate, and Al-bearing standard compounds were also studied to provide an explanation for the observed changes in the clay samples. The Si XANES results indicated the preferential dissolution of the phyllosilicates (pH ? 1.0, t ? 14 d), the persistence of quartz even at pH ? −3.0 and t ? 365 d, and the formation of an amorphous silica-like phase that was confined to the surface layer of the altered clay samples at pH ? 0.0 and t ? 90 d). Al XANES results demonstrated dissolution of Al-octahedral layers (pH ? 1.0, t ? 14 d), the persistence of four-fold relative to six-fold coordinated Al, and the precipitation of an Al-SO₄-rich phase (pH ? −1.0, t ? 90 d). An existing conceptual model of phyllosilicate dissolution under extremely acidic conditions was modified to include the results of this study.     

A possible glacially-forced tectonic mechanism for late Neogene surface uplift and subsidence around the North Atlantic

Analysis of old erosion surfaces and estimates of exhumation from apatite fission track data can be used to infer late Neogene surface uplift of Britain, Greenland, Norway and Svalbard of 1–2 km. Subsidence and sedimentation in adjacent offshore basins can be found from interpretation of seismic and well log data. Various mechanisms for surface uplift have been proposed but the underlying cause remains unexplained. Since the multiple glaciations that took place during the late Neogene were a common factor, a possible glacially-forced tectonic mechanism to thicken the crust and produce surface uplift has been investigated. This could result from the relatively slow accumulation of ice that loads the crust as an ice sheet grows during a glacial period, followed by relatively rapid retreat and unloading around its periphery at the end. Unloading could create transient stresses that induce lateral flow in a ductile lower crust to thicken it onshore and produce surface uplift, with associated thinning beneath adjacent offshore basins, producing subsidence. Simple calculations show that the proposed mechanism is feasible and indicate that crustal thickening and surface uplift accumulated from a number of glacial cycles can account for the observed surface uplift, with an acceptable flow rate in the lower crust at the end of each cycle if the viscosity of ductile flow is sufficiently low.     

Rapid incorporation of lipids into macromolecules during experimental decay of invertebrates: Initiation of geopolymer formation

Diagenetic alteration is critical for the preservation of fossil cuticles of plant and animal origin and to the formation of kerogen. The process takes place over millions of years, but the stage at which it is initiated is not known. Laboratory decay experiments were carried out on shrimps, scorpions and cockroaches to monitor changes in the chitin–protein of the arthropod cuticle and associated lipids. The cockroach and scorpion exoskeleton remained largely unaltered morphologically, but the shrimp experienced rapid decomposition within a month, which progressed through the 44 week duration of the experiment as revealed using electron microscopy. Mass spectrometry and ¹³C NMR (nuclear magnetic resonance) spectroscopy revealed the association of an n-alkyl component with labile lipids, such as fatty acids with up to 24 carbon atoms, which were incorporated into the decaying macromolecule. The scorpion and cockroach cuticle did not reveal the incorporation of additional lipids, indicating that decay is important in initiating in situ lipid association. This experiment provides evidence that lipids can become associated with carbohydrate and proteinaceous macromolecules during the very early stages of decay, representing the first stage in the transformation process that contributes to the aliphatic rich composition ubiquitous in organic fossils and kerogen.     

N-body simulations of cohesion in dense planetary rings: A study of cohesion parameters

We present results from a large suite of simulations of Saturn’s dense A and B rings using a new model of particle sticking in local simulations (Perrine, R.P., Richardson, D.C., Scheeres, D.J. [2011]. Icarus 212, 719–735). In this model, colliding particles can be incorporated into or help fragment rigid aggregations on the basis of certain user-specified parameters that can represent van der Waals forces or interlocking surface frost layers.Our investigation is motivated by laboratory results that show that interpenetration of surface layers can allow impacting frost-covered ice spheres to stick together. In these experiments, cohesion only occurs below specific impact speeds, which happen to be characteristic of impact speeds in Saturn’s rings. Our goal is to determine if weak bonding is consistent with ring observations, to constrain cohesion parameters in light of existing ring observations, to make predictions about particle populations throughout the rings, and to discover other diagnostics that may constrain bonding parameters.We considered the effects of five parameters on the equilibrium characteristics of our ring simulations: speed-based merge and fragmentation limits, bond strength, ring surface density, and patch orbital distance (i.e., the A or B ring), some with both monodisperse and polydisperse comparison cases. In total, we present data from 95 simulations.We find that weak cohesion is consistent with observations of the A and B rings (e.g., French, R.G., Nicholson, P.D. [2000]. Icarus 145, 502–523), and we present a range of simulation parameters that reproduce the observed size distribution and maximum particle size. It turns out that the parameters that match observations differ between the A and B rings, and we discuss the potential implications of this result. We also comment on other observable consequences of cohesion for the rings, such as optical depth and scale height effects, and discuss whether very large objects (e.g., “propeller” source objects) are grown bottom-up from cohesion of smaller ring particles.