Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

Discretization error for the superposition of Gaussian beams

Summary. A high-frequency asymptotic integral expansion of a time-harmonic wavefield into Gaussian beams was derived in a previous paper by Klimeš. The discretization error caused by replacing this integral superposition by a discrete summation of Gaussian beams is estimated in this paper.     

Testing the correctness of the sequential algorithm for simulating Gaussian random fields

The sequential algorithm is widely used to simulate Gaussian random fields. However, a rigorous application of this algorithm is impractical and some simplifications are required, in particular a moving neighborhood has to be defined. To examine the effect of such restriction on the quality of the realizations, a reference case is presented and several parameters are reviewed, mainly the histogram, variogram, indicator variograms, as well as the ergodic fluctuations in the first and second-order statistics. The study concludes that, even in a favorable case where the simulated domain is large with respect to the range of the model, the realizations may poorly reproduce the second-order statistics and be inconsistent with the stationarity and ergodicity assumptions. Practical tips such as the multiple-grid strategy do not overcome these impediments. Finally, extending the original algorithm by using an ordinary kriging should be avoided, unless an intrinsic random function model is sought after.     

区域集合预报系统2 m温度预报的校准技术

采用非齐次高斯回归 (NGR) 技术对国家气象中心区域集合预报系统的2 m温度预报结果开展了一阶偏差和二阶离散度的校准研究。对预报结果比较详尽的检验分析表明:校准后的2 m温度预报可靠性和预报技巧均显著提高,表现为校准后集合预报成员的均方根误差与离散度更为接近;原Talagrand直方图中的“L”形分布现象得到有效改善;Brier评分、最小连续分级概率评分 (CRPS) 明显减小,相对作用特征 (ROC) 面积增大,说明校准后的2 m温度预报表现出更好的预报技能。此外,NGR技术与自适应误差订正技术的对比试验表明,NGR在消除集合平均偏差和提高集合离散度两个方面均有优势。     

Determination of ammonia emission rates from a tunnel ventilated chicken house using passive samplers and a Gaussian dispersion model

Atmospheric deposition can provide a significant fraction of the nitrogen loading to coastal waters. The Delmarva Peninsula, on the eastern shore of the Chesapeake Bay, is a region with intense poultry production that may supply a significant source of atmospheric ammonia (NH₃). There is a need to quantify ammonia NH₃ from representative growing methods in this region in order to more accurately estimate agricultural NH₃ emissions and to develop best management practices. In this study, NH₃ emissions were determined at an 18,600-bird tunnel-ventilated chicken house using a modified sampling grid (a planar arrangement normal to the length of the house) with Ogawa passive samplers to characterize the emission plume downwind from the house. This improvement in the sampling strategy, compared to a previous study, simplified the inverse Gaussian plume analysis which improved the confidence in the emission factors. In this study, a total of four separate plume characterizations were conducted over the final 3 weeks of the 6-week broiler grow-out cycle. The mean emission factor observed at the tunnel-ventilated house, 0.13 g NH₃–N/bird-day, was an order of magnitude lower than that previously observed at a nearby side-wall ventilated house. Although not all growing variables were measured, the large difference in emission factors between the two ventilation regimes suggest that modern, tunnel-ventilated houses may result in a significant decrease in NH₃ emissions compared with traditional growing methods. This variability in emission factors underscores the need for characterizing chicken houses under various conditions and determining the factors that control these atmospheric emissions.     

蒙古国植被物候特征及其对地理要素的响应

蒙古高原是中国重要的北方生态屏障。在全球气候变化的背景下,研究蒙古国植被物候变化特征对于认识蒙古国草地生态系统对气候变化的响应和促进区域畜牧业可持续发展具有重要意义。本研究利用非对称高斯拟合法对蒙古国2001—2019年MOD13Q1产品中的归一化植被指数（Normalized Differential Vegetation Index,NDVI）数据拟合,得到较为平滑的NDVI时间序列数据;基于TIMESAT平台,采用动态阈值法分析获得蒙古国连续19a植被物候数据。研究分析了蒙古国植被物候的空间分布及年际变化趋势,发现蒙古国植被返青期（Start of growing season,SOS）主要集中在110~150d,总体呈微弱推迟趋势,植被枯黄期（End of growing season,EOS）主要集中在270~310d,总体呈提前趋势,从而导致蒙古国生长季长度（Length of growing season,LOS）呈缩短趋势,且缩短时间最长可达2d以上。采用偏相关分析方法分析了植被物候对地形、降水、地表温度等地理要素的响应,表明蒙古国植被物候具有明显的空间异质性和海拔依赖性,不同植被物候对降水、地表温度（Land Surface Temperature,LST）的响应不同,SOS与日间LST呈显著正相关,EOS与夜间LST呈显著正相关,而LOS与年均降水呈显著负相关关系。     

Probability Distribution of the Hull Motion and Mooring Line Tension of Two Floating Systems

The statistical and distribution characteristics of the responses of a floater and its mooring lines are essential in designing floating/mooring systems.In general,the dynamic responses of offshore structures obey a Gaussian distribution,assuming that the structural system,and sea loads are linear or weakly nonlinear.However,mooring systems and wave loads are considerably nonlinear,and the dynamic responses of hull/mooring systems are non-Gaussian.In this study,the dynamic responses of two types of floaters,semi-submersible and spar platforms,and their mooring lines are computed using coupled dynamic analysis in the time domain.Herein,the statistical characteristics and distributions of the hull motion and mooring line tension are discussed and compared.The statistical distributions of the dynamic responses have strong non-Gaussianity and are unreasonably fitted by a Gaussian distribution for the two floating and mooring systems.Then,the effects of water depth,wave parameters,and low-frequency and wave-frequency components on the non-Gaussianity of the hull motion,and mooring line tension are investigated and discussed.A comparison of the statistical distributions of the responses with various probability density functions,including the Gamma,Gaussian,General Extreme Value,Weibull,and Gaussian Mixture Model(GMM)distributions,shows that the GMM distribution is better than the others for characterizing the statistical distributions of the hull motion,and mooring line tension responses.Furthermore,the GMM distribution has the best accuracy of response prediction.     

知识引导的稀疏时间序列遥感数据拟合

在多云多雨的地区,光学遥感存在着获取无云数据困难的难题,这会导致时间序列应用中可用数据匮乏。因此,本文面向稀疏时间序列遥感数据,根据噪声造成遥感影像上归一化差分植被指数(NDVI)被低估的事实,提出了一种知识引导的拟合方法。首先,在遥感影像预处理的基础上,利用先验知识和时序差分法对噪声进行识别和剔除;然后,采用高斯二阶模型对原始数据进行拟合;最后,根据拟合残差更新权重,进行迭代拟合,重复上述过程直至获得稳定的结果。本文以Landsat 8 OLI作为数据源,对浙江省杭州地区的森林数据进行拟合,结果表明:在稀疏时间序列数据的情况下,本文方法与MODIS数据拟合结果的相关系数达到0.92,关键时点(如NDVI峰值点等)的时间误差在5 d;相比当前主流方法的0.88与8 d具有更高的精度。     

利用波形去噪算法的森林最大冠顶高估测研究

针对大光斑激光雷达回波信号噪声影响森林冠顶高估测精度,且回波分析法判定回波位置受限于平坦地区的问题,利用高斯低通滤波和小波去噪两种方法对GLAS波形进行去噪处理,提出了结合均方根倍差法和回波分析法来判定回波位置的有效算法。经小波去噪后信号的信噪比23.360 704,均方根误差为0.000 233 3,经均方根倍差法和回波分析法相结合来判定回波位置估测的冠顶高结果与实测结果相关性系数r值为0.864,效果均优于高斯低通滤波去噪。基于GLAS回波数据实验结果表明:小波去噪较好地实现了对回波信号的去噪处理,均方根倍差法和回波分析法相结合,实现了对坡度相对较大地区的GLAS波形的回波开始位置和地面回波位置的准确判定,对森林冠顶高的精确估算具有重要意义。     

基于全波形激光雷达数据分解方法的研究

波形分解是全波形激光雷达数据处理的关键技术之一,通过雅克比矩阵更新与否以及阻尼系数的调整,对优化LM(Levenberg-Marquard)算法进行了研究。实验表明,该方法可最优化求解足印内各波形分量的振幅、中心位置和半宽,而不必每次迭代时都需计算雅克比矩阵,简化了计算过程,提高了计算效率,最终得到的模拟波形能够非常好地拟合原始波形。