Improvement of the inversion accuracy of the Ca II line profiles using response functions

The line profile inversion is as a rule used to model the atmosphere. The algorithm for determining the temperature and electron density in the middle chromosphere has been developed based on inversion of the profiles of five Ca II chromospheric lines. Necessary information is contained in the function called by us the population function. This function is used to calculate the frequency-dependent source functions in lines and, consequently, the rates of radiative transitions. However, the accuracy of population function reconstruction in the inversion procedure is insufficient for the goal to be achieved. The accuracy is increased by using the intensity response function conception. The response functions of five Ca II lines to the population function disturbances have been calculated. The accuracy to be achieved has been estimated using one of the semiempirical models of the sunspot umbra as an example.     

基于半解析法有效和快速估计GRACE全球重力场的精度

首先基于半解析法建立了新的GRACE卫星K波段测量系统星间测速、GPS接收机轨道位置和加速度计非保守力误差联合影响累计大地水准面的误差模型;其次,基于各关键载荷精度指标的匹配关系,论证了误差模型的可靠性;最后,基于美国喷气动力实验室（JPL）公布的2006年的GRACE Level 1B实测误差数据,有效和快速地估计了120阶全球重力场的精度,在120阶处累计大地水准面的精度为18.368 cm,其结果和德国地学研究中心（GFZ）公布的EIGEN-GRACE02S全球重力场模型符合较好. 本文的研究为将来国际卫星重力测量计划（如GRACE Follow-On, 360阶）中高阶全球重力场模型精度的有效和快速估计提供了理论基础和计算保证.     

Ray-centred coordinate systems in anisotropic media

Whereas the ray-centred coordinates for isotropic media by Popov and Pšenčík are uniquely defined by the selection of the basis vectors at one point along the ray, there is considerable freedom in selecting the ray-centred coordinates for anisotropic media. We describe the properties common to all ray-centred coordinate systems for anisotropic media and general conditions, which may be imposed on the basis vectors. We then discuss six different particular choices of ray-centred coordinates in an anisotropic medium. This overview may be useful in choosing the ray-centred coordinates best suited for a particular application. The equations are derived for a general homogeneous Hamiltonian of an arbitrary degree and are thus applicable both to the anisotropic-ray-theory rays and anisotropic common S-wave rays.     

Solution of the transport equations using a moving coordinate system

A convection-diffusion equation arises from the conservation equations in miscible and immiscible flooding, thermal recovery, and water movement through desiccated soil. When the convection term dominates the diffusion term, the equations are very difficult to solve numerically. Owing to the hyperbolic character assumed for dominating convection, inaccurate, oscillating solutions result. A new solution technique minimizes the oscillations. The differential equation is transformed into a moving coordinate system which eliminates the convection term but makes the boundary location change in time. We illustrate the new method on two one-dimensional problems: the linear convection-diffusion equation and a non-linear diffusion type equation governing water movement through desiccated soil. Transforming the linear convection diffusion equation into a moving coordinate system gives a diffusion equation with time dependent boundary conditions. We apply orthogonal collocation on finite elements with a Crank-Nicholson time discretization. Comparisons are made to schemes using fixed coordinate systems. The equation describing movement of water in dry soil is a highly non-linear diffusion-type equation with coefficients varying over six orders of magnitude. We solve the equation in a coordinate system moving with a time-dependent velocity, which is determined by the location of the largest gradient of the solution. The finite difference technique with a variable grid size is applied, and a modified Crank-Nicholson technique is used for the temporal discretization. Comparisons are made to an exact solution obtained by similarity transformation, and with an ordinary finite difference scheme on a fixed coordinate system.     

On the estimation of accuracy for global models of gravitational field of the earth

The methods and techniques for estimating the accuracy of global models of the Earth’s gravity field in the form of spherical harmonic expansion of the geopotential are analyzed. Various methods for obtaining the a priori and a posteriori estimates for the accuracy are considered and classified. The application of different approaches is illustrated by numerical examples for nine models, including those recently developed using the modern methods of space geodesy. The basic requirements for the database and software for estimating the accuracy are formulated.     

新一代GRACE重力卫星反演地球重力场的预期精度

基于低低卫卫跟踪模式,本文主要探讨利用动力学法融合精密轨道数据和星间测距或距离变率数据求解地球重力场的基本原理与方法,该方法既可对两颗低低跟踪卫星的初始状态误差进行有效校正,也可充分利用低轨卫星轨道所包含的低频重力场信息.为探讨适合我国国情的低低跟踪模式下的重力卫星指标,本文以不同星载设备精度指标的组合进行模拟计算,模拟结果显示:(1)把GRACE卫星的星间距离变率指标提高一个量级,其余指标保持与GRACE卫星设计指标一致时,可使地球重力场的精度获得同量级的提高;(2)若星间距离变率为1.0×10-8 m·s-1,轨道高度为300 km,加速度计精度为3.0×10-10 m·s-2,轨道精度为0.03 m, 星间距离100 km,与利用GRACE的设计指标反演出的重力场精度相比,可提高约121倍,并建议我国未来低低跟踪重力卫星计划参考此指标.     

Evaluating the accuracy of Chinese pasture data in global historical land use datasets

Global historical land use datasets are widely used in global environmental change studies; however, uncertainties in the included pasture data have not been evaluated. In this study, using trend, relative difference ratio(RDR), and spatial comparisons, we evaluated the accuracy of China pasture data in the HYDE(versions 3.1 and 3.2), SAGE, and PJ datasets using historical document-based reconstructions, including the land use datasets for China covering 1935–1997(hereafter Ge-dataset)and Northeast China covering the 1700 s–2000(hereafter Ye-dataset), and the satellite-based China's Land-Use/cover Datasets(CLUDs) for 1980–2015. Four important results were obtained.(1) China pasture area in the HYDE, SAGE, and PJ datasets and grassland area in the Ge-dataset, Ye-dataset, and CLUDs show both disparate trends and large differences in absolute values.(2)Spatially, 50.9%, 52.8%, and 63.0% pasture/grassland grids in 2000 had RDRs greater than 60% between HYDE3.1, HYDE3.2,and SAGE datasets, and CLUDs, respectively. The percentage of grids with RDRs less than 20% were 24.9%, 26.7%, and 16.0%,respectively.(3) Based on HYDE3.2, the spatial distribution of pasture in Northeast China over the last 300 year has expanded,which is in contrast to the results from the Ye-dataset, which shows the spatial distribution of grassland shrinking because of human reclamation activities.(4) The large contrasts between the global datasets and Chinese native datasets are due to differences in pasture/grassland definitions, land use practices, and spatial reconstruction methods.     

新一代GRACE重力卫星反演地球重力场的预期精度

基于低低卫卫跟踪模式,本文主要探讨利用动力学法融合精密轨道数据和星间测距或距离变率数据求解地球重力场的基本原理与方法,该方法既可对两颗低低跟踪卫星的初始状态误差进行有效校正,也可充分利用低轨卫星轨道所包含的低频重力场信息.为探讨适合我国国情的低低跟踪模式下的重力卫星指标,本文以不同星载设备精度指标的组合进行模拟计算,模拟结果显示:(1)把GRACE卫星的星间距离变率指标提高一个量级,其余指标保持与GRACE卫星设计指标一致时,可使地球重力场的精度获得同量级的提高;(2)若星间距离变率为1.0×10^-8 m·s^-1,轨道高度为300 km,加速度计精度为3.0×10^-10 m·s^-2,轨道精度为0.03 m, 星间距离100 km,与利用GRACE的设计指标反演出的重力场精度相比,可提高约121倍,并建议我国未来低低跟踪重力卫星计划参考此指标.     

地磁与空间物理资料的组织和相关坐标系

在地磁与空间物理学实验研究和理论研究中,无论是观测资料和实验数据,还是理论计算和数值模拟结果,都需要组织在一定的坐标框架中,或者用一定的函数(函数系)来表达.坐标系或函数系的选择对于合理描述资料和发现自然规律是十分重要的环节.本文综述了地磁与空间物理研究中经常使用的20多种坐标系,着重探索这些坐标系的基本设计思想和相互联系,追踪它们的发展演化轨迹,从而加深对资料组织和坐标系选择的理解.     

2000国家GPS大地控制网的数据处理和精度评估

2000国家GPS大地控制网集合了各类GPS网近10年的观测数据, 提供GPS点在全球参考框架下2000.0历元的绝对位置.文中简要介绍GPS数据处理的方法,给出数据处理的精度统计;以部分点的实测速度为基础,建立中国大陆地壳运动速度场;基于GPS测量的动态系统状态,选择序贯卡尔曼滤波的广义测量平差方法,对GPS点做顾及地壳运动的整网平差,并给出整网平差结果的精度统计和外部检核.一系列的手段和措施保证了平差结果的合理性和可靠性.     

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ^° horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This “eddy-permitting” resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.     

Effect of ionospheric irregularities on accuracy of dual-frequency GPS systems

The residual error of the ionospheric correction of signals from dual frequency GPS systems, related to a ray deviation in random ionospheric irregularities, is studied. The formulas taking into account ray deviation from a straight line in an inhomogeneous ionosphere have been obtained based on the disturbance theory when solving the ray equations. These formulas have been used to estimate the extreme accuracy of a dual frequency GPS method. The algorithm for correcting this ionospheric error of the second order using three-frequency reception has been proposed.     

Spectral analysis for GNSS coordinate time series using chirp Fourier transform

Spectral analysis for global navigation satellite system (GNSS) coordinate time series provides a principal tool to understand the intrinsic mechanism that affects tectonic movements. Spectral analysis methods such as the fast Fourier transform, Lomb–Scargle spectrum, evolutionary power spectrum, wavelet power spectrum, etc. are used to find periodic characteristics in time series. Among spectral analysis methods, the chirp Fourier transform (CFT) with less stringent requirements is tested with synthetic and actual GNSS coordinate time series, which proves the accuracy and efficiency of the method. With the length of series only limited to even numbers, CFT provides a convenient tool for windowed spectral analysis. The results of ideal synthetic data prove CFT accurate and efficient, while the results of actual data show that CFT is usable to derive periodic information from GNSS coordinate time series.     

Ocean in peril: Reforming the management of global ocean living resources in areas beyond national jurisdiction

This article presents the outcome of research aimed at assisting governments in meeting their commitments and legal obligations for sustainable fisheries, based on increasing evidence that global fisheries are in crisis. The article assesses the effectiveness of the existing legal and institutional framework for high seas living resources. It focuses on: (1) the role of regional fisheries management organizations (RFMOs); (2) tools for compliance and enforcement to stem illegal fishing; and (3) mechanisms for habitat protection. The article further highlights a variety of options for addressing key weaknesses and gaps in current ocean governance, including United Nations General Assembly (UNGA) resolutions, reforms at the regional level, as well as a possible new legal instrument, with a view to informing international discussions on ways to ensure the sustainable use of high seas resources without compromising the health of the marine environment.     

两斜交柱坐标系间谐波函数坐标变换的辅助平面方法.

针对利用波函数展开法进行三维地震响应研究中的坐标变换问题,提出了在两斜交柱坐标系间谐波函数空间坐标变换的辅助平面方法．通过建立与空间点对应的一系列辅助平面,将地震波三维散射研究中的谐波函数表达式从柱坐标系变换至辅助平面内的极坐标系,然后在辅助平面内运用Graf加法公式,将柱坐标系下的波函数表达式变换至与该柱坐标系斜交的另一柱坐标系下,从而将三维问题转换成二维问题进行处理,最终得到了斜交柱坐标系下内域问题和外域问题的坐标变换公式．      

全球水储量变化的GRACE卫星检测

利用GRACE月尺度变化的地球重力场反演了全球水储量变化,并与陆地水文资料、卫星测高资料及海洋模式得到的结果进行了比对.通过对SOURE台站重力变化的陆地水储量变化计算结果和GRACE重力场系数截断为15阶得到的结果比较,发现两者比较接近,且年周期变化特征明显.对于亚马逊流域,当重力场系数截断为15阶且平滑半径使用106 m时,GRACE反演的区域平均水储量厚度的周年变化振幅为15.6×10-2m,小于使用平滑半径为4×105m的23.7×10-2m.在研究长江流域时,本文对水文资料做球谐系数展开,并与GRACE数据做同样的截断和平滑处理,结果发现GRACE反演的水厚度变化与水文资料结果基本上符合.对于纬度±66°之间的海洋区域,GRACE反演的海水质量变化接近于结合卫星测高和海洋模式得到的结果,但对于2°×2°网格,则在一些区域差异明显,最大超过了0.2 m,中误差为3.8×10-2m.可见,当前GRACE卫星时变重力场只能确定出上千公里及以上尺度区域的水储量变化.     

On the measurement of global radiation using black and white atmometers

Summary The difference in the evaporation from black and white Livingston atmometers was compared with the global radiation as measured by a Bellani spherical pyranometer. Results show that one unit (cm³) of difference in the evaporation corresponds to 20±3 cal/cm² of global radiation. The total of 95 daily values of radiation measured by pairs of atmometers differed from the corresponding Bellani figure by 2.3%. The operation of the atmometers and the accuracy of the results are discussed. The method is considered sufficiently accurate for general climatological purposes.Zusammenfassung Die Differenz in der Verdunstung, welche von schwarzen und weissen Livingston Atmometern angezeigt wird, wurde mit der Globalstrahlung, gemessen mit einem Bellani Kugelpyranometer, verglichen. Die Ergebnisse lehren, dass die Einheit (cm³) der Differenz in der Verdunstung 20±3 cal/cm² der Globalstrahlung entspricht. Der Unterschied der Summe von 95 Tageswerten der Strahlung, ermittelt auf der Verdunstung von weissen und schwarzen Atmometern, und jener, gemessen mit dem Bellani Kugelpyranometer, betrug 2.3%. Die Handhabung der Atmometer und die Genauigkeit der Resultate wird diskutiert. Die Methode scheint für klimatologische Zwecke hinreichend genau.     

利用半参数核估计法预报全球电离层总电子含量

本文将半参数平差模型引入电离层球谐函数系数的预报中,建立了半参数球谐函数模型（Semiparametric-Spherical Harmonic,Semi-SH）来预测全球电离层总电子含量.首先,通过快速傅里叶变换获得球谐函数系数的周期和振幅,将振幅高的主周期归入趋势函数,振幅低的剩余周期归入随机信号,建立了半参数模型,同时利用核估计方法拟合趋势函数,解算随机信号,并在时间域上进行外推,得到了预报时间的球谐函数系数,代入15阶电离层球谐函数模型,最后得出电离层总电子含量（Total Electron Content,TEC）的预报值.本文基于欧洲定轨中心（CODE）发布的球谐函数系数进行电离层TEC长期预报和短期预报分析,其中长期预报采用四年预报两年的模式对球谐函数系数进行预报,短期预报设计了三个算例,采用前30天预报后一天的模式,分别预报1天、滑动预报7天和滑动预报30天.实验结果表明：长期预报能够较好地反映全球电离层TEC的变化趋势和波动情况,Semi-SH模型对全球电离层TEC平均值（Mean TEC global,MTEC_global）的拟合值和预报值与MTEC_global实际值的相关系数分别为0.8743和0.8010,呈现出高度相关性.短期预报中,在太阳活动高年和太阳活动低年,Semi-SH模型在中纬度地区预报精度较CODE发布的电离层TEC 1天预报产品（CODE'S 1-Day Predicted GIM,C1PG）有较大提升,在高纬度与低纬度地区两种模型预报精度相当;Semi-SH模型在太阳活动高年和太阳活动低年30天滑动预报精度的均值均高于C1PG模型.实验结果说明了Semi-SH模型预报电离层TEC值的有效性.

     

首都圈地区“重复地震”及其在区域地震台网定位精度评价中的应用

根据Schaff和Richards[1]对“重复地震”的定义,利用波形相关识别出首都圈地区2002~2006年发生的“重复地震事件”.对宽频带/甚宽带和短周期记录,分别采用0.5~5.0 Hz、1.0~5.0 Hz带通滤波,识别出859例“重复地震”事件,占总数约24%.假定“重复地震对”间的距离小于1 km,地震台网测得的两个地震之间的距离系台网定位误差所致,而相应的震相走时差为人工拾取误差,给出了首都圈数字地震台网的定位精度估计.结果显示观测条件较好的首都圈东北部地区平均定位精度约为5 km;西南部地区约为13 km,是观测条件有待改善地区.     

On the possibility of using an electromagnetic ionosphere in global MHD simulations

Global magnetohydrodynamic (MHD) simulations of the Earths magnetosphere must be coupled with a dynamical ionospheric module in order to give realistic results. The usual approach is to compute the Reld-aligned current (FAC) from the magnetospheric MHD variables at the ionospheric boundary. The ionospheric potential is solved from an elliptic equation using the FAC as a source term. The plasma velocity at the boundary is the E × B velocity associated with the ionospheric potential. Contemporary global MHD simulations which include a serious ionospheric model use this method, which we call the electrostatic approach in this paper. We study the possibility of reversing the flow of information through the ionosphere: the magnetosphere gives the electric Reld to the ionosphere. The Reld is not necessarily electrostatic, thus we will call this scheme electromagnetic. The electric Reld determines the horizontal ionospheric current. The divergence of the horizontal current gives the FAC, which is used as a boundary condition for MHD equations. We derive the necessary formulas and discuss the validity of the approximations necessarily involved. It is concluded that the electromagnetic ionosphere-magnetosphere coupling scheme is a serious candidate for future global MHD simulators, although a few problem areas still remain. At minimum, it should be investigated further to discover whether there are any differences in the simulation using the electrostatic or the electromagnetic ionospheric coupling.