Joint PDF Modelling of Turbulent Flow and Dispersion in an Urban Street Canyon

The joint probability density function (PDF) of turbulent velocity and concentration of a passive scalar in an urban street canyon is computed using a newly developed particle-in-cell Monte Carlo method. Compared to moment closures, the PDF methodology provides the full one-point one-time PDF of the underlying fields containing all higher moments and correlations. The small-scale mixing of the scalar released from a concentrated source at the street level is modelled by the interaction by exchange with the conditional mean (IECM) model, with a micro-mixing time scale designed for geometrically complex settings. The boundary layer along no-slip walls (building sides and tops) is fully resolved using an elliptic relaxation technique, which captures the high anisotropy and inhomogeneity of the Reynolds stress tensor in these regions. A less computationally intensive technique based on wall functions to represent the boundary layers and its effect on the solution are also explored. The calculated statistics are compared to experimental data and large-eddy simulation. The present work can be considered as the first example of computation of the full joint PDF of velocity and a transported passive scalar in an urban setting. The methodology proves successful in providing high level statistical information on the turbulence and pollutant concentration fields in complex urban scenarios.     

Channel response to a dam-removal sediment pulse captured at high-temporal resolution using routine gage data

In this study, we captured how a river channel responds to a sediment pulse originating from a dam removal using multiple lines of evidence derived from streamflow gages along the Patapsco River, Maryland, USA. Gages captured characteristics of the sediment pulse, including travel times of its leading edge (~7.8 km yr⁻¹) and peak (~2.6 km yr⁻¹) and suggest both translation and increasing dispersion. The pulse also changed local hydraulics and energy conditions, increasing flow velocities and Froude number, due to bed fining, homogenization and/or slope adjustment. Immediately downstream of the dam, recovery to pre-pulse conditions occurred within the year, but farther downstream recovery was slower, with the tail of the sediment pulse working through the lower river by the end of the study 7 years later. The patterns and timing of channel change associated with the sediment pulse were not driven by large flow or suspended sediment-transporting events, with change mostly occurring during lower flows. This suggests pulse mobility was controlled by process-factors largely independent of high flow. In contrast, persistent changes occurred to out-of-channel flooding dynamics. Stage associated with flooding increased during the arrival of the sediment pulse, 1 to 2 years after dam removal, suggesting persistent sediment deposition at the channel margins and nearby floodplain. This resulted in National Weather Service-indicated flood stages being attained by 3–43% smaller discharges compared to earlier in the study period. This study captured a two-signal response from the sediment pulse: (1) short- to medium-term (weeks to months) translation and dispersion within the channel, resulting in aggradation and recovery of bed elevations and changing local hydraulics; and (2) dispersion and persistent longer-term (years) effects of sediment deposition on overbank surfaces. This study further demonstrated the utility of US Geological Survey gage data to quantify geomorphic change, increase temporal resolution, and provide insights into trajectories of change over varying spatial and temporal scales.     

Correction of the ionospheric distortion on the MARSIS surface sounding echoes

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft has made numerous measurements of the Martian surface and subsurface. However, all of these measurements are distorted by the ionosphere and must be compensated before any analysis. We have developed a technique to compensate for the ionospheric distortions. This technique provides a powerful tool to derive the total electron content (TEC) and other higher-order terms of the limited expansion of the plasma dispersion function that are related to overall shape of the electron column profile. The derived parameters are fitted by using a Chapman model to derive ionospheric parameters like n₀, electron density primary peak (maximum for solar zenith angle (SZA) equal 0), and the neutral height scale H.

Our estimated ionospheric parameters are in good agreement with Mars Global Surveyor (MGS) radio-occultation data. However, since MARSIS does not have the observation geometry limitations of the radio occultation measurements, our derived parameters extend over a large range of SZA for each MEX orbit.

The first results from our technique have been discussed by Safaeinili et al. [2007, Estimation of the total electron content of the Martian ionosphere using radar sounder surface echoes. Geophys. Res. Lett. 34, L23204, doi:10.1029/2007GL032154].     

Acoustic-Gravity Waves from Bolide Sources

We have developed a new approach to modeling the acoustic-gravity wave (AGW) radiation from bolide sources. This first effort involves entry modeling of bolide sources that have available satellite data through procedures developed in ReVelle (Earth Moon Planets 95, 441–476, 2004a; in: A. Milani, G. Valsecchi, D. Vokrouhlicky (eds) NEO Fireball Diversity: Energetics-based Entry Modeling and Analysis Techniques, Near-earth Objects: Our Celestial Neighbors (IAU S236), 2007b). Results from the entry modeling are directly coupled to AGW production through line source blast wave theory for the initial wave amplitude and period at (at 10 blast wave radii and perpendicular to the trajectory). The second effort involves the prediction of the formation and or dominance of the propagation of the atmospheric Lamb, edge-wave composite mode in a viscous fluid (Pierce, J. Acoust. Soc. Amer. 35, 1798–1807, 1963) as a function of the source energy, horizontal range and source altitude using the Lamb wave frequency that was deduced directly during the entry modeling and that is used as a surrogate for the source energy. We have also determined that Lamb wave production by bolides at close range decreases dramatically as either the source energy decreases or the source altitude increases. Finally using procedures in Gill (Atmospheric-Ocean Dynamics, 1982) and in Tolstoy (Wave Propagation, 1973), we have analyzed two simple dispersion relationships and have calculated the expected dispersion for the Lamb edge-wave mode and for the excited, propagating internal acoustic waves. Finally, we have used the above formalism to fully evaluate these techniques for four large bolides, namely: the Tunguska bolide of June 30, 1908; the Revelstoke bolide of March 31, 1965; the Crete bolide of June 6, 2002 and the Antarctic bolide of September 3, 2004. Due to page limitations, we will only present results in detail for the Revelstoke bolide.     

基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

胡焕庸线两侧人口的空间分异性及其变化

胡焕庸线作为人地关系的重要地理大发现之一,揭示了中国人口分布东密西疏的重要特征,是人文与经济地理学对社会的重要贡献之一。本文采用地理探测器等定量方法,以人口普查数据为基础,对胡焕庸线两侧人口分布的分异性特征进行了深入讨论,并尝试解释中国人口分布特征变化的原因及其机制。结果表明：① 胡焕庸线作为中国最重要的人口分界线仍基本保持稳定,但从空间统计学分异特征来看该分界线仍可进一步优化调整;② 胡焕庸线两侧内部人口集聚模式发生着明显变化,东南半壁人口分布由改革开放之前和初期的相对均衡状态,转变为以长三角、珠三角等少数区域为中心的集聚模式;而西北半壁人口分布的集聚程度则趋于下降;③ 从人口空间分布的影响因素看,尽管总体上,经济发展类因素对于人口空间分布影响不断增强,但三大阶梯等自然地理本底条件依然对中国人口分布有着重要影响。     

A comparison of results from geomorphological diversity evaluation methods in the Polish Lowland (Toruń Basin and Chełmno Lakeland)

Geomorphological diversity is part of geodiversity. Study and evaluation of geodiversity, including geomorphological diversity, is often conducted in uplands and mountains, despite the fact that lowland areas are of equal importance. This paper evaluates geomorphological diversity in a small area of the Polish Lowland, using a variety of methods that have been applied in recent times for evaluating geodiversity, and presents the results on maps. By comparing these maps and analyzing the correlation coefficients of the results obtained, it was possible to identify the two methods that were best suited to indicating areas with the greatest geomorphological diversity in the lowlands. These two methods are least affected by the choice of elementary fields and data classification methods applied. The study identified the two areas with the greatest relief diversity and showed that they distinctly differ from one another. They demonstrate the major influence of processes, not only on the topographic parameters and landform types, but above all on identifying and defining total geomorphological diversity. These methods, which can be used to identify the areas with the greatest total geomorphological diversity, could readily be used in applied studies relating to abiotic ecosystem services and landscape management.     

An analysis of horizontal dispersion due to the drift current with an Ekman layer

An analytical method for describing horizontal matter dispersion in shear currents is presented using a tensor expression from the point of view that matter dispersion due to the shear effect should be one of the principal mixing dilution processes. Although the behavior of horizontal dispersion is considerably more complicated than common longitudinal dispersion, the present study elucidates the vertical structure of dispersion and the dispersing process from the initial to the stationary stage, besides the usual depth-averaged dispersion coefficient at the stationary stage. As one of the typical applications of horizontal dispersion, dispersion due to the pure drift current with an Ekman layer is examined theoretically using the present method. This examination reveals that the displacement of the centroid and the major axis of dispersion are twisted in the vertical direction more than the direction of the current vector forming the Ekman spiral; that the variance increases in proportion to the third power of the elapsed time; and that the dispersion coefficient at the stationary stage remains constant, independent of the depth normalized by an Ekman layer thickness. Such dependence of the dispersion coefficient in the steady current is shown to be different from that in the oscillatory current, which is inversely proportional to the depth normalized by a Stokes layer thickness. This is considered to be induced by the difference of the vertical profiles of the first order moment in both currents, that is, the shear region of the first order moment is restricted around the floor by the alternation of the current shear in the oscillatory current while it is diffused in the whole depth in the steady current.     

瑞雷波多阶模式频散曲线稀疏正则化反演方法研究

目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型...     

储层中的导波

储层常常表现为低速,在井间地震中往往可以观察到导波的存在,其主要特点是在储层附近能量集中,而在围岩中则能量迅速衰减,同时波形具有频散特征。导波特性对储层内部结构与岩性十分敏感,它有可能作为储层描述和采油过程监测的有效手段,本文针对三层简单波导模型进行了频散方程及曲线的计算,并进行了人工合成槽波记录的研究。