PnP-AdaNet: Plug-and-Play Adversarial Domain Adaptation Network with a Benchmark at Cross-modality Cardiac Segmentation

Deep convolutional networks have demonstrated the state-of-the-art performance on various medical image computing tasks. Leveraging images from different modalities for the same analysis task holds clinical benefits. However, the generalization capability of deep models on test data with different distributions remain as a major challenge. In this paper, we propose the PnPAdaNet (plug-and-play adversarial domain adaptation network) for adapting segmentation networks between different modalities of medical images, e.g., MRI and CT. We propose to tackle the significant domain shift by aligning the feature spaces of source and target domains in an unsupervised manner. Specifically, a domain adaptation module flexibly replaces the early encoder layers of the source network, and the higher layers are shared between domains. With adversarial learning, we build two discriminators whose inputs are respectively multi-level features and predicted segmentation masks. We have validated our domain adaptation method on cardiac structure segmentation in unpaired MRI and CT. The experimental results with comprehensive ablation studies demonstrate the excellent efficacy of our proposed PnP-AdaNet. Moreover, we introduce a novel benchmark on the cardiac dataset for the task of unsupervised cross-modality domain adaptation. We will make our code and database publicly available, aiming to promote future studies on this challenging yet important research topic in medical imaging

Progress on stochastic background search codes for LIGO

One of the types of signals for which the LIGO interferometric gravitational wave detectors will search is a stochastic background of gravitational radiation. We review the technique of searching for a background using the optimally-filtered cross-correlation statistic, and describe the state of plans to perform such cross-correlations between the two LIGO interferometers as well as between LIGO and other gravitational-wave detectors, in particular the preparation of software to perform such data analysis.Comment: 7 pages, 1 encapsulated PostScript figure, uses IOP class files, submitted to the proceedings of the 4th Amaldi meeting (which will be published in Classical and Quantum Gravity

A spectral survey of an ultra-hot Jupiter: Detection of metals in the transmission spectrum of KELT-9 b

Context: KELT-9 b exemplifies a newly emerging class of short-period gaseous exoplanets that tend to orbit hot, early type stars - termed ultra-hot Jupiters. The severe stellar irradiation heats their atmospheres to temperatures of $\sim 4,000$ K, similar to the photospheres of dwarf stars. Due to the absence of aerosols and complex molecular chemistry at such temperatures, these planets offer the potential of detailed chemical characterisation through transit and day-side spectroscopy. Studies of their chemical inventories may provide crucial constraints on their formation process and evolution history. Aims: To search the optical transmission spectrum of KELT-9 b for absorption lines by metals using the cross-correlation technique. Methods: We analyse 2 transits observed with the HARPS-N spectrograph. We use an isothermal equilibrium chemistry model to predict the transmission spectrum for each of the neutral and singly-ionized atoms with atomic numbers between 3 and 78. Of these, we identify the elements that are expected to have spectral lines in the visible wavelength range and use those as cross-correlation templates. Results: We detect absorption of Na I, Cr II, Sc II and Y II, and confirm previous detections of Mg I, Fe I, Fe II and Ti II. In addition, we find evidence of Ca I, Cr I, Co I, and Sr II that will require further observations to verify. The detected absorption lines are significantly deeper than model predictions, suggesting that material is transported to higher altitudes where the density is enhanced compared to a hydrostatic profile. There appears to be no significant blue-shift of the absorption spectrum due to a net day-to-night side wind. In particular, the strong Fe II feature is shifted by $0.18 \pm 0.27$ km~s$^{-1}$, consistent with zero. Using the orbital velocity of the planet we revise the steller and planetary masses and radii.Comment: Submitted to Astronomy and Astrophysics on January 18, 2019. Accepted on May 3, 2019. 26 pages, 11 figure

Three-dimensional distribution of ejecta in Supernova 1987A at 10 000 days

Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] \lambda \lambda 7292, 7324, [O I] \lambda \lambda 6300, 6364 and Mg II \lambda \lambda 9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 \mu m and He I 2.058 \mu m. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 \mu m, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.Comment: Accepted for publication in Ap

The unstable CO2 feedback cycle on ocean planets

Ocean planets are volatile-rich planets, not present in our Solar system, which are thought to be dominated by deep, global oceans. This results in the formation of high-pressure water ice, separating the planetary crust from the liquid ocean and, thus, also from the atmosphere. Therefore, instead of a carbonate-silicate cycle like on the Earth, the atmospheric carbon dioxide concentration is governed by the capability of the ocean to dissolve carbon dioxide (CO2). In our study, we focus on the CO2 cycle between the atmosphere and the ocean which determines the atmospheric CO2 content. The atmospheric amount of CO2 is a fundamental quantity for assessing the potential habitability of the planet's surface because of its strong greenhouse effect, which determines the planetary surface temperature to a large degree. In contrast to the stabilizing carbonate-silicate cycle regulating the long-term CO2 inventory of the Earth atmosphere, we find that the CO2 cycle feedback on ocean planets is negative and has strong destabilizing effects on the planetary climate. By using a chemistry model for oceanic CO2 dissolution and an atmospheric model for exoplanets, we show that the CO2 feedback cycle can severely limit the extension of the habitable zone for ocean planet

Bayesian detection of unmodeled bursts of gravitational waves

The data analysis problem of coherently searching for unmodeled gravitational-wave bursts in the data generated by a global network of gravitational-wave observatories has been at the center of research for almost two decades. As data from these detectors is starting to be analyzed, a renewed interest in this problem has been sparked. A Bayesian approach to the problem of coherently searching for gravitational wave bursts with a network of ground-based interferometers is here presented. We demonstrate how to systematically incorporate prior information on the burst signal and its source into the analysis. This information may range from the very minimal, such as best-guess durations, bandwidths, or polarization content, to complete prior knowledge of the signal waveforms and the distribution of sources through spacetime. We show that this comprehensive Bayesian formulation contains several previously proposed detection statistics as special limiting cases, and demonstrate that it outperforms them.Comment: 18 pages, 3 figures, revisions based on referee comment

Stochastic Background Search Correlating ALLEGRO with LIGO Engineering Data

We describe the role of correlation measurements between the LIGO interferometer in Livingston, LA, and the ALLEGRO resonant bar detector in Baton Rouge, LA, in searches for a stochastic background of gravitational waves. Such measurements provide a valuable complement to correlations between interferometers at the two LIGO sites, since they are sensitive in a different, higher, frequency band. Additionally, the variable orientation of the ALLEGRO detector provides a means to distinguish gravitational wave correlations from correlated environmental noise. We describe the analysis underway to set a limit on the strength of a stochastic background at frequencies near 900 Hz using ALLEGRO data and data from LIGO's E7 Engineering Run.Comment: 8 pages, 2 encapsulated PostScript figures, uses IOP class files, submitted to the proceedings of the 7th Gravitational Wave Data Analysis Workshop (which will be published in Classical and Quantum Gravity

Three-dimensional distribution of ejecta in Supernova 1987A at 10 000 days

Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] \lambda \lambda 7292, 7324, [O I] \lambda \lambda 6300, 6364 and Mg II \lambda \lambda 9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 \mu m and He I 2.058 \mu m. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 \mu m, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.Comment: Accepted for publication in Ap

Parameter Changes from Upscaling of a Local Scale, Process-Based Erosion Model

Soil erosion affects agricultural productivity, the natural environment and infrastructure security. Soil loss and its associated impacts are important environmental problems. Consequently, model-based predictions of erosion are beneficial for a variety of applications. Process-based erosion models are used to forecast sediment transport concentration as it varies temporally and spatially. Of these, the one-dimensional Hairsine-Rose model describes multiple particle size classes, rainfall detachment, flow-driven entrainment and deposition. This model has been evaluated for different experiments, and has been shown to reliably explain experimental data in a consistent manner. It is common on both the hillslope and laboratory scales to apply one-dimensional erosion models even though the overland flow and sediment transport is two-dimensional. One-dimensional parameter determinations, which are based typically on outflow data, implicitly average the two-dimensional flow. Here we compare experimentally and numerically this averaging process for the Hairsine-Rose model. For this purpose, laboratory experiments were performed using different configurations of the 2 m × 6 m EPFL erosion flume. The flume was divided into 4 smaller flumes, with widths of 1 m, 0.5 m, and 2 × 0.25 m, but otherwise identical. A series of experiments was to provide data sets for analysis by the Hairsine-Rose model. After running the experiments, the amount of the eroded sediment in each subplot was assessed by comparing the temporal variation of eroded mass to evaluate the effect of, and sensitivity to, transverse width on erosion dynamics. The surface elevation changes due to erosion were examined to provide further understanding of the erosion data. A high resolution laser scanner provided details of the soil surface in the form of digital terrain maps before and after the experiment. This method presents a promising way for identification of spatial distribution pattern of eroded soil. In addition, we ran simulations using a fully two dimensional implementation of the Hairsine-Rose model for erosive flows with varying topography with spatially dependent flow and erosion input parameters to produce both outflow hydrographs and suspended sediment graphs. The data were integrated transversely and, as for the experimental data, the one-dimensional Hairsine-Rose model was used to fit the integrated data and so provide parameter estimates to compare with the two-dimensional input values