A six-step practical approach to semivariogram modeling

Geostatistical prediction and simulation are being increasingly used in the earth sciences and engineering to address the imperfect knowledge of attributes that fluctuate over large areas or volumes—pollutant concentration, electromagnetic fields, porosity, thickness of a geological formation. Central to the application of such techniques is the need to know the spatial continuity, knowledge that is commonly condensed in the form of covariance or semivariogram models. Their preparation is subdivided here into the following steps: (1) Data editing, (2) Exploratory data analysis, (3) Semivariogram estimation, (4) Directional investigation, (5) Simple modeling, (6) Nested modeling. I illustrate these stages practically with a real data set from a geophysical survey from Elk County, Kansas, USA. The applicability of the approach is not limited by the physical nature of the attribute of interest.     

南极GPS数据质量分析

以多路径误差和电离层延迟为评价指标,利用TEQC软件和QCVIEW软件对南极埃默里冰架、内陆冰盖、纳拉峡海湾、拉斯曼丘陵的实测GPS数据以及部分IGS站点的数据进行了数据质量分析,得到结论： 1)在同一站点上,L2上的多路径误差总比L1上的多路径误差要大;2)在南极不同环境下的多路径误差不同;3)纬度和太阳活动剧烈程度是影响南极地区GPS数据电离层延迟大小的两个重要因素。     

岷县漳县6.6级地震前区域地壳运动变形背景与断层形变异常特征

利用甘肃岷县漳县6.6级地震区所在的青藏块体东北缘地区的区域水准、GPS、流动重力和跨断层短测线等地形变监测资料,结合地质构造、动力环境和已往的研究结果,分析了不同类型资料反映的震前区域性地壳变形背景、断层形变异常特征和可能的机理.结果认为：（1）岷县漳县6.6级地震前西秦岭、六盘山等构造区不同程度地存在着GPS水平挤压闭锁高应变积累、垂直隆升异常高梯度带和重力升、降差异剧烈变化等中长期背景;（2）震前到震时发震断裂附近及其外围相关构造区域断层形变异常在空间和时间上的起伏波动变化显著,尤其是汶川地震以来的波动变化在一定程度上反映了与本区构造比邻的龙门山断裂带剧烈右旋错动对本区的影响,与本次岷县6.6级地震过程有关;（3）本区有地形变监测资料积累以来缺乏6级以上震例,虽然存在不同程度的中长期形变背景异常,但何时进入短期-短临阶段确实很难把握,需要不断积累总结和探索提高.     

由空间直角坐标计算大地坐标的简便公式

对大地测量教科书中给出的空间三维直角坐标与大地坐标的转换公式进行了大量的计算与分析，总结得出一组由空间直角坐标计算大地坐标的简便公式。     

胶东金矿金成色及其指示意义

影响金成色的因素较多,包括矿床成因类型、形成深度、变质程度、成矿时代、成矿温度和矿化阶段等,其中金矿床形成深度是最关键的因素。矿床形成深度大,成矿的物化参数变化慢,形成的金成色高;形成深度小,成矿的物化参数变化快,形成的金成色低。以胶东金矿为例,统计并定性分析了这种关系,证明胶东金矿区存在类似的规律,据此可以推广,以矿床金的成色来判断矿床的形成深度和矿体的延深情况。胶东金矿金成色较低,以５５０～７５０居多,反映了金矿成矿以中浅成为主。南张家金矿金成色最低,暗示其形成深度最浅（小于１．５ｋｍ）。     

辽宁省积雪对自动气象站观测地气温差的影响

利用1971—2016年辽宁省61个气象站气温、地表温度、积雪日数和积雪深度资料,分析了积雪的保温作用及其对地气温差的影响。结果表明：更换自动站前后地表温度观测方式的差异导致地气温差显著增大,地气温差的增大程度受所在区域积雪日数、积雪深度的影响显著。在积雪期较长、积雪较厚的地区,积雪引起反照率增大,使得雪面温度降低,导致雪气温差减小,而雪的保温作用使得地气温差显著增大。因此,更换自动站前地（雪）气温差与积雪日数呈显著负相关,而更换自动站后地气温差与积雪日数呈显著正相关。各台站之间地气温差随积雪深度的变化系数差异较大,为0.045~0.858 ℃?cm^-1,在年平均积雪日数＜40 d、年平均极端积雪深度＜10 cm的区域,积雪的保温作用随积雪深度增大而显著增大;在年平均积雪日数＞40 d、年平均极端积雪深度＞10 cm的区域,10 cm以下的积雪对土壤保温作用随积雪深度增大显著,当积雪深度＞10 cm后,其保温作用随积雪深度增大的幅度明显减小。     

近50年气候变化对滦河上游径流的影响

为了解近50年滦河上游气温和降水气候的变化特征、趋势及其对该流域径流量的影响,利用1956-2009年滦河上游的实测气温、降水量和径流深资料,分析了该流域气温、降水和径流深的年均和季度变化的时间序列,并建立了该流域气候变化对径流影响的复相关回归模型。结果表明,年径流深随着年降水量的减少而减少,随着年平均气温的升高而减少;春、夏和秋季的径流深随着同期降水量和气温的变化趋势与年际变化趋势基本一致,但冬季径流深则相反,而且其变化幅度非常小。     

整体最小二乘的迭代解法

在引入整体最小二乘平差准则的基础上,推导了整体最小二乘的迭代解法;同时,引入多元函数隐函数求导的方法以确定未知参数对观测数据的线性信息,解决了整体最小二乘下的精度评定问题。给出了运用新的解法在拟合函数确定以及坐标转换参数确定等方面的应用实例,验证了新算法的可行性。     

ANNUAL VARIATIONS OF TEMPERATURE ON FOUR URBAN UNDERLYING SURFACES AND FITTING ANALYSIS

Based on the observed 2-year temperature data for four kinds of typical urban underlying surfaces, including asphalt, cement, bare land and grass land, the annual variations and influencing factors of land surface temperature are analyzed. Then fitting equations for surface temperature are established. It is shown that the annual variation of daily average, maximum and minimum temperature and daily temperature range on the four urban underlying surfaces is consistent with the change in air temperature. The difference of temperature on different underlying surfaces in the summer half year (May to October) is much more evident than that in the winter half year (December to the following April). The daily average and maximum temperatures of asphalt, cement, bare land and grass land are higher than air temperature due to the atmospheric heating in the daytime, with that of asphalt being the highest, followed in turn by cement, bare land and grass land. Moreover, the daily average, maximum and minimum temperature on the four urban underlying surfaces are strongly impacted by total cloud amount, daily average relative humidity and sunshine hours. The land surface can be cooled (warmed) by increased total cloud amount (relative humidity). The changes in temperature on bare land and grass land are influenced by both the total cloud amount and the daily average relative humidity. The temperature parameters of the four land surfaces are significantly correlated with daily average, maximum and minimum temperature, sunshine hours, daily average relative humidity and total cloud amount, respectively. The analysis also indicates that the range of fitting parameter of a linear regression equation between the surface temperature of the four kinds of typical land surface and the air temperature is from 0.809 to 0.971, passing the F-test with a confidence level of 0.99.     

Dynamical response to supernova-induced gas removal in spiral galaxies with dark matter halo

We investigate the dynamical response, in terms of disc size and rotation velocity, to mass loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin density profile embedded in a spherical dark matter halo having a density profile proposed by Navarro, Frenk & White is considered. For the purpose of comparison, we also consider the homogeneous and   r ⁻¹  profiles for dark matter in a truncated spherical halo. Assuming for simplicity that the dark matter distribution is not affected by mass-loss from discs and the change of baryonic disc matter distribution is homologous, we evaluate the effects of dynamical response in the resulting discs. We found that the dynamical response only for an adiabatic approximation of mass-loss can simultaneously account for the rotation velocity and disc size as observed particularly in dwarf spiral galaxies, thus reproducing the Tully–Fisher relation and the size versus magnitude relation over the full range of magnitude. Furthermore, we found that the mean specific angular momentum in discs after the mass-loss becomes larger than that before the mass-loss, suggesting that the mass-loss would most likely occur from the central disc region where the specific angular momentum is low.