Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

甘肃省东南部未来地震危险区的模式识别

在综合分析甘肃省东南部地震地质特征的基础上，利用模式识别方法判定未来地震危险区的可能位置．结果表明，本区易发生地震的地点与全新世断裂强烈活动段相吻合，从而提高了对未来地震危险区判定的可信度     

Proposed support design, Kaletepe tunnel, Turkey

In this paper, preliminary support design of Kaletepe tunnel, located on Bilecik-Istanbul highway, Turkey, was analyzed by empirical and numerical methods. The rock mass rating (RMR) and rock mass quality (Q) systems were employed for empirical rock mass quality determination. Numerical analysis for the stress–strain distribution of the tunnel excavation and support systems was also carried out. The applied support performance was investigated at different sections of the tunnel route. It was seen that empirical and numerical approaches showed similar results. This indicates that when the empirical method is supported by numerical method, the preliminary support design will be more reliable.     

Temporal Characteristics of Pacific Decadal Oscillation (PDO) and ENSO and Their Relationship Analyzed with Method of Empirical Mode Decomposition (EMD)

1. IntroductionPacific Decadal Oscillation (PDO) is a long-termENSO-like variability of the North Pacific. It can becharacterized by the first principal component of EOFof the North Pacific SST (Zhu and Yang, 2003; Tren-berth, 1990; Yang and Zhang, 2003). ENSO is thestrongest signal of annular change of global climatesystem (Trenberth, 1997). The spatial pattern of PDOis a wedge similar to El Nino. In the cool (warm)phases of PDO, the central and northwest Pacific is ofwarm (co…     

南大洋混合层的时空变化特征

过去对南大洋的研究受限于长期观测的缺乏,而现在地转海洋学实时观测阵(Arrayfor Real-timeGeostrophicOceanography,Argo)项目自开始以来持续提供了高质量的温度盐度观测,使系统地研究南大洋海洋上层结构成为可能。本研究使用2000—2018年的Argo浮标观测数据,分析了南大洋混合层深度(Mixed Layer Depth, MLD)的时空分布特征。结果表明:南大洋混合层存在明显的季节变化,冬春两季MLD在副南极锋面北侧达到最高值并呈带状分布,夏秋两季由于海表加热导致混合层变浅,季节变化幅度达到400m以上;在年际尺度上,MLD受南半球环状模(Southern HemisphereAnnularMode,SAM)调制,呈现纬向不对称空间分布特征,这与前人结果一致;本文指出在所研究时段,南大洋混合层在90°E以东,180°以西有加深趋势,而在60°W以西,180°以东有变浅趋势,显示出偶极子分布特征,并且这种趋势特征主要是风场的作用。     

基于EEMD和ARIMA的海温预测模型研究

类型丰富、时空分辨率高的海洋探测数据,为信号分解和机器学习算法的应用提供了可能。本文针对如何建立有效的海温预测模型这一问题,使用高时空分辨率的海表温度(SST)融合产品,引入信号处理领域的集合经验模态分解(EEMD)和机器学习领域的自回归积分滑动平均模型(ARIMA)。首先利用最适于分解自然信号的EEMD方法,将海温数据分解成多个确定频率的序列;再利用ARIMA分别对各个频率的序列进行预测,最后将各个序列的预测结果进行组合。该方法在丰富数据的支撑下,比以往直接使用海温数据所建立的预测模型精度更高,为更好地进行海温预测提供了新方法。     

Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

北太平洋副热带模态水盐度的年代际变化

副热带模态水(Subtropical Mode Water;STMW)在气候变化中起着重要作用。本文利用全球高分辨率数值模拟结果,研究了北太平洋STMW核心层盐度(Core Layer Salinity;CLS)的年代际变化及其物理机制。结果表明,CLS存在显著的年代际变化,其空间分布则与背景流场分布特征有关。侵蚀区CLS滞后生成区CLS约1~2年,这主要是海流平流输运引起的。生成区内,STMW的季节循环一般可分为生成期(12-4月)、隔离期(5-6月)和侵蚀期(7-11月),生成期混合层盐度(Mixed Layer Salinity;MLS)决定着隔离期和侵蚀期的CLS,而MLS年代际变化则主要由同太平洋年代际涛动存在负相关性的海表面淡水通量的变化引起。     

全球百年海表面温度年际和年代际变化特征分析

利用100 a(1903—2002年)HadiSST的逐月资料,将全球海表面温度异常(SSTA)作为整体进行经验正交函数分解(EOF分解),提取了控制各大洋SSTA的主导模态和各大洋之间的联合模态,分析了相应的空间分布和时间序列。研究结果表明:SSTA变化最剧烈的海区是赤道中东太平洋、西北大西洋湾流海区和北太平洋黑潮延伸体海区。热带太平洋厄尔尼诺和南方涛动是主导模态并具有2~7 a周期的年际变化;SSTA变化第二主模态和第三主模态都是以大约70 a为周期的年代际变化为主的跨大洋联合模态。第二主模态的空间分布主要表现为中纬度北太平洋和北大西洋反位相、热带太平洋与大西洋反位相的太平洋-大西洋双偶极子型分布。SSTA变化的第三模态主要呈现南北半球海洋反位相的特征,为北太平洋-北大西洋-南大洋联合模态。第四模态基本上是反映各个不同海域特有的局地海洋-大气相互作用模态,该模态的时间序列具有1~4 a周期的年际变化和约9 a周期的年代际变化。     

不同大尺度强迫条件下考虑初始场与侧边界条件不确定性的对流尺度集合预报试验

研究了不同大尺度强迫条件下的暴雨个例中,考虑不同尺度特征的初始扰动与侧边界扰动相互作用构造对流尺度集合预报的可行性,为进一步构建“自适应”于不同强对流天气的对流尺度集合预报系统提供依据。结果表明,在大尺度强迫显著的个例1中,以大尺度扰动信息为主的动力降尺度的增长趋势较集合转换卡尔曼滤波(ETKF)更为显著,且总扰动能量在预报中后期超过集合转换卡尔曼滤波,而在大尺度强迫较弱的个例2中,集合转换卡尔曼滤波扰动能量始终高于动力降尺度。此外,当大尺度强迫显著时,初始扰动与侧边界扰动相匹配会产生相互促进的作用,而不匹配时初始扰动会在预报中后期抑制侧边界扰动的发展,当大尺度强迫较弱时,即使是互相间不匹配的初始扰动与侧边界扰动也能在大部分预报时段起到相互促进的作用,说明初始扰动与侧边界扰动的相互作用机制不仅与天气形势相关,也与二者是否匹配挂钩,另外,扰动的发展特征同样依赖于天气形势;从集合离散度的角度来看,当大尺度强迫明显时,侧边界扰动的作用会在更短的时间内取代初始扰动,从而对离散度起到主导地位;两种初始扰动方法相比,集合转换卡尔曼滤波在多数情况下对总离散度的贡献均大于动力降尺度;从降水量预报及概率预报情况来看,大尺度强迫明显的个例可预报性更高,且各集合成员间的差异较小,大尺度强迫较弱的个例则相反,且当两种初始扰动方案与侧边界扰动相结合时,较仅侧边界扰动均有一定提高。