Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion     

Eric Blanger  Alain Vincent 《Journal of Hydrology》2005,300(1-4):114-125

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.  相似文献   

Model identification for hydrological forecasting under uncertainty   总被引：2，自引：2，他引：2  

Thorsten Wagener  Hoshin V. Gupta 《Stochastic Environmental Research and Risk Assessment (SERRA)》2005,19(6):378-387

Methods for the identification of models for hydrological forecasting have to consider the specific nature of these models and the uncertainties present in the modeling process. Current approaches fail to fully incorporate these two aspects. In this paper we review the nature of hydrological models and the consequences of this nature for the task of model identification. We then continue to discuss the history (“The need for more POWER‘’), the current state (“Learning from other fields”) and the future (“Towards a general framework”) of model identification. The discussion closes with a list of desirable features for an identification framework under uncertainty and open research questions in need of answers before such a framework can be implemented.  相似文献   

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan   总被引：2，自引：0，他引：2  

KURITANI  TAKESHI; KITAGAWA  HIROSHI; NAKAMURA  EIZO 《Journal of Petrology》2005,46(7):1421-1442

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope  相似文献   

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

《Ocean Modelling》2015

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.  相似文献   

High precision analysis of an embryonic extensional fault-related fold using 3D orthorectified virtual outcrops: The viewpoint importance in structural geology     

《Journal of Structural Geology》2016

Image-based 3D modeling has recently opened the way to the use of virtual outcrop models in geology. An intriguing application of this method involves the production of orthorectified images of outcrops using almost any user-defined point of view, so that photorealistic cross-sections suitable for numerous geological purposes and measurements can be easily generated. These purposes include the accurate quantitative analysis of fault-fold relationships starting from imperfectly oriented and partly inaccessible real outcrops. We applied the method of image-based 3D modeling and orthorectification to a case study from the northern Apennines, Italy, where an incipient extensional fault affecting well-layered limestones is exposed on a 10-m-high barely accessible cliff. Through a few simple steps, we constructed a high-quality image-based 3D model of the outcrop. In the model, we made a series of measurements including fault and bedding attitudes, which allowed us to derive the bedding-fault intersection direction. We then used this direction as viewpoint to obtain a distortion-free photorealistic cross-section, on which we measured bed dips and thicknesses as well as fault stratigraphic separations. These measurements allowed us to identify a slight difference (i.e. only 0.5°) between the hangingwall and footwall cutoff angles. We show that the hangingwall strain required to compensate the upward-decreasing displacement of the fault was accommodated by this 0.5° rotation (i.e. folding) and coeval 0.8% thickening of strata in the hangingwall relatively to footwall strata. This evidence is consistent with trishear fault-propagation folding. Our results emphasize the viewpoint importance in structural geology and therefore the potential of using orthorectified virtual outcrops.  相似文献   

集合卡尔曼滤波同化中雷达位置的敏感性研究     

乔小湜  闵锦忠  王世璋 《气象学报》2016,74(5):796-814

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。   相似文献   

变分同化框架通过背景误差协方差构建动力平衡约束的研究进展   总被引：1，自引：1，他引：0  

王瑞春  龚建东 《气象》2016,42(9):1033-1044

通过背景误差协方差构建动力平衡约束是变分同化框架设计的重要环节。它不仅帮助实现变量间的协同分析,提高观测使用效率,还能改善变分极小化问题的性状。本文在系统梳理通过背景误差协方差引入动力平衡约束的方式、流程的基础上,对求解目前全球和有限区域变分同化系统普遍采用的准地转平衡和静力平衡约束的共性问题和存在的不足作了归纳总结。分析了求解准地转平衡约束的三类方案:动力平衡方程方案、统计方案和动力-统计相结合方案的优缺点。对照比较了不同垂直离散方案下求解静力平衡约束时遇到的欠定问题的表现以及解决途径。最后,展望了基于背景误差协方差构建动力平衡约束在赤道等特殊地区、高分辨率同化系统、以及集合-变分混合同化系统发展中面临的挑战和机遇。  相似文献   

自动气象观测数据的质量控制   总被引：1，自引：1，他引：0  

郁珍艳  吴利红  高大伟  李正泉  樊高峰 《气象科学》2016,36(5):703-708

精细化气候变化分析及气候预测服务亟待发展,对高分辨率优质量的气候资料需求不断增长。本文利用缺测率、对比差值、空间一致性、粗差率、一致率5个指标对2009—2014年浙江省2 209个自动站温度数据进行评估分析,运用开放式可调整的评分标准,选取数据质量稳定可靠的站点,满足气候研究与业务的发展需求。共取得91个一级可用站,499个二级可用站,363个三级可用站,一级站主要分布于宁波和嘉兴地区,二级和三级站在各地都有较均匀的分布。在全省气候统计分析中,可以根据需求利用可用站。此外,有可用自动站参与的温度分布图完美呈现了地形对温度的影响,也使城市化发展突显出来,精细化的气候分析能更有效的为高山地区耕种制度及城市建设与规划提供科学依据。  相似文献   

EnSRF雷达资料同化在一次飑线过程中的应用研究   总被引：3，自引：1，他引：2  

高士博  闵锦忠  黄丹莲 《大气科学》2016,40(6):1127-1142

本文利用包含复杂冰相微物理过程的WRF（Weather Research and Forecasting）模式,针对2007年4月23日发生在我国华南地区的一次典型飑线天气过程,分别进行了确定性预报和集合预报试验,发现确定性预报能大致捕捉到飑线系统的发生发展过程,但对飑线后部的层云区模拟效果较差。集合预报能够有效地减少模式的不确定性,大部分集合成员对飑线的模拟效果优于确定性预报。进一步将集合预报得到的40个成员作为背景场,采用EnSRF（Ensemble Square Root Filter）同化多普勒天气雷达资料,并将分析得到的集合作为初始场进行集合预报,通过与未同化雷达资料的集合对比,考察了EnSRF同化多部雷达资料对飑线系统的影响。结果表明:EnSRF雷达资料同化增加了模式初始场的中小尺度信息,大部分集合成员的分析场能够较准确地再现飑线的热力场、动力场和微物理场的细致特征,并且模拟出飑线后部的层云结构。通过对EnSRF分析的集合进行模拟发现,大部分集合成员较未同化雷达资料时模拟效果有明显改善。同化后的集合预报ETS（Equitable Threat Score）评分最高,其次是未同化的集合预报,确定性预报的最低。  相似文献   

中国东部夏季暴雨的年代际跃变及其大尺度环流背景   总被引：4，自引：0，他引：4  

陈栋  陈际龙  黄荣辉  刘永 《大气科学》2016,40(3):581-590

本文利用1960~2011年中国东部地面测站的逐日降水资料和JRA-55再分析资料探讨了夏季暴雨分布的年代际跃变及其相关联的大尺度环流异常特征。基于暴雨频数和占比(夏季暴雨占比是指5~8月暴雨降水量对总降水量的贡献百分比)的分析结果表明:中国东部夏季暴雨分布在20世纪70年代末和90年代初经历两次反相的经向"三极子"跃变。中国东部夏季暴雨的年代际演变过程可分为三个时段:1960~1979年为华南和华北暴雨偏多、江淮流域暴雨偏少的经向"三极子"分布;1980~1991年为南方和华北暴雨偏少、江淮流域暴雨偏多的经向三极子"分布;1992~2011年为南方暴雨显著偏多、华北暴雨持续偏少,逐渐形成经向"偶极子"分布,并导致近十多年我国夏季"南涝北旱"的整体格局。1970年代末(1990年代初)跃变相关联的大尺度环流异常配置:东亚夏季风的减弱(增强),西太平洋副高的增强西伸但南撤(北抬),南亚高压的减弱南缩(增强东扩),以及蒙古高原中低层的气旋式(反气旋式)环流异常。与此同时,低层局地环流也发生调整:华北和黄淮地区以及华南和江南地区均为反气旋式(气旋式)环流异常,而江淮流域和四川盆地受控于风场切变式辐合(辐散)异常;涡度场发生相应变化,南北方大部分地区的负(正)涡度异常不(有)利于低涡的发展,而江淮流域和四川盆地的正(负)涡度异常有(不)利于低涡的发展,进而引发江南和华南暴雨减少(增加)、江淮流域和四川盆地暴雨增加(减少)、黄淮和华北暴雨减少(增加)的经向"三极子"跃变。  相似文献