全球农情遥感速报系统20年

面向国家粮食安全的重大战略需求,1998年中国科学院建立了"中国农情遥感速报系统"(CropWatch),持续运行20年后,现已发展成为"全球农情遥感速报系统"(CropWatch)。本文重点论述了2013年建立参与式全球农情遥感监测云平台(CropWatch-Cloud)以来,所采用的农情监测体系、可定制的农情监测云平台理念以及CropWatch-Cloud在国内外的应用推广情况,介绍了技术方法与农情信息服务方式的创新与进步带来的国际影响力的提升。系统总结了全球农情遥感速报系统发展的农情监测指标、农情预警能力、作物长势综合监测方法以及众源数据支持的作物面积监测方法,论文进一步阐述了CropWatch未来的发展方向,借助众源地理信息、大数据技术等的发展,打通从地块—村—镇—县—市—省—国家—全球的体系化全链条监测,满足从农户到政府决策部门对农情信息的差异化需求。     

A model of present-day tectonic plate motions from 12 years of DORIS measurements

In the frame of the International DORIS Service (IDS), the Laboratoire d’Etudes en Géophysique et Océanographie Spatiales (LEGOS)/Collecte Localisation Satellites (CLS) Analysis Center (LCA) processes DORIS measurements from the SPOT, TOPEX/Poseidon and Envisat satellites and provides weekly station coordinates of the whole network to the IDS. Based on DORIS measurements, the horizontal and vertical velocities of 57 DORIS sites are computed. The 3D positions and velocities of the stations with linear motion are estimated simultaneously from the 12-year (1993–2004) combined normal equation matrix. We include 35 DORIS sites assumed to be located in the stable zones of 9 tectonic plates. For the motion of these plates, we propose a model (LCAVEL-1) of angular velocities in the ITRF2000 reference frame. Based on external comparison with the most recent global plate models (PB2002, REVEL, GSRM-1 and APKIM2000) and on internal analysis, we estimate an average velocity error of the DORIS solution of less than 3 mm/year. The LCAVEL-1 model presents new insights of the Somalia/Nubia pair of plates, as the DORIS technique has the advantage of having a few stations located on those two plates. We also computed (and provide in this article) the horizontal motion of the sites located close to plate boundaries or in the deformation zones defined in contemporary models. These computations could be used in further analysis for these particular regions of the Earth not moving as rigid plates.     

Defining a DORIS core network for Jason-1 precise orbit determination based on ITRF2000: methods and realization

In view of the future adoption of the new precise orbit determination (POD) standards for the TOPEX/Poseidon and Jason-1 satellites, we propose a method to evaluate terrestrial reference frames for POD. We applied this method to the ITRF2000 realization of the DORIS network using local geodetic ties, plate motion models, the recent DORIS IGN04D02 cumulative solution and DORIS weekly time-series of coordinates. We propose to adopt a selection of the ITRF2000 realization based on specific criteria that we define here, and to extend it with ground stations for which we propose new coordinates and velocities. Only 13 out of 131 stations were considered to be inappropriate for POD activities. The result is a robust and well-distributed DORIS core network of 118 stations (DPOD2000) suitable for POD during the 1993–2008 period considered here.     

Stochastic models in the DORIS position time series: estimates for IDS contribution to ITRF2014

This paper focuses on the investigation of the deterministic and stochastic parts of the Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) weekly time series aligned to the newest release of ITRF2014. A set of 90 stations was divided into three groups depending on when the data were collected at an individual station. To reliably describe the DORIS time series, we employed a mathematical model that included the long-term nonlinear signal, linear trend, seasonal oscillations and a stochastic part, all being estimated with maximum likelihood estimation. We proved that the values of the parameters delivered for DORIS data are strictly correlated with the time span of the observations. The quality of the most recent data has significantly improved. Not only did the seasonal amplitudes decrease over the years, but also, and most importantly, the noise level and its type changed significantly. Among several tested models, the power-law process may be chosen as the preferred one for most of the DORIS data. Moreover, the preferred noise model has changed through the years from an autoregressive process to pure power-law noise with few stations characterised by a positive spectral index. For the latest observations, the medians of the velocity errors were equal to 0.3, 0.3 and 0.4 mm/year, respectively, for the North, East and Up components. In the best cases, a velocity uncertainty of DORIS sites of 0.1 mm/year is achievable when the appropriate coloured noise model is taken into consideration.     

多普勒系统对历元地球参考框架的影响研究

针对历元地球参考框架在确定站点高频非线性运动和季节性变化方面具有国际地球参考框架不具备的优势问题,该文通过比较两种技术间组合策略建立的多源融合历元地球参考框架的精度,研究DORIS对于多源融合历元地球参考框架的影响。通过基于坐标的法方程叠加方法进行技术内法方程叠加和技术间组合,并利用并置站条件联系不同技术的法方程,从而建立多源融合历元地球参考框架。结果表明,DORIS的引入并不会对多源融合历元地球参考框架的基准定义或其他3种技术站点的点位精度产生较大的影响,同时能够体现4种技术中精度较差的DORIS对于ETRF的精度及稳定性的影响。     

Using climate-FVS to project landscape-level forest carbon stores for 100 years from field and LiDAR measures of initial conditions

Background

Forest resources supply a wide range of environmental services like mitigation of increasing levels of atmospheric carbon dioxide (CO2). As climate is changing, forest managers have added pressure to obtain forest resources by following stand management alternatives that are biologically sustainable and economically profitable. The goal of this study is to project the effect of typical forest management actions on forest C levels, given a changing climate, in the Moscow Mountain area of north-central Idaho, USA. Harvest and prescribed fire management treatments followed by plantings of one of four regionally important commercial tree species were simulated, using the climate-sensitive version of the Forest Vegetation Simulator, to estimate the biomass of four different planted species and their C sequestration response to three climate change scenarios.

Results

Results show that anticipated climate change induces a substantial decrease in C sequestration potential regardless of which of the four tree species tested are planted. It was also found that Pinus monticola has the highest capacity to sequester C by 2110, followed by Pinus ponderosa, then Pseudotsuga menziesii, and lastly Larix occidentalis.

Conclusions

Variability in the growth responses to climate change exhibited by the four planted species considered in this study points to the importance to forest managers of considering how well adapted seedlings may be to predicted climate change, before the seedlings are planted, and particularly if maximizing C sequestration is the management goal.     

Estimating the noise in space-geodetic positioning: the case of DORIS

The noise spectrum in DORIS ground- station motion is investigated by means of the Allan variance method applied to the decomposition of the 3D signal into its principal components in the time domain. Sets of weekly position time-series from 1994 to 2005 derived by three IDS Analysis Centres (IGN-JPL, INASAN, and LEGOS-CLS) for 119 stations at 69 sites are considered. The observing satellites are SPOT-2, SPOT-3, SPOT-4, and SPOT-5, TOPEX/Poseidon, and ENVISAT. Annual and semi-annual perturbations, as well as the 117.3-day term associated with the TOPEX/Poseidon orbit, are found at most stations. Their amplitudes reach up to 19.3, 23.7, and 13.3 mm, respectively, for the three analysis centres (ACs). When corrected for these components and a linear drift, the time-series dominantly show white noise (WN) at the 10–45mm level the noise level is the highest in the East direction, probably in connection with the high orbit inclinations. The noise level is minimum for the high latitude stations, mostly and intensively observed by the SPOT satellites, and the determination of the noise type is unclear; longer observation spans would be needed to decide between interannual variations and flicker noise. The improvement in positioning due to the DORIS constellation extension from three to five satellites in 2002, and the network rejuvenation program initiated in 2000, results in a decrease of the noise level by a factor of 1.7 in a WN context, both before and after the changes. One example of the benefit of studying the signal in the time eigenspace domain is the detection of anomalously large WN in the East direction for station HBKB (Hartebeesthoek, Africa) that masks the above-mentioned improvement. Studying the projection on the local frame of the second and third time-eigenspace components, a noise excess is detected in the North direction for some of the ACs. Station stability derived from our time-series analysis confirms, in general, the expected performance based on the careful technical review of the station components (antenna, pillar, etc.). The respective merits of our noise qualification method, based on direct time-series analysis in the time-eigenspace domain without any a priori statistical model, in comparison with other methods, such as the selection of a mixed-noise model by maximum likelihood estimation, are discussed.     

Future global SLR network evolution and its impact on the terrestrial reference frame

Satellite laser ranging (SLR) is an important technique that contributes to the determination of terrestrial geodetic reference frames, especially to the realization of the origin and the scale of global networks. One of the major limiting factors of SLR-derived reference frame realizations is the datum accuracy which significantly suffers from the current global SLR station distribution. In this paper, the impact of a potential future development of the SLR network on the estimated datum parameters is investigated. The current status of the SLR network is compared to a simulated potential future network featuring additional stations improving the global network geometry. In addition, possible technical advancements resulting in a higher amount of observations are taken into account as well. As a result, we find that the network improvement causes a decrease in the scatter of the network translation parameters of up to 24%, and up to 20% for the scale, whereas the technological improvement causes a reduction in the scatter of up to 27% for the translations and up to 49% for the scale. The Earth orientation parameters benefit by up to 15% from both effects.     

A corrective model for Jason-1 DORIS Doppler data in relation to the South Atlantic Anomaly

The DORIS Doppler measurements collected by Jason-1 are abnormally perturbed by the influence of the South Atlantic Anomaly (SAA). The DORIS ultra-stable oscillators on-board Jason-1 are not as stable as they should be; their frequency is sensitive both to the irradiation rate and to the total irradiation encountered in orbit. The consequence is that not only are the DORIS measurement residuals higher than they ought to be, but also large systematic positioning errors are introduced for stations located in the vicinity of the SAA. In this paper, we present a method that has been devised to obtain a continuous observation of Jason-1 frequency offsets. This method relies on the precise determination of the station frequency and troposphere parameters via the use of other DORIS satellites. More than 3 years of these observations have then been used to construct a model of response of the oscillators of Jason-1 to the SAA. The sensitivity of the Jason-1 oscillators to the SAA perturbations has evolved over time, multiplied by a factor of four between launch and mid-2004. The corrective performances of the model are discussed in terms of DORIS measurement residuals, precise orbit determination and station positioning. The average DORIS measurement residuals are decreased by more than 7 % using this model. In terms of precise orbit determination, the 3D DORIS-only orbit error decreases from 5 to 4.2 cm, but the DORIS+SLR orbit error is almost unaffected, due to the already good quality of this type of orbit. In terms of station positioning, the model brings down the average 3D mono-satellite monthly network solution discrepancy with the International Terrestrial Reference Frame ITRF2000 from 11.3 to 6.1 cm, and also decreases the scatter about that average from 11.3 to 3.7 cm. The conclusion is that, with this model, it is possible to re-incorporate Jason-1 in the multi-satellite geodetic solutions for the DORIS station network.     

河南省高速交通网通达性时空特征及演化规律

针对河南省高速公路建设对城市通达性影响问题,该文通过定义时间阈值,设计了基于人口加权时间距离倒数的通达性指数,在4个时间断面上,研究河南省高速交通网络中通达性时空特征及演化规律,时间距离倒数有效减小了多路径时间距离对分析结果的数字影响。结果表明,随着河南高速交通网络的拓展,城市通达性水平得到了显著的提升,空间收敛显著,边缘城市的提升始终高于中心区域城市;高速公路通车里程增长率与通达性水平的改善存在非协调现象;城市通达性的空间聚集现象震荡走弱,聚集区域减少;城市通达性空间分布变化规律和城市化人口变化规律逐渐趋于协调。     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

Ionospheric applications of the scintillation and tomography receiver in space (CITRIS) mission when used with the DORIS radio beacon network

The scintillation and tomography receiver in space (CITRIS) instrument will orbit the Earth near 560 km altitude to detect signals from the ground-based array of more than 50 DORIS UHF/S-band radio beacons established at sites around the world by the French Centre National d‘Etudes Spatiales (CNES) and the Institut Géographique National (IGN). The CITRIS receiver is on the US Air Force Space Test Program satellite STPSAT1, which is scheduled for launch in November 2006. CITRIS will record ionospheric total electron content (TEC) and radio scintillations with a unique ground-to-space geometry. The new instrument has been developed to study the ionosphere using data obtained with the UHF and S-band radio transmissions from the DORIS beacons because ionospheric radio scintillations can seriously degrade the performance of many space-geodetic systems, including the DORIS precise satellite orbitography system and GNSS (Global Navigation Satellite Systems). The ionospheric data will be based on radio signals sampled at a rate of 200 Hz by the CITRIS receiver. Numerical models have been used to predict that the DORIS signals measured by CITRIS may have 30 dB fluctuations in amplitude and 30 rad in phase as the satellite flies over kilometer-scale ionospheric structures. The data from the space-based CITRIS receiver will help update and validate theories on the generation and effect of ionospheric irregularities known to influence radio systems. By using simultaneous beacon transmissions from DORIS on the ground and from low-Earth-orbit beacons in space, the concept of reciprocity in a non-bilateral propagation medium like the ionosphere will be tested. Computer simulations are used to predict the magnitude of amplitude and phase scintillations that are expected to be recorded with the CITRIS instrument.     

GPS控制网起算点兼容性分析方法研究与实践

不同等级或分属不同区域的控制点间,以及由不同单位布设的控制点间可能存在较大的误差。若将这些存在较大误差的已知点作为平差计算或坐标转换的约束,其结果必然会扭曲GPS测量的原有精度,特别是当这些点误差较大或含有粗差时,将严重影响GPS成果的可靠性,使高精度的GPS定位成果失去其本来的意义。本文分别探讨了平差结果直接分析法、实测基线比较法、单位权方差假设检验法和附合路线坐标闭合差检验法等对GPS控制网起算数据进行兼容性分析的方法,并采用其中几种对GPS控制网起算数据进行具体分析和计算,剔除了含有粗差的已知点,保证了控制测量成果的精度和可靠性。     

Systematic errors in the Z-geocenter derived using satellite tracking data: a case study from SPOT-4 DORIS data in 1998

Within the scope of the Global Geodetic Observing System, Doppler Orbit Determination and Radiopositioning Integrated by Satellite – as a geodetic technique – can provide precise and continuous monitoring of the geocenter motion related to mass redistribution in the Earth, ocean and atmosphere system. We have reanalyzed 1998 DORIS/SPOT-4 (Satellite pour l’ Observation de la Terre) data that were previously generating inconsistent geocenter positions (?65 cm offset). We show here that this error is due to an incorrect phase center correction provided with the DORIS preprocessed data resulting from a +12 cm offset in the cross-track direction that has been confirmed since. We also conclude that a 1 mm error in the cross-track offset of non-yawing sun-synchronous SPOT satellites will generate a ?6.5 mm error in the derived Z-geocenter. Other non-yawing satellites would also be affected by a similar effect whose amplitude could be easily estimated from the orbit inclination     

跨进新世纪的特种精密工程测量

为了对二十一世纪的特种精密工程测量有一个新概念,在回顾特种精密工程测量发展概况的基础上,给出了如下几点展望：特种精密工程测量将走向自动化、智能化、实时化和系统化;特种精密工程测量将服务于国民经济各部门;特种精密工程测量将与其它科学技术一起发展。     

Use of multitemporal remote sensing data for mapping the Alfios River network changes from 1977 to 2000

Alfios is the biggest river of the Peloponnese and the ninth longest river in Greece. It drains an area of almost 2575 km² in Western Peloponnese and discharges at Kiparissiakos Gulf. Due to its extent, the Alfios basin presents complex physiography and geomorphology.

During the last 50 years, major direct human activities have had an important impact on the river network shape and on the valley floor morphology. The most important human activities were the channel straightening caused by meander cut-offs in 1950, the construction of two dams: the Ladonas dam in 1955 and the Flokas dam in 1968 and the illegal in-stream gravel extraction that started in the 1960s and continues until today. As a result, the Alfios riverbed has presented very rapid vertical and horizontal (lateral) erosion and significant changes to its network shape. The drainage network follows a straighter course, the number of meanders has been reduced and the seventh order branch is deeply incised. Besides the morphological impacts to the river network shape there are also serious economic damages caused by human activity. During the winter of 1999, the Flokas dam bridge was closed for a long period due to damages at its foundation caused by the heavy rain and illegal gravel extraction. Transportation between the villages in the area became extremely difficult and time consuming and the cost of the repairs was estimated at 500 000 euros. In this study, there is an effort to map the network shape transformation for the period 1977–2000 using multitemporal and multisensor satellite images. One Landsat MSS image, three Landsat TM images, two Landsat ETM images, and one Terra ASTER image have been orthorectified and processed in order to cover the specific period. Principal Component Analysis (PCA) method and Geographical Information System (GIS) techniques have been applied to map the changes in the Alfios River channel. The drainage network straightening and the cut-off of five big meanders were detected and mapped. These changes occurred between 1986 and 2000.     

GPS网平差的多元响应-非线性最小二乘模型及其解的最优化混合算法

在控制网中,用GPS测量,可以同时测得一个点的基线向量,在具有起算数据时可获得观测点的三维坐标等多个响应变量值,将这些响应变量包含的信息综合起来,可以得到参数的更精确的估计。本文针对GPS测量数据为多元响应数据的特点,建立了一个GPS网平差的多元响应-非线性最小二乘模型,针对该模型,结合拟牛顿法和信赖域算法建立了一个新非线性优化的混合算法,该算法具有全局收敛性和超线性收敛性。     

近32年塔里木盆地与准噶尔盆地湿地演化遥感分析

以1975年MSS,2000年ETM和2007年CBERS 02三期遥感数据为基础信息源,对塔里木盆地和准噶尔盆地湿地演化特征及主要驱动力进行了分析.结果表明:①32 a间,准噶尔盆地湿地面积增加了254.94 km2,塔里木盆地湿地面积减少了2 739.88 km2,沼泽湿地的严重退化引起了塔里木盆地湿地的锐减;②塔里木盆地湿地演化可分为破碎化和消亡2个发展阶段,湿地主导类型已由沼泽湿地转变为河流湿地,景观多样性指数与均匀度指数降低,说明湿地多样性减小;准噶尔盆地湿地景观多样性指数和均匀度指数持续增加,而优势度指数持续减小,说明湿地景观结构趋于均匀化;③塔里木盆地和准噶尔盆地年降水量和年均气温呈逐年上升趋势,增加量分别为0.59 mm/a和0.042℃/a,1.69 mm/a和0.047℃/a;准噶尔盆地湿地演化与气候变化密切相关;④1975-2000年和2000-2007年间,塔里木盆地因“围湿造田”直接丧失的湿地面积高达1 035.06 km2和1 030.30 km2,分别占2个时期湿地减少面积的81.74％和69.92％.人类活动是引起塔里木盆地湿地面积锐减的主要原因.     

近百年来洞庭湖区垸内湖泊时空演变分析

研究洞庭湖区垸内湖泊时空演变情况,本文利用1930年代以来10个时段的历史地形图和航天航空遥感数据,采用遥感解译、数据统计分析与历史对比方法,分析民国中期（1930年）以来各时段洞庭湖区垸内湖泊的空间分布特征、面积时序变化,利用2018年度基础性地理国情监测数据中内湖数据进行精度检验,精度达98.62%。结果表明：洞庭湖区内湖数量与面积分别从1938年的640个、994.30 km²增加到2018年的1057个、1578.33 km²,其变化过程并非呈持续增加趋势,而是经历了增加、减少、再增加的3个时序变化阶段,呈现波状起伏特征。民国中期（1930年）至建国初期（1958年）,因通江湖泊大量围垦,内湖数量和面积均呈现快速增加趋势;其后,1959年—1978年,在“以粮为纲”的农业背景下,内湖被大量围垦,数量和面积萎缩至民国中期（1930年）以来的最低值;改革开放（1979年）以来,随着“退田还湖”和水产养殖等农业结构调整,内湖面积呈现大幅增加趋势。洞庭湖区垸内湖泊总体呈增加趋势,但不同规模、不同类型、不同区域等不同特征的垸内湖泊变化情况不同。总体而言,围湖垦殖是内湖形成的主要驱动,粮食保障与水产养殖是内湖变化的主要原因。研究成果可为洞庭湖区生态修复与环境保护提供客观资料。