SURF: understanding and predicting urban convection and haze

Transcription

1 SURF: understanding and predicting urban convection and haze Article Liang, X., Miao, S., Li, J., Bornstein, R., Zhang, X., Gao, Y., Cao, X., Chen, F., Cheng, Z., Clements, C., Dabberdt, W., Ding, A., Ding, D., Dou, J. J., Dou, J. X., Dou, Y., Grimmond, C. S. B., Gonzalez Cruz, J., He, J., Huang, M., Huang, X., Ju, S., Li, Q., Niyogi, D., Quan, J., Sun, J., Sun, J. Z., Yu, M., Zhang, J., Zhang, Y., Zhao, X., Zheng, Z. and Zhou, M. (2018) SURF: understanding and predicting urban convection and haze. Bulletin of the American Meteorological Society, 99 (7). pp ISSN doi: D Available at It is advisable to refer to the publisher s version if you intend to cite from the work. To link to this article DOI: D Publisher: American Meteorological Society All outputs in CentAUR are protected by Intellectual Property Rights law, including copyright law. Copyright and IPR is retained by the creators or other copyright holders. Terms and conditions for use of this material are defined in

2 the End User Agreement. CentAUR Central Archive at the University of Reading Reading s research outputs online

3 Liang X, S Miao, J Li, R Bornstein, X Zhang, Y Gao, X Cao, F Chen, Z Cheng, C Clements, W Dabberdt, A Ding, D Ding, J Dou, J Dou, Y Dou, C Grimmond, J Gonzalez-Cruz, J He, M Huang, X Huang, S Ju, Q Li, D Niyogi, J Quan, J Sun, J Sun, M Yu, J Zhang, Y Zhang, X Zhao, Z Zheng, and M Zhou, 2018: SURF: Understanding And Predicting Urban Convection And Haze. Bull. Amer. Meteor. Soc. doi: /bams-d , in press. Comprehensive multi-scale observations and modeling are enabling an international team of scientists to understand and better predict urban effects on heavy summertime convective precipitation and wintertime fine particle pollutant episodes in the Greater Beijing Metro Area. Abstract Urbanization modifies atmospheric energy and moisture balances, forming distinct features, e.g., urban heat islands (UHIs) and enhanced or decreased precipitation. These produce significant challenges to science and society, including rapid and intense flooding, heat waves strengthened by UHIs, and air pollutant haze. The Study of Urban-Impacts on Rainfall and Fog/Haze (SURF) has brought together international expertise on observations and modeling, meteorology and atmospheric chemistry, and research and operational forecasting. The SURF overall science objective is a better understanding of urban, terrain, convection, and aerosol interactions for improved forecast accuracy. Specific objectives include: a) promoting cooperative international research to improve understanding of urban summer convective precipitation and winter particulate episodes via extensive field studies; b) improving high-resolution urban weather and air quality forecast-models; and c) enhancing urban weather forecasts for societal applications, e.g., health, energy, hydrologic, climate change, air quality, planning, and emergency-response management. Preliminary SURF observational and modeling results are shown, i.e., turbulent PBL structure, bifurcating thunderstorms, haze events, urban canopy model development, and model forecast evaluation. PROGRAM SYNOPSIS The Study of Urban-Impacts on Rainfall and Fog/Haze (SURF) Project investigates urban, terrain, convection, and aerosol interactions for improved weather and air quality forecast accuracy achieved by the combination of expertise from an international group of multidisciplinary scientists and by the employment of complementary research methodologies. Principal Chinese participants were from the Institute of Urban Meteorology (IUM), Environmental Meteorology Forecast Center (EMFC), Chengdu University of Information Technology (CUIT), Beijing Weather Modification Office (BJWMO), and National Marine Environment Forecast Center (NMERC). Principle international groups include San Jose State University (SJSU), National Center for Atmospheric Research (NCAR), Vaisala Group (VG), University of Reading (UR), City College of New York (CCNY), and Purdue University (PU). SURF consists of four symbiotic components: Observations with research-level boundary layer (BL) instrumentation to augment extensive operational surface, boundary layer, and upper air networks during SURF summer thunderstorm and winter haze field study periods. Resulting data are archived, quality controlled, and available to Chinese and international partners. Numerical model development and verification activities, which involved use of SURF observations to develop and test advanced modules to improve urban parameterizations in the operational and research-grade urbanized Weather Research and Forecasting (WRF) model (Skamarock et al. 2008) and in air quality models. New parameterizations are available for inclusion in publically available versions of the WRF model. Modeling of case studies that involves use of the newly upgraded research grade weather and air quality numerical models to simulate SURF golden-case thunderstorms and haze events, whose model input and output data are available to interested modelers. Applications with the newly upgraded operational numerical models for improved IUM urban weather and air quality forecasts and for use by health, energy, hydrologic, climate change, air quality, planning, and emergency response managers. The order is rapidly fadin And the first one now will later be last For the times they are a-changin Bob Dylan 1 Global urbanization has accelerated over the past 65 years, in 1950, 30% (746 million people) of the global population resided in cities, but by 2014 that number had increased to 54% (3.9 billion), and the United Nations projects a value of 66% (6.4 billion) by 2050 (UN 2014). Half of all urbanites currently live in small cities, while 28 mega-cities (> 10 million) already have 13% of the global population; and by 2030, 41 mega-cities will exist. Urbanization challenges society and science, overcrowding, infrastructure, and economic activities make cities vulnerable to disasters like flooding and haze, according to an UN-Habitat (2016) report. The High Impact Weather Program of the World Meteorological Organization (WMO 2013) lists urban-induced floods as well as heat waves and air pollution as two of its five major foci. US urban heat waves produce more weather-related deaths than any other weather phenomenon (NRC 2010). More than 700 deaths occurred from a 1995 five-day Chicago heat wave, while a 2015 event in India resulted in nearly deaths (Di Liberto 2015). The 14-day European August 2003 heat wave claimed lives (Robine et al. 2007), with a large number in cities (Vandentorren et al. 2004). Heavy rainfall caused flooding in July 2012, killing 79 people in Beijing and with losses of US $1.7B. 1 From The Times They Are A-Changin copyright 1963, 1964 by Warner Bros. Inc..; renewed 1991,1992 by Special Rider Music - 1 -

4 Urban areas also enhance the severity of gaseous and (solid and liquid) aerosol pollution, producing photochemical haze, smog, and acid rain episodes (Cao et al. 2013). One billion people experience excessive outdoor pollution, and during 2012 air pollution was linked to one of eight deaths globally (UNICEP 2016). Urban air pollution costs approximately 2% of GDP in developed countries and 5% in developing countries. Anthropogenic aerosols penetrate deep into human organs (Li et al. 2011), and cause severe health consequences, e.g., asthma, heart attacks, nervous system damage, and lung cancer (Tie et al. 2009). Aerosol impacts on weather and climate include alterations of: radiation fluxes (Andreae et al. 2005), horizontal visibility (Quan et al. 2015), atmospheric stability (Zhao et al. 2006), planetary boundary layer (PBL) development (Quan et al. 2013), cloud formation (Andreae et al. 2004), and convective precipitation (Rosenfeld et al. 2008). While significant knowledge concerning aerosol impacts on some aspects of weather systems has been obtained, their impacts on radiation (Ao et al. 2016; Wang et al. 2016), clouds, and rainfall is far from complete, especially in urban areas (Stevens 2008). Observation and numerical modeling are two important methods used in urban weather studies. Concurrent observations of processes in adjacent urban and rural atmospheres are fundamental to understanding interactions between underlying surfaces, aerosols, and the atmosphere. Numerical models for urban weather and air quality provide valuable information to decision makers and the public for many applications. The WMO has proposed that National Meteorological and Hydrological Services produce impact-based forecasts for environmental, energy, transportation, and health effects (Tan et al. 2015). Urban forecast model skill should increase as better a understanding of urban atmospheric processes emerges from large field measurement programs (Shepherd 2005). URBAN THUNDERSTORMS. Urban impacts on precipitation were the focus of the now-classic 1970s Metropolitan Meteorological Experiment (METROMEX) field study, which investigated modification of summer convective rainfall by the St. Louis urban area. Observations showed that the urban area increased precipitation, both within the city and km downwind of it, by 5 25% over upwind values (Huff and Vogel 1978; Changnon, 1979; Changnon et al. 1981; Braham et al. 1981; Changnon et al. 1991). Similar results were found by Shepherd et al. (2002) around Atlanta and four other large US cities (Burian and Shepherd 2005). Radar images of 91 summer thunderstorms during over Indianapolis by Niyogi et al. (2011), however, showed 60% of them with a variety of alterations (i.e., they initiated, intensified, split, or dissipated) over the city. These results agreed with the older observations by Bornstein and LeRoy (1990) for New York City (NYC) and Bornstein and Lin (2000) for Atlanta, as both studies also showed both urban initiated and bifurcated (i.e., splitting and going around the city, respectively) thunderstorms. Numerical modeling enables controlled experiments to characterize physical processes producing urban-precipitation modification. MM5 simulations by Craig and Bornstein (2002) showed an Atlanta thunderstorm arising from urban heat island (UHI)-induced convergence over the city. Simulations of two thunderstorms over Atlanta with WRF by Shem and Shepherd (2009) also showed the convergence effect, with rainfall increases downwind of the city. A case simulation with WRF and its Single Level Urban Canopy Module (SLUCM) option by Miao et al. (2011) showed that early (prior to 1980) Beijing urbanization decreased (as compared with a no-urban simulation) maximum rainfall over Beijing due to increased urban surface dryness. Further urbanization decreased rainfall, however, as storms bifurcated around the city. A case-study WRF simulation, coupled with its multi-level Building Environment Parameterization (BEP) urbanization option, by Yu and Liu (2015) showed that Beijing produced upstream and downstream precipitation increases. Yu et al. (2016) used SLUCM to simulate a heavy August Beijing storm, and results showed suburban increases of sensible heat flux, temperature, and wind speed, which facilitated storm merging and enhanced precipitation. Zhang et al. (2017b) simulated two Beijing summer thunderstorms with SLUCM, and showed contrasting urban effects on precipitation under strong and weak UHIs, i.e., initiation vs. bifurcation, respectively. Recent reviews have summarized the processes producing urban impacts on precipitation, e.g., Shepherd (2005) reviewed publications from , covering observations and simulations, and stressing both thermodynamic and microphysical impacts. Shepherd (2013) reviewed post-metromex lightning and winter precipitation studies, focusing on fluxes, aerosols, and bifurcation. Han et al. (2014) stressed a range of urban precipitation enhancers, from UHIs to surface roughness to aerosols, while Mitra and Shepherd (2016) focused on resulting storm modifications (e.g., bifurcation, merging). To address the above apparently conflicting urban precipitation impacts in the literature, a synthesis hypothesis was proposed by Bornstein (2011), in which strong UHIs initiate convective activity (and thus produce a city-center precipitation maximum) during periods with near calm regional flows. During strong wind, weak UHI conditions, however, moving thunderstorms bifurcate and move around cities because of building-barrier effects. This latter effect thus produces a precipitation minimum over and downwind of cities, and maximum around the city in lateral boundary convergence zones and in the downwind re-convergence area. Dou et al. (2015) for the first time showed UHI magnitude as key to determe urban effects on summer thunderstorms. Based on Beijing hourly precipitation observations, they showed that strong (greater than average) UHIs induced or enhanced thunderstorm formation. With weak UHI conditions, however, thunderstorms bifurcated and bypassed the city due to building-barrier effects. This study was thus the first validation of the Bornstein (2011) hypothesis concerning how UHImagnitude impacts summer convective precipitation, but an analysis of simulation results to understand bifurcation dynamics is still lacking

5 Numerical modeling of urban thunderstorms relies heavily on numerical weather prediction simulation results to provide background conditions. Such models, however, still have significant uncertainties in the estimation of the timing and amount of precipitation. Uncertainties in the numerical simulation of urbanization impacts on heavy rainfall, in particular, thus still require further investigation (Yu et al. 2018). URBAN AEROSOLS AND HAZE. Atmospheric aerosols, also called particulate matter (PM), directly influence radiation transfer processes and cause direct climate effects because of their optical scattering and absorption properties (Charlson et al.1992). As cloud condensation nuclei (CCN), they indirectly influence climate by changing cloud formation, precipitation, and radiation transfer processes (Rosenfeld et al. 1997; Liu and Daum 2002). Atmospheric aerosols have a variety of sources, as well as complex chemical compositions and physical properties, with anthropogenic fossil fuel combustion as their largest source (Andreae and Crutzen 1997). Biomass burning, wind-blown dust, sea salt, and volcanos are other important global sources. Chemical composition effects radiative transfer and cloud formation processes, e.g., black carbon (also called soot) and dust are the two most important absorption aerosols (Bond et al. 2013), while important scattering aerosols include sulfates, nitrates, and ammonium. Aerosol impacts on cloud microphysical processes, (e.g., evaporation, condensation, collision and coalescence, and droplet breakup) are better understood on larger scales than for city scales (van den Heever and Cotton 2007). Generalization of natural aerosol impacts is difficult, however, as they range from precipitation enhancement to suppression, e.g., high concentrations in winter orographic clouds yield more cloud water, but less rainwater (Rosenfeld et al. 2008). Too many CCN reduce autoconversion rates, since droplets coalesce less efficiently than giant CCN (GCCN), such as sea salts (Grabowski 1999). When droplets ascend above the freezing level, precipitation enhancement results, however, by dynamic invigoration from frozen droplets (Carrió et al. 2010), forming hail and graupel. An excellent review of all aspects of these processes appears in Fan et al. (2016). Biomass smoke also significantly alters synoptic weather (Robock 1991), and fossil fuel aerosols alter regional (or even global) climates (Jacobson 2001). Urban environments, with abundant pollutant CCN and GCCN, thus also modify cloud microphysical processes (after dynamical effects initiate it, as discussed above). Pollution aerosols have resulted in annual precipitation losses of 15 25% over topographic barriers downwind of major coastal urban areas in California and Israel (Givati and Rosenfeld 2004). Aerosol particle size distribution (PSD) also affects microphysical processes, e.g., RAMS idealized-simulations by Comarazamy et al. (2006) for Puerto Rico showed that cloud droplet production is significantly larger in polluted air than in clear skies, and that rainwater in polluted air is less than that in unpolluted air. Building from this study, Hosannah and González (2014) simulated impacts on a local NYC thunderstorm event by ingestion of observed AErosol RObotic NETwork (AERONET) PSDs ( Results showed impacted precipitation totals, rates, and spatial patterns, while month-long simulations showed smaller prediction errors with observed PSDs than with assumed model distributions. As one of the most densely populated countries, China has high aerosol loadings, because of intense anthropogenic and natural emissions. Its complex aerosol sources and geographical distribution, together with the large scale circulations of the East Asian monsoons, cause a complex transport and mixing of air masses, producing strong chemical and physical interactions (Ding et al. 2016). Long-term multi-satellite observations of rainfall, clouds, and aerosols indicate that they influence cloud properties, particularly cloud droplet size (Jin et al. 2008). Beijing is one of the largest cities in China, with a population over 20 million, and it suffers from high aerosol concentrations, and with its nearby complex terrain, is an excellent region to study aerosol-cloud-precipitation interactions. Many observational and numerical studies have investigated the role of aerosol on its radiation balance and precipitation processes. Surface radiation observations by Wang et al. (2016) at four Beijing sites showed polluted days having both attenuated direct and downwardscattered shortwave radiation fluxes. Modeling by Chen et al. (2015) found that while aerosols decreased Beijing urban-area rainfall, they increased domain-averaged values. Zhong et al. (2015) compared impacts from UHIs and aerosols on a summer precipitation event over Beijing with a WRF-Chem/UCM ensemble. Results showed aerosol impacts dominant, with enhanced downstream convection and precipitation. Direct observational evidence on the relationship between aerosols and cloud/ precipitation in Beijing is still limited, however, thus requiring time-resolving measurements of both aerosols and weather processes. The SURF project, which includes radiometer, lidar, radar, tower, and aircraft measurements, was thus initialized in 2015 to study aerosol effects on Beijing haze and precipitation processes. CRITICAL SCIENCE NEEDS. The Beijing-Tianjin-Hebei (BTH) area (Fig. 1a) includes the megacities of Beijing, Tianjin, Shijia-zhuang, Tangshan, with 90 million people, with 40 million just in the Beijing-Langfang-Tianjin belt. The BTH area incudes complex topography with mountains, plains, and coastal areas (Fig. 1a), and seven of the 10 most polluted Chinese cities, with 40% of days during 2013 (mostly in winter) having very hazardous air quality (CMEP 2014). Severe spring pollution frequently impacts Beijing and Northern China, with hourly PM 10 concentrations as high as 1000 µg m -3. Regional transport by summer monsoons subjects Beijing and surrounding cities to high concentrations of secondary aerosols from the North China Plain and southern China (Wang et al. 2010). Severe autumn and winter haze events also occur in Beijing e.g., PM 2.5 concentration were > 800 µg m -3 in January 2013 (Guo et al., 2014). Such episodes occur under stagnant, weak synoptic conditions, especially one or two days before a cold front passage. During haze episodes, the sky is colored light yellow or orange gray, and has a pungent smell

6 Given the human and environmental tolls from Beijing haze and floods (resulting from the precipitation processes discussed above), the Chinese Ministry of Science and Technology provided supports to the Institute of Urban Meteorology (IUM) for extensive field measurement and modeling campaigns to study both Beijing summer convective thunderstorm and winter haze events. A preliminary workshop (March 2014) defined study goals and methodologies. Additional planning workshops were held before the first summer thunderstorm (March 2015), first winter haze (November 2015), and second summer (May 2016) studies. All workshops were held at the IUM Beijing headquarters, with more than 40 invited Chinese experts from government agencies, laboratories, academia, and industry. An International Advisory Committee (IAC) was constituted with 10 representatives from the WMO, US, UK, Finland, and Israel. The entire group defined research needs and barriers, and workshop reports are available from the corresponding author. This paper thus summarizes the SURF effort, carried out to better understand atmospheric interactions between urban areas and aerosols that effect summer convective processes and winter haze episodes. From the first Workshop, the overall SURF goal was identified as: A better understanding of urban, terrain, convection, and aerosol interactions on both summer thunderstorm precipitation and winter haze episodes. The SURF project period is 2014 to 2019, and has the following specific goals: Conduct international collaborative field studies to improve understanding of urban impacts on convective weather systems and haze episodes. Improve high-resolution (1 km) urban numerical weather and air quality forecast models for research and operational applications. Increase use and effectiveness of urban forecasts for social, economic, and environmental applications, e.g., for health, energy, climate change mitigation and adaptation, air quality, planning, and emergency-response management. The critical science needed to achieve these goals was identified as increased understandings of: urban (i.e., UHI circulation, surface roughness, and building-barrier) and complex terrain impacts on 4-D PBL structures (i.e., height, turbulence, and inversions) by analysis of mountain-valley and urban produced wind and thermal patterns urban model parameterizations of PBL, building canopy, and water balance processes, as well as multi-scale data assimilation for observation-rich cities complex terrain, urban, and aerosol impacts and interactions on summer convective rainfall distribution, intensity, and timing, e.g., UHI-initiated vs. bifurcated thunderstorms long-range transport vs. local (urban and complex terrain) processes and pollutant emission impacts on horizontal, PBL, and diurnal structures of aerosol concentration patterns and on the duration of regional winter haze and fog episodes urban and complex-terrain 4-D meteorological and chemistry impacts on regional winter haze and fog episodes, including radiation and thermodynamic aerosolseffects on PBL structures how to improve the performance of high-resolution urban-model produced weather and air quality forecasts of summer precipitation and winter haze/fog episodes. SURF FIELD STUDIES. To gain a greater understanding and improved predictions of urban summer convective precipitation and winter haze, the greater BTH region is ideal for observational studies, as it is the largest metroplex in North China and experiences both summer flooding and winter haze. SURF summer field studies thus focused on urban impacts on rainfall initialization, movement, and intensity, while winter field studies focused on aerosol (regional and local) sources, transport, and transformation processes. Three field programs were conducted, two of which focused on thunderstorms (July - August 2015, July-September 2016) and one on haze (November 2015 January 2016). Field sites and instrumentation. Operational meteorological data are collected every 5 min at 1992 surface BTH-area Automatic Weather Stations (AWSs), 348 of which are in Beijing. Twice daily radiosondes are launched at 0800 and 2000 LST (where LST = GMT+8 h), at NAN (Fig. 1b; all site and agency acronyms are defined in Supplement A), with additional 1400 LST soundings during flood seasons (JuneAugust). Surface energy-balance fluxes have been observed operationally since 2012 by eddy-covariance techniques at one urban (IAP), three suburban (MIY, SHU, WUQ), and one rural (SDZ) tower (all PBL observational sites are shown in Fig. 1). Other operational PBL sensors (Table 1) include four S-band weather radars and six X-band radars, nine wind profilers, and a GPS/MET network, all supported by BMS, TMS, or HMS. During field study periods, existing operational PBL networks were augmented, e.g., the summer 2016 study (Table 1) added five wind profilers; one scanning Doppler lidar; a PBL aerosol Micro Pulse lidar (MPLs) each at the IAP tower and DAX sites; and 10 Vaisala CL51 ceilometers for cloud height and aerosol backscatter profiles. Continuous 10 Hz turbulent flux and mean meteorological data were also collected at multiple levels at two sites (the 325 m Beijing super site tower and the 38 m Tianjin tower). A van-mounted MPL aerosol lidar at DAX drove from the IAP fixed MPL site to IUM (Figs. 1b) on several days, while it circled the 4th Ring Road on one hazy day. During Intensive Operational Periods (IOPs), five additional radiosondes were launched daily at NAN, and eight small-balloon soundings were taken daily at each of three sites (IUM, CHA, DAX) from 28 August to 2 September 2016 under all conditions. Photos of a selection of these instruments and sites are shown in Fig

7 About 250 BTH-area Environmental Protection Bureau (EPB) stations (35 in Beijing), in addition to 14 BMS stations, operationally measure hourly concentrations of surface: PM 2.5, PM 10, SO 2, NO 2, NO, NO x, CO, and O 3. Special mobile lidar observations were also carried out during three 2016 winter IOPs. During fair weather winter-2015 and IOPs (15 and 19, respectively) King Air aircraft flights were flown over pre-approved flight paths (Fig. 1) at altitudes from 600 to m (at 300 m intervals). The aircraft is equipped with atmospheric gas and aerosol instrumentation, including cloud images (resolution of 2.3 μm), cloud drop spectra (2 to μm), aerosol spectra (0.11 to 3 μm), and meteorology probes (1 Hz sampling frequency) for temperature (T), relative humidity (RH), and horizontal wind speed (V) and direction. Data sets. The Urban Meteorological Data Sharing Platform (UMDSP) archives the 2015 and 2016 SURF datasets, includes site and instrument metadata, manufacturerprovided software, and quality controlled (QCed) observations from the wind profilers, MPLs, Doppler lidar, flux towers, microwave radiometers, and radiosondes.additional observation data (e.g., weather radar, aircraft, and lightning) will be added. Golden case SURF datasets (defined as well documented, clear signal events) are available to research partners. The status of SURF PBL observational data are summarized in Supplement B, and all SURF instrumentation sources, accuracies, sampling periods, data, etc. can be viewed and downloaded by the international research community through granted access to the UMDSP website: IOP forecasts. A daily procedure was followed by the SURF Forecast Group (SFG) to decide whether or not to call for summer 12 h IOP, with its additional rawinsonde soundings and/or mobile observations. Each day at 0830 LST during a field study period, SFG forecasters received the BMS 72 h operational synoptic forecast. They focused on evolving synoptic patterns to determine if thunderstorms were likely on the next day; if so, a tentative go (for the IOP) was issued. At 1500 LST, BMS updated its 48 h forecast, and if thunderstorms still seemed likely on the next morning, SFG forecasters evaluated the latest synoptic maps and NWP forecasts. The models considered included the European Centre for Medium Range Weather Forecast (ECMWF) system and IUM Urban (described below). Based on this evaluation, a final "go" or "no-go" decision was made. The process resulted in a five-part daily forecast summary: previous-day weather situation; current (500 and 850 hpa) upper level and surface charts; 36 h European Center forecast; 36 h IUM 3 km forecast; 24 h IUM 1-km urban forecast; and SFG IOP conclusions. The forecast accuracy for summer-2015 and 2016 thunderstorm events (Table 2) shows that correct forecasts increased (due to the year of additional experience) from 67% for IOP days in 2015 to 74% for days in Correct forecasts predicted thunderstorm activity (or not) anywhere in the Beijing study area anytime on the given IOP day. The forecast procedure for haze IOPs is simpler, as winter haze mostly occurs under stable weather condition and as severe episodes can last more than three days. The IUM operational WRF-Chem model provides 96 h forecasts, and if a haze event is forecast to occur within four days and to last for at least three days, a haze IOP (with BJWMO aircraft observations) was carried out until the event ends. SURF MODELING. IUM operational weather and air quality forecast system. A modified version of the WRF model (ARW versions 3.5.1; Skamarock et al. 2008) is used as the IUM weather model in both operational and research modes. Its operational Rapid-refresh Multi-scale Analysis and Prediction System (RMAPS) includes five off-line, one-way nested components, each with its own horizontal grid spacing ( x). The system starts with ECMWF Global forecasts (at 3 h intervals) and terminates with hourly weather and air quality forecasts. Relationships between these internal operational components, as well as the observational data used by each, are shown in the flow chart of Fig. 3 and described below: Nowcasting (NOW) Model ( x of 5 and 2.5 km): uses input from the Variational Doppler Radar Analysis System (VDRAS) for 6 h real time forecasts, with ingested data from seven S-band weather radar and all AWS sites. Short Term (ST) Model ( x = 9 and 3 km): uses WRF Data Assimilation (WRFDA), and either the WRF Noah (for rural grid points) or SLUCM (for urban grid points) land-surface modules for operational 72 h forecasts, for which regional and local data are assimilated. Urban Model ( x =1 km): for urban weather forecasts with the multi-level BEP + Building Energy Model (BEM) urban PBL modules for 24 h real time (but not yet operational) forecasts and research studies, where BEM, BEP, and their urban LU/LC input data are described below and where its AWS assimilated input weather-data were described above. Integration (IN) Model: combines observations into the NOW, ST, and Urban model forecasts for objective 12 h forecasts, where its input (QPE) are radar data and are calibrated by AWS rain-gauge observations. CHEM Model ( x = 9 km and 3 km): for 96 h chemical forecasts from the WRFChem model, where its input emission (Q*) and concentration (χ) data are discussed below. In the operational mode, the 1-km Domain-3 of the Urban model (459 by 402 grid cells) is driven by its 3-km Domain-2 (549 by 423 grid cells) output. The 3-km domain encompasses most of North China, while the 1-km domain covers most of the BTH area centered on Beijing (a somewhat larger area than in Fig. 1a). Initial conditions (ICs) for the 3-km operational ST model are obtained through the WRF preprocessor NDOWN option and a 3-D VAR approach, which uses 9-km ST model outputs, with its IC and boundary conditions (BCs) from ECMWF, the National Centers for Environmental Prediction Global Forecast System (NCEP GFS), or CMA Global and Regional Assimilation and PrEdiction System (GRAPES), and with the preferred order of usage of: ECMWF (if available), GFS, and GRAPES (if necessary). The 3-km 72-h operational simulations have their ICs adjusted by WRFDA-ingested of: S/C band weather radar, Global Navigation Satellite System-Zenith Total Delay (GNSS-ZTD), radiosonde, Aircraft Meteorological Data Relay - 5 -

8 (AMDAR), and AWS surface observations. NDOWN then provides BCs for the 1-km Domain-3 Urban model from its 3-km output, and VDRAS output is assimilated into the 1 km domain via FDDA. Additional details on all of these steps are provided by Zhang et al. (2017a). RMAPS analysis and forecast products are operationally provided to BMS for weather warning and many forecasting applications, e.g., urban flooding, air pollution, weather modification, energy consumption, and emergency-response management. Application of this system for modeling SURF UHI, thunderstorm, and haze casestudy events are discussed below. Parameterization options for each application are also given. Urban weather-model. The 1-km WRF-Urban cooling-tower and 3 km thunderstorm simulations used the sophisticated multilayer urban canopy model option in WRF-Urban (Chen et al. 2011), comprising advanced versions of BEP of Martilli et al. (2002) and (in the cooling tower study) of BEM of Salamanca et al. (2010) to estimate anthropogenic heat fluxes (Q A). The single NASA 30-m Landsat urban land-use class is scaled up to 1 or 3 -km (as appropriate) and then resampled into three new urban building (i.e., low, medium, and high) classes according to gridded urban-fraction values. Required gridded urban parameters (e.g., building-plan area fraction, mean building height) are derived from remotely-sensed building morphological data maintained by the Beijing Institute of Surveying and Mapping and scaled up according to building-planarea-weighted average; details provided in Zhang et al. (2013). Other WRF parameterization options selected for the two urban weather studies are listed in Table 3, where the Noah landsurface module has four soil moisture and temperature layers, and where the BouLac PBL scheme uses TKE prediction. The 9km domain also uses Kain-Fritsch cumulus parameterization, not needed in the finer resolution domains (Fan et al. 2013). Chemical model. The operational IUM 3 km WRF-Chem model simulates formation, transformation, and transport processes of both primary and secondary atmospheric pollutants, including gases and PM species. It used two nested grids (Table 3), with 38 vertical layers. Physical parameterizations include: SLUCM; Carbon Bond Mechanism-Z (CBM-Z) chemical mechanism, which include comprehensive reactions and alterable scenarios (Zaveri et al. 1999); MOdel for Simulating Aerosol Interactions; MOSAIC chemistry-aerosol mechanism, with its four bin sizes (Zaveri et al. 2008); and the fast photolysis computation code of Fast-J, Wild et al. (2000). Anthropogenic emission- data for 2012, with a resolution of , was developed by Multi-resolution Emission Inventory for China (MEIC) group ( and were interpolated to the 3 and 9 km model grids via a bilinear method. Meteorological ICs for the 9 km and 3km WRF-Chem model, and BCs for the 9 km grid, were obtained through the WRF Preprocessing System (WPS), via output from the 9 km operational ST model. An idealized vertical profile for each chemical species, based upon northern hemispheric, mid-latitude, clean environment conditions, is used to provide chemical ICs and BCs (WRF-Chem User Guide 2017). Meteorological and chemical BCs for the 3 km model are directly obtained from the outer 9 km domain. PRELIMINARY FINDINGS. The following preliminary results are presented only to highlight the range of SURF research activities, with additional details available in cited references and on the UMDSP website: Urban HOST. Nighttime relationships between turbulence intensity and wind speed over urban canopies were investigated to verify the HOckey-Stick Transition (HOST) hypothesis of Sun et al. (2012, 2016), previously developed from observations over a grass field. The hypothesis describes the roles of non-local large coherent turbulence eddies during turbulent mixing via the linear relationship between turbulence intensity and mean wind speed V under near neutral conditions, i.e.,, where τ is turbulent momentum stress, ρ is air density, and u is friction velocity. Turbulent flux measurements at the 47 m level (near the urban canopy top) of the 325 IAP tower (Fig. 1) were used to show that the dominant linear relationship between u and V also exists over urban canopies. The strong wind shear from the rough urban surfaces produces turbulence in near neutral urban stability conditions, shown by Bornstein 1968 to exist over NYC, as stable boundary layers over urban canopies are thus hard to maintain. The role of surface roughness on turbulent mixing is also reflected in the increasing slope (Fig. 4) of the u -V relationship with increasing surface roughness (sorted by wind direction) over the heterogeneous Beijing surface, consistent with the observations by Mahrt et al. (2013). Urban surface energy balance. The exchange of energy between urban surfaces and the overlying atmosphere involves five important processes: absorption and emission of radiation by urban surfaces, turbulent transfer of heat energy to or away from the surface, evaporation of water stored in its surfaces and/or condensation of atmosphere water vapor onto the surface, thermal conduction of heat within the ground, and anthropogenic heat fluxes. Table 4a shows the average midday ( LST) radiative and energy fluxes from both the 140-m level of the urban IAP tower and the 36-m suburban MIY tower (sites are about 50 km apart, Fig. 1) for June to August The IAP site shows slightly lower incoming shortwave radiation, due to its higher atmospheric PM 2.5 concentrations (discussed below). Reflected shortwave radiation is lower at IAP, with an albedo of only 0.10, as compared to the 0.13 suburban-value. Both sites have similar incoming longwave radiation fluxes, as its dominant control is cloudiness and its values are similar at both sites, i.e., 14 clear, 20 cloudy, and 31 rain days; the cloudiness classification procedure is described in Miao et al. (2012). Although MIY has higher outgoing longwave radiation, the net all-wave radiations are nearly equal. Significant differences existed, however, in the surface energy partitions, as IAP has smaller turbulent sensible and latent heat fluxes, and thus a larger - 6 -

9 (estimated) residual heat storage. Another contributing factor to urban heat storage is its anthropogenic heat flux source (discussed below). Daily mean Bowen ratios (ratio of turbulent sensible to latent-heat flux; Table 4b) of 0.56 and 0.49 for IAP and MIY, respectively, are lower than previous urban and suburban estimates from sites with similar or larger fractional vegetated cover. For example, Tokyo has a ratio of 1.77 in July with a urban vegetation percentage of 29% (Moriwaki and Kanda 2004), while the values for IAP and MIY are 21.7 and 37.0, respectively (Miao et al. 2012). The high Beijing urban latent heat flux Q E is likely due to roadway greenbelt-watering, air conditioning, and roadway sprinkling for dust control, while at suburban MIY, it is from crop irrigation (Miao and Chen 2014). Urban PBL depth. PBL depth is important for many applications, including weather and air-quality forecasts. Data from the Leosphere WindCube Doppler lidar (located near the IAP tower) were used to study PBL evolution on six days during the 2015 summer SURF campaign. A Doppler Beam Swinging (DBS) scan-mode was used, which takes 4-5 or 2-3 s per beam scan with elevation angles of 75 or 90, respectively. Daytime convective boundary layer (CBL) depth was estimated by a vertical velocity variance threshold (1 m 2 s -2, Barlow et al. 2011), as lidar wind profiles have shown this method as more accurate than maximum wind-shear methods (Huang et al. 2016). This method, however, has limitations with nocturnal boundary layers (NBLs), so their tops are herein defined as where the variance decreases to 10% of its nearsurface maximum (minus a background variance), following LeMone et al. (2014). These two methods have thus herein been combined to determine diurnal variations of SURF urban PBL depths. Results (Fig. 5) show clear diurnal patterns, with depths ranging from m. They also show good agreement with the elevated inversion in the 1315 LST radiosonde potential-temperature profile at the NAN site. Normalized relative backscatter (NRB) data from two MPL lidar (in a vertical to zenith scan mode) were also used to concurrently estimate midday PBL heights on 11 August 2015 from LST, both alone a mobile route (Fig. 1b) and at a fixed site (adjacent to the urban IAP tower). As MPL instruments cannot be absolutely calibrated, small differences exist in CRB values when both were co-located (at start time of Fig. 6). Results show a fairly uniform urban CBL depth (about 1.2 km), with well mixed conditions near the tower (Fig. 6a). While the mobile rural-heights (about 0.8 km, Fig. 6b) were generally lower than the concurrent IAP urban values (s expected), they increased when the instrument passed through the urban areas. In addition, both sections show some PBL-top gravity wave activity and the mobile results also show two large 10-min rural convective periods. Urban weather-model development. RMPAS-Urban forecasts suffer from overpredicted (especially nocturnal) 2-m temperatures and 10-m wind speeds. Use of the BEP scheme (also with default look-up table input parameters), did improve Beijing operational 2-m temperature and 10-m wind speed absolute-biases by 4 C and 1 m s -1, respectively (Barlage et al. 2016). Despite the use of BEP, urban models still consistently overestimate 2-m temperatures because of incorrect predicted surface energy-balance terms. Beijing anthropogenic heat emissions are comparable to other major cities, e.g. NYC (Gutiérrez et al. 2015a), and mega-city anthropogenic heat is partly due to air conditioning (A/C) systems. These systems release a mixture of sensible and latent heat to the environment, and in summer, 50 80% of the heat released by evaporativecooled A/C systems is latent energy (Sailor et al. 2007). Air conditioning heat releases in BEM, however, are assumed as 100% sensible heat (Salamanca et al. 2011), an important reason for the over-predicted 2-m urban temperatures. To improve RMAPS-Urban performance, SURF simulations used the new Evaporative Cooling (EC) scheme of Gutiérrez et al. (2015a, b), developed for NYC and incorporated into the BEM urbanization option. Model evaluation of predicted temperatures and humidities used both data from eight urban AWS sites across Beijing and fluxes from the 140-m level of the IAP tower (Fig. 1b). The modeling period was 4-11 July 2015, which included dry (6 th -11 th ) and wet (4 th -5 th ) days. Results show major improvement for daytime total sensible plus latent heat flux values (Fig. 7a and b), although the timing of its peak was about 2 h too late. The increased latent heat flux during the EC simulation (accurately) slightly increased moisture-values throughout the diurnal cycle (Fig. 7d). While the corresponding decreased nighttime latent heat flux also improved the too high temperature peak (Fig. 7c), it remained 2-h too early. The model is also much drier than the observations in the initial period, presumably due to the specified initial land surface soil moisture. These impacts from the EC parameterization on simulated urban temperature and humidity are, however, significant (Table 5). The largest RMSE improvement for temperature is for the wet days, with the reverse true for the dry days; smaller differences are seen between dry and wet days for the corresponding humidity values. The largest bias improvements with the EC formulation for both temperature and humidity, however, are on the dry days. Improvements were also seen in the PBL height and inversion-structure (not shown), also consistent with the NYC results of Gutiérrez et al. (2015b). Future efforts to obtain improved urban thermal-response data will be carried out via in-situ or satellite infrared observations. Thunderstorm-bifurcation modeling. Two simulations with the 3 km IUM Urban model (Table 3) were carried out for the 22 July 2015 SURF IOP thunderstorm, an: (a) Urban case with 2010 land use data and (b) No-Urban case, which converted all urban Beijing-area grid points to cropland, the dominant surface surrounding the city. Results reproduced the weak observed UHI before the storm bifurcated around Beijing, but its magnitude was underestimated by 0.5 o C (not shown). While the Urban case somewhat overestimated peak precipitation and total accumulation around the city, it correctly simulated its timing, location, and bifurcation pattern (Figs. 8a vs. 8c); the No- Urban results (Fig. 8b) showed no bifurcation. Results were used to gain insights (not shown) into urban dynamic - 7 -

10 mechanisms leading to bifurcation, i.e., a deep (blocking) high-pressure perturbation upwind of Beijing and deep downward urban vertical velocities; additional details can be found in Dou et al. (2018) and Bornstein et al. (2018). Haze event. A regional PM 2.5 haze episode occurred in the North China Plain (NCP, Fig. 1a) from 26 November to 1 December 2015, when observed hourly concentrations at many Beijing sites abruptly increased and decreased several times in 24-h. Average concentrations at 35 Beijing sites (red line, Fig. 9) experienced PM 2.5decreases at 1700 LST on 29 November and 0600 LT on the next day, followed 13 h later by an increase to 500 μg m -3. This value is 15 times the annual standard of 35 μg m -3 (no hourly standard exists in China). To understand the dynamic, physical, and chemical processes responsible for these high concentrations and rapid changes, a WRF-Chem simulation was carried out (Table 3). The observations show that northward-moving near-surface polluted air met southward-moving clear air and thus formed a surface convergence zone. This zone first passed northward over Beijing by 2000 LST on 29 November (Fig. 10a), producing high PM 2.5 values (up to 500 μg m -3 ) over southern urban Beijing. When the cold air strengthened and pushed the boundary southward 12 h later, concentrations over the city fell abruptly to 100 μg m -3 (Fig. 10b). A re-strengthened southerly flow then pushed the boundary northward over Beijing for a second time at 1400 LST on the 30 th (Fig. 10c), and concentration in the southern part of the city thus rose to 700 μg m -3. These high concentrations continued to move northward, so that 6 h later, an UHI-induced convergence (Fig. 10d) drew maximum concentrations into central Beijing from the south. Data from the IUM Doppler lidar at the 325-m IAP tower and from the mini-mpl lidar at HAI (7 km NW of the tower) were also used in the analysis of PM 2.5 concentrations on the 30th. The HAI data (Fig. 11a) show increasing surface PM 2.5 starting at 0900 LST, peaking at 600 μg m -3 by 1700 LST. The near-surface flow from the northeast at IAP (up to 370 m) until about 0500 LST corresponds to decreasing PM 2.5 at HAI and to the low Carrier-to-Noise Ratio (CNR) values at IAP (Fig. 11b). As the PM 2.5 front pushed into Beijing, concentrations (as reflected in the larger CNR values) steadily rose throughout the lowest 600 m of the IAP PBL. The shape of the IAP (fixed point) internal boundary layer is consistent with urban-barrier-induced surface frontalretardation of a moving cold-front. As the IAP tower was always north of the concentration front (Fig. 10), winds in its lowest 200 m remained northerly, although above 300 m they were consistently southerly. Chemical modeling (Table 3) correctly simulated the maximum PM 2.5 south of Beijing at 2000 LST on 29 November (Figs. 12a vs. 10a), but incorrectly had the max southeast of the city (and not to its east) at 0800 LST on the 30th (Figs. 12c vs. 10c). It did, however, correctly place a narrow maximum over Beijing 12 h later (Figs. 12d vs. 10d). All simulated maximum concentrations are, however, underestimated by a factor of two. Possible reasons for this include underestimated heating-season pollutant emissions, failure to correctly simulate the complex urban meteorology conditions, and missing modeled heterologous reactions during high humidity conditions. Li et al. (2017) also showed that shifting wind directions (from northerly to southwesterly) above the stable boundary layer was the key process that transported the heavily polluted air over Beijing. Morning convective turbulent mixing then transported upper polluted air downward, leading to the observed dramatic increases of PM 2.5 on November 30. Such local scale processes (e.g., dilution due to northwesterly downslope winds from the mountains around Beijing, and the rapid increases of pollutants due to shifting wind directions and vertical mixing) are difficult to model and can lead to serious underestimation of PM 2.5 peak concentrations. SUMMARY. SURF is truly a multidiscipline effort, in which scientists (i.e. observationalists, modelers, meteorologists, engineers, and chemists) and forecasters collaborate across disciplines to create knowledge and to develop modeling tools to improve urban weather and air quality predictions. SURF has four components: experimental (field) studies, numerical model development, modeled case studies of thunderstorms and haze events, and applications (improved forecasts). The structure of its lead agency (IUM), which comprises both urban PBL researchers and operational weather and air quality urban numericalforecasters, facilitates active collaboration with stakeholders, who provided useful insights for hypothesis testing and field study planning. New scientific insights from SURF include use of flux data from a 325 m urban tower to confirm the HOckey-Stick Transition (HOST) turbulence-wind relationships over urban canopies, and a new composite diagnostic-method of Doppler lidar turbulence intensities to produce more accurate estimates of diurnal variations of urban PBL-depth. These new findings will be applied to reduce systemic differences between simultaneous fixed and mobile lidar observations. They will also be incorporated into IUM operational WRF models to evaluate their roles in forecast improvement, as the IUM urbanization schemes require additional development to reduce predicted urban temperature and wind biases. A new anthropogenic cooling tower latent-energy module was added to the BEM urban weather model; while it reduced model biases, additional development is also needed. Analyses of energy balance fluxes from the 325 m Beijing tower and at a nearby suburban tower showed reasonable urban-induced differences, which will also be useful in this effort. A SURF golden case bifurcating summer convective thunderstorm event was analyzed and correctly simulated with the IUM urbanized WRF model. Simulated outputs were useful in gaining insight into the bifurcation mechanisms (i.e. from deep blocking high-pressure perturbations and intense downward vertical velocities); other SURF storms will be analyzed for additional insights. Regional and urban aspects of a SURF haze episode were analyzed and correctly simulated with the IUM chemical model to better understand 4-D transport mechanisms. Urban land surface, anthropogenic heat, and emission impacts on such events will be further investigated

11 Many challenges still need investigation within the scope of SURF, including development of advanced data assimilation methods for the dense and diverse observations at the urban scale, (e.g., aerosol PSDs). While the BEP and BEM urban parameterizations represent current state-of-the-art for high resolution urban forecast models, problems still exist in the partition between their latent and sensible anthropogenic-heat components and in their vertical distributions. The comprehensive SURF data set thus presents unique opportunities to test, calibrate, and develop improvements for such models, as they allow for model performance assessment, parameterization development, and golden-case analyses. An example of this is the evaporative cooling (cooling tower) SURF study, as its large errors in sensible heat flux need further investigation. Comparative studies of modeled and observational energy and radiation budgets should be carried out via model sensitivity experiments (e.g., for PBL height values.) SURF looks to explore new boundaries for high-resolution urban modeling (to subkilometer scales) for a variety of summer and winter conditions. Urban precipitation processes involve large-scale and urban forcings, as well as cloud microphysics modifications from natural and urban CCN, but aerosol impacts on urban precipitation are still largely unknown. SURF may answer some of these questions by its planned urbanization of WRF-Chem and by application of the new model to its golden case data sets. SURF has now entered an intensive data processing and modeling phase. These theoretical, numerical, and operational developments, as well as the massive datasets and planned journal publications, will aid future research in a wide variety of cities into the remaining unknowns of urban meteorology, climate, climate change, and air quality. ACKNOWLEDGEMENTS: SURF was funded by the Ministry of Science and Technology of China (Grant No. 2015DFA20870), with additional funding from the Beijing Municipal Science and Technology Commission (Grant Nos. Z and D ) and the National Natural Science Foundation of China (Grant Nos , , , , and ). Without the invaluable support of the observational staffs at BMS, TMS, HMS, IAP, and the Beijing Airda Electronic Equipment Co., Ltd., the project could not have evolved. REFERENCES Andreae, M. O., C. D. Jones, and C. P. Cox, 2005: Strong present-day aerosol cooling implies a hot future. Nature, 435, , D. Rosenfeld, P. Artaxo, A. A. Costa, G. P. Frank, K. M. Longo, and M. A. F. Silva-Dias, 2004: Smoking rain clouds over the Amazon. Science, 303, , M. O., and P. J. Crutzen, 1997: Atmospheric Aerosols: Biogeochemical Sources and Role in Atmospheric Chemistry. Science, 276, Ao, X., C. Grimmond, D. Liu, Z. Han, P. Hu, Y. Wang, X. Zhen, and J. Tan, 2016: Radiation fluxes in a business district of Shanghai, China. J. Appl. Meteor. Climatol., 55, Barlage, M., S. Miao, and F. Chen, 2016: Impact of physics parameterizations on high-resolution weather prediction over two Chinese megacities. J. Geophys. Res., 121, , doi: /2015jd Barlow, J. F., T. M. Dunbar, E. G. Nemitz, C. R. Wood, M. W. Gallagher, F. Davies, E. O'Connor, and R. M. Harrison, 2011: Boundary layer dynamics over London, UK, as observed using Doppler lidar during REPARTEE-II. Atmos. Chem. Phys., 11, Beljaars, A. C. M., 1994: The parameterization of surface fluxes in large-scale models under free convection. Quart. J. Royal Meteor. Soc., 121, Bond, T. C., and co-authors, 2013: Bounding the role of black carbon in the climate system: A scientific assessment. J. Geo. Res., 118, , doi: / jgrd Bornstein, R., and G. LeRoy, 1990: Urban barrier effects on convective and frontal thunderstorms. Preprints, 4th AMS Conf. on Mesoscale Processes, Boulder, CO, , doi: / _14., and Q. Lin, 2000: Urban heat islands and summertime convective thunderstorms in Atlanta: three case studies. Atmos. Environ., 34, , 2011: Establishment of meso-met modeling case studies to evaluate the relative roles of urban dynamics and aerosols on summer thunderstorms: A proposal. Presented 18th Conf. on Planned and Inadvertent Weather Modification and Third Symp, On Aerosol-Cloud-Climate Interactions, Seattle, WA, Amer. Meteor. Soc., J5.2. (Available online at recordingredirect.cgi/id/17137)., J. Dou, and S. Miao, 2018: A Beijing urban-thunderstorm bifurcation case study. Part 2: explanation of dynamic mechanisms. Mon. Wea. Rev. (submitted) Bougeault, P., and P. Lacarrere, 1989: Parameterization of orography-induced turbulence in a mesobeta--scale model. Mon. Wea. Rev., 117, , doi: / (1989)117<1872:pooiti> 2.0.CO;2. Braham, R. R., R. G. Semonin, A. H. Auer, S. A. Changnon Jr., and J. M. Hales, 1981: Summary of urban effects on clouds and rain, METROMEX: A review and summary. Meteor. Monogr., No. 40, Amer. Meteor. Soc., , doi: / _7. Burian, S. J., and J. M. Shepherd, 2005: Effect of urbanization on the diurnal rainfall pattern in Houston. Hydrol. Processes, 19, , doi: /hyp Cao, J. J., X. Tie, W. Dabberdt, Z. Z. Zhao, and T. Jie, 2013: On potential acid rain enhancement in eastern China. J. Geophys. Res., 118, Carrió, G. G., W. R. Cotton, and W. Y. Y. Cheng, 2010: Urban growth and aerosol effects on convection over Houston. Part I: The August 2000 case. Atmos. Res., 96, Changnon, S. A., 1979: Rainfall changes in summer caused by St. Louis. Science, 205, 402 4, doi: /science , R. G. Semonin, A. H. Auer, R. R. Braham, and J. Hales, 1981: METROMEX: A review & summary. Meteor.Monogr., 40, Amer. Meteor. Soc., 81 pp

12 , R. T. Shealy, and R. W. Scott, 1991: Precipitation changes in fall, winter, and spring caused by St. Louis. J. Appl. Meteor., 30, , doi: / (1991)030<0126:pcifwa>2.0.co;2. Chen, F., and J. Dudhia, 2001: Coupling an advanced land surface hydrology model with the Penn State NCAR MM5 modeling system. Part I: Model implementation. Mon. Wea. Rev., 129, Chen, F., H. Kusaka, R. Bornstein, J. Ching, C. S. B. Grimmond, S. Grossman-Clarke, T. Loridan, K. W. Manning, A. Martilli, S. Miao, D. Sailor, F. P. Salamanca, H. Taha, M. Tewari, X. Wang, A. A. Wyszogrodzki, and C. Zhang, 2011: The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems. Int. J. Climatol., 31, Chen W., D. Fu, S. Miao, and Y. Zhang, 2015: Impacts of aerosols from Beijing and the surrounding areas on urban precipitation. Chinese Sci. Bull., 60, , doi: /N (in Chinese with English abstract). CMEP, 2014: China Ministry of Environmental Protection issues air quality statement for important regions and 74 cities in China for gkml/hbb/qt/201403/t _ htm. Charlson, R. J., S. E. Schwartz, J. M. Hales, R. D. Cess, J. A. Coakley JR., J. E. Hansen, and D. J. Hofmann, 1992: Climate forcing by anthropogenic aerosols. Science, 255, Chou M.-D., and M. J. Suarez, 1994: An efficient thermal infrared radiation parameterization for use in general circulation models. NASA Tech. Memo., , 3, 85pp. Comarazamy, D. E., J. E. Gonzalez, C. A. Tepley, S. Raizada, and R. V. R. Pandya, 2006: Effects of atmospheric particle concentration on cloud microphysics over Arecibo. J. Geophys. Res, 111, D09205, doi: /2005JD Craig, K., and R. Bornstein, 2002: MM5 simulation of urban induced convective precipitation over Atlanta. Preprints, Fourth AMS Conf. on the Urban Environment, Norfolk, VA, Amer. Meteor. Soc., 5 6. Demuzere M., K. De Ridder, and N. P. M Van Lipzig., Modeling the energy balance in Marseille: Sensitivity to roughness length parameterizations and thermal admittance. J. Geophys. Res., 113, D doi: /2007JD Di Liberto T., 2015: India heat wave kills thousands. NOAA Climate.gov, June 9, Ding, A., and co-authors, 2016: Long-term observation of air pollution-weather/ climate interactions at the SORPES station: a review and outlook. Frontiers of Environmental Science & Engineering, 10, 15, doi: /s Dou, J., Y. Wang, R. Bornstein, and S. Miao, 2015: Observed spatial characteristics of Beijing urban climate impacts on summer thunderstorms. J. Appl. Meteor. Climatol., 54, , S. Miao, R. Bornstein, 2018: A Beijing urban thunderstorm-bifurcation case study. Part 1: observations and simulations. Mon. Wea. Rev. (submitted) Dudhia, J., 1989: Numerical study of convection observed during the winter monsoon experiment using a mesoscale twodimensional model. J. Atmos. Sci., 46, Fan, J., Y. Wang, D. Rosenfeld, and X. Liu, 2016: Review of aerosol cloud interactions: Mechanisms, significance, and challenges. J. Atmos. Sci., 73, doi: /JAS-D Fan, S., H. Wang, M. Chen., and H. Gao, 2013: Study of the data assimilation of radar reflectivity with the WRF 3DVar. Acta Meteorol. Sinica, 71, , doi: /qxxb (in Chinese with English abstract). Givati, A., and D. Rosenfeld, 2004: Quantifying precipitation suppression due to air pollution. J. Appl. Meteor. Climatol., 43, Grabowski, W. W., 1999: A parameterization of cloud microphysics for long-term cloud-resolving modeling of tropical convection. Atmos. Res., 52, Grell, G. A., and D. Devenyi, 2002: A generalized approach to parameterizing convection combining ensemble and data assimilation techniques. Geophys. Res. Lett., 29, Article Guo, S., and co-authors, 2014: Elucidating severe urban haze formation in China. Proceedings US National Academy of Sciences, 111, 17, Gutierrez, E., J.E. González, M. Arend, R. Bornstein, and A. Martilli, 2015a: Simulations of a heat wave event in New York City using a multi-layer urban parameterization. J. Appl. Meteor. Climatol., 54, doi: Gutierrez E., J.E. González, M. Arend, R. Bornstein, and A. Martilli, 2015b: On the anthropogenic heat fluxes using an air conditioning evaporative cooling parameterization for mesoscale urban canopy models. J. of Solar Energy Engineering, 137, doi: / Han, J. Y., J. J. Baik, and H. Lee, 2014: Urban Impacts on Precipitation, Asia-Pac. J. Atmos. Sci., 50, 17-30, Hosannah, N., and J. E. González, 2014: Impacts of aerosol particle size distribution and land cover land use on precipitation in a coastal urban environment using a cloud-resolving Mesoscale model. Adv. in Meteor., 2014, Article ID , doi: /2014/ Hong, S.-Y., Y. Noh, and J. Dudhia, 2006: A new vertical diffusion package with an explicit treatment of entrainment processes. Mon. Wea. Rev., 134, Huang, M., Z. Gao, S. Miao, F. Chen, M. A. LeMone, J. Li, F. Hu, and L. Wang, 2016: Estimate of boundary-layer depths over Beijing, China, using Doppler lidar data during SURF Bound. Layer Meteor., 162, doi: /s Huff, F. A., and S. A. Changnon, 1972: Climatological assessment of urban effects on precipitation at St. Louis. J. Appl. Meteor., 11, , and, 1973: Precipitation modification by major urban areas. Bull. Amer. Meteor. Soc., 54, , and J. L. Vogel, 1978: Urban, topographic and diurnal effects on rainfall in the St. Louis region. J. Appl. Meteor., 17, , doi: / (1978) 017<0565:UTADEO>2.0.CO;2. Jacobson, M. Z., 2001: Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols. Nature, 409,

13 Janjic, Z. I., 2002: Nonsingular Implementation of the Mellor Yamada Level 2.5 Scheme in the NCEP Meso model, NCEP Office Note, No. 437, 61 pp. Jin, M. L., and J. M. Shepherd, 2008: Aerosol relationships to warm season clouds and rainfall at monthly scales over east China: Urban land versus ocean. J. Geophys. Res. Atmos., 113, doi: /2008 JD Kain, J. S., 2004: The Kain Fritsch convective parameterization: an update. J. Appl. Meteor., 43, , doi: / (2004)043<0170:TKCPAU> 2.0.CO;2. Kusaka, H., H. Kondo, Y. Kikegawa, and F. Kimura, 2001: A simple single-layer urban canopy model for atmospheric models: Comparison with multi-layer and slab models. Bound.-Layer Meteor., 101, LeMone, M. A., M. Tewari, F. Chen, and J. Dudhia, 2014: Objectively determined fair-weather NBL features in ARW-WRF and their comparison to CASES-97 observations. Mon. Wea. Rev., 142, Li, J., J. Sun, M. Zhou, Z. Cheng, Q. Li, X. Cao, and J. Zhang, 2017: Observational analyses of dramatic developments of a severe air pollution event in the Beijing area, Atmos. Chem. Phys. Discuss., in review. Li, Z., F. Niu, J. Fan, Y. Liu, D. Rosenfeld, and Y. Ding, 2011: Long-term impacts of aerosols on the vertical development of clouds and precipitation. Nature Geoscience, 4, Liu, Y. G., and P. H. Daum, 2002: Anthropogenic aerosols - Indirect warming effect from dispersion forcing. Nature, 419, , doi: /419580a. Mahrt, L., C. Thomas, S. Richardson, N. Seaman, D. Stauffer, and M. Zeeman, 2013: Non-stationary generation of weak turbulence for very stable and weak-wind conditions. Bound. Layer Meteor., 147, Martilli, A., A. Clappier, and M. W. Rotach, 2002: An urban surface exchange parameterization for mesoscale models. Bound. Layer Meteor., 104, Miao, S., and F. Chen, 2014: Enhanced modeling of latent heat flux from urban surfaces in the Noah/single-layer urban canopy coupled model. Sci. China Earth Sciences, 57, doi: /s ,, Q. Li, and S. Fan, 2011: Impacts of urban processes and urbanization on summer precipitation: a case study of heavy rainfall in Beijing on 1 August J. Appl. Meteor. Climatol., 50, , J. Dou, F. Chen, J. Li, A. Li, 2012: Analysis of observations on the urban surface energy balance in Beijing. Sci. China Earth Sciences, 55, doi: /s Mitra, C., and J. M. Shepherd, Urban Precipitation: A global perspective. pp Mlawer, E. J., S. J. Taubman, P. D. Brown, M. J. Iacono, and S. A. Clough, 1997: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated model for the longwave. J. Geophys. Res., 102, , doi: /97JD Moriwaki, R., and M. Kanda, 2004: Seasonal and diurnal fluxes of radiation, heat, water vapor, and carbon dioxide over a suburban area. J. Appl. Meteorol., 43, Munck, C., G. Pigeon, V. Masson, F. Meunier, P. Bousquet, B. Tremeac, M. Merchat, P. Poeuf, and C. Marchadier, 2012: How much can air conditioning increase air temperatures for a city like Paris, France? Int. J. Climatol., 33, Niyogi, D., P. Pyle, M. Lei, S. Aalarya, M. Kishtawal, M. Shepherd, F. Chen, and B. Wolfe, 2011: Urban modification of thunderstorms: An observational storm climatology and model case study for the Indianapolis urban region. J. Appl. Meteor. Climatol., 50, NRC, 2010, When Weather Matters: Science and Services to Meet Critical Societal Needs. Wash, DC: National Research Council Report, National Academies Press. ISBN 13: and 10: Pleim, J. E., 2007: A combined local and non-local closure model for the atmospheric boundary layer. Part 1: Model description and testing. J. Appl. Meteor. Clim., 46, Quan, J., Q. Liu, X. Li, Y. Gao, X. Jia, J. Sheng, and Y. Liu, 2015: Effect of heterogeneous aqueous reactions on the secondary formation of inorganic aerosols during haze events. Atmos. Environ., 122, , and co-authors, 2013: Evolution of planetary boundary layer under different weather conditions, and its impact on aerosol concentrations. Particuology, 11, Robine, J-M, and co-authors, 2007: Report on excess mortality in Europe during summer EU Community Action Programme for Public Health, Grant Agreement , 2003 Heat Wave Project, 28 February. Robock, A., 1991: Surface cooling due to forest fire smoke. J. Geophys. Res., 96, Rosenfeld, D., and co-authors, 2008: Flood or drought: how do aerosols affect precipitation? Science, 321, Rosenfeld, D., J. Dai, X. Yu, Z. Yao, X. Xu, X. Yang, and C. Du, 2007: Inverse relations between amounts of air pollution and orographic precipitation, Science, 315, Sailor D. J., A. Brooks, M. Hart, and S. Heiple, 2007: A bottom-up approach for estimating latent and sensible heat emissions from anthropogenic sources, In Proc.7th Symp. on Urban Environment, American Meteorological Society, San Diego, CA. Salamanca, F., A. Krpo, A. Martilli, and A. Clappier, 2010: A new building energy model coupled with an urban canopy parameterization for urban climate simulations:part I. formulation, verification, and sensitivity analysis of the model. Theor. Appl. Climatol., 99, Shem, W., and J. M. Shepherd, 2009: On the impact of urbanization on summertime thunderstorms in Atlanta: two numerical model case studies. Atmos. Res., 92, Shepherd, J. M., H. Pierce, and A. J. Negri, 2002: Rainfall modification by major urban areas: Observations from space borne rain radar on the TRMM satellite. J. Appl. Meteor., 41, , doi: / (2002)041<0689: RMBMUA>2.0.CO;2., J. M., 2005: A review of current investigations of urban-induced rainfall and recommendations for the future. Earth Interact., 9, 1-27, doi: /EI156.1., 2013: Impacts of Urbanization on Precipitation and Storms: Physical Insights and Vulnerabilities. In: Climate Vulnerability: Understanding and Addressing Threats to Essential Resources , Elsevier Inc., Academic Press. Skamarock, W. C., and co-authors, 2008: A description of the advanced research WRF version 3. NCAR/TN-475+STR, 125 pp, wrf/users/docs/arw_v3.pdf

14 Stevens, A. S., 2008: ileaps Newsletter Issue No.5 April ileaps.org., D. H. Lenschow, M. A. LeMone, and L. Mahrt, 2016: The role of large- coherent-eddy transport in the atmospheric surface layer based on CASES-99 observations. J. Atmos. Sci., 73, , L. Mahrt, R. M. Banta, and Y. L. Pichugina, 2012: Turbulence regimes and turbulence intermittency in the stable boundary layer during CASES-99. J. Atmos. Sci., 69, Tan, J., and co-authors, 2015: Urban integrated meteorological observations: Practice and experience in Shanghai. China. Bull. Amer. Meteor. Soc., 96, Thompson, G., R. Rasmussen, and K. Manning, 2004: Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part I: Description & sensitivity analysis. Mon. Wea. Rev., 132, , doi: /2008MWR Tie, X., D. Wu, and G. Brasseur, 2009: Lung cancer mortality and exposure to atmospheric aerosol particles in Guangzhou. China. Atmos. Environ., 43, UN, 2014: World Urbanization Prospects: The 2014 Revision Highlights. United Nations Report ST/ESA/SER.A/352. Dept. Econ. Soc. Affairs, ISBN UN-Habitat, 2016: United Nations Human Settlements Programme: Urbanization and development: emerging futures, World Cities Report 2016, UN-Habitat. UNICEP, 2016: Clear the air for children. UNICEP Report UN037345, 100 pp, available at unicep.org/environment. Van den Heever, S. C., and W. R. Cotton, 2007: Urban aerosol impacts on downwind convective storms. J. Appl. Meteor. Climatol., 46, , doi: / jam Vandentorren, S., and co-authors, 2004: Mortality in 13 French cities during the August 2003 heat wave. Amer. J. Pub. Health, 94, Wang, T., and co-authors, 2010: Air quality during the 2008 Beijing Olympics: secondary pollutants and regional impact. Atmos. Chem. Phys., 10, Wang, X., S. Miao, J. Dou, F. Dong, and J. Wang, 2016: Observation and analysis of the air pollution impacts on radiation balance of urban and suburb areas in Beijing. Chinese J. Geophys., 59: , doi: /cjg (in Chinese with English abstract). Wild, O., X. Zhu, and M. J. Prather, 2000: Fast-J: accurate simulation of in and below cloud photolysis in tropospheric chemical models. J. Atmos. Chem., 37, WMO, 2013: The HIWeather program of the World Meteorological Organization, ml. WRF-Chem user s guide, 2017: Yu, M., and Y. Liu, 2015: The possible impact of urbanization on a heavy rainfall event in Beijing. J. Geophys. Res., 120, , doi: /2015jd , S. Miao, and Q. Li, 2017: Synoptic analysis and urban signatures of a heavy rainfall on 7 August 2015 in Beijing. J. Geophys. Res., 122, ,, and H. Zhang, 2018: Uncertainties in the numerical simulation of the urbanization impact on heavy rainfall: case study of a rainfall event in Beijing on 7 August 2015, J. Geophys. Res. (submitted) Zaveri, R.A., R. C. Easter, J. D. Fast, and L. K. Peters, 2008: Model for simulating aerosol interactions and chemistry (MOSAIC). J. Geo. Res., 113, D doi: /2007jd Zaveri, R. A., and L. K. Peters, 1999: A new lumped structure photochemical mechanism for long-scale applications. J. Geo. Res., 104, , doi: /1999JD Zhang, Y., S. Miao, and M. Chen, 2017a: Simulation experiments on high-resolution assimilation of radar thermodynamic data. Acta Meteorol. Sinica, 75, (in Chinese with English abstract).,, Y. Dai, and R. Bornstein, 2017b: Numerical simulation of urban land surface effects on summer convective rainfall under different UHI intensity in Beijing. J. Geophys. Res., 122, doi: /2017JD ,,, and Y. Liu, 2013: Numerical simulation of characteristics of summer clear day boundary layer in Beijing and the impact of urban underlying surface on sea breeze. Chinese J. Geophys., 56, , doi: /cjg (in Chinese with English abstract). Zhao, C., X. Tie, and Y. Lin, 2006: A possible positive feedback of reduction of precipitation and increase in aerosols over eastern central China. Geophys. Res. Lett., 33, L11814, doi: /2006gl Zheng Y., S. Miao, Q. Zhang, and Y. Bao, 2015: Improvements of building energy model and anthropogenic heat release from cooling system. Plateau Meteorolo., 34, doi: /j.issn (in Chinese w/english abstract). Zhong, S., Y. Qian, C. Zhao, R. Leung, and X. Q. Yang, 2015: A case study of urbanization impact on summer precipitation in the Greater Beijing Metropolitan Area: Urban heat island versus aerosol effects. J. Geophys. Res., 120, , doi: /2015JD

15 Fig. 1. Study area (a) full and (b) Beijing zoom-in (area of black box in a); with topography (m, color bar); urban areas (pink); province names (full words); Tian an Square (red star); PBL instrument sites (dots; three letter abbreviations defined in Supplement A); King Air flight tracks (dashes); mobile lidar route (purple), which starts from IAP ends at IUM; and Beijing Ring Roads 2 to 6 (light blue circles)

16 Fig. 2. Photographs of some SURF instrumentation (a) 325-m tower at IAP, (b) small balloon radiosonde at IUM, (c) GPS radiosonde at NAN, (d) X-band radar at SHU, (e) ceilometer at LIA, (f) Mini-MPL at DAX, (g) wind profiler at SID, (h) mobile wind profiler at CHA, (i) wind profiler at SHA, (j) wind profiler at YAN, (k) wind profiler at GUC, (l) King Air airplane, (m) Doppler lidar and Mini-MPL at IAP, and (n) mobile MPL, where three letter abbreviations are defined in Supplement A

17 NOWCASTING MODEL ( NOW ) RADAR, AWS INTEGRATION MODEL ( IN ) EC Global Forecast SHORT TERM MODEL (ST) URBAN MODEL AWS, QPE Objective Forecast s GTS, GPS, LU/LC, AWS RADAR CHEM MODEL Q*, χ Fig. 3. IUM RMAPS modeling system flow-chart, where input and outputs (blue boxes), component models (green boxes), and input data (yellow boxes) are discussed in the text, and where EC is the ECMWF model, Q* emission inventory, and χ concentration. Fig. 4. Wind speed V vs. friction velocity u * (each dot represents 30-min average value) at 47-m Beijing IAP tower level on nights following sunny days in 2015 for rough, medium, and smooth surfaces, i.e., wind direction ranges of: , , and o, respectively

18 Fig. 5. Doppler lidar Beam Swinging (DBS) scans of 30-min mean vertical-velocity variances (colors, m 2 s -2 ) on 6 July 2015, with CBL ( ) and NBL (+) depths from threshold and fractional methods, respectively (described in text); potential temperature (K) profile (solid line) from the 1315 LST NAN radiosonde; and sunrise and sunset (triangles); from Huang et al. (2016). Fig. 6. Concurrent cross sections of MPL normalized relative backscatter (NRB, colors; note different scales) and PBL height (magenta lines) from wavelet covariance transforms on 11 August 2015: (a) at IAP and (b) on the mobile route of Fig. 1b, where the red triangles delineate period in rural areas and rectangular bars indicate periods on Ring Roads 5 (green) and 6 (blue)

19 Fig. 7. Diurnal variations of average-summer observed 140-m level IAP-tower observed (OBS, black); and WRF Control (CTL, blue) and Cooling Tower (EC, red) simulated (a) sensible and (b) latent heat-fluxes, as well as EC anthropologic heatflux component (AH, green); and 2-m average (for eight urban stations): (c) temperature and, (d) specific humidity. Fig. 8. Precipitation centers on 22 July 2015 every 15 min (colors) for the: (a) Urban and (b) No Urban simulations (min outline values 3-8 mm) and (c) observations (min outline values 1-5 mm), where grays are topographic height (m), black rectangle is Ring Road 4, short arrow is storm-movement direction, and long arrows are storm-center trajectories

20 Fig. 9. Average hourly observed PM 2.5 concentration (μg m -3 ) from 26 November to 1 December 2015 from 35 Beijing, 28 Tianjin, and 27 Shijia-zhuang urban sites. Fig. 10. Observed Beijing-area 10-m surface wind (1 full barb is 1 m s -1 ) and PM concentration (colors, μg m -3 ) during November 2015 at: (a) 20 LST on the 29th; and at (b) 08, (c) 14, and (d) 20 LST on the 30 th, with Ring Roads 4 and 6 (white circles)

21 Fig. 11. Hourly observations on 30 November 2015 of: (a) PM 2.5 concentration (μg m -3 ) at HAI, and (b) PBL horizontal wind vectors (1 full barb is 2 m s -1 ) and dbcarrier-to-noise Ratios (CNR, colors, db) from IAP Doppler lidar. The db-units version of CNR used herein is CNR(db) = 10 log 10 (CNR), and as CNR(db) is less than unity, the values in the figure are negative powers of 10. PM 2.5 ( μg m - 3 ) Fig. 12. Simulated Beijing-area 10-m wind (1 full barb is 2 m s -1 ) and PM 2.5 concentration (shades, μg m -3 ) during November 2015 at: (a) 20 LST on the 29 th, and (b) 02, (c) 08, and (d) 20 LST on the 30th; with Ring Roads 4 and 6 (white circles)