淮河流域夏季降水异常与若干气候因子的关系

基于旋转经验正交函数分解 (REOF) 方法探讨淮河流域1961—2010年夏季降水与厄尔尼诺/南方涛动 (ENSO)、北大西洋涛动 (NAO)、印度洋偶极子 (IOD)、太平洋年代际振荡 (PDO) 之间的关系,并进一步分析各气候因子不同位相单独以及联合对淮河流域夏季降水的影响。结果表明:淮河流域夏季降水与ENSO,PDO,NAO,IOD等气候因子具有较稳定的相关性,其中,PDO和IOD是影响淮河流域夏季降水的关键因子,且PDO与夏季降水呈显著负相关关系;各气候因子的冷暖位相单独及联合对淮河流域夏季降水的影响不同,PDO的冷期以及NAO,IOD冷位相使流域北部的夏季降水量呈显著增加趋势,PDO分别联合ENSO,NAO和IOD的冷、暖位相对流域北部地区和淮河上游地区的夏季降水影响显著。     

高性能计算技术在气象领域的应用

高性能计算通过应用超级计算机与并行处理技术解决复杂的计算问题,是信息技术发展比较迅猛的领域之一。气象应用始终是高性能计算的重要领域,高性能计算技术有效地解决了高分辨率、高精度气象数值预报模式发展限制,在气象预报预测业务中发挥着核心支撑作用。数十年以来,由于数值天气模式研究和业务运行对计算资源的强烈需求,国内外气象领域高性能计算机及应用迅速发展起来。气象领域对高性能计算能力及系统的可靠性需求日益提升。高性能计算技术将与气象预报应用日益融合,相互影响促进,不断创新发展。为满足气象预报预测业务和科研工作需求,中国气象部门将进一步提升高性能计算能力,并致力于优化集约高性能计算系统布局,高效管理计算资源,发挥最大效益。     

干侵入对我国东部一次强飑线过程的作用分析

利用NCEP/NCAR再分析资料、雷达资料、FY2E卫星资料等,对2014年7月30日至31日发生在我国东部的一次强飑线天气过程进行综合分析,并着重分析了干侵入对此次飑线过程的作用。结果表明：中层阶梯槽引导高层干冷空气向下入侵后叠加在低层暖湿气流之上,增强了大气不稳定,不稳定能量累积;低空切变线促使不稳定能量释放,是此次飑线的重要触发机制。对干侵入的分析结果表明：此次飑线过程中,干侵入来源于高空槽后下沉气流,干侵入底部风场调整使得飑线在移动过程中发生了转向。在飑线发展过程中,干侵入范围扩大、强度增强;干侵入使飑线发生区域上空空气变冷、变干,增强了大气不稳定。干侵入底部的低相当位温空气与飑线内部高相当位温空气混合是飑线长时间维持的重要原因。     

防雷箱残压“二次感应”的研究试验与仿真

防雷箱广泛应用在各类低压配电系统的雷电防护中,为研究防雷箱在8/20μs冲击电流作用下常出现"二次感应"的现象,基于Haefely DIAS 733测试系统通过大量试验,结合电磁感应理论对感应电压进行计算,并利用PSPICE建立感应回路。结果表明:1防雷箱内部平行布线且输入输出同端口时,"二次感应"电压幅值服从电磁感应定律,且可进行较精确计算。2采用双绞线布线、地线与输出线非平行、输入与输出线异端口且与地线非平行3种方法可有效地降低因电磁感应而增加的残压,从而获得较低的输出残压。3试验表明,当内部采用输入与输出线异端口且与地线非平行的综合布线方案后,效果最明显,在60kA冲击电流下,相对未改进的防雷箱,采样点残压下降约45%。     

贵州西北部基性侵入岩锆石U-Pb年代学、岩石地球化学特征 及其与铅锌成矿关系初探

贵州西北部发育多处基性侵入岩,总体规模较小(约0.25 Km2),多沿深大断裂侵位,出露于铅锌矿点外围,本文对猫猫厂、凉山两处矿点附近的儿马冲和白岩庆两地小型基性侵入岩进行了重点研究。侵入岩主要岩性为细粒辉长岩,造岩矿物主要为拉长石、普通辉石。SiO2范围为49.60-51.09 wt%,MgO从3.88-4.27 wt %,TiO2为3.69-3.85 wt %。LA-ICP-MS锆石U-Pb定年结果为268.3±7.4 Ma,显示岩浆侵位于二叠纪。基性侵入岩的微量元素蛛网图呈OIB型特征,富集大离子亲石元素（LILE）、轻稀土元素（LREE）,亏损重稀土元素（HFEE）,相对亏损高场强元素（Nb,Ta）,有Sr、Y亏损,Pb富集。(87Sr/86Sr)i范围0.706749～0.707069,(143Nd/144Nd)i范围0.512313～0.512363,εNd(t)范围0.2～1.2;源区熔融深度处于石榴石橄榄岩相深度,可能经历了1-3 %的部分熔融,处于亏损石榴石二辉橄榄岩相向原始石榴石二辉橄榄岩相的过渡区。成岩过程中发生了单斜辉石、斜长石等矿物分离结晶,受到了有限的地壳混染作用,未经历明显的AFC过程。地壳物质在地幔源区中的加入可能是造成地幔富集的主要原因。侵入岩与成矿作用之间的关系,主要通过两方面所表现。一方面是二者间构造活动上的耦合性。另一方面是基性岩在成矿过程中可能发挥了重要的化学屏障层作用。     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

An Immune Genetic Algorithm to Buildings Displacement in Cartographic Generalization

Spatial conflicts may occur when map data are displayed at a scale smaller than that of the source map. This study applies the displacement operator in cartographic generalization to resolve such spatial conflicts and to improve the clarity and legibility of map. The immune genetic algorithm (IGA) is used in this study for buildings displacement to solve conflicts. IGA is based on the genetic algorithm (GA) and employs the self‐adjusting mechanism of antibody concentration to enhance population diversity. Meanwhile, the elitism retention strategy is adopted in IGA to guarantee that the best individual (antibody) is not lost and destroyed in the next generation to strengthen convergence efficiency. The compared experiment between IGA and GA shows that the displacement result produced by IGA performs better than GA. Finally, in order to make the displaced map more attractive to cartographers, two constraints – the building alignment constraint and building tangent relation constraint – are applied in IGA to restrict the buildings’ displacement. The same experimental data are adopted to prove that the improved IGA is useful for maintaining the two constraints.     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.     

Optimization and Parametrical Investigation to Assess the Reconstitution of Different Types of Indian Sand Using Portable Travelling Pluviator

By using the air pluviation technique, it is aimed to achieve the desired relative density with uniform void ratio throughout the specimen in order to maintain the homogeneity and to avoid the spatial variability. Further, in order to achieve the maximum deposition intensity, a systematic optimization study has been carried out rigorously in a test tank to determine the diameter of the orifice to be employed for the sieve plates of different porosity and the number of sieve plates to be installed in the diffuser sieve sets. The study has been conducted with four different patterns of sieves with different porosity to achieve a wide range of relative densities for four different uniformly graded Indian sands. The dynamic penetrometer which is considered to be one of the cost effective instruments has been efficiently used to determine the soil resistance at various locations of the test tank for every given height of fall in order to check the uniformity of placement density throughout the sand bed. The study reveals that the sand beds of different relative densities could be achieved using different patterns of diffuser sieves at optimum sand flow rate without compromising the uniformity. The effect of height of fall as well as porosity of diffuser sieves on the relative density of different sands has been studied in detail. The deposition intensity and the relative density obtained from the present study are compared with the values available in the literature.