Crowdsourcing EO datasets to improve cloud detection algorithms and land cover change

Involving citizens in science is gaining considerable traction of late. With positive examples (e.g. Geo-Wiki, FotoQuest Austria), a number of projects are exploring the options to engage the public in contributing to scientific research, often by asking participants to collect some data or validate some results. The International Institute for Applied Systems Analysis (IIASA), with extensive experience in crowdsourcing and gamification, has joined Sinergise, Copernicus Masters 2016 winners, to engage the public in an initiative involving ESA’s Sentinel-2 satellite imagery. Sentinel-2 imagery offers high revisit times and sufficient resolution for land change detection applications. Unfortunately, simple (but fast) algorithms often fail due to many false-positives: changes in clouds are perceived as land changes. The ability to discriminate of cloudy pixels is thus crucial for any automatic or semi-automatic solutions that detect land change. A plethora of algorithms to distinguish clouds in Sentinel-2 data are available. However, there is a need for better data on where and when clouds occur to help improve these algorithms. To overcome this current gap in the data, we are engaging the public in this task. Using a number of tools, developed at IIASA, and Sentinel Hub services, which provide fast access to the entire global archive of Sentinel-2 data, the aim is to obtain a large data resource of curated cloud classifications. The resulting dataset will be published as open data and made available through Geopedia platform. The gamified process will start by asking users if there are clouds on a small image (e.g. 8x8 pixels at the highest Sentinel-2 resolution of 10 m/px), which will provide us with a screening process to pinpoint cloudy areas, employing Picture Pile crowdsourcing game from IIASA. The next step will involve a more detailed workflow, as users will get a slightly larger image (e.g. 64x64 pixels) and will then be asked to delineate different types of clouds: opaque clouds (nothing is seen through the clouds), thick clouds (where the surface is still discernible through the clouds), and thin clouds (where the surface is unequivocally covered by a cloud); the rest of the image will be implicitly cloud-free. The resulting data will be made available through the Geopedia portal, both for exploring and downloading. This paper will demonstrate this process and show some results from a crowdsourcing campaign. The approach will also allow us to collect other datasets in a rapid and efficient manner. For example, using a slightly modified configuration, a similar workflow could be used to obtain a manually curated land cover classification data set, which could be used as training data for machine learning algorithms

HADES RV Programme with HARPS-N at TNG VI. GJ 3942 b behind dominant activity signals

Short- to mid-term magnetic phenomena on the stellar surface of M-type stars cannot only resemble the effects of planets in radial velocity data, but also may hide them. We analyze 145 spectroscopic HARPS-N observations of GJ 3942 taken over the past five years and additional photometry to disentangle stellar activity effects from genuine Doppler signals as a result of the orbital motion of the star around the common barycenter with its planet. To achieve this, we use the common methods of pre-whitening, and treat the correlated red noise by a first-order moving average term and by Gaussian-process regression following an MCMC analysis. We identify the rotational period of the star at 16.3 days and discover a new super-Earth, GJ 3942 b, with an orbital period of 6.9 days and a minimum mass of 7.1 Me. An additional signal in the periodogram of the residuals is present but we cannot claim it to be related to a second planet with sufficient significance at this point. If confirmed, such planet candidate would have a minimum mass of 6.3 Me and a period of 10.4 days, which might indicate a 3:2 mean-motion resonance with the inner planet

The HADES RV Programme with HARPS-N@TNG IV. Time resolved analysis of the Ca ii H&K and H{\alpha} chromospheric emission of low-activity early-type M dwarfs

M dwarfs are prime targets for planet search programs, particularly of those focused on the detection and characterization of rocky planets in the habitable zone. Understanding their magnetic activity is important because it affects our ability to detect small planets, and it plays a key role in the characterization of the stellar environment. We analyze observations of the Ca II H&K and H{\alpha} lines as diagnostics of chromospheric activity for low-activity early-type M dwarfs. We analyze the time series of spectra of 71 early-type M dwarfs collected for the HADES project for planet search purposes. The HARPS-N spectra provide simultaneously the H&K doublet and the H{\alpha} line. We develop a reduction scheme able to correct the HARPS-N spectra for instrumental and atmospheric effects, and to provide flux-calibrated spectra in units of flux at the stellar surface. The H&K and H{\alpha} fluxes are compared with each other, and their variability is analyzed. We find that the H and K flux excesses are strongly correlated with each other, while the H{\alpha} flux excess is generally less correlated with the H&K doublet. We also find that H{\alpha} emission does not increase monotonically with the H&K line flux, showing some absorption before being filled in by chromospheric emission when H&K activity increases. Analyzing the time variability of the emission fluxes, we derive a tentative estimate of the rotation period (of the order of a few tens of days) for some of the program stars, and the typical lifetime of chromospheric active regions (a few stellar rotations). Our results are in good agreement with previous studies. In particular, we find evidence that the chromospheres of early-type M dwarfs could be characterized by different filaments coverage, affecting the formation mechanism of the H{\alpha} line. We also show that chromospheric structure is likely related to spectral type

The HADES RV Programme with HARPS-N at TNG XI. GJ 685 b: a warm super-Earth around an active M dwarf

Small rocky planets seem to be very abundant around low-mass M-type stars. Their actual planetary population is however not yet precisely understood. Currently several surveys aim to expand the statistics with intensive detection campaigns, both photometric and spectroscopic. We analyse 106 spectroscopic HARPS-N observations of the active M0-type star GJ 685 taken over the past five years. We combine these data with photometric measurements from different observatories to accurately model the stellar rotation and disentangle its signals from genuine Doppler planetary signals in the RV data. We run an MCMC analysis on the RV and activity indexes time series to model the planetary and stellar signals present in the data, applying Gaussian Process regression technique to deal with the stellar activity signals. We identify three periodic signals in the RV time series, with periods of 9, 24, and 18 d. Combining the analyses of the photometry of the star with the activity indexes derived from the HARPS-N spectra, we identify the 18 d and 9 d signals as activity-related, corresponding to the stellar rotation period and its first harmonic respectively. The 24 d signals shows no relations with any activity proxy, so we identify it as a genuine planetary signal. We find the best-fit model describing the Doppler signal of the newly-found planet, GJ 685\,b, corresponding to an orbital period $P_b = 24.160^{+0.061}_{-0.047}$ d and a minimum mass $M_P \sin i = 9.0^{+1.7}_{-1.8}$ M$_\oplus$. We also study a sample of 70 RV-detected M-dwarf planets, and present new statistical evidence of a difference in mass distribution between the populations of single- and multi-planet systems, which can shed new light on the formation mechanisms of low-mass planets around late-type stars.Comment: 18 pages, 13 figures, accepted for publication in A&

The HADES RV Programme with HARPS-N@TNG II. Data treatment and simulations

The distribution of exoplanets around low-mass stars is still not well understood. Such stars, however, present an excellent opportunity of reaching down to the rocky and habitable planet domains. The number of current detections used for statistical purposes is still quite modest and different surveys, using both photometry and precise radial velocities, are searching for planets around M dwarfs. Our HARPS-N red dwarf exoplanet survey is aimed at the detection of new planets around a sample of 78 selected stars, together with the subsequent characterization of their activity properties. Here we investigate the survey performance and strategy. From 2700 observed spectra, we compare the radial velocity determinations of the HARPS-N DRS pipeline and the HARPS-TERRA code, we calculate the mean activity jitter level, we evaluate the planet detection expectations, and we address the general question of how to define the strategy of spectroscopic surveys in order to be most efficient in the detection of planets. We find that the HARPS-TERRA radial velocities show less scatter and we calculate a mean activity jitter of 2.3 m/s for our sample. For a general radial velocity survey with limited observing time, the number of observations per star is key for the detection efficiency. In the case of an early M-type target sample, we conclude that approximately 50 observations per star with exposure times of 900 s and precisions of about 1 m/s maximizes the number of planet detections

Elastic and vibrational properties of alpha and beta-PbO

The structure, electronic and dynamic properties of the two layered alpha (litharge) and beta (massicot) phases of PbO have been studied by density functional methods. The role of London dispersion interactions as leading component of the total interaction energy between layers has been addressed by using the Grimme's approach, in which new parameters for Pb and O atoms have been developed. Both gradient corrected and hybrid functionals have been adopted using Gaussian-type basis sets of polarized triple zeta quality for O atoms and small core pseudo-potential for the Pb atoms. Basis set superposition error (BSSE) has been accounted for by the Boys-Bernardi correction to compute the interlayer separation. Cross check with calculations adopting plane waves that are BSSE free have also been performed for both structures and vibrational frequencies. With the new set of proposed Grimme's type parameters structures and dynamical parameters for both PbO phases are in good agreement with experimental data.Comment: 8 pages, 5 figure

HADES RV Programme with HARPS-N at TNG. VII. Rotation and activity of M-Dwarfs from time-series high-resolution spectroscopy of chromospheric indicators

We aim to investigate the presence of signatures of magnetic cycles and rotation on a sample of 71 early M-dwarfs from the HADES RV programme using high-resolution time-series spectroscopy of the Ca II H & K and Halpha chromospheric activity indicators, the radial velocity series, the parameters of the cross correlation function and the V-band photometry. We used mainly HARPS-N spectra, acquired over four years, and add HARPS spectra from the public ESO database and ASAS photometry light-curves as support data, extending the baseline of the observations of some stars up to 12 years. We provide log(R'hk) measurements for all the stars in the sample, cycle length measurements for 13 stars, rotation periods for 33 stars and we are able to measure the semi-amplitude of the radial velocity signal induced by rotation in 16 stars. We complement our work with previous results and confirm and refine the previously reported relationships between the mean level of chromospheric emission, measured by the log(R'hk), with the rotation period, and with the measured semi-amplitude of the activity induced radial velocity signal for early M-dwarfs. We searched for a possible relation between the measured rotation periods and the lengths of the magnetic cycle, finding a weak correlation between both quantities. Using previous v sin i measurements we estimated the inclinations of the star's poles to the line of sight for all the stars in the sample, and estimate the range of masses of the planets GJ 3998 b and c (2.5 - 4.9 Mearth and 6.3 - 12.5 Mearth), GJ 625 b (2.82 Mearth), GJ 3942 b (7.1 - 10.0 Mearth) and GJ 15A b (3.1 - 3.3 Mearth), assuming their orbits are coplanar with the stellar rotation.Comment: 19 pages, 16 figures, 10 table

Acceleration of generalized hypergeometric functions through precise remainder asymptotics

We express the asymptotics of the remainders of the partial sums {s_n} of the generalized hypergeometric function q+1_F_q through an inverse power series z^n n^l \sum_k c_k/n^k, where the exponent l and the asymptotic coefficients {c_k} may be recursively computed to any desired order from the hypergeometric parameters and argument. From this we derive a new series acceleration technique that can be applied to any such function, even with complex parameters and at the branch point z=1. For moderate parameters (up to approximately ten) a C implementation at fixed precision is very effective at computing these functions; for larger parameters an implementation in higher than machine precision would be needed. Even for larger parameters, however, our C implementation is able to correctly determine whether or not it has converged; and when it converges, its estimate of its error is accurate.Comment: 36 pages, 6 figures, LaTeX2e. Fixed sign error in Eq. (2.28), added several references, added comparison to other methods, and added discussion of recursion stabilit

Morphological Correlates of Mechanotransduction in Acousticolateral Hair Cells

The development of ideas on mechanotransduction in acousticolateral hair cells is described, leading to the current idea that transduction depends on deflection of the bundle of stereocilia by a force parallel to the plane of the sensory epithelium. Electrophysiological experiments are summarised, suggesting that transduction depends on a shear between the different rows of stereocilia, and that the transducer channels are situated towards the tips of the stereocilia. Analysis of the ways that shear between the rows of stereocilia could be detected suggests that tip links are the structures which are most likely to transmit the stimulus-induced forces to the transducer channels on the membrane. The directional selectivity of mechanotransduction is associated with the position of the kinocilium and gradation in heights of the stereocilia; evidence is presented suggesting that in the lateral line these are partly determined by the mitosis giving rise to the hair cell. Tip links differentiate out of links which initially join the stereocilia in all directions, with their final spatial organisation, which sets the directional selectivity of mechanotransduction, probably being determined by the gradient in growth of the stereocilia