Observational Field Campaign Research Articles

Integrating High-Resolution Modeling with Diverse Field Campaign Observations through an Open Science Summer School

Integrating High-Resolution Modeling with Diverse Field Campaign Observations through an Open Science Summer School

Evaluation of near‐surface and boundary‐layer meteorological conditions that support cold‐fog formation using Cold Fog Amongst Complex Terrain field campaign observations

AbstractCold fog refers to a type of fog that forms when the temperature is below 0°C. It can be composed of liquid, ice, and mixed‐phase fog particles. Cold fog happens frequently over mountainous terrain in the cold season, but it is difficult to predict. Using observations from the Cold Fog Amongst Complex Terrain (CFACT) field campaign conducted in Heber Valley, Utah, in the western United States during January and February of 2022, this study investigates the meteorological conditions in the surface and boundary layers that support the formation of wintertime ephemeral cold fog in a local area of small‐scale mountain valleys. It is found that fog formation is susceptible to subtleties in forcing conditions and is supported by several factors: (1) established high pressure over the Great Basin with associated local clear skies, calm winds, and a stable boundary layer; (2) near‐surface inversion with saturation near the surface and strong moisture gradient in the boundary layer; (3) warm (above‐freezing) daytime air temperature with a large diurnal range, accompanied with warm soil temperatures during the daytime; (4) a period of increased turbulence kinetic energy (above 0.5 m2·s−2), followed by calm conditions throughout the fog's duration; and (5) supersaturation with respect to ice. Then, the field observations and identified supporting factors for fog formation were utilized to evaluate high‐resolution (˜400 m horizontal grid spacing) Weather Research and Forecasting (WRF) model simulations. Results show that the WRF model accurately simulates the mesoscale conditions facilitating cold‐fog formation but misses some critical surface and atmospheric boundary conditions. The overall results from this paper indicate that these identified factors that support fog formation are vital to accurately forecasting cold‐fog events. At the same time, they are also critical fields for the NWP model validation.

Importance of CCN activation for fog forecasting and its representation in the two‐moment microphysical scheme LIMA

AbstractThe work presented in this article studies the impact of cloud condensation nuclei (CCN) activation for fog forecasting and improves its parameterization in the LIMA (Liquid, Ice, Multiple Aerosols) two‐moment microphysical scheme, building upon the Local And Non‐local Fog EXperiment (LANFEX) field campaign observations, specifically the intensive observation period (IOP) 1 and the DEMISTIFY intercomparison. Large‐eddy simulations were performed with the Meso‐NH model, first using a prognostic supersaturation allowing us to compute the number of activated CCN at each time step, and then with the usual saturation adjustment hypothesis and a diagnostic maximum supersaturation. The prognostic supersaturation method provided very good results, similar to those from earlier simulations of this case using a bin scheme, and was thus used as a reference. In contrast, the diagnostic maximum supersaturation method strongly overestimated droplet numbers and produced a too‐thick fog. Thus, improvements to the maximum supersaturation diagnostic were proposed, by (1) revising the temperature tendency and (2) accounting for pre‐existing cloud droplets in the activation parameterization. These improvements resulted in a simulation in good agreement with observations and the reference simulation, and are promising for use in numerical weather prediction systems with a lower resolution and/or a longer time step.

CALIPSO Observed Low Clouds Over the Subtropical Northeast Pacific: Summer Climatology and Interannual Variability

AbstractThis study examines the spatial distribution and temporal variations of low‐cloud top over the subtropical northeast Pacific (NEP), using observations from the Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). A comparison between CALIPSO observations and in‐situ soundings reveals a distinct bimodal distribution of low‐cloud tops, characterized by peaks near the inversion base and transition layer, corresponding to shallow cumulus clouds below the stratocumulus deck and boundary‐layer decoupling phenomena. Climatologically, while traversing along a great circle from Los Angeles to Honolulu, the frequency of low cloud occurrence experiences a gradual reduction from 80% to 30%. Meanwhile, cloud top height increases from 0.6 to 1.6 km, and surface roughness varies from 50 to 300 m. This study marks the first presentation of long‐term characteristics of low clouds over a vast, remote area, expanding beyond discrete field campaign observations. Accurately characterizing the marine atmospheric boundary layer (MABL) including its depth and degree of decoupling, enables further investigation into the interannual variations of low clouds, using 14 years of CALIPSO observations. Significant interannual variability in areal fraction and vertical bimodality of low clouds emerges downstream of the anomalously cold water; increased cloud cover with weakened bimodality is associated with a shallow weakly decoupled MABL. This research emphasizes the need to consider the downstream effects in studying low cloud‐sea surface temperature feedback within climate models.

Validation of the Geostationary Lightning Mapper With a Lightning Mapping Array in Argentina: Implications for Current and Future Spaceborne Lightning Observations

AbstractA validation study of the Geostationary Lightning Mapper (GLM) on board the Geostationary Operational Environmental Satellite 16 (GOES‐16) was done using a ground‐based lightning mapping array (LMA) deployed as part of the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign in Argentina. GLM detected lightning with 74.6% efficiency over 61 thunderstorm days in December 2018 through April 2019. However, GLM detection efficiency (DE) was negatively correlated (r = −0.49) with LMA flash rate. GLM DE also was negatively correlated with LMA flash altitude (r = −0.24), reflecting the influence of multiple competing trends. GLM DE was positively correlated (r = 0.27) with number of LMA sources in a flash, indicating improved DE for larger flashes. During periods with anomalously electrified storms, GLM DE was reduced to 50.9%. Statistics were found to be sensitive to analysis criteria, but most of the above trends remained consistent regardless of specific criteria. Because the methodology allowed a GLM flash to match more than one LMA flash, actual GLM flash rate was a factor of 2.9 lower than the LMA flash rate, and this ratio grew larger as LMA flash rate increased. A sensitivity study examined the impact of improved DE for smaller flashes; that is, an improved sensor (or algorithm) that was better able to detect and distinguish between separate small lightning flashes. The results showed improved correlation with LMA flash rates, as well as improved ability to identify lightning jumps associated with intensifying convection.

Characterization of wind speed and directional shear at the AWAKEN field campaign site

The American wake experiment (AWAKEN) is taking place in northern Oklahoma, USA, close to the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) atmospheric observatory. The planning for the deployment of the instruments in this observational field campaign required an assessment of the wind characteristics of the site. This paper analyzes long-term data collected by instruments at the ARM SGP observatory to characterize the winds near the AWAKEN site. The analysis shows that this site experiences high wind shear and veer events with a large number of nocturnal low-level jets. A total of 7086 low-level jet wind profiles over 6 years are examined and found to be dominant from the south and southeast. Significant nocturnal wind veer is observed, which causes southerly wind near the surface to become westerly wind aloft. By identifying a strong relationship between atmospheric stability and wind shear, the wind shear at the site is predicted using the Monin–Obukhov similarity theory (MOST) and validated with the observational data collected by a scanning Doppler lidar. The results show that wind speed at a height of 91 m, a proxy hub height for wind turbines in this area, can be predicted from data collected at a height of 10 m with a bias of −0.35 and 0.65 m s−1 in unstable and stable atmospheric boundary layers, respectively. The bias of the predicted wind speed is mostly in the region of low wind speed, and wind speed above 5 m s−1 at a height of 91 m can be predicted with a bias of less than 0.2 m s−1, and the limitations of the MOST in predicting winds during the stably stratified boundary layer is well-observed.

The Uncertainty of IMERG Over the Western Edge of the Eastern Pacific Fresh Pool: An Error Model Based on SPURS-2 Field Campaign Observations

The Uncertainty of IMERG Over the Western Edge of the Eastern Pacific Fresh Pool: An Error Model Based on SPURS-2 Field Campaign Observations

Hourly assimilation of different sources of observations including satellite radiances in a mesoscale convective system case during RELAMPAGO campaign

This paper evaluates the impact of assimilating high-resolution surface networks and satellite observations using the WRF-GSI-LETKF over central and north eastern Argentina where the surface and upper air observing networks are relatively coarse. A case study corresponding to a huge mesoscale convective system (MCS) that developed during November 22, 2018 was used. The accumulated precipitation associated with this MCS was quite high, exceeding 200 mm over northern Argentina and Paraguay. The MCS developed during the Intense Observing Period (IOP) of the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign. The GSI-4DLETKF data assimilation package is used to produce analyses by assimilating observations every hour with 10-km horizontal grid spacing and a 60-member multiphysics ensemble. Four assimilation experiments are conducted using different sets of observations: CONV, consisting of conventional observations from NCEP’s prepBUFR files; AWS, combining CONV and dense automatic surface weather station networks (AWS), SATWND, combining AWS with satellite-derived winds, and RAD, including SATWND; and satellite radiances from different microwave and infrared sensors. The assimilation of observations with high temporal and spatial frequency generates an important impact on the PBL, primarily on the precipitable water content, that leads to the development of deep convection and heavy precipitation closer to the observed in this case study. The assimilation of radiance observations produces a better development of the convection mainly during the mature state of the MCS leading to an increase in the accumulated precipitation. Ensemble forecasts initialized from each experiment were also simulated to evaluate their skill to predict precipitation. The hourly assimilation of the observations in AWS, SATWND, and RAD helped to improve the precipitation forecast.

Measurements of Size and Electrical Charges Carried by Precipitation Particles During RELAMPAGO Field Campaign

AbstractThe electrical charge carried by raindrops provides significant information about thunderstorm electrification mechanisms, since the charge acquired by hydrometeors is closely related to the microphysical processes that they undergo within clouds. Investigation of charges on raindrops was conducted during the Remote sensing of Electrification, Lightning, And Meso‐scale/micro‐scale Processes with Adaptive Ground Observations field campaign. A newly designed instrument was used to determine simultaneously the fall velocity and charge for precipitating particles. Hydrometeor size and charge were measured in Córdoba city, Argentina, during electrified storms. Temporal series of size‐charge of single raindrops were recorded for two storms, which were also monitored with a Parsivel disdrometer and Lightning Mapping Array. The results show that the magnitude of the electric charges range between 1 and 50 pC and more than 90% of the charges are mainly carried by raindrops >1 mm, even though most of the raindrops are smaller than 1 mm. Furthermore, the measurement series show charged hydrometeors of both signs all the time. A correlation between the sizes and the charges carried by the raindrops was found in both storms.

Elevated CO2 moderates the impact of climate change on future bamboo distribution in Madagascar

The distribution of bamboo is sensitive to climate change and is also potentially affected by increasing atmospheric CO2 concentrations due to its C3 photosynthetic pathway. Yet the effect of CO2 in climate impact assessments of potential changes in bamboo distribution has to date been overlooked. In this study, we proposed a simple and quantitative method to incorporate the impact of atmospheric CO2 concentration into a species distribution modeling framework. To do so, we implemented 10 niche modeling algorithms with regionally downscaled climatic variables and combined field campaign observations. We assessed future climate impacts on the distribution of an economically and ecologically important and widely distributed bamboo species in Madagascar, and examined the effect of increasing CO2 on future projections. Our results suggested that future climatic changes negatively impact potential bamboo distribution in Madagascar, leading to a decline of 34.8% of climatic suitability and a decline of 63.6 ± 3.2% in suitable areas towards 2100 under RCP 8.5. However, increasing atmosphere CO2 offsets the climate impact for bamboo, and led to a smaller reduction of 19.8% in suitability and a potential distribution expansion of +111.6 ± 9.8% in newly suitable areas. We also found that the decline in climatic suitability for bamboo was related to increasing monthly potential evapotranspiration of the warmest quarter and minimum temperature of the warmest month. Conversely, the decreasing isothermality and increasing precipitation of the warmest quarter contributed to projected increase in bamboo-suitable areas. Our study suggested that elevated CO2 may mitigate the decrease in climatic suitability and increase bamboo-suitable areas, through enhancing water use efficiency and decreasing potential evapotranspiration. Our results highlight the importance of accounting for the CO2 effect on future plant species distributions, and provide a mechanistic approach to do so for ecosystems constrained by water.

Easterly Waves in the East Pacific during the OTREC 2019 Field Campaign

Abstract Easterly waves (EWs) are off-equatorial tropical synoptic disturbances with a westward phase speed between 11 and 14 m s−1. Over the east Pacific in boreal summer, the combination of EWs and other synoptic disturbances, plus local mechanisms associated with sea surface temperature (SST) gradients, define the climatological structure of the intertropical convergence zone (ITCZ). The east Pacific ITCZ has both deep and shallow convection that is linked to deep and shallow meridional circulations, respectively. The deep convection is located around 9°N over warm SSTs. The shallow convection is located around 6°N and is driven by the meridional SST gradient south of the ITCZ. This study aims to document the interaction between east Pacific EWs and the deep and shallow meridional circulations during the Organization of Tropical East Pacific Convection (OTREC) field campaign in 2019 using field campaign observations, ERA5, and satellite precipitation. We identified three EWs during the OTREC period using precipitation and dynamical fields. Composite analysis shows that the convectively active part of the EW enhances ITCZ deep convection and is associated with an export of column-integrated moist static energy (MSE) by vertical advection. The subsequent convectively suppressed, anticyclonic part of the EW produces an increase of moisture and column-integrated MSE by horizontal advection that likely enhances shallow convection and the shallow overturning flow at 850 hPa over the southern part of the ITCZ. Therefore, EWs appear to strongly modulate shallow and deep circulations in the east Pacific ITCZ.

Enlightenment Strikes! Broadening Graduate School Training through Field Campaign Participation

AbstractThis article provides an overview of the Advanced Study Institute: Field Studies of Convection in Argentina (ASI-FSCA) program, a 3-week dynamic and collaborative hands-on experience that allowed 16 highly motivated and diverse graduate students from the United States to participate in the 2018–19 Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign. This program is unique as it represents the first effort to integrate an intensive Advanced Study Institute with a field campaign in atmospheric science. ASI-FSCA activities and successful program outcomes for five key elements are described: 1) intensive field research with field campaign instrumentation platforms; 2) recruitment of diverse graduate students who would not otherwise have opportunities to participate in intensive field research; 3) tailored curriculum focused on scientific understanding of cloud and mesoscale processes and professional/academic development topics; 4) outreach to local K–12 schools and the general public; and 5) building a collaborative international research network to promote weather and climate research. These five elements served to increase motivation and improve confidence and self-efficacy of students to participate in scientific research and field work with goals of increasing retention and a sense of belonging in STEM graduate programs and advancing the careers of students from underrepresented groups as evidenced by a formal program evaluation effort. Given the success of the ASI-FSCA program, our team strongly recommends considering this model for expanding the opportunities for a broader and more diverse student community to participate in dynamic and intensive field work in atmospheric science.

A rapid refresh ensemble based data assimilation and forecast system for the RELAMPAGO field campaign

This paper describes the lessons learned from the implementation of a regional ensemble data assimilation and forecast system during the intensive observing period of the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign (central Argentina, November–December 2018). This system is based on the coupling of the Weather Research and Forecasting (WRF) model and the Local Ensemble Transform Kalman Filter (LETKF). It combines multiple data sources both global and locally available like high-resolution surface networks, AMDAR data from local aircraft flights, soundings, AIRS retrievals, high-resolution GOES-16 wind estimates, and local radar data. Hourly analyses with grid spacing of 10 km are generated along with warm-start 36-h ensemble-forecasts, which are initialized from the rapid refresh analyses every three hours. A preliminary evaluation shows that a forecast error reduction is achieved due to the assimilated observations. However, cold-start forecasts initialized from the Global Forecasting System Analysis slightly outperform the ones initialized from the regional assimilation system discussed in this paper. The system uses a multi-physics approach, focused on the use of different cumulus and planetary boundary layer schemes allowing us to conduct an evaluation of different model configurations over central Argentina. We found that the best combinations for forecasting surface variables differ from the best ones for forecasting precipitation, and that differences among the schemes tend to dominate the forecast ensemble spread for variables like precipitation. Lessons learned from this experimental system are part of the legacy of the RELAMPAGO field campaign for the development of advanced operational data assimilation systems in South America.

A Storm Safari in Subtropical South America: Proyecto RELAMPAGO

AbstractThis article provides an overview of the experimental design, execution, education and public outreach, data collection, and initial scientific results from the Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign. RELAMPAGO was a major field campaign conducted in the Córdoba and Mendoza provinces in Argentina and western Rio Grande do Sul State in Brazil in 2018–19 that involved more than 200 scientists and students from the United States, Argentina, and Brazil. This campaign was motivated by the physical processes and societal impacts of deep convection that frequently initiates in this region, often along the complex terrain of the Sierras de Córdoba and Andes, and often grows rapidly upscale into dangerous storms that impact society. Observed storms during the experiment produced copious hail, intense flash flooding, extreme lightning flash rates, and other unusual lightning phenomena, but few tornadoes. The five distinct scientific foci of RELAMPAGO—convection initiation, severe weather, upscale growth, hydrometeorology, and lightning and electrification—are described, as are the deployment strategies to observe physical processes relevant to these foci. The campaign’s international cooperation, forecasting efforts, and mission planning strategies enabled a successful data collection effort. In addition, the legacy of RELAMPAGO in South America, including extensive multinational education, public outreach, and social media data gathering associated with the campaign, is summarized.

Characterizing Charge Structure in Central Argentina Thunderstorms During RELAMPAGO Utilizing a New Charge Layer Polarity Identification Method.

A new automated method to retrieve charge layer polarity from flashes, named Chargepol, is presented in this paper. Using data from the NASA Lightning Mapping Array (LMA) deployed during the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign in Cordoba, Argentina, from November 2018 to April 2019, this method estimates the polarity of vertical charge distributions and their altitudes and thicknesses (or vertical depth) using the very‐high frequency (VHF) source emissions detected by LMAs. When this method is applied to LMA data for extended periods of time, it is capable of inferring a storm's bulk electrical charge structure throughout its life cycle. This method reliably predicted the polarity of charge within which lightning flashes propagated and was validated in comparison to methods that require manual assignment of polarities via visual inspection of VHF lightning sources. Examples of normal and anomalous charge structures retrieved using Chargepol for storms in Central Argentina during RELAMPAGO are presented for the first time. Application of Chargepol to five months of LMA data in Central Argentina and several locations in the United States allowed for the characterization of the charge structure in these regions and for a reliable comparison using the same methodology. About 13.3% of Cordoba thunderstorms were defined by an anomalous charge structure, slightly higher than in Oklahoma (12.5%) and West Texas (11.1%), higher than Alabama (7.3%), and considerably lower than in Colorado (82.6%). Some of the Cordoba anomalous thunderstorms presented enhanced low‐level positive charge, a feature rarely if ever observed in Colorado thunderstorms.

Short-Term Response of Chlorophyll a Concentration Due to Intense Wind and Freshwater Peak Episodes in Estuaries: The Case of Fangar Bay (Ebro Delta)

Estuaries and coastal bays are areas of large spatio-temporal variability in physical and biological variables due to environmental factors such as local wind, light availability, freshwater inputs or tides. This study focuses on the effect of strong wind events and freshwater peaks on short-term chlorophyll a (Chl a) concentration distribution in the small-scale and microtidal, Fangar Bay (Ebro Delta, northwestern Mediterranean). The hydrodynamics of this bay are primarily driven by local wind episodes modulated by stratification in the water column. Results based on field-campaign observations and Sentinel-2 images revealed that intense wind episodes from both NW (offshore) and NE-E (onshore) caused an increase in the concentration of surface Chl a. The mechanisms responsible were horizontal mixing and the bottom resuspension (also linked to the breakage of the stratification) that presumably resuspended Chl a containing biomass (i.e., micropyhtobentos) and/or incorporated nutrients into the water column. On the other hand, sea-breeze was not capable of breaking up the stratification, so the chlorophyll a concentration did not change significantly during these episodes. It was concluded that the mixing produced by the strong winds favoured an accumulation of Chl a concentration, while the stratification that causes a positive estuarine circulation reduced this accumulation. However, the spatial-temporal variability of the Chl a concentration in small-scale estuaries and coastal bays is quite complex due to the many factors involved and deserve further intensive field campaigns and additional numerical modelling efforts.

Cross Validation of the Network of Ground-Based Radar with GPM during the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) Field Campaign

Cross Validation of the Network of Ground-Based Radar with GPM during the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) Field Campaign

The Influence of Terrain on the Convective Environment and Associated Convective Morphology from an Idealized Modeling Perspective

AbstractOrographic deep convection (DC) initiation and rapid evolution from supercells to mesoscale convective systems (MCSs) are common near the Sierras de Córdoba, Argentina, which was the focal point of the Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign. This study used an idealized numerical model with elongated north–south terrain similar to that of the Sierras de Córdoba to address how variations in terrain height affected the environment and convective morphology. Simulations used a thermodynamic profile from a RELAMPAGO event that featured both supercell and MCS storm modes. Results revealed that DC initiated earlier in simulations with higher terrain, owing both to stronger upslope flows and standing mountain waves. All simulations resulted in supercell formation, with higher-terrain supercells initiating closer to the terrain peak and moving slower off the terrain. Higher-terrain simulations displayed increases in both low-level and deep-layer wind shear along the eastern slopes of the terrain that were related to the enhanced upslope flows, supporting stronger and wider supercell updrafts/downdrafts and a wider swath of heavy rainfall. Deeper and stronger cold pools from these wider and stronger higher-terrain supercells led to surging outflow that reduced convective available potential energy accessible to deep convective updrafts, resulting in quicker supercell demise off the terrain. Lower-terrain supercells moved quickly off the terrain, merged with weaker convective cells, and resulted in a quasi-organized MCS. These results demonstrate that terrain-induced flow modification may lead to substantial local variations in convective morphology.

Analysis of Back-Building Convection in Simulations with a Strong Low-Level Stable Layer

Abstract In a mesoscale convective system (MCS), convection that redevelops over (i.e., back-builds), and/or repeatedly passes over (i.e., trains) a region for an extended period of time can contribute to extreme rainfall and flash flooding. Past studies have indicated that both mesoscale ascent and lifting of the inflow layer by a cold pool or bore are important when this back-building/training convection is displaced from the leading line [sometimes called rearward off-boundary development (ROD)]. However, Plains Elevated Convection At Night (PECAN) field campaign observations suggest that the stability of the nocturnal boundary layer is highly variable and some MCSs with ROD have only a weak surface cold pool. Numerical simulations presented in this study suggest that in an environment with strong boundary layer stability, ROD can be supported by mechanisms other than those mentioned above. Simulations were initialized using a sounding from ahead of a PECAN MCS with a strong stable layer and ROD, and the three-dimensional simulation produced an MCS similar to that observed despite the homogeneous initial conditions. Some of the findings presented herein challenge existing understanding of nocturnal MCSs, and especially how downdrafts interact with a stable boundary layer. Notably, downdrafts can reach the surface, and different regions of the MCS may have different propagation mechanisms and different relevant inflow layers. Unlike previous studies of ROD, parcel lifting may be supported by an intrusion (an elevated layer of downdraft air) modified by the three-dimensional vertical wind shear.

Understanding the formation of high-ozone episodes at Raoyang, a rural site in the north China plain

Abstract An observational field campaign was carried out from 20 June to July 31, 2016 at Raoyang, a rural site in the North China Plain (NCP). In-situ measurements of O3 and related gases were obtained and analyzed in this work. Heavy photochemical pollution episodes lasting for a few days were observed in the campaign. The episodes could be divided into two categories, with one being influenced by plume from northeast (NE) via Beijing-Tianjin cities and the other from southeast (SE) via Shandong province. The SE episodes were more frequent with higher O3 concentrations and lagged daily peaks of O3 despite of lower concentrations of precursors. Such episodes could last for several days or even a week. In contrast, the NE episode was associated with relatively lower O3 levels and earlier daily O3 peaks though the levels of precursors were higher. The WRF-Chem model was used to study the differences between typical SE and NE episodes in regional distributions of O3 concentrations, chemical regimes of O3 formation, ozone production efficiency (OPE), photochemical age of plume, and evolution of O3 concentration within the boundary layer. Our results suggest that the O3 production regimes were different for the Raoyang site, with the NE episode being VOCs-sensitive and the SE episode in the VOC-sensitive to NOx-sensitive transition regime. In both episodes, a major part of the NCP was under the VOCs-sensitive regime. Plumes in the SE episodes were more aged than those in the NE episode. The OPE values were higher in more aged SE plumes, resulting in higher O3 concentrations despite of lower precursor concentrations. Our results also indicate that downward mixing of O3 in the nocturnal residual layer played an important role in rapid buildup of surface O3 in the next morning in both NE and SE episodes. The formation of photochemical smog seems to be particularly effective in the SE episode since O3 is rapidly produced during daytime, stored aloft at night, and mixed downward to the ground-level, and the process can circulate for several consecutive days. Our study highlights the complicated role of chemical reaction, transport and boundary layer processes in the formation of severe photochemical smog in the region.