2000国家GPS大地控制网数据处理方法与结果

介绍了我国GPS大地控制网的基本情况,阐述了结合地壳运动观测网络工程建立全国统一的“2000国家GPS大地控制网”的意义,并介绍了整网平差中采用的技术方法和最后平差结果及精度。     

Extensions to Spatial Factor Methods with an Illustration in Geochemistry

Many applications involving spatial data require several layers of information to be simultaneously analyzed in relation to underlying geography and topographic detail. This in turn generates a need for forms of multivariate analysis particularly oriented to spatial problems and designed to handle spatial structure and dependency both within and between spatially indexed multivariate responses. In this paper we focus on one group of such methods sometimes referred to as spatial factor analysis. Use of these techniques has so far been mostly restricted to applications in the geosciences and in some forms of image processing, but the methods have potential for wider use outside these fields. They are concerned with identifying components of a multivariate data set with a spatial covariance structure that predominantly acts over a particular spatial range or zone of influence. We review the various forms of spatial factor analysis that have been proposed and emphasize links between them and with the linear model of coregionalization employed in geostatistics. We then introduce extensions to such methods that may prove useful in exploratory spatial analysis, both generally and more specifically in the context of multivariate spatial prediction. Application of our proposed exploratory techniques is demonstrated on a small but illustrative geochemical data set involving multielement measurements from stream sediments.     

Short and long memory unobserved components in hydrological time series

In this paper a semiparametric approach is introduced to decompose an ARFIMA model in the long memory and short memory unobserved components. The procedure is based on the DECOMEL method which produces a statistical decomposition by minimizing the Euclidean distance between the spectrum of the aggregated series and the sum of the parametric spectra of the components. The extension to long memory stationary models is achieved defining an approximate model where the fractional operator is replaced by the ratio of two polynomials of order one. The feasibility and performance of the proposed procedure are discussed through a case study.     

Study on the processing method of nighttime CO_2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the "true value" of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u., we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u.c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.     

Net ecosystem CO_2 exchange and controlling factors in a steppe——Kobresia meadow on the Tibetan Plateau

Knowledge of seasonal variation of net ecosystem CO2 exchange (NEE) and its biotic and abiotic controllers will further our understanding of carbon cycling process, mechanism and large-scale modelling. Eddy covariance technique was used to measure NEE, biotic and abiotic factors for nearly 3 years in the hinterland alpine steppe--Korbresia meadow grassland on the Tibetan Plateau, the present highest fluxnet station in the world. The main objectives are to investigate dynamics of NEE and its components and to determine the major controlling factors. Maximum carbon assimilation took place in August and maximum carbon loss occurred in November. In June, rainfall amount due to monsoon climate played a great role in grass greening and consequently influenced interannual variation of ecosystem carbon gain. From July through September, monthly NEE presented net carbon assimilation. In other months, ecosystem exhibited carbon loss. In growing season, daytime NEE was mainly controlled by photosynthetically active radiation (PAR). In addition, leaf area index (LAI) interacted with PAR and together modulated NEE rates. Ecosystem respiration was controlled mainly by soil temperature and simultaneously by soil moisture. Q10 was negatively correlated with soil temperature but positively correlated with soil moisture. Large daily range of air temperature is not necessary to enhance carbon gain. Standard respiration rate at referenced 10℃(R10) was positively correlated with soil moisture, soil temperature, LAI and aboveground biomass. Rainfall patterns in growing season markedly influenced soil moisture and therefore soil moisture controlled seasonal change of ecosystem respiration. Pulse rainfall in the beginning and at the end of growing season induced great ecosystem respiration and consequently a great amount of carbon was lost. Short growing season and relative low temperature restrained alpine grass vegetation development. The results suggested that LAI be usually in a low level and carbon uptake be relatively low. Rainfall patterns in the growing season and pulse rainfall in the beginning and at end of growing season control ecosystem respiration and consequently influence carbon balance of ecosystem.     

全球资料同化中误差协方差三维结构的准确估计与应用Ⅰ:观测空间协方差的准确估计

观测误差与背景误差协方差在四维资料同化和业务资料分析系统中起到决定性作用,它决定着观测信息与背景初猜值信息的相对重要性以及这些信息在空间及不同变量间的扩展方式。由于实际大气的真值并不知道,需要发展不同的技巧来估计观测误差与背景误差协方差,其中在观测空间利用观测与背景初猜值之差来分离观测误差与背景误差协方差的方法估计出的结果较为准确,其估计出的观测误差可直接用于资料分析系统中,背景误差可作为标尺来度量其他方法估计结果的可靠性。文章采用国家气象中心T213L31全球中期分析预报系统的6 h预报作为背景初猜场及同时段冬夏两个季节的北半球探空,利用贝塞尔函数拟合方法来分离观测误差与背景误差协方差,并比较了东亚区、北美区、欧洲区3个探空资料均匀密集区的区域与季节变化结果。结果表明,观测空间拟合方法所要求的水平均质、各向同性在欧洲区和北美区成立程度较好,在东亚区略差,使用时需要斟酌。此外均方差区域间差别较大,在冬季明显大于夏季,温度场偏大0.2 K,风场偏大0.9 m/s。温度场在400 hPa以下与150 hPa以上,背景误差略小于观测误差,而在200—300 hPa,背景误差略大一些。风场的特点与温度场比较一致。温度与风场背景误差主要集中在前40波,并在20波左右达到最大,水平相关季节区域差别不大,而温度垂直相关比风场窄,两者相关范围比较大的波数主要集中在前20波。此外利用贝塞尔函数拟合方法获得结果的分析表明,在质量场中不同区域季节间温度误差的稳定性要明显好于高度场,涡度散度的稳定性要明显好于流函数和势函数,特别是对于特征长度更为明显。这表明利用贝塞尔函数拟合方法获得的结果对校准在全球资料同化中采用温度、涡度散度作为资料同化的分析变量具有一定的优势。     

全球资料同化中误差协方差三维结构的准确估计与应用Ⅱ:背景误差协方差调整与数值试验分析

利用NMC方法针对背景误差协方差的方差、三维相关与特征长度来揭示T213L31模式的误差主要特点,并与传统更新矢量方法的计算结果进行了对比与调整。结果表明NMC方法结果与更新矢量方法结果在大体特征上基本吻合,但细节上的差异不可忽视,特别是对背景误差方差与特征长度的估计存在显著的差异,其主要原因是NMC方法倾向于高估天气尺度波的背景误差,而低估次天气尺度到中尺度波的背景误差。通过对背景误差方差、特征长度的调整,显著改善了背景误差功率谱的分布特点,使得NMC方法结果与更新矢量方法结果更为吻合。通过三维变分同化与最优插值中观测与背景误差相对重要性的比较,发现两者结果基本一致,但三维变分同化在850 hPa以下的温度场和300 hPa以上的风场统计结果都表现出背景误差相对于观测误差偏小的特点。背景误差相对于观测误差偏小有助于保证分析场中质量场与风场平衡,消除了大气底层和高层质量场与风场不匹配现象。在数值试验中,针对不同的背景误差均方差与特征长度的特点,分析了分析增量和预报效果的差异,结果表明,准确的背景误差估计与优化工作改善了预报效果,使得北半球三维变分同化的120 h预报效果整体好于现有最优插值。     

Seismic Moment Tensor Resolution on a Local Scale: Simulated Rockburst and Mine-induced Seismic Events in the Kopanang Gold Mine, South Africa

Seismic records contain information about the effect of the source as well as the effect of wave propagation through the rock mass. The effect of wave propagation is usually not well known as only simplified models of geological structures are available. Therefore, the information about the source retrieved by inverting seismograms may include errors due to incomplete knowledge of the rock mass along the propagation path, which in turn cause a distortion in the calculated moment tensor (MT). The distortion of the MT on a local scale was observed by inverting records of a simulated rockburst conducted at the Kopanang gold mine in South Africa. A dominant isotropic component of the explosive characteristics was found from the inversion. The deviatoric components retrieved from the blast are spurious. A test of their stability indicated that they are not significant, assuming an uncertainty above 5% for velocities and 10% for attenuation within the homogeneous model available for the mine. Thus, the retrieval of the MT from records of local networks in mines using a homogeneous model of the rock mass seems to be feasible. However, the homogeneous model of the rock mass can only be applied to close stations, within a few kilometers of the source. The seismic records from distant stations were too complex to be modelled by a homogeneous rock mass. Records of six mine-induced seismic events recorded at the Kopanang gold mine were also inverted. A vertical linear dipole along the pressure (P) axis was found for three of the events, suggesting a pillar burst. The mechanism of two events contains an isotropic implosion together with a nearly vertical dip-slip, and seems to indicate a combination of a cavity collapse with a down dip-slip along a nearly vertical fault. One event corresponds to a dipole along the tensional (T) axis. However, it is vertical, thus its association with tensile faulting of the hangingwall is uncertain.     

Processing and quality control of flux data during LITFASS-2003

Different aspects of the quality assurance and quality control (QA/QC) of micrometeorological measurements were combined to create a comprehensive algorithm which was then applied to experimental data from LITFASS-2003 (Lindenberg Inhomogeneous Terrain—Fluxes between Atmosphere and Surface: a long term Study). Eddy-covariance measurements of the latent heat flux were the main focus of the QA/QC efforts. The results of a turbulence sensor intercomparison experiment showed deviations between the different eddy-covariance systems on the order of 15%, or less than 30 W m⁻², for the latent heat flux and 5%, or less than 10 W m⁻², for the sensible heat flux. In order to avoid uncertainties due to the post-processing of turbulence data, a comprehensive software package was used for the analysis of experimental data from LITFASS-2003, including all necessary procedures for corrections and quality control. An overview of the quality test results shows that for most of the days more than 80% of the available latent heat flux data are of high quality so long as there are no instrumental problems. The representativeness of a flux value for the target land-use type was analysed using a stochastic footprint model. Different methods to calculate soil heat fluxes at the surface are discussed and a sensitivity analysis is conducted to select the most robust method for LITFASS-2003. The lack of energy balance closure, which was found for LITFASS-2003, can probably be attributed to the presence of low-frequency flux contributions that cannot be resolved with an averaging time of 30 min. Though the QA/QC system has been developed for the requirements of LITFASS-2003, it can also be applied to other experiments dealing with similar objectives.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.