Training Images from Process-Imitating Methods

The lack of a suitable training image is one of the main limitations of the application of multiple-point statistics (MPS) for the characterization of heterogeneity in real case studies. Process-imitating facies modeling techniques can potentially provide training images. However, the parameterization of these process-imitating techniques is not straightforward. Moreover, reproducing the resulting heterogeneous patterns with standard MPS can be challenging. Here the statistical properties of the paleoclimatic data set are used to select the best parameter sets for the process-imitating methods. The data set is composed of 278 lithological logs drilled in the lower Namoi catchment, New South Wales, Australia. A good understanding of the hydrogeological connectivity of this aquifer is needed to tackle groundwater management issues. The spatial variability of the facies within the lithological logs and calculated models is measured using fractal dimension, transition probability, and vertical facies proportion. To accommodate the vertical proportions trend of the data set, four different training images are simulated. The grain size is simulated alongside the lithological codes and used as an auxiliary variable in the direct sampling implementation of MPS. In this way, one can obtain conditional MPS simulations that preserve the quality and the realism of the training images simulated with the process-imitating method. The main outcome of this study is the possibility of obtaining MPS simulations that respect the statistical properties observed in the real data set and honor the observed conditioning data, while preserving the complex heterogeneity generated by the process-imitating method. In addition, it is demonstrated that an equilibrium of good fit among all the statistical properties of the data set should be considered when selecting a suitable set of parameters for the process-imitating simulations.     

Petrological imaging of an active pluton beneath Cerro Uturuncu,Bolivia

Uturuncu is a dormant volcano in the Altiplano of SW Bolivia. A present day ~70 km diameter interferometric synthetic aperture radar (InSAR) anomaly roughly centred on Uturuncu’s edifice is believed to be a result of magma intrusion into an active crustal pluton. Past activity at the volcano, spanning 0.89 to 0.27 Ma, is exclusively effusive and almost all lavas and domes are dacitic with phenocrysts of plagioclase, orthopyroxene, biotite, ilmenite and Ti-magnetite plus or minus quartz, and microlites of plagioclase and orthopyroxene set in rhyolitic groundmass glass. Plagioclase-hosted melt inclusions (MI) are rhyolitic with major element compositions that are similar to groundmass glasses. H₂O concentrations plotted versus incompatible elements for individual samples describe a trend typical of near-isobaric, volatile-saturated crystallisation. At 870 °C, the average magma temperature calculated from Fe–Ti oxides, the average H₂O of 3.2 ± 0.7 wt% and CO₂ typically <160 ppm equate to MI trapping pressures of 50–120 MPa, approximately 2–4.5 km below surface. Such shallow storage precludes the role of dacite magma emplacement into pre-eruptive storage regions as being the cause of the observed InSAR anomaly. Storage pressures, whole-rock (WR) chemistry and phase assemblage are remarkably consistent across the eruptive history of the volcano, although magmatic temperatures calculated from Fe–Ti oxide geothermometry, zircon saturation thermometry using MI and orthopyroxene-melt thermometry range from 760 to 925 °C at NNO ± 1 log. This large temperature range is similar to that of saturation temperatures of observed phases in experimental data on Uturuncu dacites. The variation in calculated temperatures is attributed to piecemeal construction of the active pluton by successive inputs of new magma into a growing volume of plutonic mush. Fluctuating temperatures within the mush can account for sieve-textured cores and complex zoning in plagioclase phenocrysts, resorption of quartz and biotite phenocrysts and apatite microlites. That Fe–Ti oxide temperatures vary by ~50–100 °C in a single thin section indicates that magmas were not homogenised effectively prior to eruption. Phenocryst contents do not correlate with calculated magmatic temperatures, consistent with crystal entrainment from the mush during magma ascent and eruption. Microlites grew during ascent from the magma storage region. Variability in the proportion of microlites is attributed to differing ascent and effusion rates with faster rates in general for lavas >0.5 Ma compared to those <0.5 Ma. High microlite contents of domes indicate that effusion rates were probably slowest in dome-forming eruptions. Linear trends in WR major and trace element chemistries, highly variable, bimodal mineral compositions, and the presence of mafic enclaves in lavas demonstrate that intrusion of more mafic magmas into the evolving, shallow plutonic mush also occurred further amplifying local temperature fluctuations. Crystallisation and resorption of accessory phases, particularly ilmenite and apatite, can be detected in MI and groundmass glass trace element covariation trends, which are oblique to WRs. Marked variability of Ba, Sr and La in MI can be attributed to temperature-controlled, localised crystallisation of plagioclase, orthopyroxene and biotite within the evolving mush.     

What is a flux tube? On the magnetic field topology of buoyant flux structures

We study the topology of field lines threading buoyant magnetic flux structures. The magnetic structures, visually resembling idealized magnetic flux tubes, are generated self-consistently by numerical simulation of the interaction of magnetic buoyancy and a localized velocity shear in a stably stratified atmosphere. Depending on the parameters, the system exhibits varying degrees of symmetry. By integrating along magnetic field lines and constructing return maps, we show that, depending on the type of underlying behaviour, the stages of the evolution, and therefore the degree of symmetry, the resulting magnetic structures can have field lines with one of three distinct topologies. When the x -translational and y -reflectional symmetries remain intact, magnetic field lines lie on surfaces but individual lines do not cover the surface. When the y symmetry is broken, magnetic field lines lie on surfaces and individual lines do cover the surface. When both x and y symmetries are broken, magnetic field lines wander chaotically over a large volume of the magnetically active region. We discuss how these results impact our simple ideas of a magnetic flux tube as an object with an inside and an outside, and introduce the concept of 'leaky' tubes.     

Variable Grid Method: An Intuitive Approach for Simultaneously Quantifying and Visualizing Spatial Data and Uncertainty

Efforts to develop applications and methods that effectively quantify and communicate uncertainty associated with spatial data remains a focus within many scientific communities. However, the inherent complexity of uncertainty makes it difficult to define, characterize, and represent. Frequently, the products of spatial and spatio‐temporal data are presented without a clear explanation of the inherent uncertainty underlying the data. As uses and applications for spatial data and their products continues to increase, so does the importance for utilizing reliable approaches to effectively communicate spatial data along with their inherent uncertainties. To address this need, the Variable Grid Method (VGM) was developed as an intuitive approach that simultaneously communicates both spatial patterns and trends and the uncertainty associated with data or their analyses. This article details the VGM approach and demonstrates the utility of the VGM to provide critical information about the relationship between uncertainty and spatial data, necessary to support the increasing utilization of spatial information for a wide range of research and other needs.     

The High-Resolution Coronal Imager (Hi-C)

The High-Resolution Coronal Imager (Hi-C) was flown on a NASA sounding rocket on 11 July 2012. The goal of the Hi-C mission was to obtain high-resolution (≈?0.3?–?0.4′′), high-cadence (≈?5 seconds) images of a solar active region to investigate the dynamics of solar coronal structures at small spatial scales. The instrument consists of a normal-incidence telescope with the optics coated with multilayers to reflect a narrow wavelength range around 19.3 nm (including the Fe xii 19.5-nm spectral line) and a 4096×4096 camera with a plate scale of 0.1′′?pixel^?1. The target of the Hi-C rocket flight was Active Region 11520. Hi-C obtained 37 full-frame images and 86 partial-frame images during the rocket flight. Analysis of the Hi-C data indicates the corona is structured on scales smaller than currently resolved by existing satellite missions.     

H2 in molecular shocks: observations from 2 to 12 microns

We present spectroscopy of the molecular shock in the supernova remnant IC 443. The H₂ level populations in IC 443 are nearly identical to Orion Peak 1 (Brand et al. 1988; Parmar et al. 1994) and are consistent with a partially dissociative J-shock. The similarity with Orion suggests the shock processes which dominate the emission in both regions are essentially identical.