High Resolution Simulation of the Variability of Surface Energy Balance Fluxes Across Central London with Urban Zones for Energy Partitioning

The parameterization of surface heat-flux variability in urban areas relies on adequate representation of surface characteristics. Given the horizontal resolutions (e.g. $\approx $ 0.1–1 km) currently used in numerical weather prediction (NWP) models, properties of the urban surface (e.g. vegetated/built surfaces, street-canyon geometries) often have large spatial variability. Here, a new approach based on Urban Zones to characterize Energy partitioning (UZE) is tested within a NWP model (Weather Research and Forecasting model; WRF v3.2.1) for Greater London. The urban land-surface scheme is the Noah/Single-Layer Urban Canopy Model (SLUCM). Detailed surface information (horizontal resolution 1 km) in central London shows that the UZE offers better characterization of surface properties and their variability compared to default WRF-SLUCM input parameters. In situ observations of the surface energy fluxes and near-surface meteorological variables are used to select the radiation and turbulence parameterization schemes and to evaluate the land-surface scheme and choice of surface parameters. For radiative fluxes, improved performance (e.g. $>$ 25 W m $^{-2}$ root-mean-square error reduction for the net radiation) is attained with UZE parameters compared to the WRF v3.2.1 default for all three methods from the simplest to the most detailed. The UZE-based spatial fluxes reproduce a priori expectations of greater energy storage and less evaporation in the dense city centre compared to the residential surroundings. Problems in Noah/SLUCM partitioning of energy between the daytime turbulent fluxes are identified with the overestimation of the turbulent sensible heat and underestimation of the turbulent latent heat fluxes.     

The Effects of Heat Advection on UK Weather and Climate Observations in the Vicinity of Small Urbanized Areas

Weather and climate networks traditionally follow rigorous siting guidelines, with individual stations located away from frost hollows, trees or urban areas. However, the diverse nature of the UK landscape suggests that the feasibility of siting stations that are truly representative of regional climate and free from distorting local effects is increasingly difficult. Whilst the urban heat island is a well-studied phenomenon and usually accounted for, the effect of warm urban air advected downwind is rarely considered, particularly at rural stations adjacent to urban areas. Until recently, urban heat advection (UHA) was viewed as an urban boundary-layer process through the formation of an urban plume that rises above the surface as it is advected. However, these dynamic UHA effects are shown to also have an impact on surface observations. Results show a significant difference in temperatures anomalies (\(p\,< \,0.001\)) between observations taken downwind of urban and rural areas. For example, urban heat advection from small urbanized areas (\(\sim \)1\(\,\hbox {km}^{2}\)) under low cloud cover and wind speeds of 2–3\(\,\hbox {m}\,\hbox {s}^{-1}\) is found to increase mean nocturnal air temperatures by 0.6\(\,^{\circ }\hbox {C}\) at a horizontal distance of 0.5 km. Fundamentally, these UHA results highlight the importance of careful interpretation of long-term temperature data taken near small urban areas.     

Using Temperature Fluctuation Measurements to Estimate Meteorological Inputs for Modelling Dispersion During Convective Conditions in Urban Areas

We examine the performance of several methods to estimate meteorological inputs for modelling dispersion in urban areas during convective conditions. Sensible heat flux, surface friction velocity and turbulent velocities are estimated from measurements of mean wind speed and the standard deviation of temperature fluctuations at a single level on a tower at two suburban sites and at one urban site in Riverside, California. These estimates are compared with observations made at these sites during a field study conducted in 2007. The sensible heat flux is overestimated in the urban area, while it is underestimated at a suburban site when temperature fluctuations are used in the free convection formulation to estimate heat flux. The bias in heat flux estimates can be reduced through a correction that depends on stability. It turns out that the bias in heat flux estimates has a minor effect on the prediction of surface friction velocity and turbulent velocities. Estimates of sensible heat flux, surface friction velocity and turbulent velocities are sensitive to estimates of aerodynamic roughness length, and we suggest estimating the aerodynamic roughness length through detailed micrometeorological measurements made during a limited field study. An examination of the impact of the uncertainty in estimating surface micrometeorology on concentrations indicates that, at small distances from a surface release, ground-level concentrations computed using estimates of heat flux and surface friction compare well with the those based on observed values: the bias is small and the 95% confidence interval of the ratio of the two concentrations is 1.7. However, at distances much larger than the Obukhov length, this confidence interval is close to 2.3 because errors in both friction velocity and heat flux affect plume spread. Finally, we show that using measurements of temperature fluctuations in estimating heat flux is an improvement on that based on the surface energy balance, even when net radiation measurements are available.     

A Backward-Lagrangian-Stochastic Footprint Model for the Urban Environment

Built terrains, with their complexity in morphology, high heterogeneity, and anthropogenic impact, impose substantial challenges in Earth-system modelling. In particular, estimation of the source areas and footprints of atmospheric measurements in cities requires realistic representation of the landscape characteristics and flow physics in urban areas, but has hitherto been heavily reliant on large-eddy simulations. In this study, we developed physical parametrization schemes for estimating urban footprints based on the backward-Lagrangian-stochastic algorithm, with the built environment represented by street canyons. The vertical profile of mean streamwise velocity is parametrized for the urban canopy and boundary layer. Flux footprints estimated by the proposed model show reasonable agreement with analytical predictions over flat surfaces without roughness elements, and with experimental observations over sparse plant canopies. Furthermore, comparisons of canyon flow and turbulence profiles and the subsequent footprints were made between the proposed model and large-eddy simulation data. The results suggest that the parametrized canyon wind and turbulence statistics, based on the simple similarity theory used, need to be further improved to yield more realistic urban footprint modelling.     

Coupling Mesoscale Modelling with a Simple Urban Model: The Lisbon Case Study

The ongoing trend of urbanisation worldwide is leading to a growing requirement for detailed flow and transport parameterisations to be included within numerical weather prediction (NWP) models. Such models often employ a simple roughness parameterisation for urban areas, which is not particularly accurate in predicting or assessing the flow and dispersion at street scale. Moreover, this kind of parameterisation offers too poor a representation of the mechanical and thermal forcing exerted by urban areas on the larger scale flow. At present, high computational costs and long simulation running times are among the constraints for the implementation of more detailed urban sub-models within NWP models. To overcome such limitations, a downscaling procedure from the atmospheric flow at the synoptic scale to the neighbourhood scale and below, is presented in this study. This is achieved by means of a simple urban model based on a parameterised formulation of the drag exerted by the building on the airflow. Application of the urban model for estimating spatially-averaged mean wind speed and the urban heat island over a selected neighbourhood area in Lisbon, Portugal, is presented. The results show the capability of the urban model to provide more accurate mean wind and temperature profiles. Moreover, the urban model has the advantage of being cost effective, as it requires small computational resources, and thus is suitable to be adopted in an operational context. The model is simple enough to be also used to assess how the resolving of urban surface processes may affect those at the larger scales.     

The surface energy balance and the mixing height in urban areas—activities and recommendations of COST-Action 715

The specific problems of determining and simulating the surface energy balance (SEB) and the mixing height (MH) over urban areas are examined. The SEB and MH are critical components of algorithms and numerical models for the urban boundary layer, though the constituent parts of the SEB and the MH are not routinely measured by national weather services. Parameterisations are thus needed in applications. In this investigation, several recently developed algorithms and models for estimating the SEB and MH were applied to new datasets and assessed. Results are discussed in terms of the need for spatial resolution and the parameters needed to describe the urban atmosphere. Limitations of models are identified and recommendations for further development and observations are given. Having identified gaps in knowledge, key findings from new urban experiments and numerical modelling for the SEB and MH are given. The diurnal cycle for the SEB is significantly different from rural conditions—urban heat storage is needed in urban parameterisations. The urban MH is increased over the rural MH, as shown by several numerical schemes and careful sodar analyses. This work has been carried out within the COST-715 Action “Meteorology applied to urban air pollution problems (1998–2004). COST 715 reached a consensus proposing representatively sited measurements of meteorological parameters and turbulent fluxes above roof-tops, and recognised that such data are needed to improve numerical models of the urban surface processes.     

基于Hadoop的数值预报产品服务平台设计与实现

数值预报产品数据与日俱增,采用传统的关系型数据库对其进行存储和管理存在效率低和存储能力不足的问题。另外,基于文件的存储方式在数据存储处理、数据读取和算法计算等方面存在性能瓶颈。针对这一问题,基于Hadoop技术体系设计了分布式的数据存储模型,实现了数值预报产品数据的分布式存储和处理,开发了数值预报产品数据接入处理模块;并实现了基于Rest Web Service的获取数值预报产品要素场数据访问接口、时间序列数据访问接口、数据下载接口等业务应用接口。多业务用户的实际业务测试表明, 该平台在诸如数值预报产品气象数据处理和业务应用方面较传统技术架构具有一定优势。     

Rapid scale-up of negative emissions technologies: social barriers and social implications

We describe a new approach that allows for systematic causal attribution of weather and climate-related events, in near-real time. The method is designed so as to facilitate its implementation at meteorological centers by relying on data and methods that are routinely available when numerically forecasting the weather. We thus show that causal attribution can be obtained as a by-product of data assimilation procedures run on a daily basis to update numerical weather prediction (NWP) models with new atmospheric observations; hence, the proposed methodology can take advantage of the powerful computational and observational capacity of weather forecasting centers. We explain the theoretical rationale of this approach and sketch the most prominent features of a “data assimilation–based detection and attribution” (DADA) procedure. The proposal is illustrated in the context of the classical three-variable Lorenz model with additional forcing. The paper concludes by raising several theoretical and practical questions that need to be addressed to make the proposal operational within NWP centers.     

The Atmospheric Boundary Layer Over Baltic Sea Ice

A new parametrization for the surface energy balance of urban areas is presented. It is shown that this new method can represent some of the important urban phenomena, such as an urban heat island and the occurrence of a near-neutral nocturnal boundary layer with associated positive turbulent heat fluxes, unlike the traditional method for representing urban areas within operational numerical weather prediction (NWP) models. The basis of the new parametrization is simple and can be applied easily within an operational NWP model. Also, it has no additional computational expense compared to the traditional scheme and is hence applicable for operational forecasting requirements. The results show that the errors for London within the Met Office operational mesoscale model have been significantly reduced since the new scheme was introduced. The bias and root-mean-square (rms) errors have been approximately halved, with the rms error now similar to the model as a whole. The results also show that a seasonal cycle still exists in the model errors, but it is suggested that this may be caused by anthropogenic heat sources that are neglected in the urban scheme.The British Crowns right to retain a non-exclusive royalty-free license in and to any copyright is acknowledged.     

Optimizing the Determination of Roughness Parameters for Model Urban Canopies

We present an objective optimization procedure to determine the roughness parameters for very rough boundary-layer flow over model urban canopies. For neutral stratification the mean velocity profile above a model urban canopy is described by the logarithmic law together with the set of roughness parameters of displacement height d, roughness length \(z_0\), and friction velocity \(u_*\). Traditionally, values of these roughness parameters are obtained by fitting the logarithmic law through (all) the data points comprising the velocity profile. The new procedure generates unique velocity profiles from subsets or combinations of the data points of the original velocity profile, after which all possible profiles are examined. Each of the generated profiles is fitted to the logarithmic law for a sequence of values of d, with the representative value of d obtained from the minima of the summed least-squares errors for all the generated profiles. The representative values of \(z_0\) and \(u_*\) are identified by the peak in the bivariate histogram of \(z_0\) and \(u_*\). The methodology has been verified against laboratory datasets of flow above model urban canopies.     

应用MODIS资料对城乡地表热通量的估算

为了揭示城市热岛形成机制,基于MODIS资料,结合自动气象站实测的气象资料,利用地表能量参数化方法估算了地表热通量,分析了城乡地表热通量的空间分布及变化特征.结果表明城乡地气热交换差异明显,与相关文献对比证明该方法是可行、有效的.     

Simulation of a Summer Urban Breeze Over Paris

Numerical simulations for an anticyclonic summer episode in the Paris area have been performed at the meso- scale for a 48-hour period, and compared to observations from a dense operational observational network. The meteorological stations have been classified, according to the extent of urbanization of their surroundings, into four classes (central Paris, urban, suburban, and rural). The atmospheric model, coupled with an urban surface scheme, correctly reproduces the temperature (within 1 K from the observations) and humidity. The intense urban heat island during the night is also well represented.Following the validation, the model is used to quantify atmospheric effects of Paris on the boundary layer, through a comparison with a purely rural simulation. At night, the model simulates a neutral or even slightly unstable boundary layer to a depth of 200 m over the city. In contrast, a very stable layer formed in the countryside. During the day, the boundary layer was more turbulent and 500 m deeper over Paris; vertical velocities of up to 1 m s^-1 were created over the city. This leads to an urban breeze with convergence at low levels (with winds around 5 to 7 m s^-1), and divergence at the boundary-layer top (with similar wind speeds). The horizontal extent of the breeze reaches for more than 50 km from the city centre, and could have an important impact on pollutant diffusion in the area for calm days.Finally, three other spring cases are presented briefly. These show that an urban breeze develops if the synoptic wind is weak enough or disorganized; an urban plume develops otherwise.     

FY-3E: The First Operational Meteorological Satellite Mission in an Early Morning Orbit

Fengyun-3 E(FY-3E),the world’s first early-morning-orbit meteorological satellite for civil use,was launched successfully at the Jiuquan Satellite Launch Center on 5 July 2021.The FY-3E satellite will fill the vacancy of the global early-morning-orbit satellite observation,working together with the FY-3C and FY-3D satellites to achieve the data coverage of early morning,morning,and afternoon orbits.The combination of these three satellites will provide global data coverage for numerical weather prediction(NWP)at 6-hour intervals,effectively improving the accuracy and time efficiency of global NWP,which is of great significance to perfect the global earth observing system.In this article,the background and meteorological requirements for the early-morning-orbit satellite are reviewed,and the specifications of the FY-3E satellite,as well as the characteristics of the onboard instrumentation for earth observations,are also introduced.In addition,the ground segment and the retrieved geophysical products are also presented.It is believed that the NWP communities will significantly benefit from an optimal temporal distribution of observations provided by the early morning,mid-morning,and afternoon satellite missions.Further benefits are expected in numerous applications such as the monitoring of severe weather/climate events,the development of improved sampling designs of the diurnal cycle for accurate climate data records,more efficient monitoring of air quality by thermal infrared remote sensing,and the quasicontinuous monitoring of the sun for space weather and climate.     

Short-term prediction of local wind conditions

Using Numerical Weather Prediction (NWP) models it has been shown that they, combined with models (either physical or statistical) taking local effects into account, can be used to predict the wind locally better than the models commonly used today (as eg persistence). By local is meant at one distinct spot, as eg the location of a meteorological mast. The physical model of local effects takes the following into account: shelter from near-by obstacles, the effect of roughness changes and the effect of the local orography. The large-scale flow is linked to the surface flow by the geostrophic drag law, and the logarithmic wind profile. The predictions are made up to 36 hours ahead. The model is tested on data from 50 meteorological stations scattered all over Europe.     

Characteristics of the Drag Coefficient in the Roughness Sublayer over a Complex Urban Surface

The statistics of momentum exchange in the urban roughness sublayer are investigated. The analysis focuses on the characteristics of the dimensionless friction velocity, \({u_{*}}/U\) , which is defined as the square root of the drag coefficient. The turbulence observations were made at a height of 47 m above the ground on the 325-m meteorological tower, which is located in a very inhomogeneous urban area in Beijing. Under neutral conditions, the dependence of the drag coefficient on wind speed varies with wind direction. When the airflow is from the area of densely built-up buildings, the drag coefficient does not vary with wind speed, while when the airflow is from the area covered by vegetation, the drag coefficient appears to decrease with increasing wind speed. Also, the drag coefficient does not vary monotonically with the atmospheric stability. Both increasing stability and increasing instability lead to the decrease of the drag coefficient, implying that the roughness length and zero-plane displacement may vary in urban areas.     

Pollutant dispersion close to an urban surface – the BUBBLE tracer experiment

Summary In this publication first results of an urban tracer experiment are reported. This experiment was realized in the framework of the Basel UrBan Boundary-Layer Experiment (BUBBLE) in an area with abundant information on turbulence and flow conditions available. Release height was close to roof level and so was the height of the concentration samplers. The meteorological conditions during the experiments were mainly convective, but due to the rough character of the underlying surface also the mechanical turbulence was substantial.The concentration distribution is found to be essentially Gaussian in the horizontal plane and some commonly used methods to estimate the plume widths in applied dispersion models are compared to the observations. From measurements at one site downwind of the source it is found that for a near-roof level source, only an insignificant vertical gradient in tracer concentration is present within a street canyon. Using a Lagrangian Particle Dispersion Model the tracer experiments are simulated. It is shown that the exact form of the parameterization for the flow and turbulence structure within the urban roughness sublayer is of great importance for the simulation results. Also the numerical simulation results underline the necessity (and difficulty) to describe the vertical profile of the dissipation rate of turbulent kinetic energy close to an urban surface.     

A Simple Model for Spatially-averaged Wind Profiles Within and Above an Urban Canopy

This paper deals with the modelling of the flow in the urban canopy layer. It critically reviews a well-known formula for the spatially-averaged wind profile, originally proposed by Cionco in 1965, and provides a new interpretation for it. This opens up a number of new applications for modelling mean wind flow over the neighbourhood scale. The model is based on a balance equation between the obstacle drag force and the local shear stress as proposed by Cionco for a vegetative canopy. The buildings within the canopy are represented as a canopy element drag formulated in terms of morphological parameters such as λ _f and λ _p (the ratios of plan area and frontal area of buildings to the lot area). These parameters can be obtained from the analysis of urban digital elevation models. The shear stress is parameterised using a mixing length approach. Spatially-averaged velocity profiles for different values of building packing density corresponding to different flow regimes are obtained and analysed. The computed solutions are compared with published data from wind-tunnel and water-tunnel experiments over arrays of cubes. The model is used to estimate the spatially-averaged velocity profile within and above neighbourhood areas of real cities by using vertical profiles of λ _f .     

Verification of a One-Dimensional Model of $$\hbox {CO}_{2}$$ Atmospheric Transport Inside and Above a Forest Canopy Using Observations at the Norunda Research Station

A model of \(\hbox {CO}_{2}\) atmospheric transport in vegetated canopies is tested against measurements of the flow, as well as \(\hbox {CO}_{2}\) concentrations at the Norunda research station located inside a mixed pine–spruce forest. We present the results of simulations of wind-speed profiles and \(\hbox {CO}_{2}\) concentrations inside and above the forest canopy with a one-dimensional model of profiles of the turbulent diffusion coefficient above the canopy accounting for the influence of the roughness sub-layer on turbulent mixing according to Harman and Finnigan (Boundary-Layer Meteorol 129:323–351, 2008; hereafter HF08). Different modelling approaches are used to define the turbulent exchange coefficients for momentum and concentration inside the canopy: (1) the modified HF08 theory—numerical solution of the momentum and concentration equations with a non-constant distribution of leaf area per unit volume; (2) empirical parametrization of the turbulent diffusion coefficient using empirical data concerning the vertical profiles of the Lagrangian time scale and root-mean-square deviation of the vertical velocity component. For neutral, daytime conditions, the second-order turbulence model is also used. The flexibility of the empirical model enables the best fit of the simulated \(\hbox {CO}_{2}\) concentrations inside the canopy to the observations, with the results of simulations for daytime conditions inside the canopy layer only successful provided the respiration fluxes are properly considered. The application of the developed model for radiocarbon atmospheric transport released in the form of \(^{14}\hbox {CO}_{2}\) is presented and discussed.     

Simulation of Meteorological Fields Within and Above Urban and Rural Canopies with a Mesoscale Model

Accurate simulation of air quality at neighbourhood scales (on order of 1-km horizontal grid spacing) requires detailed meteorological fields inside the roughness sub-layer (RSL). Since the assumptions of the roughness approach, used by most of the mesoscale models, are unsatisfactory at this scale, a detailed urban and rural canopy parameterisation, called DA-SM2-U, is developed inside the Penn State/NCAR Mesoscale Model (MM5) to simulate the meteorological fields within and above the urban and rural canopies. DA-SM2-U uses the drag-force approach to represent the dynamic and turbulent effects of the buildings and vegetation, and a modified version of the soil model SM2-U, called SM2-U(3D), to represent the thermodynamic effects of the canopy elements. The turbulence length scale is also modified inside the canopies. SM2-U(3D) assesses the sensible and latent heat fluxes from rural and urban surfaces in each of the computational layers inside the canopies by considering the shadowing effect, the radiative trapping by the street canyons, and the storage heat flux by the artificial surfaces. DA-SM2-U is tested during one simulated day above the city of Philadelphia, U.S.A. It is shown that DA-SM2-U is capable of simulating the important features observed in the urban and rural RSL, as seen in the vertical profiles of the shear stress, turbulent kinetic energy budget components, eddy diffusivity, potential air temperature, and specific humidity. Within the canopies, DA-SM2-U simulates the decrease of the wind speed inside the dense canopies, the skirting of the flow around the canopy blocks, warmer air inside the vegetation canopy than above open areas during the night and conversely during the day, and constantly warmer air inside the urban canopy. The comparison with measurements shows that the surface air temperature above rural and urban areas is better simulated by DA-SM2-U than by the `standard version' of MM5.     

The Impact of Land-Surface Parameter Properties and Resolution on the Simulated Cloud-Topped Atmospheric Boundary Layer

A method to simulate characteristics of wind speed in the boundary layer of tropical cyclones in an idealized manner is developed and evaluated. The method can be used in a single-column modelling set-up with a planetary boundary-layer parametrization, or within large-eddy simulations (LES). The key step is to include terms in the horizontal velocity equations representing advection and centrifugal acceleration in tropical cyclones that occurs on scales larger than the domain size. Compared to other recently developed methods, which require two input parameters (a reference wind speed, and radius from the centre of a tropical cyclone) this new method also requires a third input parameter: the radial gradient of reference wind speed. With the new method, simulated wind profiles are similar to composite profiles from dropsonde observations; in contrast, a classic Ekman-type method tends to overpredict inflow-layer depth and magnitude, and two recently developed methods for tropical cyclone environments tend to overpredict near-surface wind speed. When used in LES, the new technique produces vertical profiles of total turbulent stress and estimated eddy viscosity that are similar to values determined from low-level aircraft flights in tropical cyclones. Temporal spectra from LES produce an inertial subrange for frequencies \(\gtrsim \)0.1 Hz, but only when the horizontal grid spacing \(\lesssim \)20 m.