Sensitivity of boundary-layer variables to PBL schemes in the WRF model based on surface meteorological observations, lidar, and radiosondes during the HygrA-CD campaign

Air quality forecast systems need reliable and accurate representations of the planetary boundary layer (PBL) to perform well. An important question is how accurately numerical weather prediction models can reproduce conditions in diverse synoptic flow types. Here, observations from the summer 2014 HygrA-CD (Hygroscopic Aerosols to Cloud Droplets) experimental campaign are used to validate simulations from the Weather Research and Forecasting (WRF) model over the complex, urban terrain of the Greater Athens Area. Three typical atmospheric flow types were identified during the 39-day campaign based on 2-day backward trajectories: Continental, Etesians, and Saharan. It is shown that the numerical model simulations differ dramatically depending on the PBL scheme, atmospheric dynamics, and meteorological parameter (e.g., 2-m air temperature). Eight PBL schemes from WRF version 3.4 are tested with daily simulations on an inner domain at 1-km grid spacing. Near-surface observations of 2-m air temperature and relative humidity and 10-m wind speed are collected from multiple meteorological stations. Estimates of the PBL height come from measurements using a multiwavelength Raman lidar, with an adaptive extended Kalman filter technique. Vertical profiles of atmospheric variables are obtained from radiosonde launches, along with PBL heights calculated using bulk Richardson number. Daytime maximum PBL heights ranged from 2.57 km during Etesian flows, to as low as 0.37 km during Saharan flows. The largest differences between model and observations are found with simulated PBL height during Saharan synoptic flows. During the daytime, campaign-averaged near-surface variables show WRF tended to have a cool, moist bias with higher simulated wind speeds than the observations, especially near the coast. It is determined that non-local PBL schemes give the most agreeable solutions when compared with observations.

[1]  E. Akylas,et al.  Model evaluation of the atmospheric boundary layer and mixed-layer evolution , 2007 .

[2]  I. Troen,et al.  A simple model of the atmospheric boundary layer; sensitivity to surface evaporation , 1986 .

[3]  Jordan G. Powers,et al.  A Description of the Advanced Research WRF Version 2 , 2005 .

[4]  Gerhard Wotawa,et al.  The Lagrangian particle dispersion model FLEXPART-WRF version 3.1 , 2013 .

[5]  Alexandros Papayannis,et al.  Systematic lidar observations of Saharan dust layers over Athens, Greece in the frame of EARLINET project (2004–2006) , 2009 .

[6]  G. Rybicki.  Radiative transfer , 2019, Climate Change and Terrestrial Ecosystem Modeling.

[7]  Alexandros Papayannis,et al.  Inter-comparison of lidar and ceilometer retrievals for aerosol and Planetary Boundary Layer profiling over Athens, Greece , 2011 .

[8]  J. Baldasano,et al.  Performance Evaluation of the Boundary-Layer Height from Lidar and the Weather Research and Forecasting Model at an Urban Coastal Site in the North-East Iberian Peninsula , 2015, Boundary-Layer Meteorology.

[9]  P. Chazette,et al.  Evaluation of the Weather Research and Forecast/Urban Model Over Greater Paris , 2013, Boundary-Layer Meteorology.

[10]  Stuart A. McKeen,et al.  Numerical uncertainty at mesoscale in a Lagrangian model in complex terrain , 2012 .

[11]  A. Holtslag,et al.  A High Resolution Air Mass Transformation Model for Short-Range Weather Forecasting , 1990 .

[12]  J. Pleim.  A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part I: Model Description and Testing , 2007 .

[13]  E. Mlawer,et al.  Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave , 1997 .

[14]  C. Draxl,et al.  Evaluating winds and vertical wind shear from Weather Research and Forecasting model forecasts using seven planetary boundary layer schemes , 2014 .

[15]  Thorsten Mauritsen,et al.  Performance of an Eddy Diffusivity-Mass Flux Scheme for Shallow Cumulus Boundary Layers , 2010 .

[16]  Francesc Rocadenbosch,et al.  Atmospheric Boundary Layer Height Monitoring Using a Kalman Filter and Backscatter Lidar Returns , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[17]  B. Galperin,et al.  ‘Application of a New Spectral Theory of Stably Stratified Turbulence to the Atmospheric Boundary Layer over Sea Ice’ , 2005 .

[18]  G. Danabasoglu,et al.  The Community Climate System Model Version 4 , 2011 .

[19]  Francesc Rocadenbosch,et al.  Retrieval of boundary layer height from lidar using extended Kalman filter approach, classic methods, and backtrajectory cluster analysis , 2014, Remote Sensing.

[20]  J. Dudhia,et al.  A New Vertical Diffusion Package with an Explicit Treatment of Entrainment Processes , 2006 .

[21]  J. Dudhia,et al.  2 A IMPLEMENTATION AND VERIFICATION OF THE UNIFIED NOAH LAND SURFACE MODEL IN THE WRF MODEL , 2003 .

[22]  Laser remote sensing of the atmosphere , 1986 .

[23]  Gelsomina Pappalardo,et al.  EARLINET: the European aerosol research lidar network , 2008 .

[24]  Zavisa Janjic,et al.  The Step-Mountain Coordinate: Physical Package , 1990 .

[25]  A. Holtslag,et al.  Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes, Boundary Conditions and Spin-Up , 2014, Boundary-Layer Meteorology.

[26]  F Rocadenbosch,et al.  Adaptive filter solution for processing lidar returns: optical parameter estimation. , 1998, Applied optics.

[27]  J. Baldasano,et al.  Interactive dust‐radiation modeling: A step to improve weather forecasts , 2006 .

[28]  P. Lacarrére,et al.  Parameterization of Orography-Induced Turbulence in a Mesobeta--Scale Model , 1989 .

[29]  F Rocadenbosch,et al.  Lidar inversion of atmospheric backscatter and extinction-to-backscatter ratios by use of a Kalman filter. , 1999, Applied optics.

[30]  Philip H. Ramsey.  Statistical Methods in the Atmospheric Sciences , 2005 .

[31]  R. F. Brown,et al.  PERFORMANCE EVALUATION , 2019, ISO 22301:2019 and business continuity management – Understand how to plan, implement and enhance a business continuity management system (BCMS).

[32]  Fuqing Zhang,et al.  Evaluation of Three Planetary Boundary Layer Schemes in the WRF Model , 2010 .

[33]  J. Dudhia.  A Nonhydrostatic Version of the Penn State–NCAR Mesoscale Model: Validation Tests and Simulation of an Atlantic Cyclone and Cold Front , 1993 .

[34]  Song-You Hong,et al.  A new stable boundary‐layer mixing scheme and its impact on the simulated East Asian summer monsoon , 2010 .

[35]  J. Dudhia.  Numerical Study of Convection Observed during the Winter Monsoon Experiment Using a Mesoscale Two-Dimensional Model , 1989 .

[36]  O. Jorba,et al.  Influence of the PBL scheme on high-resolution photochemical simulations in an urban coastal area over the Western Mediterranean , 2006 .

[37]  Jimy Dudhia,et al.  Objectively Determined Fair-Weather CBL Depths in the ARW-WRF Model and Their Comparison to CASES-97 Observations , 2013 .

[38]  J. Dudhia,et al.  An evaluation of WRF's ability to reproduce the surface wind over complex terrain based on typical circulation patterns , 2013 .

[39]  C. Zerefos,et al.  From hygroscopic aerosols to cloud droplets: The HygrA-CD campaign in the Athens basin - An overview. , 2017, The Science of the total environment.

[40]  J. Dudhia,et al.  A Revised Approach to Ice Microphysical Processes for the Bulk Parameterization of Clouds and Precipitation , 2004 .

[41]  C. Bretherton,et al.  A New Moist Turbulence Parameterization in the Community Atmosphere Model , 2009 .

[42]  A. Clappier,et al.  An Urban Surface Exchange Parameterisation for Mesoscale Models , 2002 .

[43]  H. Niino,et al.  An Improved Mellor–Yamada Level-3 Model: Its Numerical Stability and Application to a Regional Prediction of Advection Fog , 2006 .

[44]  John S. Kain,et al.  The Kain–Fritsch Convective Parameterization: An Update , 2004 .

[45]  James Correia,et al.  Verification of Convection-Allowing WRF Model Forecasts of the Planetary Boundary Layer Using Sounding Observations , 2013 .

[46]  Richard A. Brown,et al.  Introduction to random signals and applied kalman filtering (3rd ed , 2012 .

[47]  Christos Zerefos,et al.  Boundary layer dynamics in an urban coastal environment under sea breeze conditions , 1995 .

[48]  E. Athanasopoulou,et al.  The Role of Planetary Boundary-Layer Parameterizations in the Air Quality of an Urban Area with Complex Topography , 2009 .

[49]  A. Holtslag,et al.  Improving Stable Boundary-Layer Height Estimation Using a Stability-Dependent Critical Bulk Richardson Number , 2013, Boundary-Layer Meteorology.

[50]  William C. Skamarock,et al.  A time-split nonhydrostatic atmospheric model for weather research and forecasting applications , 2008, J. Comput. Phys..

[51]  Zaviša I. Janić.  Nonsingular implementation of the Mellor-Yamada level 2.5 scheme in the NCEP Meso model , 2001 .

[52]  H. Pan,et al.  Nonlocal Boundary Layer Vertical Diffusion in a Medium-Range Forecast Model , 1996 .

[53]  C. Srinivas,et al.  Mesoscale atmospheric flow-field simulations for air quality modeling over complex terrain region of Ranchi in eastern India using WRF , 2015 .

[54]  Francesc Rocadenbosch,et al.  Atmospheric Boundary Layer Height Estimation Using a Kalman Filter and a Frequency‐Modulated Continuous‐Wave Radar , 2015, IEEE Transactions on Geoscience and Remote Sensing.