LiDAR-based digital terrain analysis of an area exposed to the risk of lava flow invasion: the Zafferana Etnea territory, Mt. Etna (Italy)

The town of Zafferana Etnea, located on the southeastern slope of Mt. Etna volcano (Italy), has been repeatedly threatened by lava flows in recent centuries. The last serious threat occurred during the 1991–1993 eruption, when the lava front came to a halt only 1.7 km from the centre of town. Morphostructural data derived from light detection and ranging (LiDAR) surveys carried out on Etna in 2005 have enabled us to evaluate the risk of lava invasion in a section (16 km²) of the Zafferana Etnea territory. Qualitative and quantitative results are obtained combining the information derived from LiDAR analysis with geological, morphological and structural data using geographic information systems technology (GIS). The study quantifies in unprecedented detail the areal extent and volume of forested and urban areas and its degree of exposure to different levels of hazard from future lava invasion. Nearly 52% of the urban texture fall into areas of moderate to high risk from lava invasion. Future land use planning should take these findings into account and promote new development preferentially in areas of lower risk.     

基于探空资料的江西典型冻雨天气过程垂直结构分析

利用探空资料和地面观测资料,对1959—2008年发生在江西的典型冻雨天气过程的大气垂直结构进行了分析。结果表明:冻雨的形成与大气饱和层、冰晶层、暖层、逆温层、冷层及特征物理量等因素的关系密切。对流层中下层的大气呈饱和或准饱和的逆温状态,饱和层顶位于暖层之上是冻雨发生的主要层结特征。冰晶层底部的平均高度为675hPa,平均厚度为2309m;暖层底部的平均高度为834hPa,平均温度为4.1℃,平均厚度为1765m;逆温层的平均温差为7.5℃;冷层的平均厚度为1668m,平均最低温度为-3.9℃;平均地面最低温度为-1.2℃;冻雨发生时大气整层比湿积分指数须增大为1000以上,干暖盖强度指数须降到-5℃以下。     

流域坡面汇流研究现状述评

对当前坡面汇流计算方法的研究进展进行了较为系统的总结与分析,并对坡面汇流的非线性效应以及城市低影响开发中的雨水入渗与蓄集对坡面汇流的控制作用进行了简要分析。从模型简单实用的角度出发,认为以流域时间-面积关系与线性水库相串联的ModClark法等为代表的概念性分布式坡面汇流模型具有良好的发展前景;考虑到基于等流时单元的变动等流时线法在反映雨强非线性影响中存在的问题,认为根据水文响应单元在不同雨强条件下汇流时间的变化,调整其汇流参数以反映坡面汇流的非线性效应,对于流域坡面汇流的分布式模拟更具有实际意义;针对目前低影响开发设施长时间序列大空间尺度的室外降雨径流监测资料普遍较为缺乏的现状,给出了后期应积极选择合适的技术以加强低影响开发性能监测工作的建议。     

西藏羊湖二厂电站厂房区泥石流危险性评价

泥石流的活动特征及其危险性评价对电站的施工及运行有重要意义。根据羊湖电站厂房区勘查所获取的资料,对厂房区后山5条泥石流沟的成因进行了分析,在此基础上将层次分析法(AHP)应用到泥石流危险性评价研究中。对该厂房区5条泥石流沟的危险度评价表明,5号沟泥石流危险度最高,4号沟次之,1号沟、2号和3号沟危险度相对较低。泥石流危险度与流域面积的大小呈正比关系。     

Spatial and seasonal variability of sediment accumulation potential through controlled flooding of the beels located in the polders of the Ganges-Brahmaputra-Meghna delta of Southwest Bangladesh

The Ganges-Brahmaputra-Meghna (GBM) delta plain within Bangladesh is one of the most vulnerable to relative sea level rise (RSLR) in the world especially under current anthropogenically modified (i.e., embanked) conditions. Tidal river management (TRM) as practiced in coastal regions of Bangladesh may provide an opportunity to combat RSLR by raising the land level through controlled sedimentation inside beels (depression within embanked polders) with re-opening of polders. To date, TRM has been applied to tide-dominated coastal regions, but the potential applicability of TRM for the beels within the polders of river-dominated and mixed flow (MF) regimes remains to be assessed. We apply a calibrated 2D numerical hydromorphodynamic model to quantify sediment deposition in a beel flooded through breaching of the polder dike under conditions of river-dominated, tide-dominated and MF regimes for different seasons and applying different regulation schemes for the flow into the beel. Simulation results show considerable seasonality in sediment deposition with largest deposition during the monsoon season. The potential of controlled flooding is highest in the tide-dominated region, where sediment accumulation can be up to 28 times higher than in the river-dominated region. Regulating flow into a beel increases trapping efficiency, but results in slightly lower total deposition than without regulation. We conclude that re-establishing flooding of the beel within the polder without regulating the flow into the beel through breaching of the polder dike is a promising strategy for the mixed and tide-dominated flow regions in the delta as the sediment accumulation can raise the land surface at a higher rate than RSLR and effective SLR (ESLR). In the more upstream river-dominated section of the delta, accumulation rates would be much lower, but the pressure of sea level rise on these areas is lower as well. Owing to the abundant availability of sediment, application of controlled flooding like TRM therefore provides an opportunity to counteract the impact of RSLR and ESLR by means of land raising, particularly along the tidal river reaches in the GBM delta.     

荆江三口分流变化及三峡水库蓄水影响

长江中游荆江的水沙通过三口洪道分流入洞庭湖,三口分流是荆江-洞庭湖关系调整的驱动因子,为揭示其变化特征及三峡工程运用的影响,基于大量的原观数据,系统研究了近60年三口分流比的变化过程,提出其显著调整特征及诱发因素,引入径流还原计算方法,量化了三峡水库不同运用方式对三口分流量的影响幅度。结果表明:①特大洪水及重大人类活动等诱发因素作用下,三口分流比出现4~5年持续性减小的趋势调整期,之后进入分流比稳定恢复的平衡调整期;②2003-2014年,三峡水库汛前枯水补偿调度使得三口分流量年均增加8.000亿m3,汛后蓄水使得三口分流量年均减小29.00亿m3,对三口分流综合影响量为年均减少21.00亿m3,占同期三口年均分流量的4.29%。     

太行山地形影响下的极端短时强降水分析

2015年8月2日午夜和2011年8月9日前半夜,在两种不同天气系统背景下太行山东麓都出现了小时雨量超过50 mm的极端短时强降水天气,两次过程都是雷暴先在太行山区触发加强,经过下山2 h先后在丘陵站平山和山前平原站石家庄市区产生极端短时强降水。利用常规探测资料、地面加密观测资料、石家庄SA多普勒天气雷达资料,对不同天气系统背景下太行山特殊地形影响的极端短时强降水成因进行分析。结果表明:偏东气流被南北向的太行山地形强迫抬升,且与下山雷暴出流形成中尺度辐合线触发新的雷暴,雷达回波呈现后向传播特征和列车效应造成局地极端短时强降水。太行山地形通过增强辐合上升运动、增大垂直风切变使雷暴下山加强。不同天气系统强迫下,太行山特殊地形对雷暴发展作用不同。在偏西气流引导下,暖区极端短时强降水由阵风锋触发,具有突发性、降水时间短、伴随风力大的特点,下山雷暴出流加快且与山前偏东风的辐合加强,陆续在丘陵区和山前平原触发对流与下山雷暴合并加强造成极端短时强降水;而在东北气流引导下,回流冷锋和阵风锋共同触发的极端短时强降水具有持续时间较长、降雨总量较大、伴随风力较小的特点,太行山东坡对东北冷湿回流有阻挡积聚作用,东北偏北来的雷暴出流边界西端在迎风坡上强迫抬升使雷暴触发并加强,东北气流遇山后发生气旋性偏转使雷暴出流转向东南下山,与平原的偏东风辐合加强,造成丘陵区和山前平原的总降雨时间更长、降雨总量更大。     

Modeled Sverdrup flow in the North Atlantic from 11 different wind stress climatologies

In studies of large-scale ocean dynamics, often quoted values of Sverdrup transport are computed using the Hellerman–Rosenstein wind stress climatology. The Sverdrup solution varies, however, depending on the wind set used. We examine the differences in the large-scale upper ocean response to different surface momentum forcing fields for the North Atlantic Ocean by comparing the different Sverdrup interior/Munk western boundary layer solutions produced by a 1/16° linear numerical ocean model forced by 11 different wind stress climatologies. Significant differences in the results underscore the importance of careful selection of a wind set for Sverdrup transport calculation and for driving nonlinear models. This high-resolution modeling approach to solving the linear wind-driven ocean circulation problem is a convenient way to discern details of the Sverdrup flow and Munk western boundary layers in areas of complicated geometry such as the Caribbean and Bahamas. In addition, the linear solutions from a large number of wind sets provide a well-understood baseline oceanic response to wind stress forcing and thus, (1) insight into the dynamics of observed circulation features, by themselves and in conjunction with nonlinear models, and (2) insight into nonlinear model sensitivity to the choice of wind-forcing product.The wind stress products are evaluated and insight into the linear dynamics of specific ocean features is obtained by examining wind stress curl patterns in relation to the corresponding high-resolution linear solutions in conjunction with observational knowledge of the ocean circulation. In the Sverdrup/Munk solutions, the Gulf Stream pathway consists of two branches. One separates from the coast at the observed separation point, but penetrates due east in an unrealistic manner. The other, which overshoots the separation point at Cape Hatteras and continues to flow northward along the continental boundary, is required to balance the Sverdrup interior transport. A similar depiction of the Gulf Stream is commonly seen in the mean flow of nonlinear, eddy-resolving basin-scale models of the North Atlantic Ocean. An O(1) change from linear dynamics is required for realistic simulation of the Gulf Stream pathway. Nine of the eleven Sverdrup solutions have a C-shaped subtropical gyre, similar to what is seen in dynamic height contours derived from observations. Three mechanisms are identified that can contribute to this pattern in the Sverdrup transport contours. Along 27°N, several wind sets drive realistic total western boundary current transport (within 10% of observed) when a 14 Sv global thermohaline contribution is added (COADS, ECMWF 10 m re-analysis and operational, Hellerman–Rosenstein and National Centers for Environmental Prediction (NCEP) surface stress re-analysis), a few drive transport that is substantially too high (ECMWF 1000 mb re-analysis and operational and Isemer–Hasse) and Fleet Numerical Meteorology and Oceanography Center (FNMOC) surface stresses give linear transport that is slightly weaker than observed. However, higher order dynamics are required to explain the partitioning of this transport between the Florida Straits and just east of the Bahamas (minimal in the linear solutions vs. 5 Sv observed east of the Bahamas). Part of the Azores Current transport is explained by Sverdrup dynamics. So are the basic path of the North Atlantic Current (NAC) and the circulation features within the Intra-Americas Sea (IAS), when a linear rendition of the northward upper ocean return flow of the global thermohaline circulation is added in the form of a Munk western boundary layer.     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---