Institut national de l’information géographique et forestière

Abstract For the first time in the International Terrestrial Reference Frame (ITRF) history, the ITRF2014 is generated with an enhanced modeling of nonlinear station motions, including seasonal (annual and semiannual) signals of station positions and postseismic deformation for sites that were subject to major earthquakes. Using the full observation history of the four space geodetic techniques (very long baseline interferometry (VLBI), satellite laser ranging (SLR), Global Navigation Satellite Systems (GNSS), and Doppler orbitography and radiopositioning integrated by satellite (DORIS)), the corresponding international services provided reprocessed time series (weekly from SLR and DORIS, daily from GNSS, and 24 h session‐wise from VLBI) of station positions and daily Earth Orientation Parameters. ITRF2014 is demonstrated to be superior to past ITRF releases, as it precisely models the actual station trajectories leading to a more robust secular frame and site velocities. The ITRF2014 long‐term origin coincides with the Earth system center of mass as sensed by SLR observations collected on the two LAGEOS satellites over the time span between 1993.0 and 2015.0. The estimated accuracy of the ITRF2014 origin, as reflected by the level of agreement with the ITRF2008 (both origins are defined by SLR), is at the level of less than 3 mm at epoch 2010.0 and less than 0.2 mm/yr in time evolution. The ITRF2014 scale is defined by the arithmetic average of the implicit scales of SLR and VLBI solutions as obtained by the stacking of their respective time series. The resulting scale and scale rate differences between the two solutions are 1.37 (±0.10) ppb at epoch 2010.0 and 0.02 (±0.02) ppb/yr. While the postseismic deformation models were estimated using GNSS/GPS data, the resulting parametric models at earthquake colocation sites were applied to the station position time series of the three other techniques, showing a very high level of consistency which enforces more the link between techniques within the ITRF2014 frame. The users should be aware that the postseismic deformation models are part of the ITRF2014 products, unlike the annual and semiannual signals, which were estimated internally with the only purpose of enhancing the velocity field estimation of the secular frame.

ITRF2008 is a refined version of the International Terrestrial Reference Frame based on reprocessed solutions of the four space geodetic techniques: VLBI, SLR, GPS and DORIS, spanning 29, 26, 12.5 and 16 years of observations, respectively. The input data used in its elaboration are time series (weekly from satellite techniques and 24-h session-wise from VLBI) of station positions and daily Earth Orientation Parameters (EOPs). The ITRF2008 origin is defined in such a way that it has zero translations and translation rates with respect to the mean Earth center of mass, averaged by the SLR time series. Its scale is defined by nullifying the scale factor and its rate with respect to the mean of VLBI and SLR long-term solutions as obtained by stacking their respective time series. The scale agreement between these two technique solutions is estimated to be 1.05 ± 0.13 ppb at epoch 2005.0 and 0.049 ± 0.010 ppb/yr. The ITRF2008 orientation (at epoch 2005.0) and its rate are aligned to the ITRF2005 using 179 stations of high geodetic quality. An estimate of the origin components from ITRF2008 to ITRF2005 (both origins are defined by SLR) indicates differences at epoch 2005.0, namely: −0.5, −0.9 and −4.7 mm along X, Y and Z-axis, respectively. The translation rate differences between the two frames are zero for Y and Z, while we observe an X-translation rate of 0.3 mm/yr. The estimated formal errors of these parameters are 0.2 mm and 0.2 mm/yr, respectively. The high level of origin agreement between ITRF2008 and ITRF2005 is an indication of an imprecise ITRF2000 origin that exhibits a Z-translation drift of 1.8 mm/yr with respect to ITRF2005. An evaluation of the ITRF2008 origin accuracy based on the level of its agreement with ITRF2005 is believed to be at the level of 1 cm over the time-span of the SLR observations. Considering the level of scale consistency between VLBI and SLR, the ITRF2008 scale accuracy is evaluated to be at the level of 1.2 ppb (8 mm at the equator) over the common time-span of the observations of both techniques. Although the performance of the ITRF2008 is demonstrated to be higher than ITRF2005, future ITRF improvement resides in improving the consistency between local ties in co-location sites and space geodesy estimates.

Unlike the past International Terrestrial Reference Frame (ITRF) versions where global long‐term solutions were combined, the ITRF2005 uses as input data time series (weekly from satellite techniques and 24‐h session‐wise from Very Long Baseline Interferometry) of station positions and daily Earth Orientation Parameters (EOPs). The advantage of using time series of station positions is that it allows to monitor station non‐linear motion and discontinuities and to examine the temporal behavior of the frame physical parameters, namely the origin and the scale. The ITRF2005 origin is defined in such a way that it has zero translations and translation rates with respect to the Earth center of mass, averaged by the Satellite Laser Ranging (SLR) time series spanning 13 years of observations. Its scale is defined by nullifying the scale and its rate with respect to the Very Long Baseline Interferometry (VLBI) time series spanning 26 years of observations. The ITRF2005 orientation (at epoch 2000.0) and its rate are aligned to the ITRF2000 using 70 stations of high geodetic quality. The estimated level of consistency of the ITRF2005 origin (at epoch 2000.0) and its rate with respect to the ITRF2000 is respectively 0.1, 0.8, 5.8 mm and 0.2, 0.1, 1.8 mm/yr along the X , Y and Z ‐axis. We estimate the formal errors on these components to be 0.3 mm and 0.3 mm/yr. We believe that this low level of agreement between the two frame origins is most probably due to the poor SLR network geometry and its degradation over time. The ITRF2005 combination involving 84 co‐location sites revealed a scale inconsistency of 1 ppb (6.3 mm at the equator), at epoch 2000.0, and 0.08 ppb/yr between the SLR and VLBI long‐term solutions as obtained by the stacking of their respective time series. Possible causes of this inconsistency may include the poor SLR and VLBI networks and their co‐locations, local tie uncertainties, systematic effects and possible inconsistent model corrections used in the data analysis of both techniques. For the first time of the ITRF history, the ITRF2005 rigorous combination provides self‐consistent series of EOPs, including Polar Motion from VLBI and satellite techniques and Universal Time and Length of Day from VLBI only. A velocity field of 152 sites with an error less than 1.5 mm/yr is used to estimate absolute rotation poles of 15 tectonic plates that are consistent with the ITRF2005 frame. This new absolute plate motion model supersedes and significantly improves that of the ITRF2000 which involved six major tectonic plates.

For the first time in the history of the International Terrestrial Reference Frame, the ITRF2000 combines unconstrained space geodesy solutions that are free from any tectonic plate motion model. Minimum constraints are applied to these solutions solely in order to define the underlying terrestrial reference frame (TRF). The ITRF2000 origin is defined by the Earth center of mass sensed by satellite laser ranging (SLR) and its scale by SLR and very long baseline interferometry. Its orientation is aligned to the ITRF97 at epoch 1997.0, and its orientation time evolution follows, conventionally, that of the no‐net‐rotation NNR‐NUVEL‐1A model. The ITRF2000 orientation and its rate are implemented using a consistent geodetic method, anchored over a selection of ITRF sites of high geodetic quality, ensuring a datum definition at the 1 mm level. This new frame is the most extensive and accurate one ever developed, containing about 800 stations located at about 500 sites, with better distribution over the globe compared to past ITRF versions but still with more site concentration in western Europe and North America. About 50% of station positions are determined to better than 1 cm, and about 100 sites have their velocity estimated to at (or better than) 1 mm/yr level. The ITRF2000 velocity field was used to estimate relative rotation poles for six major tectonic plates that are independent of the TRF orientation rate. A comparison to relative rotation poles of the NUVEL‐1A plate motion model shows vector differences ranging between 0.03° and 0.08°/m.y. (equivalent to approximately 1–7 mm/yr over the Earth's surface). ITRF2000 angular velocities for four plates, relative to the Pacific plate, appear to be faster than those predicted by the NUVEL‐1A model. The two most populated plates in terms of space geodetic sites, North America and Eurasia, exhibit a relative Euler rotation pole of about 0.056 (±0.005)°/m.y. faster than the pole predicted by NUVEL‐1A and located about (10°N, 7°E) more to the northwest, compared to that model.

Global pressure and temperature 2 wet (GPT2w) is an empirical troposphere delay model providing the mean values plus annual and semiannual amplitudes of pressure, temperature and its lapse rate, water vapor pressure and its decrease factor, weighted mean temperature, as well as hydrostatic and wet mapping function coefficients of the Vienna mapping function 1. All climatological parameters have been derived consistently from monthly mean pressure level data of ERA-Interim fields (European Centre for Medium-Range Weather Forecasts Re-Analysis) with a horizontal resolution of 1°, and the model is suitable to calculate slant hydrostatic and wet delays down to 3° elevation at sites in the vicinity of the earth surface using the date and approximate station coordinates as input. The wet delay estimation builds upon gridded values of the water vapor pressure, the weighted mean temperature, and the water vapor decrease factor, with the latter being tuned to ray-traced zenith wet delays. Comparisons with zenith delays at 341 globally distributed global navigation satellite systems stations show that the mean bias over all stations is below 1 mm and the mean standard deviation is about 3.6 cm. The GPT2w model with the gridded input file is provided at http://ggosatm.hg.tuwien.ac.at/DELAY/SOURCE/GPT2w/ .

As part of the Copernicus programme of the European Commission (EC), the European Space Agency (ESA) has developed and is currently operating the Sentinel-2 mission that is acquiring high spatial resolution optical imagery. This article provides a description of the calibration activities and the status of the mission products validation activities after one year in orbit. Measured performances, from the validation activities, cover both Top-Of-Atmosphere (TOA) and Bottom-Of-Atmosphere (BOA) products. The presented results show the good quality of the mission products both in terms of radiometry and geometry and provide an overview on next mission steps related to data quality aspects.

Deep-sea ecosystems represent the largest biome of the global biosphere, but knowledge of their biodiversity is still scant. The Mediterranean basin has been proposed as a hot spot of terrestrial and coastal marine biodiversity but has been supposed to be impoverished of deep-sea species richness. We summarized all available information on benthic biodiversity (Prokaryotes, Foraminifera, Meiofauna, Macrofauna, and Megafauna) in different deep-sea ecosystems of the Mediterranean Sea (200 to more than 4,000 m depth), including open slopes, deep basins, canyons, cold seeps, seamounts, deep-water corals and deep-hypersaline anoxic basins and analyzed overall longitudinal and bathymetric patterns. We show that in contrast to what was expected from the sharp decrease in organic carbon fluxes and reduced faunal abundance, the deep-sea biodiversity of both the eastern and the western basins of the Mediterranean Sea is similarly high. All of the biodiversity components, except Bacteria and Archaea, displayed a decreasing pattern with increasing water depth, but to a different extent for each component. Unlike patterns observed for faunal abundance, highest negative values of the slopes of the biodiversity patterns were observed for Meiofauna, followed by Macrofauna and Megafauna. Comparison of the biodiversity associated with open slopes, deep basins, canyons, and deep-water corals showed that the deep basins were the least diverse. Rarefaction curves allowed us to estimate the expected number of species for each benthic component in different bathymetric ranges. A large fraction of exclusive species was associated with each specific habitat or ecosystem. Thus, each deep-sea ecosystem contributes significantly to overall biodiversity. From theoretical extrapolations we estimate that the overall deep-sea Mediterranean biodiversity (excluding prokaryotes) reaches approximately 2805 species of which about 66% is still undiscovered. Among the biotic components investigated (Prokaryotes excluded), most of the unknown species are within the phylum Nematoda, followed by Foraminifera, but an important fraction of macrofaunal and megafaunal species also remains unknown. Data reported here provide new insights into the patterns of biodiversity in the deep-sea Mediterranean and new clues for future investigations aimed at identifying the factors controlling and threatening deep-sea biodiversity.

GipsyX/RTGx is the Jet Propulsion Laboratory’s (JPL) next generation software package for positioning, navigation, timing, and Earth science using measurements from three geodetic techniques: Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS); with Very Long Baseline Interferometry (VLBI) under development. The software facilitates combined estimation of geodetic and geophysical parameters using a Kalman filter approach on real or simulated data in both post-processing and in real-time. The estimated parameters include station coordinates and velocities, satellite orbits and clocks, Earth orientation, ionospheric and tropospheric delays. The software is also capable of full realization of a dynamic terrestrial reference through analysis and combination of time series of ground station coordinates. Applying lessons learned from its predecessors, GIPSY-OASIS and Real Time GIPSY (RTG), GipsyX/RTGx was re-designed from the ground up to offer improved precision, accuracy, usability, and operational flexibility. We present some key aspects of its new architecture, and describe some of its major applications, including Real-time orbit determination and ephemeris predictions in the U.S. Air Force Next Generation GPS Operational Control Segment (OCX), as well as in JPL’s Global Differential GPS (GDGPS) System, supporting User Range Error (URE) of <5 cm RMS; precision post-processing GNSS orbit determination, including JPL’s contributions to the International GNSS Service (IGS) with URE in the 2 cm RMS range; Precise point positioning (PPP) with ambiguity resolution, both statically and kinematically, for geodetic applications with 2 mm horizontal, and 6.5 mm vertical repeatability for static positioning; Operational orbit and clock determination for Low Earth Orbiting (LEO) satellites, such as NASA’s Gravity Recovery and Climate Experiment (GRACE) mission with GRACE relative clock alignment at the 20 ps level; calibration of radio occultation data from LEO satellites for weather forecasting and climate studies; Satellite Laser Ranging (SLR) to GNSS and LEO satellites, DORIS-based and multi-technique orbit determination for LEO; production of terrestrial reference frames and Earth rotation parameters in support of JPL’s contribution to the International Terrestrial Reference Frame (ITRF).

The recent development of operational small unmanned aerial systems (UASs) opens the door for their extensive use in forest mapping, as both the spatial and temporal resolution of UAS imagery better suit local-scale investigation than traditional remote sensing tools. This article focuses on the use of combined photogrammetry and “Structure from Motion” approaches in order to model the forest canopy surface from low-altitude aerial images. An original workflow, using the open source and free photogrammetric toolbox, MICMAC (acronym for Multi Image Matches for Auto Correlation Methods), was set up to create a digital canopy surface model of deciduous stands. In combination with a co-registered light detection and ranging (LiDAR) digital terrain model, the elevation of vegetation was determined, and the resulting hybrid photo/LiDAR canopy height model was compared to data from a LiDAR canopy height model and from forest inventory data. Linear regressions predicting dominant height and individual height from plot metrics and crown metrics showed that the photogrammetric canopy height model was of good quality for deciduous stands. Although photogrammetric reconstruction significantly smooths the canopy surface, the use of this workflow has the potential to take full advantage of the flexible revisit period of drones in order to refresh the LiDAR canopy height model and to collect dense multitemporal canopy height series.

The publication familiarizes the reader with MicMac - a free, open-source photogrammetric software for 3D reconstruction. A brief history of the tool, its organisation and unique features vis-à-vis other software tools are in the highlight. The essential algorithmic aspects of the structure from motion and image dense matching problems are discussed from the implementation and the user’s viewpoints.

The ITRF2008 velocity field is demonstrated to be of higher quality and more precise than past ITRF solutions. We estimated an absolute tectonic plate motion model made up of 14 major plates, using velocities of 206 sites of high geodetic quality (far from plate boundaries, deformation zones and Glacial Isostatic Adjustment (GIA) regions), derived from and consistent with ITRF2008. The precision of the estimated model is evaluated to be at the level of 0.3 mm/a WRMS. No GIA corrections were applied to site velocities prior to estimating plate rotation poles, as our selected sites are outside the Fennoscandia regions where the GIA models we tested are performing reasonably well, and far from GIA areas where the models would degrade the fit (Antarctica and North America). Our selected velocity field has small origin rate bias components following the three axis (X, Y, Z), respectively 0.41 ± 0.54, 0.22 ± 0.64 and 0.41 ± 0.60 (95 per cent confidence limits). Comparing our model to NNR‐NUVEL‐1A and the newly available NNR‐MORVEL56, we found better agreement with NNR‐MORVEL56 than with NNR‐NUVEL‐1A for all plates, except for Australia where we observe an average residual rotation rate of 4 mm/a. Using our selection of sites, we found large global X‐rotation rates between the two models (0.016°/Ma) and between our model and NNR‐MORVEL56 of 0.023°/Ma, equivalent to 2.5 mm/a at the Earth surface.

According to the weak equivalence principle, all bodies should fall at the same rate in a gravitational field. The MICROSCOPE satellite, launched in April 2016, aims to test its validity at the 10^{-15} precision level, by measuring the force required to maintain two test masses (of titanium and platinum alloys) exactly in the same orbit. A nonvanishing result would correspond to a violation of the equivalence principle, or to the discovery of a new long-range force. Analysis of the first data gives δ(Ti,Pt)=[-1±9(stat)±9(syst)]×10^{-15} (1σ statistical uncertainty) for the titanium-platinum Eötvös parameter characterizing the relative difference in their free-fall accelerations.

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p ≤ 5 × 10⁻⁷). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10⁻⁸) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p =2 × 10⁻⁸) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 × 10⁻⁸); rs1229984-ADH1B, p = 7 × 10⁻⁹; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

Abdel Omran's 1971 theory of "Epidemiologic Transition" was the first attempt to account for the extraordinary advances in health care made in industrialized countries since the 18th century. In the framework of the Demographic Transition, it implied a

International audience

For various geodetic and geophysical applications, users need to have access to a plate motion model (PMM) that is consistent with the ITRF2014 frame. This paper describes the approach used for determining a PMM from the horizontal velocities of a subset of the ITRF2014 sites away from plate boundaries, Glacial Isostatic Adjustment regions and other deforming zones. In theory it would be necessary to include in the inversion model a translational motion vector (called in this paper origin rate bias, ORB) that would represent the relative motion between the ITRF2014 origin (long-term averaged centre of mass of the Earth as sensed by SLR) and the centre of tectonic plate motion. We show that in practice, the magnitude of the estimated ORB is strongly dependent on the selection of ITRF2014 sites used for the PMM adjustment. Its Z-component can in particular range between 0 and more than 1 mm yr−1 depending on the station network used, preventing any geophysical interpretation of the estimated value. Relying on rigorous statistical criteria, the site selection finally adopted for the ITRF2014-PMM adjustment leads to a relatively small ORB (0.30 ± 0.18 mm yr−1 in the Z-component), which is statistically insignificant at the 2-sigma level, but also according to an F-ratio test. Therefore we opted for an ITRF2014-PMM without estimating the ORB, which in turn accommodates geodetic applications that require access to the ITRF2014 frame through pure plate rotation poles.

Global Positioning System (GPS) measurement campaigns in Myanmar, conducted in 1998 and 2000, allow quantifying the present‐day crustal deformation around the Sagaing fault system in central Myanmar. Both a regional network installed at four points within the country and a local 18‐station network centered on the city of Mandalay across the Sagaing fault demonstrate that active deformation related to the northward motion of India is distributed across Myanmar in a platelet that extends from the western edge of the Shan Plateau in the east to the Andaman Trench in the west. In this platelet, deformation is rather diffuse and distributed over distinct fault systems. In the east, the Sagaing/Shan Scarp fault system absorbs &lt;20 mm/yr of the 35 mm/yr India/Sundaland strike‐slip motion. Along this major plate boundary, strain is partitioned along the N‐S trending Sagaing fault and the transtensile N160°E trending Shan Scarp fault. Shortening and wrenching within the inverted central Myanmar basins, strike‐slip faults affecting the Arakan Yoma fold‐and‐thrust belt, and oblique subduction along the Andaman trench should absorb the remaining India/Sundaland motion (&gt;10 mm/yr). This GPS study combined with an on land geotectonic survey demonstrate that oblique slip of India along the rigid Sundaland block is accommodated by a partitioned system characterized by distribution of deformation over a wide zone.

Abstract. Global navigation satellite systems (GNSSs) have revolutionised positioning, navigation, and timing, becoming a common part of our everyday life. Aside from these well-known civilian and commercial applications, GNSS is now an established atmospheric observing system, which can accurately sense water vapour, the most abundant greenhouse gas, accounting for 60–70 % of atmospheric warming. In Europe, the application of GNSS in meteorology started roughly two decades ago, and today it is a well-established field in both research and operation. This review covers the state of the art in GNSS meteorology in Europe. The advances in GNSS processing for derivation of tropospheric products, application of GNSS tropospheric products in operational weather prediction and application of GNSS tropospheric products for climate monitoring are discussed. The GNSS processing techniques and tropospheric products are reviewed. A summary of the use of the products for validation and impact studies with operational numerical weather prediction (NWP) models as well as very short weather prediction (nowcasting) case studies is given. Climate research with GNSSs is an emerging field of research, but the studies so far have been limited to comparison with climate models and derivation of trends. More than 15 years of GNSS meteorology in Europe has already achieved outstanding cooperation between the atmospheric and geodetic communities. It is now feasible to develop next-generation GNSS tropospheric products and applications that can enhance the quality of weather forecasts and climate monitoring. This work is carried out within COST Action ES1206 advanced global navigation satellite systems tropospheric products for monitoring severe weather events and climate (GNSS4SWEC, http://gnss4swec.knmi.nl).

Abstract. IGN has developed a set of photogrammetric tools, APERO and MICMAC, for computing 3D models from set of images. This software, developed initially for its internal needs are now delivered as open source code. This paper focuses on the presentation of APERO the orientation software. Compared to some other free software initiatives, it is probably more complex but also more complete, its targeted user is rather professionals (architects, archaeologist, geomophologist) than people. APERO uses both computer vision approach for estimation of initial solution and photogrammetry for a rigorous compensation of the total error; it has a large library of parametric model of distortion allowing a precise modelization of all the kind of pinhole camera we know, including several model of fish-eye; there is also several tools for geo-referencing the result. The results are illustrated on various application, including the data-set of 3D-Arch workshop.

Abstract To better describe the shape of the constantly deforming Earth’s surface, the ITRF2020 is provided as an augmented terrestrial reference frame that precisely models nonlinear station motions for both seasonal (annual and semi-annual) signals present in the station position time series and Post-Seismic Deformation (PSD) for sites impacted by major earthquakes. Reprocessed solutions in the form of station position time series and Earth Orientation Parameters using the full observation history provided by the four space geodetic techniques (DORIS, GNSS, SLR and VLBI) were used as input data, spanning 28, 27, 38 and 41 years of observations, respectively. The ITRF2020 long-term origin follows linearly with time the Earth’s Center of Mass (CM) as sensed by SLR, based on observations collected over the time span 1993.0–2021.0. We evaluate the accuracy of the ITRF2020 long-term origin position and time evolution by comparison to previous solutions, namely ITRF2014, ITRF2008 and ITRF2005, to be at the level of or better than 5 mm and 0.5 mm/yr, respectively. The ITRF2020 long-term scale is defined by a rigorous weighted average of selected VLBI sessions up to 2013.75 and SLR weekly solutions covering the 1997.75–2021.0 time span. For the first time of the ITRF history, the scale agreement between SLR and VLBI long-term solutions is at the level of 0.15 ppb (1 mm at the equator) at epoch 2015.0, with no drift. To accommodate most of ITRF2020 users, the seasonal station coordinate variations are provided in the CM as well as in the Center of Figure frames, together with a seasonal geocenter motion model. While the PSD parametric models were determined by fitting GNSS data only, they also fit the station position time series of the three other techniques that are colocated with GNSS, demonstrating their high performance in describing site post-seismic trajectories.