Fusion of infrared polarization and intensity images using support value transform and fuzzy combination rules

Infrared polarization and intensity imagery provide complementary and discriminative information in image understanding and interpretation. In this paper, a novel fusion method is proposed by effectively merging the information with various combination rules. It makes use of both low-frequency and highfrequency images components from support value transform (SVT), and applies fuzzy logic in the combination process. Images (both infrared polarization and intensity images) to be fused are firstly decomposed into low-frequency component images and support value image sequences by the SVT. Then the low-frequency component images are combined using a fuzzy combination rule blending three sub-combination methods of (1) region feature maximum, (2) region feature weighting average, and (3) pixel value maximum; and the support value image sequences are merged using a fuzzy combination rule fusing two sub-combination methods of (1) pixel energy maximum and (2) region feature weighting. With the variables of two newly defined features, i.e. the low-frequency difference feature for low-frequency component images and the support-value difference feature for support value image sequences, trapezoidal membership functions are proposed and developed in tuning the fuzzy fusion process. Finally the fused image is obtained by inverse SVT operations. Experimental results of visual inspection and quantitative evaluation both indicate the superiority of the proposed method to its counterparts in image fusion of infrared polarization and intensity images.

Jet characterization in the upper troposphere/lower stratosphere (UTLS): applications to climatology and transport studies

A method of classifying the upper tropospheric/lower stratospheric (UTLS) jets has been developed that allows satellite and aircraft trace gas data and meteorological fields to be efficiently mapped in a jet coordinate view. A detailed characterization of multiple tropopauses accompanies the jet characterization. Jet climatologies show the well-known high altitude subtropical and lower altitude polar jets in the upper troposphere, as well as a pattern of concentric polar and subtropical jets in the Southern Hemisphere, and shifts of the primary jet to high latitudes associated with blocking ridges in Northern Hemisphere winter. The jet-coordinate view segregates air masses differently than the commonly-used equivalent latitude (EqL) coordinate throughout the lowermost stratosphere and in the upper troposphere. Mapping O3 data from the Aura Microwave Limb Sounder (MLS) satellite and the Winter Storms aircraft datasets in jet coordinates thus emphasizes different aspects of the circulation compared to an EqL-coordinate framework: the jet coordinate reorders the data geometrically, thus highlighting the strong PV, tropopause height and trace gas gradients across the subtropical jet, whereas EqL is a dynamical coordinate that may blur these spatial relationships but provides information on irreversible transport. The jet coordinate view identifies the concentration of stratospheric ozone well below the tropopause in the region poleward of and below the jet core, as well as other transport features associated with the upper tropospheric jets. Using the jet information in EqL coordinates allows us to study trace gas distributions in regions of weak versus strong jets, and demonstrates weaker transport barriers in regions with less jet influence. MLS and Atmospheric Chemistry Experiment-Fourier Transform Spectrometer trace gas fields for spring 2008 in jet coordinates show very strong, closely correlated, PV, tropopause height and trace gas gradients across the jet, and evidence of intrusions of stratospheric air below the tropopause below and poleward of the subtropical jet; these features are consistent between instruments and among multiple trace gases. Our characterization of the jets is facilitating studies that will improve our understanding of upper tropospheric trace gas evolution.

A Gaussian-mixture ensemble transform filter

We generalize the popular ensemble Kalman filter to an ensemble transform filter, in which the prior distribution can take the form of a Gaussian mixture or a Gaussian kernel density estimator. The design of the filter is based on a continuous formulation of the Bayesian filter analysis step. We call the new filter algorithm the ensemble Gaussian-mixture filter (EGMF). The EGMF is implemented for three simple test problems (Brownian dynamics in one dimension, Langevin dynamics in two dimensions and the three-dimensional Lorenz-63 model). It is demonstrated that the EGMF is capable of tracking systems with non-Gaussian uni- and multimodal ensemble distributions. Copyright © 2011 Royal Meteorological Society

A nonparametric ensemble transform method for Bayesian inference

Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step, while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a priori assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables. © 2013, Society for Industrial and Applied Mathematics

Using paint to investigate fires: an ATR-IR study of the degradation of paint samples upon heating

Fire investigation is a challenging area for the forensic investigator. The aim of this work was to use spectral changes to paint samples to estimate the temperatures to which a paint has been heated. Five paint samples (one clay paint, two car paints, one metallic paint, and one matt emulsion) have been fully characterized by a combination of attenuated total reflectance Fourier transform infrared (ATR-IR), Raman, X-ray fluorescence spectroscopy and powder X-ray diffraction. The thermal decomposition of these paints has been investigated by means of ATR-IR and thermal gravimetric analysis. Clear temperature markers are observed in the ATR-IR spectra namely: loss of m(C = O) band, >300°C; appearance of water bands on cooling, >500°C; alterations to m(Si–O) bands due to dehydration of silicate clays, >700°C; diminution of m(CO3) and d(CO3) modes of CaCO3, >950°C. We suggest the possible use of portable ATR-IR for nondestructive, in situ analysis of paints.

Ensemble transform Kalman-Bucy filters

Two recent works have adapted the Kalman–Bucy filter into an ensemble setting. In the first formulation, the ensemble of perturbations is updated by the solution of an ordinary differential equation (ODE) in pseudo-time, while the mean is updated as in the standard Kalman filter. In the second formulation, the full ensemble is updated in the analysis step as the solution of single set of ODEs in pseudo-time. Neither requires matrix inversions except for the frequently diagonal observation error covariance. We analyse the behaviour of the ODEs involved in these formulations. We demonstrate that they stiffen for large magnitudes of the ratio of background error to observational error variance, and that using the integration scheme proposed in both formulations can lead to failure. A numerical integration scheme that is both stable and is not computationally expensive is proposed. We develop transform-based alternatives for these Bucy-type approaches so that the integrations are computed in ensemble space where the variables are weights (of dimension equal to the ensemble size) rather than model variables. Finally, the performance of our ensemble transform Kalman–Bucy implementations is evaluated using three models: the 3-variable Lorenz 1963 model, the 40-variable Lorenz 1996 model, and a medium complexity atmospheric general circulation model known as SPEEDY. The results from all three models are encouraging and warrant further exploration of these assimilation techniques.

The 22-Year Hale Cycle in cosmic ray flux: evidence for direct heliospheric modulation

The ability to predict times of greater galactic cosmic ray (GCR) fluxes is important for reducing the hazards caused by these particles to satellite communications, aviation, or astronauts. The 11-year solar-cycle variation in cosmic rays is highly correlated with the strength of the heliospheric magnetic field. Differences in GCR flux during alternate solar cycles yield a 22-year cycle, known as the Hale Cycle, which is thought to be due to different particle drift patterns when the northern solar pole has predominantly positive (denoted as qA>0 cycle) or negative (qA<0) polarities. This results in the onset of the peak cosmic-ray flux at Earth occurring earlier during qA>0 cycles than for qA<0 cycles, which in turn causes the peak to be more dome-shaped for qA>0 and more sharply peaked for qA<0. In this study, we demonstrate that properties of the large-scale heliospheric magnetic field are different during the declining phase of the qA<0 and qA>0 solar cycles, when the difference in GCR flux is most apparent. This suggests that particle drifts may not be the sole mechanism responsible for the Hale Cycle in GCR flux at Earth. However, we also demonstrate that these polarity-dependent heliospheric differences are evident during the space-age but are much less clear in earlier data: using geomagnetic reconstructions, we show that for the period of 1905 - 1965, alternate polarities do not give as significant a difference during the declining phase of the solar cycle. Thus we suggest that the 22-year cycle in cosmic-ray flux is at least partly the result of direct modulation by the heliospheric magnetic field and that this effect may be primarily limited to the grand solar maximum of the space-age.

X‑ray crystal structure of rac-[Ru(phen)2dppz]2+ with d(ATGCAT)2 shows enantiomer orientations and water ordering

We report an atomic resolution X-ray crystal structure containing both enantiomers of rac-[Ru(phen)2dppz]2+ with the d-(ATGCAT)2 DNA duplex (phen = phenanthroline; dppz = dipyridophenazine). The first example of any enantiomeric pair crystallized with a DNA duplex shows different orientations of the Λ and Δ binding sites, separated by a clearly defined structured water monolayer. Job plots show that the same species is present in solution. Each enantiomer is bound at a TG/CA step and shows intercalation from the minor groove. One water molecule is directly located on one phenazine N atom in the Δ-enantiomer only.

A synchrotron powder X-ray diffraction study of the skutterudite-related phases AB1.5Te1.5 (A=Co, Rh, Ir; B= Ge, Sn)

X-ray resonant scattering has been exploited to investigate the crystal structure of the AB1.5Te1.5 phases (A = Co, Rh, Ir; B = Ge, Sn). Analysis of the diffraction data reveals that CoGe1.5Te1.5 and ASn1.5Te1.5 adopt a rhombohedral skutterudite-related structure, containing diamond-shape B2Te2 rings, in which the B and Te atoms are ordered and trans to each other. Anion ordering is however incomplete, and with increasing the size of both cations and anions, the degree of anion ordering decreases. By contrast, the diffraction data of IrGe1.5Te1.5 are consistent with an almost statistical distribution of the anions over the available sites, although some ordered domains may be present. The thermoelectric properties of these materials are discussed in the light of these results.

Multiwavelength Observations of EXO 0748-676. I. Reprocessing of X-Ray Bursts

We present high time-resolution multiwavelength observations of X-ray bursts in the low-mass X-ray binary UY Vol. Strong reprocessed signals are present in the ultraviolet and optical, lagged and smeared with respect to the X-rays. The addition of far-ultraviolet coverage for one burst allows much tighter constraints on the temperature and geometry of the reprocessing region than previously possible. A blackbody reprocessing model for this burst suggests a rise in temperatures during the burst from 18,000 to 35,000 K and an emitting area comparable to that expected for the disk and/or irradiated companion star. The lags are consistent with those expected. The single-zone blackbody model cannot reproduce the ratio of optical to ultraviolet flux during the burst, however. The discrepancy seems too large to explain with deviations from a local blackbody spectrum and more likely indicates that a range of reprocessing temperatures are required. Comparable results are derived from other bursts, and in particular the lag and smearing both appear shorter when the companion star is on the near side of the disk as predicted. The burst observed by HST also yielded a spectrum of the reprocessed light. It is dominated by continuum, with a spectral shape consistent with the temperatures derived from lightcurve modeling. Taken as a whole, our observations confirm the standard paradigm of prompt reprocessing distributed across the disk and companion star, with the response dominated by a thermalized continuum rather than by emission lines.

Control and analysis of oriented thin films of lipid inverse bicontinuous cubic phases using grazing incidence small-angle X‑ray scattering

Lipid cubic phases are complex nanostructures that form naturally in a variety of biological systems, with applications including drug delivery and nanotemplating. Most X-ray scattering studies on lipid cubic phases have used unoriented polydomain samples as either bulk gels or suspensions of micrometer-sized cubosomes. We present a method of investigating cubic phases in a new form, as supported thin films that can be analyzed using grazing incidence small-angle X-ray scattering (GISAXS). We present GISAXS data on three lipid systems: phytantriol and two grades of monoolein (research and industrial). The use of thin films brings a number of advantages. First, the samples exhibit a high degree of uniaxial orientation about the substrate normal. Second, the new morphology allows precise control of the substrate mesophase geometry and lattice parameter using a controlled temperature and humidity environment, and we demonstrate the controllable formation of oriented diamond and gyroid inverse bicontinuous cubic along with lamellar phases. Finally, the thin film morphology allows the induction of reversible phase transitions between these mesophase structures by changes in humidity on subminute time scales, and we present timeresolved GISAXS data monitoring these transformations.

Galactic cosmic ray modulation near the heliospheric current sheet

Galactic cosmic rays (GCRs) are modulated by the heliospheric magnetic field (HMF) both over decadal time scales (due to long-term, global HMF variations), and over time scales of a few hours (associated with solar wind structures such as coronal mass ejections or the heliospheric current sheet, HCS). Due to the close association between the HCS, the streamer belt, and the band of slow solar wind, HCS crossings are often associated with corotating interaction regions where fast solar wind catches up and compresses slow solar wind ahead of it. However, not all HCS crossings are associated with strong compressions. In this study we categorize HCS crossings in two ways: Firstly, using the change in magnetic polarity, as either away-to-toward (AT) or toward-to-away (TA) magnetic field directions relative to the Sun and, secondly, using the strength of the associated solar wind compression, determined from the observed plasma density enhancement. For each category, we use superposed epoch analyses to show differences in both solar wind parameters and GCR flux inferred from neutron monitors. For strong-compression HCS crossings, we observe a peak in neutron counts preceding the HCS crossing, followed by a large drop after the crossing, attributable to the so-called ‘snow-plough’ effect. For weak-compression HCS crossings, where magnetic field polarity effects are more readily observable, we instead observe that the neutron counts have a tendency to peak in the away magnetic field sector. By splitting the data by the dominant polarity at each solar polar region, we find that the increase in GCR flux prior to the HCS crossing is primarily from strong compressions in cycles with negative north polar fields due to GCR drift effects. Finally, we report on unexpected differences in GCR behavior between TA weak compressions during opposing polarity cycles.

Estimating correlated observation error statistics using an ensemble transform Kalman filter

For certain observing types, such as those that are remotely sensed, the observation errors are correlated and these correlations are state- and time-dependent. In this work, we develop a method for diagnosing and incorporating spatially correlated and time-dependent observation error in an ensemble data assimilation system. The method combines an ensemble transform Kalman filter with a method that uses statistical averages of background and analysis innovations to provide an estimate of the observation error covariance matrix. To evaluate the performance of the method, we perform identical twin experiments using the Lorenz ’96 and Kuramoto-Sivashinsky models. Using our approach, a good approximation to the true observation error covariance can be recovered in cases where the initial estimate of the error covariance is incorrect. Spatial observation error covariances where the length scale of the true covariance changes slowly in time can also be captured. We find that using the estimated correlated observation error in the assimilation improves the analysis.

Observing the in situ chiral modification of Ni nanoparticles using scanning transmission X-ray microspectroscopy

Enantioselective heterogeneous hydrogenation of Cdouble bond; length as m-dashO bonds is of great potential importance in the synthesis of chirally pure products for the pharmaceutical and fine chemical industries. One of the most widely studied examples of such a reaction is the hydrogenation of β-ketoesters and β-diketoesters over Ni-based catalysts in the presence of a chiral modifier. Here we use scanning transmission X-ray microscopy combined with near-edge X-ray absorption fine structure spectroscopy (STXM/NEXAFS) to investigate the adsorption of the chiral modifier, namely (R,R)-tartaric acid, onto individual nickel nanoparticles. The C K-edge spectra strongly suggest that tartaric acid deposited onto the nanoparticle surfaces from aqueous solutions undergoes a keto-enol tautomerisation. Furthermore, we are able to interrogate the Ni L2,3-edge resonances of individual metal nanoparticles which, combined with X-ray diffraction (XRD) patterns showed them to consist of a pure nickel phase rather than the more thermodynamically stable bulk nickel oxide. Importantly, there appears to be no “particle size effect” on the adsorption mode of the tartaric acid in the particle size range ~ 90–~ 300 nm.