The HARPS search for southern extra-solar planets XXXVII. Five new long-period giant planets and a system update

We describe radial-velocity time series obtained by HARPS on the 3.60 m telescope in La Silla (ESO, Chile) over ten years and report the discovery of five new giant exoplanets in distant orbits; these new planets orbit the stars HD 564, HD 30669, HD 108341, and BD -114672. Their periods range from 492 to 1684 days, semi-major axes range from 1.2 to 2.69 AU, and eccentricities range from 0 to 0.85. Their minimum mass ranges from 0.33 to 3.5 MJup. We also refine the parameters of two planets announced previously around HD 113538, based on a longer series of measurements. The planets have a period of 663 ± 8 and 1818 ± 25 days, orbital eccentricities of 0.14 ± 0.08 and 0.20 ± 0.04, and minimum masses of 0.36 ± 0.04 and 0.93 ± 0.06 MJup. Finally, we report the discovery of a new hot-Jupiter planet around an active star, HD 103720; the planet has a period of 4.5557 ± 0.0001 days and a minimum mass of 0.62 ± 0.025 MJup. We discuss the fundamental parameters of these systems and limitations due to stellar activity in quiet stars with typical 2 m s-1 radial velocity precision.

The STEREO/PLASTIC response to solar wind ions (Flight measurements and models)

The Plasma and Supra-Thermal Ion Composition (PLASTIC) instrument is one of four experiment packages on board of the two identical STEREO spacecraft A and B, which were successfully launched from Cape Canaveral on 26 October 2006. During the two years of the nominal STEREO mission, PLASTIC is providing us with the plasma characteristics of protons, alpha particles, and heavy ions. PLASTIC will also provide key diagnostic measurements in the form of the mass and charge state composition of heavy ions. Three measurements (E/qk, time of flight, ESSD) from the pulse height raw data are used to characterize the solar wind ions from the solar wind sector, and part of the suprathermal particles from the wide-angle partition with respect to mass, atomic number and charge state. In this paper, we present a new method for flight data analysis based on simulations of the PLASTIC response to solar wind ions. We present the response of the entrance system / energy analyzer in an analytical form. Based on stopping power theory, we use an analytical expression for the energy loss of the ions when they pass through a thin carbon foil. This allows us to model analytically the response of the time of flight mass spectrometer to solar wind ions. Thus we present a new version of the analytical response of the solid state detectors to solar wind ions. Various important parameters needed for our models were derived, based on calibration data and on the first flight measurements obtained from STEREO-A. We used information from each measured event that is registered in full resolution in the Pulse Height Analysis words and we derived a new algorithm for the analysis of both existing and future data sets of a similar nature which was tested and works well. This algorithm allows us to obtain, for each measured event, the mass, atomic number and charge state in the correct physical units. Finally, an important criterion was developed for filtering our Fe raw flight data set from the pulse height data without discriminating charge states.

Holocene flood frequency across the Central Alps – solar forcing and evidence for variations in North Atlantic atmospheric circulation

The frequency of large-scale heavy precipitation events in the European Alps is expected to undergo substantial changes with current climate change. Hence, knowledge about the past natural variability of floods caused by heavy precipitation constitutes important input for climate projections. We present a comprehensive Holocene (10,000 years) reconstruction of the flood frequency in the Central European Alps combining 15 lacustrine sediment records. These records provide an extensive catalog of flood deposits, which were generated by flood-induced underflows delivering terrestrial material to the lake floors. The multi-archive approach allows suppressing local weather patterns, such as thunderstorms, from the obtained climate signal. We reconstructed mainly late spring to fall events since ice cover and precipitation in form of snow in winter at high-altitude study sites do inhibit the generation of flood layers. We found that flood frequency was higher during cool periods, coinciding with lows in solar activity. In addition, flood occurrence shows periodicities that are also observed in reconstructions of solar activity from C-14 and Be-10 records (2500-3000, 900-1200, as well as of about 710, 500, 350, 208 (Suess cycle), 150, 104 and 87 (Gleissberg cycle) years). As atmospheric mechanism, we propose an expansion/shrinking of the Hadley cell with increasing/decreasing air temperature, causing dry/wet conditions in Central Europe during phases of high/low solar activity. Furthermore, differences between the flood patterns from the Northern Alps and the Southern Alps indicate changes in North Atlantic circulation. Enhanced flood occurrence in the South compared to the North suggests a pronounced southward position of the Westerlies and/or blocking over the northern North Atlantic, hence resembling a negative NAO state (most distinct from 4.2 to 2.4 kyr BP and during the Little Ice Age). South-Alpine flood activity therefore provides a qualitative record of variations in a paleo-NAO pattern during the Holocene. Additionally, increased South Alpine flood activity contrasts to low precipitation in tropical Central America (Cariaco Basin) on the Holocene and centennial time scale. This observation is consistent with a Holocene southward migration of the Atlantic circulation system, and hence of the ITCZ, driven by decreasing summer insolation in the Northern hemisphere, as well as with shorter-term fluctuations probably driven by solar activity. (C) 2013 Elsevier Ltd. All rights reserved.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

The HARPS search for southern extra-solar planets

The meter-per-second precision achieved by today’s velocimeters enables us to search for 1−10 M⊕ planets in the habitable zone of cool stars. This paper reports on the detection of three planets orbiting GJ 163 (HIP 19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ 163 over a period of eight years. We searched the radial-velocity time series for coherent signals and found five distinct periodic variabilities. We investigated the stellar activity and called into question the planetary interpretation for two signals. Before more data can be acquired we concluded that at least three planets are orbiting GJ 163. They have orbital periods of Pb = 8.632 ± 0.002, Pc = 25.63 ± 0.03, and Pd = 604 ± 8 days and minimum masses msini = 10.6 ± 0.6, 6.8 ± 0.9, and 29 ± 3 M⊕, respectively. We hold our interpretations for the two additional signals with periods P(e) = 19.4 and P(f) = 108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. The planet GJ 163c, in particular, is a super-Earth with an equilibrium temperature of Teq = (302 ± 10)(1 − A)1/4 K and may lie in the so-called habitable zone for albedo values (A = 0.34 − 0.89) moderately higher than that of Earth (A⊕ = 0.2−0.3).

Solar forcing of North Atlantic surface temperature and salinity over the past millennium

There were several centennial-scale fluctuations in the climate and oceanography of the North Atlantic region over the past 1,000 years, including a period of relative cooling from about AD 1450 to 1850 known as the Little Ice Age1. These variations may be linked to changes in solar irradiance, amplified through feedbacks including the Atlantic meridional overturning circulation2. Changes in the return limb of the Atlantic meridional overturning circulation are reflected in water properties at the base of the mixed layer south of Iceland. Here we reconstruct thermocline temperature and salinity in this region from AD 818 to 1780 using paired δ18O and Mg/Ca ratio measurements of foraminifer shells from a subdecadally resolved marine sediment core. The reconstructed centennial-scale variations in hydrography correlate with variability in total solar irradiance. We find a similar correlation in a simulation of climate over the past 1,000 years. We infer that the hydrographic changes probably reflect variability in the strength of the subpolar gyre associated with changes in atmospheric circulation. Specifically, in the simulation, low solar irradiance promotes the development of frequent and persistent atmospheric blocking events, in which a quasi-stationary high-pressure system in the eastern North Atlantic modifies the flow of the westerly winds. We conclude that this process could have contributed to the consecutive cold winters documented in Europe during the Little Ice Age.

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

Cosmogenic 180W variations in meteorites and re-assessment of a possible 184Os–180W decay system

We measured tungsten (W) isotopes in 23 iron meteorites and the metal phase of the CB chondrite Gujba in order to ascertain if there is evidence for a large-scale nucleosynthetic heterogeneity in the p-process isotope 180W in the solar nebula as recently suggested by Schulz et al. (2013). We observed large excesses in 180W (up to ≈ 6 ε) in some irons. However, significant within-group variations in magmatic IIAB and IVB irons are not consistent with a nucleosynthetic origin, and the collateral effects on 180W from an s-deficit in IVB irons cannot explain the total variation. We present a new model for the combined effects of spallation and neutron capture reactions on 180W in iron meteorites and show that at least some of the observed within-group variability is explained by cosmic ray effects. Neutron capture causes burnout of 180W, whereas spallation reactions lead to positive shifts in 180W. These effects depend on the target composition and cosmic-ray exposure duration; spallation effects increase with Re/W and Os/W ratios in the target and with exposure age. The correlation of 180W/184W with Os/W ratios in iron meteorites results in part from spallogenic production of 180W rather than from 184Os decay, contrary to a recent study by Peters et al. (2014). Residual ε180W excesses after correction for an s-deficit and for cosmic ray effects may be due to ingrowth of 180W from 184Os decay, but the magnitude of this ingrowth is at least a factor of ≈2 smaller than previously suggested. These much smaller effects strongly limit the applicability of the putative 184Os-180W system to investigate geological problems.