Comparison of parameter uncertainty analysis techniques for a TOPMODEL application

Parameter uncertainty in hydrologic modeling is crucial to the flood simulation and forecasting. The Bayesian approach allows one to estimate parameters according to prior expert knowledge as well as observational data about model parameter values. This study assesses the performance of two popular uncertainty analysis (UA) techniques, i.e., generalized likelihood uncertainty estimation (GLUE) and Bayesian method implemented with the Markov chain Monte Carlo sampling algorithm, in evaluating model parameter uncertainty in flood simulations. These two methods were applied to the semi-distributed Topographic hydrologic model (TOPMODEL) that includes five parameters. A case study was carried out for a small humid catchment in the southeastern China. The performance assessment of the GLUE and Bayesian methods were conducted with advanced tools suited for probabilistic simulations of continuous variables such as streamflow. Graphical tools and scalar metrics were used to test several attributes of the simulation quality of selected flood events: deterministic accuracy and the accuracy of 95 % prediction probability uncertainty band (95PPU). Sensitivity analysis was conducted to identify sensitive parameters that largely affect the model output results. Subsequently, the GLUE and Bayesian methods were used to analyze the uncertainty of sensitive parameters and further to produce their posterior distributions. Based on their posterior parameter samples, TOPMODEL’s simulations and the corresponding UA results were conducted. Results show that the form of exponential decline in conductivity and the overland flow routing velocity were sensitive parameters in TOPMODEL in our case. Small changes in these two parameters would lead to large differences in flood simulation results. Results also suggest that, for both UA techniques, most of streamflow observations were bracketed by 95PPU with the containing ratio value larger than 80 %. In comparison, GLUE gave narrower prediction uncertainty bands than the Bayesian method. It was found that the mode estimates of parameter posterior distributions are suitable to result in better performance of deterministic outputs than the 50 % percentiles for both the GLUE and Bayesian analyses. In addition, the simulation results calibrated with Rosenbrock optimization algorithm show a better agreement with the observations than the UA’s 50 % percentiles but slightly worse than the hydrographs from the mode estimates. The results clearly emphasize the importance of using model uncertainty diagnostic approaches in flood simulations.     

东北地区水稻障碍型低温冷害变化对区域气候增暖的响应

利用东北地区153个气象站1961—2010年逐日气温资料,采用统计学方法分析了水稻障碍型低温冷害的气候变化特征及其对区域气候变暖的响应情况。结果表明,东北大部地区水稻障碍型低温冷害事件呈减少趋势,但区域性较为明显;障碍型低温冷害对关键发育期气温变化响应较为敏感,二者呈显著的负相关关系,表现为气温每升高1 ℃,东北地区冷害减少35个站次。东北地区关键发育期气温均呈上升趋势,但吉林西部地区障碍型冷害却随之增加,分析了关键发育期气温变率和气候变率,将其解释为局地障碍型冷害增加主要受气候变率增大的影响,逐日气温变率对其影响不大。     

跨越断层埋地管线地震反应研究述评

本文论述了跨越断层埋地管线地震反应研究工作概况,包括理论研究和实验研究的进展,并介绍了不同情况下采用的研究方法和建立的相应的研究模型,通过对数据结果进行整理分析,得出跨越断层埋地管线地震反应规律。同时,本文还对埋地管线研究中的一些重要参数对研究结果的影响作了详尽的阐述,提出了需进一步研究的问题和今后可能的研究发展方向。     

Three-dimensional analysis of ozone and PM<Subscript>2.5</Subscript> distributions obtained by observations of tethered balloon and unmanned aerial vehicle in Shanghai,China

A lightweight unmanned aerial vehicle (UAV) and a tethered balloon platform were jointly used to investigate three-dimensional distributions of ozone and PM_2.5 concentrations within the lower troposphere (1000 m) at a localized coastal area in Shanghai, China. Eight tethered balloon soundings and three UAV flights were conducted on May 25, 2016. Generalized additive models (GAMs) were used to quantitatively describe the relationships between air pollutants and other obtained parameters. Field observations showed that large variations were captured both in the vertical and horizontal distributions of ozone and PM_2.5 concentrations. Significant stratified layers of ozone and PM_2.5 concentrations as well as wind directions were observed throughout the day. Estimated bulk Richardson numbers indicate that the vertical mixing of air masses within the lower troposphere were heavily suppressed throughout the day, leading to much higher concentrations of ozone and PM_2.5 in the planetary boundary layer (PBL). The NO and NO₂ concentrations in the experimental field were much lower than that in the urban area of Shanghai and demonstrated totally different vertical distribution patterns from that of ozone and PM_2.5. This indicates that aged air masses of different sources were transported to the experimental field at different heights. Results derived from the GAMs showed that the aggregate impact of the selected variables for the vertical variations can explain 94.3% of the variance in ozone and 94.5% in PM_2.5. Air temperature, relative humidity and atmospheric pressure had the strongest effects on the variations of ozone and PM_2.5. As for the horizontal variations, the GAMs can explain 56.3% of the variance in ozone and 57.6% in PM_2.5. The strongest effect on ozone was related to air temperature, while PM_2.5 was related to relative humidity. The output of GAMs also implied that fine aerosol particles were in the stage of growth in the experimental field, which is different from ozone (aged air parcels of ozone). Geographical parameters influenced the horizontal variations of ozone and PM_2.5 concentrations by changing underlying surface types. The differences of thermodynamic properties between land and sea resulted in quick changes of PBL height, air temperature and dew point over the coastal area, which was linked to the extent of vertical mixing at different locations. The results of GAMs can be used to analyze the sources and formation mechanisms of ozone and PM_2.5 pollutions at a localized area.     

木本植被覆盖岸坡上波浪爬升过程的数值模拟研究

近岸木本植物构成的生态缓冲带作为新型的海岸软防护结构,兼具功能性和生态友好性,在沿海工程建设中愈发受到关注,如何深入开展其防护效果的机理研究是目前亟待解决的问题。本文采用数值模拟方法,在N-S方程中分别考虑树枝和树干的拖曳力影响,提出了木本植被作用下波浪沿斜坡爬升的表面波衰减的连续介质等效模型,并采用MAC法来跟踪自由曲面上的水颗粒轨迹。本文以波浪沿1/30的斜坡爬升为算例,对比讨论了有无植被作用下波浪的传播过程,并将算例结果与以往试验结果规律进行对照,验证了数值模型的有效性。最后,分别讨论了植物枝干的高度、密度、树枝倾斜角度等植被特性和波浪因素对植被消浪效果的影响,得到植被消浪的基本规律。文中的计算结果也可为实际的护岸工程和生态景观设计提供参考。     

北京天文台的WET三通道高速光度计

对北京天文台新近改造成的三通道高速光度计的结构和性能进行介绍,并给出了在兴隆站85公分望远镜的一些实测结果.     

A bio-optical inversion model to retrieve absorption contributions and phytoplankton size structure from total minus water spectral absorption using genetic algorithm

We propose a bio-optical inversion model that retrieves the absorption contributions of phytoplankton and colored detrital matter(CDM),as well as the phytoplankton size classes(PSCs),from total minus water absorption spectra.The model is based on three-component separation of phytoplankton size structure and a genetic algorithm.The model performance was tested on two independent datasets(the NASA bio-Optical Marine Algorithm Dataset(NOMAD) and the northern South China Sea(NSCS) dataset).The relationships between the estimated and measured values were strongly linear,especially for aCDM(412),and the Root Mean Square Error(RMSE) of the CDM exponential slope(SCDM) was relatively low.Next,the inversion model was directly applied to in-situ total minus water absorption spectra determined by an underwater meter during a cruise in September 2008,to retrieve the phytoplankton size structure in the seawater.By comparing the measured and retrieved chlorophyll a concentrations,we demonstrated that total and size-specific chlorophyll a concentrations could be retrieved by the model with relatively high accuracy.Finally,we applied the bio-optical inversion model to investigate changes in phytoplankton size structure induced by an anti-cyclonic eddy in the NSCS.     

拱形波纹钢屋盖动力特性测试及分析

拱形波纹钢屋盖是一种新颖的金属拱形薄壳结构，在我国已经得到了广泛的应用。由于结构的复杂性，目前有关的理论研究尚不够深入，动力研究还很少有人涉及。以18m和30m两种跨度的W666型拱形波纹钢屋盖为研究对象，测试了拱在竖向平面内振动时的频率、振型及阻尼。在所提简化计算模型的基础上，采用有限元方法计算了结构的动力特性，其结果与试验数据符合较好。测试和分析结果为该结构的动力反应研究提供了方便。     

线翼截断方式对大气辐射计算的影响

在大气辐射传输计算方法中,有3种基本方法,即,逐线积分方法,k-分布方法和带模式方法。其中,逐线积分方法是最精确的计算大气透过率的方法,本文根据透过率计算方式的不同,将逐线积分方法分为追线积分法和追点积分法。由于逐线积分计算需要耗费大量的计算时间,在大气遥感和大气探测业务中使用时,必须减少计算成本,提高计算速度。本文在追线积分法的基础上,给出了简化的逐线积分的基本方法,在保证同样计算精度的同时,大大提高了计算速度。对在精确的和简化的逐线积分下,不同线翼截断方式(CUTOFF)对吸收系数、大气透过率和冷却率的影响进行了更详细的探讨。通过数值试验发现,对谱线线翼的截断方式是影响辐射计算精度和计算速度的重要因子。在不同压力下,用CUTOFF=2计算的吸收系数误差最大;对CUTOFF=1,在大多数取样点上误差都小于2%;对CUTOFF=3或4,对绝大多数取样点上计算的吸收系数误差都在5%以内,但所用的计算时间却明显减少。大气低层的透过率对不同的计算方法和不同的线翼截断方式不敏感;对大气高层,无论是对精确的还是简化的逐线积分方法,当CUTOFF=2时的透过率结果与其他线翼截断方式的结果差别较大。通过比较,本文给出线翼截断的优选方案。