Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach

The analysis step of the (ensemble) Kalman filter is optimal when (1) the distribution of the background is Gaussian, (2) state variables and observations are related via a linear operator, and (3) the observational error is of additive nature and has Gaussian distribution. When these conditions are largely violated, a pre-processing step known as Gaussian anamorphosis (GA) can be applied. The objective of this procedure is to obtain state variables and observations that better fulfil the Gaussianity conditions in some sense. In this work we analyse GA from a joint perspective, paying attention to the effects of transformations in the joint state variable/observation space. First, we study transformations for state variables and observations that are independent from each other. Then, we introduce a targeted joint transformation with the objective to obtain joint Gaussianity in the transformed space. We focus primarily in the univariate case, and briefly comment on the multivariate one. A key point of this paper is that, when (1)-(3) are violated, using the analysis step of the EnKF will not recover the exact posterior density in spite of any transformations one may perform. These transformations, however, provide approximations of different quality to the Bayesian solution of the problem. Using an example in which the Bayesian posterior can be analytically computed, we assess the quality of the analysis distributions generated after applying the EnKF analysis step in conjunction with different GA options. The value of the targeted joint transformation is particularly clear for the case when the prior is Gaussian, the marginal density for the observations is close to Gaussian, and the likelihood is a Gaussian mixture.

The high-resolution infrared spectrum of HOF near 2700 cm-1: the ground and 2v2 states

The a/b hybrid-type ν1 fundamental and 2ν2 overtone bands of HOF were investigated by FTIR spectroscopy with a resolution close to 0.008 cm−1. Improved ground state parameters of HOF were determined from a merge of more than 3000 ground state combination differences formed from ν1 and previously measured ν2 transitions with the reported pure rotational lines. Excited state parameters of the v2 = 2 state, ν0 = 2686.924 6(1) and χ22 = −9.942 4(1) cm−1, were determined employing Watson's A-reduced Hamiltonian up to sixth order in I′ representation. The 2ν2 state was found to be unperturbed, the excited state parameters being closely related to those of ν2.

Have aerosols caused the observed Atlantic multidecadal variability?

Identifying the prime drivers of the twentieth-century multidecadal variability in the Atlantic Ocean is crucial for predicting how the Atlantic will evolve in the coming decades and the resulting broad impacts on weather and precipitation patterns around the globe. Recently, Booth et al. showed that the Hadley Centre Global Environmental Model, version 2, Earth system configuration (HadGEM2-ES) closely reproduces the observed multidecadal variations of area-averaged North Atlantic sea surface temperature in the twentieth century. The multidecadal variations simulated in HadGEM2-ES are primarily driven by aerosol indirect effects that modify net surface shortwave radiation. On the basis of these results, Booth et al. concluded that aerosols are a prime driver of twentieth-century North Atlantic climate variability. However, here it is shown that there are major discrepancies between the HadGEM2-ES simulations and observations in the North Atlantic upper-ocean heat content, in the spatial pattern of multidecadal SST changes within and outside the North Atlantic, and in the subpolar North Atlantic sea surface salinity. These discrepancies may be strongly influenced by, and indeed in large part caused by, aerosol effects. It is also shown that the aerosol effects simulated in HadGEM2-ES cannot account for the observed anticorrelation between detrended multidecadal surface and subsurface temperature variations in the tropical North Atlantic. These discrepancies cast considerable doubt on the claim that aerosol forcing drives the bulk of this multidecadal variability.

Decays of new nuclides and isomers beyond the proton drip line — the influence of neutron configurations

The energy of the vh9/2 orbital in nuclei above N = 82 drops rapidly in energy relative to the vf7/2 orbital as the occupancy of the πh11/2 orbital increases. These two neutron orbitals become nearly degenerate as the proton drip line is approached. In this work, we have discovered the new nuclides 161Os and 157W, and studied the decays of the proton emitter 160Re in detail. The 161Os and 160Re nuclei were produced in reactions of 290, 300 and 310 MeV 58Ni ions with an isotopically enriched 106Cd target, separated in‐flight using the RITU separator and implanted into the GREAT spectrometer. The 161Os α a decays populated the new nuclide 157W, which decayed by β‐particle emission. The β decay fed the known α‐decaying 1/2+ and 11/2− states in 157Ta, which is consistent with a vf7/2 ground state in 157W. The measured α‐decay energy and half‐life for 161Os correspond to a reduced α‐decay width that is compatible with s‐wave α‐particle emission, implying that its ground state is also a vf7/2 state. Over 7000 160Re nuclei were produced and the γ decays of a new isomeric state feeding the πd3/2 level in 160Re were discovered, but no evidence for the proton or a decay of the expected πh11/2 state could be found. The isomer decays offer a natural explanation for this non‐observation and provides a striking example of the influence of the near degeneracy of the vh9/2 and vf7/2 orbitals on the properties of nuclei in this region.

Climatic impacts of land-use change due to crop yield increases and a universal carbon tax from a scenario model*

Future land cover will have a significant impact on climate and is strongly influenced by the extent of agricultural land use. Differing assumptions of crop yield increase and carbon pricing mitigation strategies affect projected expansion of agricultural land in future scenarios. In the representative concentration pathway 4.5 (RCP4.5) from phase 5 of the Coupled Model Intercomparison Project (CMIP5), the carbon effects of these land cover changes are included, although the biogeophysical effects are not. The afforestation in RCP4.5 has important biogeophysical impacts on climate, in addition to the land carbon changes, which are directly related to the assumption of crop yield increase and the universal carbon tax. To investigate the biogeophysical climatic impact of combinations of agricultural crop yield increases and carbon pricing mitigation, five scenarios of land-use change based on RCP4.5 are used as inputs to an earth system model [Hadley Centre Global Environment Model, version 2-Earth System (HadGEM2-ES)]. In the scenario with the greatest increase in agricultural land (as a result of no increase in crop yield and no climate mitigation) there is a significant -0.49 K worldwide cooling by 2100 compared to a control scenario with no land-use change. Regional cooling is up to -2.2 K annually in northeastern Asia. Including carbon feedbacks from the land-use change gives a small global cooling of -0.067 K. This work shows that there are significant impacts from biogeophysical land-use changes caused by assumptions of crop yield and carbon mitigation, which mean that land carbon is not the whole story. It also elucidates the potential conflict between cooling from biogeophysical climate effects of land-use change and wider environmental aims.