Review of the Research Progress on Static Earth Gravity Field and Vertical Datum in China during 2019—2023

The contribution presents the representative research progress on global static gravity field modeling,regional geoid/quasigeoid determination,vertical datum study,as well as the theory,algorithm and software for gravity field study in China from 2019 to 2023,which are the highlights of the chapter 6“Progress in Earth Gravity Model and Vertical Datum”in the“2019—2023 China National Report on Geodesy”that submitted to the International Association of Geodesy(IAG).In addition,suggestions are proposed to promote the research in the fields of earth gravity field,geoid/quasigeoid and vertical datumin China according to trends of international geodesy and related disciplines.

The Method and Practice of Constructing 3D Geological Model from Coalfield Exploration 2D Maps

3D geological modeling is an inevitable choice for coal exploration to adapt to the transformation of coal mining for green, fine, transparent and Intelligent mining. In the traditional Coalfield exploration geological reports, the spatial expression form for the coal seams and their surrounding rocks are 2D maps. These 2D maps are excellent data sources for constructing 3D geological models of coal field exploration areas. How to construct 3D models from these 2D maps has been studying in coal exploration industry for a long time, and still no breakthrough has been achieved so far. This paper discusses the principle, method and software design idea of constructing 3D geological model of an exploration area with 2D maps made by AutoCAD/MapGIS. At first, the paper analyzes 3D geological surface expression mode in 3D geological modeling software. It is pointed out that although contour method has unique advantages in coal field exploration, TIN (Triangular Irregular Network) is still the standard configuration of 3D modeling software for coal field. Then, the paper discusses the method of 2D line features obtaining elevation and upgrading 2D curve to 3D curve. Next, the method of semi-automatic partition is introduced to build the boundary ring of the surface patch, that is, the user clicks and selects the line feature to build the outer boundary ring of the surface patch. Then, Auto-process method for fault line inside of the outer boundary ring is discussed, it including construction of fault ring, determining fault ring being normal fault ring or reverse fault ring and an algorithm of dealing with normal fault ring. An algorithm of dealing with reverse fault ring is discussed detailly, the method of expanding reverse fault ring and dividing the duplicate area in reverse fault into two portions is introduced. The paper also discusses the method of extraction ridge line/valley line, the construction of fault plane, the construction of stratum and coal body. The above ideas and methods have been initially implemented in the “3D modeling platform for coal field exploration” software, and applied to the 3D modeling practice of data from several coal field exploration areas in Ningxia, Shanxi, Qinghai, etc.

Realization of an Optimal Dynamic Geodetic Reference Frame in China:Methodology and Applications

China Geodetic Coordinate System 2000(CGCS2000)has been used for several years as a formal published reference frame.The coordinates of all global navigation satellite system(GNSS)stations in China need to be corrected to align with the CGCS2000 frame.Different strategies can be adopted for the realization of an optimal reference frame.However,different strategies lead to different results,with differences as great as several decimeters when GNSS station coordinates are transformed into CGCS2000-defined coordinates.The two common methods for the coordinate correction of a GNSS station are quasi-stable adjustment under CGCS2000 and plate movement correction,and the differences between their results can be greater than 10 cm.In this study,a statistic method called"supervised clustering"is applied to the selection of GNSS reference stations;a new scheme named"partition spacing"for the grouping of all processed GNSS stations is proposed;and the plate movement correction method is used to correct the coordinates of all GNSS stations from the GNSS epoch to the CGCS2000 epoch.The results from the new partitioning method were found to be significantly better than those from the conventional station-blocking approach.When coordinates from the stations without grouping were used as the standard,the accuracy of all the three-dimensional coordinate components from the new partitioning method was better than 2 mm.The root mean squares(RMSs)of the velocities in the x,y,and z directions resulting from the supervised clustering method were 0.19,0.45,and 0.32 mm∙a1,respectively,which were much smaller than the values of 0.92,0.72,and 0.97 mm∙a1 that resulted from the conventional approach.In addition,singular spectrum analysis(SSA)was used to model and predict the position nonlinear movements.The modeling accuracies of SSA were better than 3,2,and 5 mm in the east(E),north(N),and up(U)directions,respectively;and its prediction accuracies were better than 5 mm and 1 cm for the horizontal and vertical domains,respectively.

Progresses and Prospects of Marine Geodetic Datum and Marine Navigation in China

Territorial water is a significant part of national sovereignty of China,thus the infrastructures of national space datum and location services should cover the sea areas except for the land areas.China has established relatively complete geodetic coordinate frame in land areas over the past decades,including the lastest developed China Geodetic Coordinate System 2000(CGCS 2000)with its reference frame and the national gravity datum 2000.However,the currently used geodetic infrastructures have not well covered the sea areas of China.The marine geodetic datum and marine navigation technologies need to be further developed and extended to satisfy the national demands of marine environment and resources detection,scientific investigation as well as marine economy development in new era of China.This paper mainly reviews the development and the progress of Chinese marine geodetic datum and marine navigation,analyses related key technologies in establishing our national marine geodetic datum.Some current trends and future directions for independently developing our national marine geodetic datum and marine navigation technologies are discussed.

Adjustment Model and Colored Noise Compensation of Continuous Observation System

The affection caused by the colored noises should be taken into account to the adjustment model.As useful signals,these colored noises should be accurately identified and extracted by Fourier analysis.A continuous adjustment model is introduced with respect to the colored noises,and then it can be generalized from the finite space to the infinite space so called as Hilbert space.This extension is to provide a new technique to perform the continuous observational system design,Fourier analysis as well as the parameter estimation.It shows that the Gramer’s determinant provides maximization criteria in the system optimization design as well as a rule in diagnosing the adjustment model.Related with the definition of the integral,the least squares solution of the continuous adjustment model becomes the limit of the traditional least squares solution in finite space.Moreover,the influence caused by the colored noises is systematic,but it can be eliminated or compensated by optimally designing the observational system.

Maintenance of Millimeter-level Geodetic Reference Framework

The high-precision terrestrial reference frame,as the spatial benchmark for geodesy,is an important national infrastructure.However,due to the influence of nonlinear factors related to geophysical phenomena,the overall maintenance accuracy of the ITRF framework is still at the centimeter level.Therefore,accurately characterizing the true trajectories of linear motion,nonlinear motion,and geocentric motion of the reference station is the key to achieve the construction and maintenance technology of a millimeter level terrestrial reference framework.Based on long-term global and regional GNSS observation data,more Chinese geodesy scientists devoted much efforts to the maintenance of millimeter-level geodetic reference framework.The main contributions of this work included the followings:①Dynamic maintenance of millimeter-level terrestrial reference frame;②Research progress on the method of maintenance of regional reference frame based on GNSS;③The progress of CGCS2000 frame maintenance in millimeter level accuracy;④Reprocessing and reanalysis of two-decade GNSS observation in continental China;⑤Research on current GNSS velocity field model and deformation in Chinese mainland;⑥The preliminary realization and evaluation of CTRF2020.

Precise Positioning Method for Seafloor Geodetic Stations Based on the Temporal Variation of Sound Speed Structure

At present,GNSS-Acoustic(GNSS-A)combined technology is widely used in positioning for seafloor geodetic stations.Based on Sound Velocity Profiles(SVPs)data,the equal gradient acoustic ray-tracing method is applied in high-precision position inversion.However,because of the discreteness of the SVPs used in the forementioned method,it ignores the continuous variation of sound velocity structure in time domain,which worsens the positioning accuracy.In this study,the time-domain variation of Sound Speed Structure(SSS)has been considered,and the cubic B-spline function is applied to characterize the perturbed sound velocity.Based on the ray-tracing theory,an inversion model of“stepwise iteration&progressive corrections”for both positioning and sound speed information is proposed,which conducts the gradual correction of seafloor geodetic station coordinates and disturbed sound velocity.The practical data was used to test the effectiveness of our method.The results show that the Root Mean Square(RMS)errors of the residual values of the traditional methods without sound velocity correction,based on quadratic polynomial correction and based on cubic B-spline function correction are 1.43 ms,0.44 ms and 0.21 ms,respectively.The inversion model with sound velocity correction can effectively eliminate the systematic error caused by the change of SSS,and significantly improve the positioning accuracy of the seafloor geodetic stations.

General Structure Physics of an Aerial Remote Sensing Platform and Its Systemic Accuracy Criterion

Accuracy is a key factor in high-resolution remote sensing and photogrammetry. The factors that affect accuracy are imaging system errors and data processing errors. Due to the complexity of aerial camera errors, this paper focuses on the design of digital aerial camera systems and the means to reduce system error and data processing inefficiencies. There are many kinds of digital aerial camera systems at present;however, these systems lack a unified physical model, which ultimately leads to more complicated designs and multi-camera modes. Such a system is complex and costly, as it is easily affected by factors such as vibration and temperature. Thus, the installed accuracy can only reach the millimeter level. Here, we describe a unified physical structure for a digital aerial camera that imitates an out-of-field multi-charge-coupled device (CCD), an in-field multi-CCD, and once-imaging and twice-imaging digital camera systems. This model is referred to as the variable baseline-height ratio spatiotemporal model. The variable ratio allows the opto-mechanical spatial parameters to be linked with height accuracy, thus providing a connection to the surface elevation. The twice-imaging digital camera prototype system and the wideband limb imaging spectrometer provide a transformation prototype from the current multi-rigid once-imaging aerial camera to a single rigid structure. Thus, our research lays a theoretical foundation and prototype references for the construction and industrialization of digital aerial systems.

Recent Advances in Marine Geodesy of China

The ocean accounts for approximately 71%of the total area of the Earth.Whether it is studying the shape of the Earth itself through geodesy or the future development of earth system science,strengthening the construction of ocean geodesy disciplines and innovating ocean geodetic observation technologies have evident theoretical and practical significance.In recent years,the discipline of ocean geodesy in China has been continuously developing and growing,and notable breakthroughs have been made in ocean satellite geodesy and seafloor geodetic observation technology.Research on ocean geodetic observation models and algorithms has also made great progress.

Classification of CBERS-2 Imagery with Fuzzy ARTMAP Classifier

A fuzzy ARTMAP classifier is adopted for a classification experiment of CBERS-2 imagery. The fundamental theory and processing about the algorithm are first introduced, followed with a land-use classification experiment in Shihezi County on CBERS-2 high resolution imagery. Three classifiers are compared: maximum likelihood classifier (MLC), error back propagation (BP) classifier, and fuzzy ARTMAP classifier. The comparison shows comparably better results for the fuzzy ARTMAP classifier, with overall classification accuracy of 9.9% and 4.6% higher than that of MLC and BP. The results also prove that the fuzzy ARTMAP classifier has better discernment in identifying bare soil on CBERS-2 imagery.

Treatment of discontinuities inside Earth models:Effects on computed coseismic deformations

In this paper,we study how coseismic deformations calculated in 1066 Earth models are affected by how the models treat Earth discontinuities.From the results of applying models 1066A(continuous)and 1066B(discontinuous),we find that the difference in Love numbers of strike-slip and horizontal tensile sources are bigger than dip-slip and vertical tensile sources.Taken collectively,discontinuities have major effects on Green’s functions of four independent sources.For the near-field coseismic deformations of the 2013 Okhotsk earthquake(Mw 8.3),the overall differences between theoretical calculations in vertical displacement,geoid,and gravity changes caused by discontinuities are 10.52 percent,9.07 percent and 6.19 percent,with RMS errors of 0.624 mm,0.029 mm,and 0.063μGal,respectively.The difference in far-field displacements is small,compared with GPS data,and we can neglect this effect.For the shallow earthquake,2011 Tohoku-Oki earthquake(Mw 9.0),the differences in near-field displacements are 0.030 m(N-S),0.093 m(E-W),and 0.025 m(up-down)in our study area with the ARIA slip model,which gives results closer to GPS data than those from the USGS model.The difference in vertical displacements and gravity changes on the Earth’s surface caused by discontinuities are larger than 10 percent.The difference in the theoretical gravity changes at spatially fixed points truncated to degrees 60,as required by GRACE data,is 0.0016μGal and the discrepancy is 11 percent,with the theoretical spatial gravity changes from 1066B closer to observations than from 1066A.The results show that an Earth model with discontinuities in the medium has a large effect on the calculated coseismic deformations.

On the application of DInSAR to deformation monitoring in desert areas

DInSAR 技术被用于监视荒芜的高度变化学习沙暴。作为测试区域， Hunshandake 沙的土地是在北京的沙暴的主要来源之一。为了学习沙暴来源和它的影响，，一双 EnviSat ASAR 2005 年 10 月 11 日，和 2004 年 10 月 26 日想象，根据六 ERS-2 的分析被处理， EnviSat 雷达想象。在图象配置，扁平的地球效果修正，过滤的数据，阶段 unwrapping，和编码 geo 以后，在 Hunshandake 沙漠上的一个变丑模型被造。根据结果，高度在大多数区域减少了并且在一些区域增加了，它基本上与一致，强壮的沙暴在 2005 的春天出现在北京。结果证明 DInSAR 在荒芜的表面变丑监视有一个重要角色。

Earth rotation parameter and variation during 2005—2010 solved with LAGEOS SLR data

Time series of Earth rotation parameters were estimated from range data measured by the satellite laser ranging technique to the Laser Geodynamics Satellites（LAGEOS）-1/2 through 2005 to 2010 using the dynamic method. Compared with Earth orientation parameter（EOP）C04, released by the International Earth Rotation and Reference Systems Service, the root mean square errors for the measured X and Y of polar motion（PM） and length of day（LOD）were 0.24 and 0.25 milliarcseconds（mas）, and 0.068 milliseconds（ms）, respectively.Compared with ILRSA EOP, the X and Y of PM and LOD were 0.27 and 0.30 mas, and 0.054 ms, respectively. The time series were analyzed using the wavelet transformation and least squares methods. Wavelet analysis showed obvious seasonal and interannual variations of LOD, and both annual and Chandler variations of PM; however, the annual variation could not be distinguished from the Chandler variation because the two frequencies were very close. The trends and periodic variations of LOD and PM were obtained in the least squares sense, and PM showed semi-annual, annual, and Chandler periods.Semi-annual, annual, and quasi-biennial cycles for LOD were also detected. The trend rates of PM in the X and Y directions were 3.17 and 1.60 mas per year, respectively, and the North Pole moved to 26.8E relative to the crust during 2005—2010. The trend rate of the LOD change was 0.028 ms per year.

Accuracy Analysis of Low Altitude Photogrammetry with Wide-angle Camera

Firstly,the relationship between the accuracy of low altitude aerial photogrammetry and the field angle of camera is made by a quantitative analysis from the theory.The conclusion that the low altitude photogrammetry should use wide-angle camera as much as possible is done.Then,the limitation of the single lens camera to expand field angle and the combined wide-angle camera existing on the market not suitable for light load of low altitude UAV(Unmanned Aerial Vehicle)due to excessive weight are pointed out.The characteristics of combined wide-angle low altitude light camera with self-calibration and self-stabilization developed by the author are described,especially the principle of self-calibration for the combination of static error and dynamic error.Based on the practice of large scale mapping,a technical procedure in aerial photography by taking with wide-angle camera and large overlap simultaneously for improving the accuracy of low altitude photogrammetry is proposed.The typical engineering produced data is used to verity the above theoretical analysis.A technical route for increasing accuracy of low altitude photogrammetry with combined wide-angle camera is expounded.

Assimilating ASAR Data for Estimating Soil Moisture Profile Using an Ensemble Kalman Filter

Active microwave remote sensing data were used to calculate the near-surface soil moisture in the vegetated areas.In this study,Advanced Synthetic Aperture Radar(ASAR)observations of surface soil moisture content were used in a data assimilation framework to improve the estimation of the soil moisture profile at the middle reaches of the Heihe River Basin,Northwest China.A one-dimensional soil moisture assimilation system based on the ensemble Kalman filter(EnKF),the forward radiative transfer model,crop model,and the Distributed Hydrology-Soil-Vegetation Model(DHSVM)was developed.The crop model,as a semi-empirical model,was used to estimate the surface backscattering of vegetated areas.The DHSVM is a distributed hydrology-vegetation model that explicitly represents the effects of topography and vegetation on water fluxes through the landscape.Numerical experiments were conducted to assimilate the ASAR data into the DHSVM and in situ soil moisture at the middle reaches of the Heihe River Basin from June20 to July 15,2008.The results indicated that EnKF is effective for assimilating ASAR observations into the hydrological model.Compared with the simulation and in situ observations,the assimilated results were significantly improved in the surface layer and root layer,and the soil moisture varied slightly in the deep layer.Additionally,EnKF is an efficient approach to handle the strongly nonlinear problem which is practical and effective for soil moisture estimation by assimilation of remote sensing data.Moreover,to improve the assimilation results,further studies on obtaining more reliable forcing data and model parameters and increasing the efficiency and accuracy of the remote sensing observations are needed,also improving estimation accuracy of model operator is important.

Assessment Model of Atmosphere Transmitting Influence on High-resolution Airborne SAR Stereo Positioning

The influence derived from atmosphere transmitting of radar wave, in the application of high-resolution airborne Synthetic Aperture Radar (SAR) stereo positioning, may produce some phase errors, and eventually be introduced into positioning model. This paper described the principle of airborne SAR stereo positioning and the error sources of stereo positioning accuracy that arose from atmosphere transmitting, established a corresponding assess- ment model of atmosphere transmitting influence, and testified the model and the assessment principle taking the 1-m resolution airborne SAR images of Zigong City, Sichuan Province in China, as the test dataset. The test result has proved that the assessment model is reliable and reasonable. And, it has shown that the phase error arisen from time delay is the main error source during the atmosphere transmitting, which has much more influences on cross-track di- rection and introduces a stereo positioning error of about eight meters, but less on the along-track direction.

Potential application of particle based simulations in reservoir security management

In order to model the movement progress in case of risks such as dam collapse and coastal inundation, particle-based simulation methods, including the discrete-element method and smoothed particle hydrodynamics, which have specific advantages in modeling complex three-dimensional environmental fluid and particulate flows, are adopted as an effective way to illustrate environmental applications possibly happening in the real world. The theory of these methods and their relative advantages compared with traditional methods are discussed. Examples of 3-D flows on realistic topography including the flooding of a river valley as a result of a dam collapse and coastal inundation by a tsunami are introduced. Issues related to validation and quality data availability are also discussed. The results show that the simulations provide a valuable insight in a given situation for the security management of reservoir dams. Validation can only be performed where both the initial and final states can be very well characterized.

A Spatialization-based Method for Checking and Updating Metadata

In this paper the application of spatialization technology on metadata quality check and updating was dis-cussed. A new method based on spatialization was proposed for checking and updating metadata to overcome the defi-ciency of text based methods with the powerful functions of spatial query and analysis provided by GIS software. Thismethod employs the technology of spatialization to transform metadata into a coordinate space and the functions ofspatial analysis in GIS to check and update spatial metadata in a visual environment. The basic principle and technicalflow of this method were explained in detail, and an example of implementation using ArcMap of GIS software wasillustrated with a metadata set of digital raster maps. The result shows the new method with the support of interactionof graph and text is much more intuitive and convenient than the ordinary text based method, and can fully utilize thefunctions of GIS spatial query and analysis with more accuracy and efficiency.

Research of Directory Service Model for Geographic Information

With the deepening application of geographic information web services, as a significant part, geographic information directory service sets up a service transaction platform between the service providers and the service requestors, providing one-stop registration, access and management functions of geographic information web services, so it is very important to carry out the research on geographic information directory service. According to the domestic and overseas metadata standards, the service expression-oriented metadata information model for geographic information web service is established, and the discovery mechanism dominated by geographic information service registration and active discovery is realized accordingly. After that, overall design of directory service model is made, and finally the research results are integrated into the geographic information service software NewMapServer and deployment test is made under the Google Cloud Environment Google App Engine, providing a certain reference for sharing the Smart City construction achievements.

Thermodynamic driving mechanisms for the formation of global precipitation extremes and ecohydrological effects

Global warming has altered the thermodynamic and dynamic environments of climate systems,affecting the biogeochemical processes between the geosphere and atmosphere,which has significant impacts on precipitation extremes and the terrestrial carbon budget of ecosystems.Existing studies have reported a hook structure for precipitation extreme-temperature relationships but have rarely examined the underlying physical mechanisms.Previous studies have also failed to quantify the impact of precipitation on ecosystem productivity,hindering the assessment of future extreme climatic hazards and potential ecosystem risks.To reveal the thermodynamic driving mechanisms for the formation of global precipitation extremes and ecohydrological effects,this study utilizes over ten multisource datasets(i.e.,satellite,reanalysis,climate model,land surface model,machine learning reconstruction,and flux tower measurements).We first assess the response of water-heat-carbon flux to precipitation extremes and explain the underlying physical mechanisms behind the hook structures in terms of atmospheric thermodynamics and dynamics.Based on outputs from five global climate models(GCMs)under ISIMIP3b,we project future changes in the hook structures as well as their impacts on precipitation extremes.Finally,we discuss the impact of precipitation on the terrestrial carbon budget by using outputs from the CLM4.5 model.The results show that precipitation extremes are usually accompanied by strong exchanges of water and heat and demonstrate a nonlinear relationship between precipitation and ecosystem productivity.The intensity(duration)of extreme precipitation is intensifying(decreasing)over most areas of the globe,whereas three-dimensional precipitation events are becoming more concentrated.Atmospheric dynamics play a key role in shaping the hook structure.The structure is not stable;it shifts under climate change and is projected to result in a 10–40%intensification in precipitation by the end of this century.Moderate levels of precipitation contribute to carbon assimilation in ecosystems,and the response of the carbon budget to precipitation is relatively stable under climate change.