Land Water Storage Changes from Ground and Space Geodesy: First Results from the GHYRAF (Gravity and Hydrology in Africa) Experiment

This paper is devoted to the first results from the GHYRAF (Gravity and Hydrology in Africa) experiment conducted since 2008 in West Africa and is aimed at investigating the changes in water storage in different regions sampling a strong rainfall gradient from the Sahara to the monsoon zone. The analysis of GPS vertical displacement in Niamey (Niger) and Djougou (Benin) shows that there is a clear annual signature of the hydrological load in agreement with global hydrology models like GLDAS. The comparison of GRACE solutions in West Africa, and more specifically in the Niger and Lake Chad basins, reveals a good agreement for the large scale annual water storage changes between global hydrology models and space gravity observations. Ground gravity observations done with an FG5 absolute gravimeter also show signals which can be well related to measured changes in soil and ground water. We present the first results for two sites in the Sahelian band (Wankama and Diffa in Niger) and one (Djougou in Benin) in the Sudanian monsoon region related to the recharge–discharge processes due to the monsoonal event in summer 2008 and the following dry season. It is confirmed that ground gravimetry is a useful tool to constrain local water storage changes when associated to hydrological and subsurface geophysical in situ measurements.     

Hydro-gravimetry in West-Africa: First results from the Djougou (Benin) superconducting gravimeter

The increasing number of hydro-gravimetry studies proves the rising interest of the hydrology community toward this monitoring method. The accuracy of superconducting gravimeters (SG) potentially allows the retrieval of small water storage changes (WSC) down to a few millimeters of equivalent water thickness. However, the importance of corrections applied to SG data to achieve such a precision in gravity residuals should be recalled. The Djougou permanent gravity station presented in this paper and located in northern Benin, West-Africa, provides a good opportunity to review these considerations. This station is equipped since July 2010 with the superconducting gravimeter SG-060 aimed at deriving WSC at different time-scales, daily to inter-annual. In this area, WSC are (1) part of the control system for evapotranspiration (ET) process, a key variable of the West-African monsoon cycle and (2) the state variable for resource management, a critical issue in storage-poor hard rock basement contexts such as in northern Benin. The potential for deriving WSC from time-lapse gravity data partly depends on environmental features such as topography and the instrument shelter. Therefore, this issue is addressed first, with the background idea that such sensitivity analysis should be undertaken before setting up any new instrument. In Djougou, local topography is quite flat leading to a theoretical straightforward relationship between gravity changes and WSC, close to the standard Bouguer value. However, the shelter plays a significant masking role, which is the principal limitation to the retrieval of fast hydrological processes such as ET following a rain event. Several issues concerning classical gravity corrections are also addressed in the paper. These include gap-filling procedures during rain-events and drift estimates for short time series. Special attention is provided to atmospheric corrections, and different approaches are tested: a simple scalar admittance, a filtered scalar admittance, a frequency-dependent admittance and direct atmospheric loading calculations. It is shown that the physically based approach of direct loading calculations performs better in both residual minimization and ET retrieval. Moreover, non-local hydrological effects are investigated and account for about 20% of the gravity residuals. Finally, gravity residuals are briefly analyzed at two distinct time scales: rapid (up to a few days) and seasonal. At the rapid time-scale, it is shown that ET retrieval is hardly achievable given shelter size and state-of-the-art atmospheric corrections. Still, mean values retrieved from this study are in accordance with known values of potential ET and lateral flow. Direct comparison of gravity changes with hydrological data (neutron probe monitoring and water table levels) show some discrepancies, particularly for the hydrological year of 2011, for which all hydrological data show a deficit, but SG and FG5 data do not. This preliminary analysis both provides a basis and call for further hydro-gravity modeling, to comprehensively investigate the water-cycle at the Djougou station.     

A search for atmospheric effects on gravity at different time and space scales

We investigate in this paper various approaches to correct gravity changes for the effect of atmospheric pressure changes. Two specific locations are considered: Strasbourg (France) as mid-latitude station, where regular pressure fronts occur and Djougou (Benin) as equatorial station with large thermally driven S₁ and S₂ waves of planetary extent. We first review the classical approaches based on a constant or frequency-dependent admittance using only local pressure and gravity data. We consider then a model of atmospheric loading and show the barometric admittance in terms of elastic, Newtonian and total load, as a function of the distance from the station. We consider both a 2D pressure model (surface loading) and a 2.5D model, where the density decreases with height (standard atmosphere). Assuming horizontal advection in the atmospheric dynamics, we convert this spatially dependent admittance into a frequency-dependent admittance. Using global pressure data from European Centre for Medium-Range Weather Forecasts (ECMWF) at about 12 km spatial resolution and 3 h sampling, we compute the model-predicted pressure admittance for Djougou and Strasbourg and we simulate the frequency dependence inferred from gravity and pressure observations below 4 cycle per day. A long gravity and pressure data set (1996–2013) from Strasbourg is used to investigate the low frequency part of the pressure admittance while a common 2.5 year data set (August 2010–February 2013) for Strasbourg and Djougou is then analyzed to investigate the high frequency part of the admittance. In both cases, our results are in close agreement with the predictions inferred from an atmospheric 2.5D loading model with a distance–time relationship due to horizontal advection. The frequency dependence of the barometric admittance is explained by the competing contributions of Newtonian attraction and elastic surface deformation according to the distance from the gravimeter. In the far field (low frequencies), the magnitude of the admittance decreases with frequency because of the combined elasticity effect and Newtonian attraction (when the atmosphere is below the horizon) while, on the contrary, in the near field (high frequencies), elasticity becomes negligible and the pressure admittance mainly decreases with increasing frequency because of the decreasing attraction effect of the atmospheric masses inside the cylindrical pressure cell centered on the sensor location of decreasing radius. In the last part, we show that there is variability in time in the pressure admittance for both stations.     

Long-term and seasonal gravity changes at the Strasbourg station and their relation to crustal deformation and hydrology

The time-varying gravity is observed at the Strasbourg station with super-conducting gravimeters (SG) since 1987, with a first record from 1987 to 1996 (GWR T005) and a second one in continuity from 1996 till now (GWR C026). The long-term behaviour of the SG is constrained by regular absolute gravity (AG) measurements, which are performed in parallel since 1989, first with the JILAg-5 instrument and later on with the FG5#206. Moreover, a permanent GPS station, which belongs to the French geodetic network, has been installed at the end of 1999. We will show that the AG measurements suggest that the gravity is slowly increasing in time at a rate of about 1.6 μGal/year, and besides, exhibits a quasi-annual component of several μGal variable amplitude. We present an analysis of the GPS data obtained at the SG station in the Rhine graben as well as at another regional station in the Vosges mountains distant by about 40 km in order to constrain the gravity contribution due to the vertical displacement of the station in the earth's gravity field (geometrical part). The tectonic context of the region is described and our first results from our two data sets of limited duration suggest a small subsidence of our station in the graben; however, this fact needs further confirmation when more geodetic data will be available. We also analyze the water table changes beneath the station (local scale) and in the Alsatian Plain (regional scale) to estimate the hydrological contributions from ground water to the gravity residual signal, and we show that some similarities exist between the gravity residuals and the hydrological contributions, especially for the seasonal terms. Other missing contributions of annual period (air mass motion in the atmosphere, ocean circulation, continental hydrology) have to be considered.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

Regional gravity variations in Europe from superconducting gravimeters

Recent satellite missions (CHAMP, GRACE) are now returning data on the time variation of the gravity field with harmonic coefficients computed every 4 weeks. The promise is to achieve a sub-microgal accuracy that will define continental mass variations involving large-scale hydrology. With this in mind, we examine the time varying gravity field over central Europe using a limited number of high quality ground-based superconducting gravimeter stations within the Global Geodynamics Project (GGP). Our purpose is to see whether there are coherent signals between the individual stations and to compare the regional component with that predicted from models of continental hydrology. The results are encouraging. We have found, using empirical orthogonal eigenfunctions of the gravity data that a clear annual signal is present that is consistent in phase (low amplitudes in summer) and amplitude (1–3 microgal) with that determined from a large-scale model of land water in connection with global climate modeling. More work is required to define how the gravity field is related to large-scale soil moisture and other mass variations, and we have yet to compare our results to the latest satellite-derived data.     

Assessment of GPS data for meteorological applications over Africa: Study of error sources and analysis of positioning accuracy

The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L₂ data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L₂ data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.     

三峡水库的建立对库区地球重力场及地壳形变的影响

作为人造工程的三峡水库的建立，使得库区水体聚集，地球质量重新分布，将使地壳的物理结构以及局部地球重力场发生变化．针对上述问题，对库区的水准面、点绝对重力值、垂线方向和高程、高程基准面等的重力场变化以及地壳形变进行了研究．主要结果是：当蓄水水位达到峰值１７５ｍ（坝高）时，大地水准面有２．２３－１１．２ｍｍ的变化，点重力值有（０．８３－４．６）×１０^－５ｍｓ^－２的变化，垂线偏差的变化分别为－０．６２＂－６．５０＂（南北方向）和－６．４２＂－１．３４＂（东西方向）地壳的形变量为１．３２－６．６５ｍｍ，这些变化将引起测区的高程产生３．５５－１７．８４ｍｍ的变化．因此，原有库区及其附近的测量资料（包括水准、天文、重力等）必须审慎使用，并应建立库区形变监测网，对地壳形变、地震、滑坡、大坝变形等进行监测．     

联合星载GPS和KBRR星间测速数据反演GRACE时变重力场模型

高精度GRACE卫星时变重力场反演一直是卫星重力测量中的难题.为了恢复高精度的时变地球重力场模型,本文联合GRACE卫星的星载GPS和KBR星间测速观测数据,在对GRACE卫星进行精密定轨的同时,解算出60阶月平均地球重力场模型.通过对GRACE卫星的定轨精度、星载GPS相位和KBR星间测速数据的拟合残差以及时变地球重力场模型解算精度等分析,表明:(1)与美国宇航局喷气推进实验室(JPL)发布的约化动力学精密轨道相比,本文确定GRACE卫星轨道三维位置误差小于5 cm.(2)星载GPS相位数据拟合残差为5~8 mm,KBR星间测速数据拟合残差为0.18~0.30μm·s~(-1).(3)解算的月平均重力场模型与美国德克萨斯大学空间研究中心(CSR)、德国地学研究中心(GFZ)和JPL发布的RL05模型精度接近,时变信号在全球范围内具有很好的空间分布一致性.通过计算亚马逊流域和长江流域的水储量变化,本文与上述三个机构的计算结果无明显差异,且相关系数均达0.9以上.可见,本文建立的卫星轨道与重力场同解算法具有反演高精度GRACE时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

Tidal loading along a profile Europe-East Africa-South Asia-Australia and the Pacific Ocean

A knowledge of the vertical component of the oceanic tidal load to a precision of at least one microgal is essential for the geophysical exploitation of the high-precision absolute and differential gravity measurements which are being made at ground level and in deep boreholes. On the other hand the ocean load and attraction signal contained in Earth tide gravity measurements can be extracted with a precision which is sufficient to characterize the behaviour of the oceanic tides in different basins and this provides a check of the validity of the presently proposed cotidal maps. The tidal gravity profiles made since 1971 from Europe to Polynesia, through East Africa, Asia and Australia, with correctly intercalibrated gravimeters, comprise information from 91 tidal gravity stations which is used in this paper with this goal in mind.A discussion of all possible sources of error is presented which shows that at the level of 0.5 μgal the observed effects cannot be ascribed to computational or instrumental errors. Cotidal maps which generate computed loads in agreement with the Earth tide gravity measurements over a sufficiently broad area can be used with confidence as a working standard to apply tidal corrections to high-precision measurements made by using new techniques in geodesy, geophysics and geodynamics, satellite altimetry, very long baseline interferometry, Moon and satellite laser ranging and absolute gravity. The recent cotidal maps calculated by Schwiderski for satellite altimetry reductions agree very well with land-based gravimeter observations of the diurnal components of the tides (O₁, K₁ and P₁ waves) but his semi-diurnal component maps (M₂, S₂ and N₂ waves) strangely appear less satisfactory in some large areas. The maps of Hendershott and Parke give good results in several large areas but not everywhere. More detailed investigations are needed not only for several coastal stations but mainly in the Himalayas.     

Seasonal water storage change of the Yangtze River basin detected by GRACE

1 Introduction Large-scale mass redistribution, or temporal varia- tion of mass within the Earth system, the driving force of interactions between solid Earth and geophysical fluids envelope (i.e., atmosphere, ocean, and hydro- sphere), is an important geophysical process critical to human life. Most of the interactions between solid Earth and the atmosphere/oceans happen at seasonal and inter-annual time scales. One important contribu- tor of mass redistribution at seasonal and inter-annual …     

Contributions of Satellite Laser Ranging to Past and Future Radar Altimetry Missions

Satellite laser ranging (SLR) has proven avery efficient method for contributingto the tracking of altimetric satellites anddetermining accurately their orbitalthough hampered by the non-worldwide coverageand the meteorologicalconditions. Indeed, in some cases it is the onlymethod available to determinethe satellite orbit (e.g., the orbits of the ERS-1and Geosat-Follow-On missions).Moreover, any operational and non-weather dependenttechniques, like GPS,DORIS, PRARE, can exhibit systematic errors inpositioning and orbitography. Acomparison with SLR results allows to evidence sucherrors and vice versa. Fordoing that, two different approaches for determiningprecise orbits can beconsidered: one based on global orbit determination,the other on a short-arctechnique used to locally improve a global orbitdetermined by another trackingtechniques, such as DORIS or GPS. We can thusvalidate a global orbit andachieve orbit quality control to a level of2 to 3 centimeters at present and expectto achieve a level of 1 to 2 centimeters inthe near future. Errors induced bystation coordinates or by the gravity field(geographically correlated errors, forexample) can be estimated from SLR tracking data.Colocation experiments withdifferent techniques in the same geodetic siteplay also a key role to ensure preciserelationships between the geodetic referenceframes linked to each technique. Inparticular, the role of the SLR technique is tostrengthen the vertical component(including velocity) of the positioning, whichis crucial for altimetry missions.The role of SLR data in the modelling of the firstterms of the gravity field has finally to be emphasized,which is of primary importance in orbitography,whatever the tracking technique used.Another application of SLR technology is thesatellite altimeter calibration. Examples of past calibrationand future experiments are given, including theaccuracy we can expect from the Jason-1 and EnviSatspace oceanography missions.     

中国电磁监测试验卫星工程研制进展

中国电磁监测试验卫星工程项目于2013年正式批准立项,卫星预计于2017年8月发射入轨,设计在轨运行5年,这是我国地震立体观测体系的第一个专用天基平台. 本文简要介绍了电磁卫星工程设计和卫星工程六大系统的功能与任务,并对卫星系统和应用系统的研制情况予以阐述. 目前,已基本完成卫星平台和有效载荷设计初样阶段的研制,将于2016年6月转入正样研制阶段. 应用系统基础平台的建设工作按计划进行,将在卫星发射前半年开始试运行,以确保电磁监测试验卫星数据的有效应用.      

Gravity variation before Kunlun mountain pass western M_s=8.1 earthquake

AbstractThe relation between the gravity variation features and M_S=8.1 earthquake in Qinghai-Xizang monitoring area isanalyzed preliminarily,by using spatial dynamic variation results of regional gravity field from absolute gravityand relative gravity observation in 1998 and 2000.The results show that:1)M_S=8.1 earthquake in Kulun mountainpass western occurred in the gravity variation high gradient near gravity's high negative variation; 2)The maintectonic deformation and energy accumulation before M_S=8.1 earthquake are distributed at south side of theepicenter;3)The range of gravity’s high negative variation at east of the M_s=8.1 earthquake epicenter relativelycoincides with that rupture region according to field geology investigation; 4)Gravity variation distribution in highnegative value region is just consistent with the second shear strain’s high value region of strain field obtainedfrom GPS observation.     

A revised fit of South America and South Central Africa

A new set of parameters for the total plate tectonic reconstruction of South America and south central Africa is presented: euler pole 46.75°N, 32.65°W; rotation angle 56.40°. This fit is constrained by at least three pre-drift tectonic features crossing from one continent to the other: (1) the geophysically defined eastern and western boundaries of the submarine Jurassic Outeniqua Basin (South Africa) and the Falkland Plateau Basin; (2) the Late Precambrian transcurrent fault and mylonite belts of Pernambuco (Brazil) and Foumban (West Africa); and (3) the Triassic northern tectonic front of the Cape Fold Belt and the major morphological feature on the Falkland Plateau with which it is closely lined up. Isotopic ages of Falkland Plateau gneisses correspond to Cape Pluton and Cape Fold Belt ages, suggesting their palaeoposition was within the realm of the Cape Fold Belt.In addition, the bathymetrically and geophysically defined northeastern apex of the Falkland Plateau fits into the re-entrant angle defined on the South African margin by the steep southeast-facing sheared Agulhas margin and the southern face of the Tugela Cone. Simultaneously known Precambrian outcrops in northeastern Brazil and in the Gulf of Benin area of West Africa are juxtaposed rather than overlapped. Reconstructions producing a closer fit of these cratonic areas are considered untenable.     

GRACE和地面重力测量监测到的中国大陆长期重力变化

自2002年以来,GRACE卫星探测计划可提供高精度的时变地球重力场,用以探测地球系统的物质分布.自1998年中国大陆重力监测网建立以来,利用FG5绝对重力仪和LCR-G型相对重力仪每2年对该网进行重复测量获取重力场时变信息.基于此,本文利用GRACE和地面重力测量获得了中国大陆重力场的长期年变率,利用位错理论根据USGS发布的断层模型计算了2008年汶川M_s8.0级地震的同震重力变化并进行了300 km高斯滤波.GRACE卫星重力和地面重力结果均表明华北地区地下水流失严重,在绝对重力基准站上,GRACE卫星重力与绝对重力变化率较为一致,汶川区域的地面重力变化结果可视为大地震前兆信息.     

What Can be Expected from the GRACE-FO Laser Ranging Interferometer for Earth Science Applications?

The primary objective of the gravity recovery and climate experiment follow-on (GRACE-FO) satellite mission, due for launch in August 2017, is to continue the GRACE time series of global monthly gravity field models. For this, evolved versions of the GRACE microwave instrument, GPS receiver, and accelerometer will be used. A secondary objective is to demonstrate the effectiveness of a laser ranging interferometer (LRI) in improving the satellite-to-satellite tracking measurement performance. In order to investigate the expected enhancement for Earth science applications, we have performed a full-scale simulation over the nominal mission lifetime of 5 years using a realistic orbit scenario and error assumptions both for instrument and background model errors. Unfiltered differences between the synthetic input and the finally recovered time-variable monthly gravity models show notable improvements with the LRI, on a global scale, of the order of 23 %. The gain is realized for wavelengths smaller than 240 km in case of Gaussian filtering but decreases to just a few percent when anisotropic filtering is applied. This is also confirmed for some typical regional Earth science applications which show randomly distributed patterns of small improvements but also degradations when using DDK4-filtered LRI-based models. Analysis of applied error models indicates that accelerometer noise followed by ocean tide and non-tidal mass variation errors are the main contributors to the overall GRACE-FO gravity model error. Improvements in these fields are therefore necessary, besides optimized constellations, to make use of the increased LRI accuracy and to significantly improve gravity field models from next-generation gravity missions.     

基于GRACE卫星重力数据确定地球重力场模型WHU-GM-05

基于卫星轨道运动的能量积分方程,可导出利用卫星跟踪卫星数据求解地球重力场的实用公式.本文在Jekeli给出的公式基础上导出了基于能量守恒方程利用两颗低-低卫星跟踪的扰动位差求解重力位系数的严密关系式.基于两颗GRACE卫星的观测数据,采用本文导出的严密能量积分方法求解得到120阶的GRACE地球重力场模型,命名为WHU-GM-05;将WHU-GM-05模型与国际上同类重力场模型EIGEN-GRACE系列和GGM02S分别在阶方差和大地水准面高等方面作了比较,并与美国和中国的部分地区GPS水准观测值进行了精度分析.结果表明基于本文推导的严密双星能量守恒方程得到的WHU-GM-05重力场模型精度与国际上同类重力场模型的精度相当.     

Contributions of terrestrial and GRACE data to the study of the secular geoid changes in North America

This paper tests and discusses different statistical methods for modelling secular rates of change of the geoid in North America. In particular, we use the method of principal component/empirical orthogonal functions (PC/EOF) analysis to model the geoid rates from Gravity Recovery and Climate Experiment (GRACE) satellite data. As demonstrated, the PC/EOF analysis is useful for studying the contributions from different signals (mainly residual hydrology signals and leakage effects) to the GRACE-derived geoid rates. The PC/EOF analysis leads to smaller geoid rates compared to the conventional least-squares fitting of a trend and annual and semi-annual cycles to the time series of the spherical harmonic coefficients. This is because we filter out particular spatiotemporal modes of the regional geoid changes.We apply the method of least-squares collocation with parameters to combine terrestrial data (GPS vertical velocities from the Canadian Base Network and terrestrial gravity rates from the Canadian Gravity Standardization Net) with the GRACE-derived vertical motion to obtain again the geoid rates. The combined model has a peak geoid rate of 1.4 mm/year in the southeastern area of Hudson Bay contrary to the GRACE-derived geoid rates that show a large peak of 1.6–1.7 mm/year west of Hudson Bay. We demonstrate that the terrestrial data, which have a longer time span than the GRACE data, are important for constraining the GRACE-derived secular signal in the areas that are well sampled by the data.