污染天气下成都城市热岛环流结构的数值模拟

本文利用中尺度模式WRF V3.9对2016年7月16日成都一次污染天气下的城市热岛环流个例进行了数值模拟,分析了城市热岛环流的三维结构及演变特征.并通过设计减少气溶胶光学厚度的敏感性试验,研究了污染对城市热岛环流的可能影响.结果表明:当地时间17∶00城市热岛环流开始形成.随着热岛强度增强,环流增强,城乡边界处的城市风锋不断向城市中心推进.19∶00热岛环流结构最显著.21∶00环流结构被破坏,仅低层存在微弱的乡村风.在热岛环流结构最显著时刻,近地面风场由郊区向城市辐合,地面以上2.0km处风场由城市向乡村辐散,辐散中心与辐合中心位置大致对应.此时城市风锋在城市处合并,环流的水平尺度约为城市尺度的2～2.5倍.当气溶胶光学厚度减小后,城市热岛环流尺度和强度以减小为主.特别是在热岛环流最显著时刻,低层乡村风风速减小,城市风锋最大上升速度降低,环流的水平尺度在西、南、北三个方向均减小,且高空回流高度降低.可能的影响机制是,气溶胶光学厚度减小后,净辐射通量增大,城乡地表能量通量差异增大,城市边界层高度升高.但城市边界层高度升高对城市热岛的抑制作用超过了城乡地表能量通量差异增加对城市热岛的增...     

通过地震成像分辨率研究近地表复杂介质散射对成像质量的影响

近地表广泛分布的小尺度非均匀介质严重影响地震数据和成像质量.本文引入参数化的随机介质来描述近地表非均匀介质;借助功率谱、相关长度、均方根扰动等随机统计特征量研究近地表非均匀介质层对地震波传播、散射及成像的影响;通过成像过程的点弥散函数表述成像分辨率.利用数值方法建立模型的随机统计特征量与成像分辨率特征量之间的关系,以定量或半定量的方式研究二者之间的相关性.数值计算结果表明,近地表非均匀介质层的厚度、速度扰动的幅度以及非均匀体的尺度等都对地震数据质量和成像品质有明显的影响.以全数值方法建立起关于近地表非均匀介质复杂性和地下成像质量之间的关系,有望成为研究近地表非均匀介质对地震数据采集和成像分辨率影响的有效工具.     

长白山天池火山区介质非均匀性

高频S波随着传播距离的增大其均方根(RMS)包络逐渐变宽,我们把这种现象称为S波包络展宽现象.S波传播路径上随机分布的非均匀体对S波的多次前向散射和绕射作用是导致S波包络展宽现象的主要原因,因此可用S波包络展宽现象来研究介质非均匀性.本文采用S波包络峰值延时来对S波包络展宽现象进行量化.S波包络峰值延时定义为S波初至与其均方根包络峰值最大值之间的时间差.本文选用长白山天池火山区的小震记录,运用S波包络峰值延时对长白山天池火山口地区的介质非均匀性进行了评价.结果发现长白山天池火山区呈现强烈的介质非均匀性,在0~2 km深度范围内介质非均匀性表现出南部强北部弱的空间分布特征;在2~5 km深度范围内介质非均匀性的空间分布特征与频率具有相关性;天池火山区介质非均匀性具有明显的多尺度特性;强弱非均匀性接触带,往往是地震频发地带.根据地震与非均匀体在空间分布的相对位置,我们认为火山区介质非均匀性可能反映了火山早期喷发堆积物介质结构的差异.     

基于尺度分析的断层热信息遥感图像增强方法——以江山—绍兴断裂金衢段为例

根据热红外遥感影像上断层热信息具有特征几何尺寸的特点,以断层系统热信息分析为目标,提出了一种基于尺度分析的断层热信息遥感图像增强方法.在江山—绍兴断裂金衢段的实际工作中,通过断层两侧地表高温区域尺度分析、特征尺度网格抽样和样本插值成图等步骤,有效地降低了背景干扰,客观地描述了研究区与断层相关的热信息的空间分布形态及特征.多种尺度分析结果表明,在9 km2特征尺度上,断层热信息特征规律明显：地表高温区域沿北东走向的江山—绍兴断裂带和常山—漓渚断裂带两侧分布,呈线性特征;在淳安—温州断裂带与衢州—天台断裂带交叉位置地表温度较高.研究结果经实测资料验证,基本特征与实测资料相符.     

长江三角洲地区春季低空大气臭氧垂直分布特征

介绍分析了2001年3月3日~4月13日浙江临安臭氧探空观测5 km以下臭氧垂直分布特征. 结果表明, 臭氧浓度垂直分布与湿球位温、风场有密切的关系. 臭氧浓度在2 km以下变化幅度很大, 明显的东风分量伴随臭氧高值. 5 km以下臭氧垂直分布可以分为峰值型、均匀型、分层结构型、低空污染型和线性增长型5个基本类型. 此外, 还分析了3种情形下区域尺度输送对低空污染型臭氧分布的影响.     

非均匀灌溉棉田能量平衡特征研究

运用国际能量平衡实验(EBEX-2000)的湍流、净辐射和土壤观测资料,运用涡动相关法分析了非均匀灌溉引起的热内边界层发展条件下近地层感热、潜热通量特征,并对有无灌溉两种条件下的能量闭合度进行了对比分析.在计算感热、潜热通量过程中,分别将Schotanus订正和Webb订正纳入了考虑范围,研究了两种订正方法对计算湍流热通量的影响.研究结果发现,由于非均匀灌溉生成的热内边界层使得近地层感热通量受到抑制,潜热通量出现波动,该现象在8.7 m比2.7 m 更为显著.非均匀灌溉导致的热内边界层的存在使得近地层能量闭合度偏低,能量平衡比率约为0.65;而没有热内边界层存在时,近地层能量平衡比率约为0.70.本实验中,Schotanus订正使得感热通量显著减小,其订正量日平均值约为-8 W/m²,占净辐射的近4%;Webb订正量日平均值约为2 W/m²,对能量平衡的影响较小.     

一次冷锋过境过程的中尺度通量观测

根据Taylor假定 ,通过对铁塔定点观测冷锋过境湍流资料的谱分析 ,分离出其中尺度过程 ,从而计算出其中尺度通量 .计算结果的分析表明 ,在强背景风场条件下 ,湍流动量输送占据绝对优势 ;当背景风较弱时 ,中尺度动量通量不能被忽视 ,它甚至强于湍流动量通量 .而中尺度感热通量的强弱受多种因素的制约 .总的说来 ,冷锋后的中尺度感热通量大大强于湍流感热通量 .对于湍流通量参数化计算的理论分析表明 ,在较长的时间尺度进行湍流通量参数化时 ,有必要考虑因中尺度扰动而产生的修正 .合理的通量计算参数化方案需要全面包含湍流通量、中尺度通量以及中尺度过程对湍流通量的修正     

京津唐张地区地壳上地幔三维P波速度结构

采用作者最近新提出的频谱参数化方法,对北京电信传输台网2447个远震P波到时残差数据进行了反演,说明新的参数化方法能有效地应用于实际数据.反演结果表明,京津唐张地区深部100多km内存在明显的速度横向不均匀性;50 km以上部分的结构与地表的地质构造有较好的对应关系:在燕山隆起区下部表现为高速异常区,而凹陷区下部则呈现低速异常区;深度100 km以下,地表构造的痕迹逐步消失,异常蝠度也相应降低.研究区域范围内存在三个大的低速异常块体,其中以天津、唐山附近的水平尺度为100 km左右的低速异常区幅度最大,向下一直延伸到约100 km的深度,而且在100 km以下存在一低速异常构造与之对应.研究区域内发生的强震大部分都位于低速异常区边缘向高速区过渡的梯度带上.     

泉州盆地地震效应的二维非线性分析

基于ABAQUS显式有限元并行计算平台,采用大尺度二维精细化有限元非线性分析方法,研究泉州盆地地震效应特征.结果表明:(1)与基岩输入地震动相比,地表峰值加速度总体呈放大效应,且盆地不同位置的放大效应存在明显差异;(2)地表或基岩剧烈起伏及土层横向分布极不均匀处,地震动易产生局部聚集效应,相应地表地震动呈显著放大或缩小效应;地表加速度反应谱产生双峰甚至多峰现象.(3)在土层竖向、横向分布不均匀处,峰值加速度沿深度方向呈非单调递减,加速度反应谱沿深度存在较大差异.场地地震效应的二维非线性分析法能在一定程度上反映特殊场地条件对地震动效应的影响,该结果可为泉州盆地及类似盆地的场地设计地震动参数的确定提供参考依据.     

鄱阳湖地区大气边界层特征的数值模拟

采用WRFV2.2中尺度数值模式对鄱阳湖地区200 km×200 km范围内,2009年11月5日00∶00至2009年11月6日12∶00不同高度的气象要素进行了数值模拟,得到了水平分辨率为1 km的鄱阳湖地区大气边界层风、温、湿度场和廓线分布的大气边界层物理特征.模拟结果发现:白天鄱阳湖面上空存在着冷岛效应并伴随湖风,而夜间湖面上空存在着热岛效应并伴随陆风,湖面与陆地之间最大温差可达6 ℃;同时地形以及下垫面类型对鄱阳湖区风场的分布具有很大影响,夜间存在一条东北西南走向的低空辐合带,白天逐渐消失;此外受风场和地形作用湖面上空的湿度分布也不均匀,白天湿度层厚度低而夜晚湿度层厚度高,湖中心右侧湿度值大于左侧湿度值.模拟结果能较好地反映鄱阳湖的大气边界层物理特征,有助于了解鄱阳湖地区区域气候的特点,以及由于地形、地理环境、地表特征所形成的不同高度上的风、温、湿的分布规律和大气边界层物理特征,为鄱阳湖地区局地天气预报、风能资源开发、环境保护等提供了科学依据.     

Scaling effects on modeled surface energy-balance components using the NOAH-OSU land surface model

As surface exchange processes are highly non-linear and heterogeneous in space and time, it is important to know the appropriate scale for the reasonable prediction of these exchange processes. For example, the explicit representation of surface variability has been vital in predicting mesoscale weather events such as late-afternoon thunderstorms initiated by latent heat exchanges in mid-latitude regions of the continental United States. This study was undertaken to examine the effects of different spatial scales of input data on modeled fluxes, so as to better understand the resolution needed for accurate modeling. A statistical procedure was followed to select two cells from the Southern Great Plains 1997 hydrology experiment region, each 20 km×20 km, representing the most homogeneous and the most heterogeneous surface conditions (based on soil and vegetation) within the study region. The NOAH-OSU (Oregon State University) Land Surface Model (LSM) was employed to estimate surface energy fluxes. Three scales of study (200 m, 2 and 20 km) were considered in order to investigate the impacts of the aggregation of input data, especially soil and vegetation inputs, on the model output. Model results of net radiation and latent, sensible and ground heat fluxes were compared for the three scales. For the heterogeneous area, the model output at the 20-km resolution showed some differences when compared with the 200-m and 2-km resolutions. This was more pronounced in latent heat (12% decrease), sensible heat (22% increase), and ground heat flux (44% increase) estimation than in net radiation. The scaling effects were much less for the relatively homogeneous land area with 5% increase in sensible heat and 4% decrease in ground heat flux estimation. All of the model outputs for the 2- and 20-km resolutions were in close agreement. The results suggested that, for this study region, soils and vegetation input resolution of about 2 km should be chosen for realistic modeling of surface exchange processes. This resolution was sufficient to capture the effects of sub-grid scale heterogeneity, while avoiding the data and computational difficulties associated with higher spatial resolutions.     

Natural integration of scalar fluxes from complex terrain

Large eddy simulations of turbulent flow and transport in the atmospheric boundary layer were conducted over heterogeneous sources of heat and water vapor to identify the blending properties of the turbulent mixing in an unstably stratified boundary layer. The numerical simulations show that the concept of blending in the ABL is in fact a useful one, even under convective conditions, for a range of surface conditions. Since the transport eddies that are responsible for the blending have sizes that are constrained by the boundary layer depth, and since the vertical motion is so important under the unstable density stratification studied here, we see that when the length scales of the source variability on the land surface become significantly greater than the ABL depth the blending is lost. In this case the source fields remain relatively uncoupled by the important eddy motion. However, for smaller surface length scales, the dynamic eddy motion couples the surface patches. Hence, there is good reason that the land surface exchange phenomenon would not be scale invariant over the entire range of scales. Because of the active role of temperature the effects of inhomogeneous surface sources of sensible heat persist higher into the ABL than do the effects of surface sources from more passive scalars, such as water vapor. Moreover, the mean fields of potential temperature and specific humidity blend at much lower heights than do the vertical turbulent flux fields of these two scalars. We propose a useful measure of blending efficiency for simulation studies and show how this bridges from the dynamics responsible for the blending to the horizontal homogeneity of scalar flux fields at measurement heights in the ABL.     

Patch scale aggregation of heterogeneous land surface cover for mesoscale meteorological models

Theoretical studies dealing with aggregation of surface parameters at small scale are reviewed. Finding effective parameters for surface resistance is possible for most cases by taking simple geometric or arithmetic averages of the component resistances. The use of more sophisticated techniques such as the blending height improves the calculations. Resistances for heat and water vapour behave differently in heterogeneous terrain. A simple surface energy balance model is adapted to show the behaviour of the roughness length of heat and water vapour in heterogeneous terrain. It is suggested that this simple parameterization can adequately take into account the effect of variation in surface cover on the fluxes of heat and water vapour.     

Multiscale analysis of vegetation surface fluxes: from seconds to years

The variability in land surface heat (H), water vapor (LE), and CO₂ (or net ecosystem exchange, NEE) fluxes was investigated at scales ranging from fractions of seconds to years using eddy-covariance flux measurements above a pine forest. Because these fluxes significantly vary at all these time scales and because large gaps in the record are unavoidable in such experiments, standard Fourier expansion methods for computing the spectral and cospectral statistical properties were not possible. Instead, orthonormal wavelet transformations are proposed and used. The are ideal at resolving process variability with respect to both scale and time and are able to isolate and remove the effects of missing data (or gaps) from spectral and cospectral calculations. Using the spectra, we demonstrated unique aspects in three appropriate ranges of time scales: turbulent time scales (fractions of seconds to minutes), meteorological time scales (hour to weeks), and seasonal to interannual time scales corresponding to climate and vegetation dynamics. We have shown that: (1) existing turbulence theories describe the short time scales well, (2) coupled physiological and transport models (e.g. CANVEG) reproduce the wavelet spectral characteristics of all three land surface fluxes for meteorological time scales, and (3) seasonal dynamics in vegetation physiology and structure inject strong correlations between land surface fluxes and forcing variables at monthly to seasonal time scales. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the flux variability is dominated by a narrow band or bands of modes, and these modes “resonate” with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE.     

Intercomparison of oceanic and atmospheric forced and coupled mesoscale simulations. Part I: Surface fluxes

A mesoscale non-hydrostatic atmospheric model has been coupled with a mesoscale oceanic model. The case study is a four-day simulation of a strong storm event observed during the SEMAPHORE experiment over a 500 × 500 km² domain. This domain encompasses a thermohaline front associated with the Azores current. In order to analyze the effect of mesoscale coupling, three simulations are compared: the first one with the atmospheric model forced by realistic sea surface temperature analyses; the second one with the ocean model forced by atmospheric fields, derived from weather forecast re-analyses; the third one with the models being coupled. For these three simulations the surface fluxes were computed with the same bulk parametrization. All three simulations succeed well in representing the main oceanic or atmospheric features observed during the storm. Comparison of surface fields with in situ observations reveals that the winds of the fine mesh atmospheric model are more realistic than those of the weather forecast re-analyses. The low-level winds simulated with the atmospheric model in the forced and coupled simulations are appreciably stronger than the re-analyzed winds. They also generate stronger fluxes. The coupled simulation has the strongest surface heat fluxes: the difference in the net heat budget with the oceanic forced simulation reaches on average 50 Wm⁻² over the simulation period. Sea surface-temperature cooling is too weak in both simulations, but is improved in the coupled run and matches better the cooling observed with drifters. The spatial distributions of sea surface-temperature cooling and surface fluxes are strongly inhomogeneous over the simulation domain. The amplitude of the flux variation is maximum in the coupled run. Moreover the weak correlation between the cooling and heat flux patterns indicates that the surface fluxes are not responsible for the whole cooling and suggests that the response of the ocean mixed layer to the atmosphere is highly non-local and enhanced in the coupled simulation.     

Mesoscale fluxes over heterogeneous flat landscapes for use in larger scale models

We have developed a framework to parameterize mesoscale fluxes of heat, moisture, and momentum, in a cloud-free environment, which result from heterogeneous heating of flat land surfaces. The importance and parameterizability of these mesoscale fluxes is demonstrated using the mathematical concept of predictability. This methodology is used to estimate the relative importance of mesoscale, as contrasted with turbulent fluxes, in the Konza Prairie of Kansas during the FIFE field experiment.     

Mapping mesoscale currents by optimal interpolation of satellite radiometer and altimeter data

 Mapping the mesoscale surface velocity stream function by combining estimates of surface height from satellite altimetry and surface currents from sequential infrared (sea-surface temperature) imagery using optimal interpolation is described. Surface currents are computed from infrared images by the method of maximum cross-correlations (MCC) and are combined with altimeter sea-level anomaly data from the TOPEX/Poseidon and ERS satellites. The analysis method was applied to 6 years of data from the East Australian Current region. The covariance of velocity and sea-level data is consistent with the statistical assumptions of homogeneous, isotropic turbulence, with typical length scales of order 220 km and time scales of 10 days in this region. Augmenting the analysis of altimeter data with MCC velocity observations improves the resolution of the surface currents, especially near the Australian coast, and demonstrates that the two data sources provide consistent and complementary observations of the surface mesoscale circulation. The volume of MCC data is comparable to that from a satellite altimeter, but with a more variable distribution of spatial and temporal resolution. In concert with altimetry, satellite radiometer velocimetry represents a technique useful for retrospective analysis of currents from high-resolution satellite radiometer data-sets. Received: 3 July 2001 / Accepted: 16 November 2001     

Hydrodynamic modelling of mesoscale eddies in the Black Sea

A three dimensional structure of mesoscale circulation in the Black Sea is simulated using the Proudman Oceanographic Laboratory Coastal Ocean Modelling System. A number of sensitivity tests reveal the response of the model to changes in the horizontal resolution, time steps, and diffusion coefficients. Three numerical grids are examined with x-fine (3.2 km), fine (6.7 km) and coarse (25 km) resolution. It is found that the coarse grid significantly overestimates the energy of the currents and is not adequate even for the study of basin-scale circulation. The x-fine grid, on the other hand, does not give significant advantages compared to the fine grid, and the latter is used for the bulk of simulations. The most adequate parameters are chosen from the sensitivity study and used to model both the basin-scale circulation and day-to-day variability of mesoscale currents for the months of May and June of 2000. The model is forced with actual wind data every 6 h and monthly climatic data for evaporation, precipitation, heat fluxes and river run-off. The results of the fine grid model are compared favourably against the satellite imagery. The model adequately reproduces the general circulation and many mesoscale features including cyclonic and anticyclonic eddies, jets and filaments in different parts of the Black Sea. The model gives a realistic geographical distribution and parameters of mesoscale currents, such as size, shape and evolution of the eddies.     

Spatio‐temporal variability of water and energy fluxes – a case study for a mesoscale catchment in pre‐alpine environment

Water and energy fluxes at and between the land surface, the subsurface and the atmosphere are inextricably linked over all spatio‐temporal scales. Our research focuses on the joint analysis of both water and energy fluxes in a pre‐alpine catchment (55 km²) in southern Germany, which is part of the Terrestrial Environmental Observatories (TERENO). We use a novel three‐dimensional, physically based and distributed modelling approach to reproduce both observed streamflow as an integral measure for water fluxes and heat flux and soil temperature measurements at an observation location over a period of 2 years. While heat fluxes are often used for comparison of the simulations of one‐dimensional land surface models, they are rarely used for additional validation of physically based and distributed hydrological modelling approaches. The spatio‐temporal variability of the water and energy balance components and their partitioning for dominant land use types of the study region are investigated. The model shows good performance for simulating daily streamflow (Nash–Sutcliffe efficiency > 0.75). Albeit only streamflow measurements are used for calibration, the simulations of hourly heat fluxes and soil temperatures at the observation site also show a good performance, particularly during summer. A limitation of the model is the simulation of temperature‐driven heat fluxes during winter, when the soil is covered by snow. An analysis of the simulated spatial fields reveals heat flux patterns that reflect the distribution of the land use and soil types of the catchment. The water and energy partitioning is characterized by a strong seasonal cycle and shows clear differences between the selected land use types. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Impact of soil moisture heterogeneity length scale and gradients on daytime coupled land‐cloudy boundary layer interactions

Using a coupled large‐eddy simulation–land surface model framework, the impact of two‐dimensional soil moisture heterogeneity on the cloudy boundary layer under varied free‐atmosphere stabilities is investigated. Specifically, the impacts of soil moisture heterogeneity length scale and heterogeneity in terms of soil moisture gradients on micrometeorological states, surface fluxes, boundary layer characteristics, and cloud development are examined. The results show that mesoscale circulations due to surface heterogeneity in soil moisture play an important role in transferring water vapour within the boundary layer and in regulating cloud distribution at the entrainment zone, which, in turn, provides feedbacks on boundary layer/surface energy budgets. The initial domain‐averaged soil moisture is identical for all homogenous and heterogeneous cases; however, the soil moisture heterogeneity in gradient and length scale between dry and wet regions has a significant impact on the estimates of near‐surface micrometeorological properties and surface fluxes, which further affect the boundary layer states and characteristics. Both liquid water potential temperature and liquid water mixing ratio increase with an increasing soil moisture gradient, whereas the amount of specific humidity decreases. Heterogeneity length scale and free atmosphere stability also amplify these impacts on the boundary layer structure and cloud formation. In a low atmospheric stability condition that potentially allows for a deeper boundary layer and a higher entrainment rate, cloud base height and cloud thickness significantly increase as the soil moisture gradient and length scale increase. Analysis to differentiate the influences of surface heterogeneity type (i.e. length scale vs gradient) shows that in general soil moisture gradient provides a larger impact than heterogeneity length scale, although the heterogeneity length scale is large enough to initiate circulation features responsible for differences in the coupled system between homogeneous and heterogeneous soil moisture cases. Copyright © 2012 John Wiley & Sons, Ltd.