Improved Tropical Cyclone Forecasts with Increased Vertical Resolution in the TRAMS Model

The numerical simulation of typhoons has been found to be very sensitive to the vertical resolution of the model. During the updating of the TRAMS model from version 1.0 to 3.0, the horizontal resolution has been increased from 36 km to 9 km, while the vertical layer number only increased from 55 to 65 layers. The lack of high vertical resolution limits the performance of the TRAMS model in typhoon forecasting to a certain extent. In order to study the potential improvement of typhoon forecasting by increasing the vertical resolution, this paper increases the vertical resolution of the TRAMS model from 65 to 125 layers for the first time for a comparative simulation test. The results of the case study with Typhoon Hato (2017) show that the model with high vertical resolution can significantly enhance the warm structure caused by water vapor flux convergence and vertical transport, thus accurately simulating the rapid strengthening process of the typhoon. Meanwhile, the model with 125-layer vertical resolution can simulate the asymmetric structural characteristics of the wind field, which are closer to the observations and can help to reduce the bias in typhoon track forecasting. The improvement of vertical resolution is also trialed by using the batch test results of several landfalling typhoons in 2016-2017. The experimental results show that the typhoon forecast of the model becomes consistent with the observations only when the number of vertical layers of the model increases to about 125 layers, which in turn causes a large computational burden. In the next step, we will try to solve the computational burden problem caused by ultra-high vertical resolution with the top boundary nesting technique, and realize the application of high vertical resolution in the actual operation of the TRAMS model.     

Simulation of the Effect of Water-vapor Increase on Temperature in the Stratosphere

To analyze the mechanism by which water vapor increase leads to cooling in the stratosphere, the effects of water-vapor increases on temperature in the stratosphere were simulated using the two-dimensional, interactive chemical dynamical radiative model (SOCRATES) of NCAR. The results indicate that increases in stratospheric water vapor lead to stratospheric cooling, with the extent of cooling increasing with height, and that cooling in the middle stratosphere is stronger in Arctic regions. Analysis of the radiation process showed that infrared radiative cooling by water vapor is a pivotal factor in middle-lower stratospheric cooling. However, in the upper stratosphere (above 45 km), infrared radiation is not a factor in cooling; there, cooling is caused by the decreased solar radiative heating rate resulting from ozone decrease due to increased stratospheric water vapor. Dynamical cooling is important in the middle-upper stratosphere, and dynamical feedback to temperature change is more distinct in the Northern Hemisphere middle-high latitudes than in other regions and signiffcantly affects temperature and ozone in winter over Arctic regions. Increasing stratospheric water vapor will strengthen ozone depletion through the chemical process. However, ozone will increase in the middle stratosphere. The change in ozone due to increasing water vapor has an important effect on the stratospheric temperature change.     

Variations of sea ice temperature from CHINARE 2003 and its application on sea ice model evaluation

Variation of vertical profiles of sea ice temperature and adjacent atmosphere and ocean temperatures were measured by ice drifting buoys deployed in the northeast Chukchi Sea as part of the 2003 Chinese Arctic Research Expedition.The buoy observations (September 2003 to February 2005) show that the cooling of the ice began in late September,propagated down through the ice,reaching the bottom of the ice in December,and continued throughout the winter.In winter 2003/04,some obvious warmings were observed in the upper portion of the ice in response to major warmings in the overlying atmosphere associated with the periodicity of storms in the northeast Chukchi Sea.It is found that the melt season at the buoy site in 2004 was about 15% longer than normal.The buoy observed vertical ice temperature profiles were used as a diagnostic for sea ice model evaluation.The results show that the simulated ice temperature profiles have large discrepancies as compared with the observations.     

RETRIEVING ATMOSPHERIC SOUNDING PROFILES AROUND TYPHOON YUNNA USING INFRARED HYPERSPECTRAL MEASUREMENTS AIRS

In this study,we derived atmospheric profiles of temperature,moisture,and ozone,along with surface emissivity,skin temperature,and surface pressure,from infrared-sounder radiances under clear sky (cloudless) condition.Clouds were detected objectively using the Atmospheric Infrared Sounder under a relatively low spatial resolution and cloud-mask information from the Moderate Resolution Imaging Spectroradiometer under a high horizontal resolution;this detection was conducted using space matching.Newton’s nonlinear physical iterative solution technique is applied to the radiative transfer equation (RTE) to retrieve temperature profiles,relative humidity profiles,and surface variables simultaneously.This technique is carried out by using the results of an eigenvector regression retrieval as the background profile and using corresponding iterative forms for the weighting functions of temperature and water-vapor mixing ratio.The iterative forms are obtained by applying the variational principle to the RTE.We also compared the retrievals obtained with different types of observations.The results show that the retrieved atmospheric sounding profile has great superiority over other observations by accuracy and resolution.Retrieved profiles can be used to improve the initial conditions of numerical models and used in areas where conventional observations are sparse,such as plateaus,deserts,and seas.     

MODIS brightness temperature data assimilation under cloudy conditions Ⅱ: Impacts on rainstorm forecasting

Although satellite observations provide large amount of information of clouds and precipitation and play an important role in the forecast of heavy rainfall, they have not been fully taken advantage of in data assimilation of numerical weather predictions, especially those in infrared channels. Assimilating radiances is common under clear-sky conditions since it is extremely difficult to simulate infrared transmittance in cloudy sky. Based on the Global and Regional Assimilation and Prediction Enhanced System 3-dimensional variance (GRAPES-3DVar), cloud liquid water content, ice-water content and cloud cover are employed as governing variables in the assimilation system. This scheme can improve the simulation of infrared transmittance by a fast radiative transfer model for TOVS (RTTOV) and adjust the atmospheric and cloud parameters based on infrared radiance observations. In this paper, we investigate a heavy rainfall over Guangdong province on May 26, 2007, which is right after the onset of a South China Sea monsoon. In this case, channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) for observing water vapor (Channel 27) and cloud top altitude (Channel 36) are selected for the assimilation. The process of heavy rainfall is simulated by the Weather Research and Forecasting (WRF) model. Our results show that the assimilated MODIS data can improve the distribution of water vapor and temperature in the first guess field and indirectly adjust the upper-level wind field. The tendency of adjustment agrees well with the satellite observations. The assimilation scheme has positive impacts on the short-range forecasting of rainstorm.     

Application of the Conditional Nonlinear Optimal Perturbation Method to the Predictability Study of the Kuroshio Large Meander

A reduced-gravity barotropic shallow-water model was used to simulate the Kuroshio path variations.The results show that the model was able to capture the essential features of these path variations.We used one simulation of the model as the reference state and investigated the effects of errors in model parameters on the prediction of the transition to the Kuroshio large meander (KLM) state using the conditional nonlinear optimal parameter perturbation (CNOP-P) method.Because of their relatively large uncertainties,three model parameters were considered:the interfacial friction coefficient,the wind-stress amplitude,and the lateral friction coefficient.We determined the CNOP-Ps optimized for each of these three parameters independently,and we optimized all three parameters simultaneously using the Spectral Projected Gradient 2 (SPG2) algorithm.Similarly,the impacts caused by errors in initial conditions were examined using the conditional nonlinear optimal initial perturbation (CNOP-I) method.Both the CNOP-I and CNOP-Ps can result in significant prediction errors of the KLM over a lead time of 240 days.But the prediction error caused by CNOP-I is greater than that caused by CNOP-P.The results of this study indicate not only that initial condition errors have greater effects on the prediction of the KLM than errors in model parameters but also that the latter cannot be ignored.Hence,to enhance the forecast skill of the KLM in this model,the initial conditions should first be improved,the model parameters should use the best possible estimates.     

The Effects of Assimilating Satellite Brightness Temperature and Bogus Data on the Simulation of Typhoon Kalmaegi (2008)

Observational and bogus satellite data are directly assimilated into the Weather Research and Forecasting (WRF) model in simulations of Typhoon Kalmaegi (2008). The data assimilation is performed using the Radiative Transfer for TIROS-N Operational Vertical Sounder (RTTOV) model and the three-dimensional variational data assimilation (3DVAR) technique, with satellite observations taken from the National Oceanic and Atmospheric Administration-16 (NOAA-16) Advanced TIROS Vertical Sounder (ATOVS) system composed of the High-resolution Infrared Radiation Sounder (HIRS), the Advanced Microwave Sounding Unit-A (AMSU-A), and the Advanced Microwave Sounding Unit-B (AMSU-B). Data assimilation experiments are initialized at three different times. Improvements in the numerical simulation of the typhoon are discussed in the context of wind, temperature, pressure, and geopotential fields. The results indicate that assimilation of satellite data can improve both the representation of the initial conditions and the subsequent simulation of the typhoon. Different satellite data have different impacts on the typhoon track. In these simulations, data from AMSU-A play a greater role in improving the simulation of the typhoon than data from AMSU-B or HIRS. Assimilation of satellite data significantly affects the simulation of the subtropical high and the steering of the typhoon by the environmental flow. The subtropical high is enhanced and extends westward in the data assimilation experiments. The background flow therefore steers the typhoon more westward, improving the simulated typhoon track. Although direct assimilation of satellite brightness temperature improves the simulated environmental conditions, it does not significantly improve the simulated intensity of the typhoon. By contrast, initializing the typhoon simulation using bogus data in tandem with satellite data improves not only the environmental conditions but also the simulated inner-core structure of the typhoon. Assimilation of both types of data therefore improves the simulation of both the typhoon track and the typhoon intensity. The results of these experiments offer new insight into improving numerical simulations of typhoons.     

Impacts of Initial Conditions on Cloud-Resolving Model Simulations

Impacts of initial conditions on cloud-resolving model simulations are investigated using a series of sensitivity experiments. Five experiments with perturbed initial temperature, moisture, and cloud conditions are conducted and compared to the control experiment. The model is forced by the large-scale vertical velocity and zonal wind observed and derived from NCEP/Global Data Assimilation System （GDAS）. The results indicate that model predictions of rainfall are much more sensitive to the initial conditions than those of temperature and moisture. Further analyses of the surface rainfall equation and the moisture and cloud hydrometeor budgets reveal that the calculations of vapor condensation and deposition rates in the model account for the large sensitivities in rainfall simulations.     

FY-3A SATELLITE MICROWAVE DATA ASSIMILATION EXPERIMENTS IN TROPICAL CYCLONE FORECAST

China’s new generation of polar-orbiting meteorological satellite FY-3A was successfully launched on May 26,2008,carrying microwave sounding devices which had similar performance to ATOVS of NOAA series.In order to study the application of microwave sounding data in numerical prediction of typhoons and to improve typhoon forecasting,we assimilated data directly for numerical forecasting of the track and intensity of the 2009 typhoon Morakot(0908)based on the WRF-3DVar system.Results showed that the initial fields of the numerical model due to direct assimilation of FY-3A microwave sounding data was improved much more than that due to assimilation of conventional observations alone,and the improvement was especially significant over the ocean,which is always without conventional observations.The model initial fields were more reasonable in reflecting the initial situation of typhoon circulation as well as temperature and humidity conditions,and typhoon central position at sea was also adjusted.Through direct 3DVar assimilation of FY-3A microwave data,the regional mesoscale model improves the forecasting of typhoon track.Therefore,the FY-3A microwave data could efficiently improve the numerical prediction of typhoons.     

Radar data assimilation of the GRAPES model and experimental results in a typhoon case

Constructing β-mesoscale weather systems in initial fields remains a challenging problem in a mesoscale numerical weather prediction (NWP) model. Without vertical velocity matching the β-mesoscale weather system, convection activities would be suppressed by downdraft and cooling caused by precipitating hydrometeors. In this study, a method, basing on the three-dimensional variational (3DVAR) assimilation technique, was developed to obtain reasonable structures of β-mesoscale weather systems by assimilating radar data in a next-generation NWP system named GRAPES (the Global and Regional Assimilation and Prediction System) of China. Single-point testing indicated that assimilating radial wind significantly improved the horizontal wind but had little effect on the vertical velocity, while assimilating the retrieved vertical velocity (taking Richardson's equation as the observational operator) can greatly improve the vertical motion. Experiments on a typhoon show that assimilation of the radial wind data can greatly improve the prediction of the typhoon track, and can ameliorate precipitation to some extent. Assimilating the retrieved vertical velocity and rainwater mixing ratio, and adjusting water vapor and cloud water mixing ratio in the initial fields simultaneously, can significantly improve the tropical cyclone rainfall forecast but has little effect on typhoon path. Joint assimilating these three kinds of radar data gets the best results. Taking into account the scale of different weather systems and representation of observational data, data quality control, error setting of background field and observation data are still requiring further in-depth study.     

Impacts of XBT, TAO, Altimetry and ARGO Observations on the Tropical Pacific Ocean Data Assimilation

This study aims at assessing the relative impacts of four major components of the tropical Pacific Ocean observing system on assimilation of temperature and salinity fields. Observations were collected over a period between January 2001 through June 2003 including temperature data from the expendable bathythermographs （XBT）, thermistor data from the Tropical Ocean Global Atmosphere Tropical Atmosphere-Ocean （TOGA-TAO） mooring array, sea level anomalies from the Topex/Poseidon and Jason-1 altimetry （T/P-J）, and temperature and salinity profiles from the Array for Real-time Geostrophic Oceanography （ARGO） floats. An efficient three-dimensional variational analysis-based method was introduced to assimilate the above data into the tropical-Pacific circulation model. To evaluate the impact of the individual component of the observing system, four observation system experiments were carried out. The experiment that assimilated all four components of the observing system was taken as the reference. The other three experiments were implemented by withholding one of the four components. Results show that the spatial distribution of the data influences its relative contribution. XBT observations produce the most distinguished effects on temperature analyses in the off-equatorial region due to the large amount of measurements and high quality. Similarly, the impact of TAO is dominant in the equatorial region due to the focus of the spatial distribution. The Topex/Poseidon-Jason-1 can be highly complementary where the XBT and TAO observations are sparse. The contribution of XBT or TAO on the assimilated salinity is made by the model dynamics because no salinity observations from them are assimilated. Therefore, T/P-J, as a main source for providing salinity data, has been shown to have greater impacts than either XBT or TAO on the salinity analysis. Although ARGO includes the subsurface observations, the relatively smaller number of observation makes it have the smallest contribution to the assimilation syst     

Three-Dimensional Wind Field Retrieved from Dual-Doppler Radar Based on a Variational Method: Refinement of Vertical Velocity Estimates

In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.     

How Are El Nio and La Nia Events Improved in an Eddy-Resolving Ocean General Circulation Model?

The present study compares the performance of two versions of the LASG/IAP(State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics) Climate System Ocean Model(LICOM) in reproducing the interannual variability associated with El Nio and La Nia events in the tropical Pacific. Both versions are forced with the identical boundary conditions from observed or reanalysis data, in which one version has a finer spatial resolution of(1/10)° in the horizontal domain and 55 vertical layers, and the other version has a coarse resolution of 1° in the horizontal domain and 30 vertical layers. ENSO simulations form the two versions are compared with observations and, in particular, the improvements with regard to ENSO by the finer resolution ocean model are emphasized. As a result of the finer spatial resolution, both the vertical temperature gradient and vertical velocity are better represented in the equatorial Pacific than they are by the coarse resolution model; and thus, the corresponding vertical advections of temperature are more reasonable. Besides the mean climatology, simulated ENSO events and relevant feedbacks are much improved in the finer resolution model. A heat budget analysis suggests that both thermocline feedback and Ekman feedback are mainly responsible for the rapid increase in temperature anomalies during the developing and mature phases of ENSO events.     

Dynamical and Microphysical Retrieval from Simulated Doppler Radar Observations Using the 4DVAR Assimilation Technique

Based on a cloud model and the four-dimensional variational (4DVAR) data assimilation method developed by Sun and Crook (1997), simulated experiments of dynamical and microphysical retrieval from Doppler radar data were performed. The 4DVAR data assimilation technique was applied to a cloud scale model with a warm rain parameterization scheme. The 3D wind, thermodynamical, and microphysical fields were determined by minimizing a cost function, defined by the difference between both radar observed radial velocities and reflectivities and their model predictions. The adjoint of the numerical model was used to provide the gradient of the cost function with respect to the control variables. Experiments have demonstrated that the 4DVAR assimilation method is able to retrieve the detailed structure of wind, thermodynamics, and microphysics by using either dual-Doppler or single-Doppler information. The quality of retrieval depends strongly on the magnitude of constraint with respect to the variables. Retrieving the temperature field, cloud water and water vapor is more difficult than the recovery of the wind field and rainwater. Accurate thermodynamic retrieval requires a longer assimilation period. The inclusion of a background term, even mean fields from a single sounding, helped reduce the retrieval errors. Less accurate velocity fields were obtained when single-Doppler data were used. It was found that the retrieved velocity is sensitive to the location of the retrieval domain relative to the radars while the other fields have very little changes. Two radar volumetric scans are generally adequate for providing the evolution, although the use of additional volumes improves the retrieval. As the amount of the observations decreases, the performance of the retrieval is degraded. However, the missing observations can be compensated by adding a background term to the cost function. The technique is robust to random errors in radial velocity and calibration errors in reflectivity. The boundary conditions from the dual-Doppler synthesized winds are sufficient for the retrieval. When the retrieval is mainly controlled by the observations in the regions away from the boundaries, the simple boundary conditions from velocity azimuth display (VAD) analysis are also available. The microphysical retrieval is sensitive to model errors.     

ON TEMPERATURE INITIALIZATION IN PRIMITIVE EQUATION FORECAST MODELS

Four procedures of specifying model initial temperature were described and tested in the present study. It was found that the use of observed temperatures along with a proper vertical interpolation scheme was not only acceptable, but produced less error than the use of temperatures derived from geopotential height through the hydrostatic equation did. Use of the difference form of the hydrostatic equation would produce unacceptable errors in the initial temperatures, unrealistic horizontal and vertical distribution of temperature, and these errors would influence the calculalion of the pressure gradient force, resulting in substantial, artificial disturbances within the model domain.In addition, an approach to check the initial data was described. Taking advantage of the fact that the geostrophic wind in sigma coordinates should be nondivergent, geopotential height and temperature were used to calculate the pressure gradient force terms and an initial divergence of the geostrophic wind. This approach ca     

Impact of Assimilation of Radiosonde and UAV Observations from the Southern Ocean in the Polar WRF Model

Weather forecasting in the Southern Ocean and Antarctica is a challenge above all due to the rarity of observations to be assimilated in numerical weather prediction(NWP)models.As observations are expensive and logistically challenging,it is important to evaluate the benefit that additional observations could bring to NWP.Atmospheric soundings applying unmanned aerial vehicles(UAVs)have a large potential to supplement conventional radiosonde sounding observations.Here,we applied UAV and radiosonde sounding observations from an RV Polarstern cruise in the ice-covered Weddell Sea in austral winter 2013 to evaluate the impact of their assimilation in the Polar version of the Weather Research and Forecasting(Polar WRF)model.Our experiments revealed small to moderate impacts of radiosonde and UAV data assimilation.In any case,the assimilation of sounding data from both radiosondes and UAVs improved the analyses of air temperature,wind speed,and humidity at the observation site for most of the time.Further,the impact on the results of 5-day-long Polar WRF experiments was often felt over distances of at least 300 km from the observation site.All experiments succeeded in capturing the main features of the evolution of near-surface variables,but the effects of data assimilation varied between different cases.Due to the limited vertical extent of the UAV observations,the impact of their assimilation was limited to the lowermost 1?2-km layer,and assimilation of radiosonde data was more beneficial for modeled sea level pressure and near-surface wind speed.     

Simulating the effects of the 1991 Mount Pinatubo volcanic eruption using the ARPEGE atmosphere general circulation model

The climate changes that occured following the volcanic eruption of Mount Pinatubo in the Phillippines on 15 June 1991 have been simulated using the ARPEGE atmosphere general circulation model （AGCM）. The model was forced by a reconstructed spatial-time distribution of stratospheric aerosols intended for use in long climate simulations. Four statistical ensembles of the AGCM simulations with and without volcanic aerosols over a period of 5 years following the eruption have been made, and the calculated fields have been compared to available observations. The model is able to reproduce some of the observed features after the eruption, such as the winter warming pattern that was observed over the Northern Hemisphere （NH） during the following winters. This pattern was caused by an enhanced Equator-to-pole temperature gradient in the stratosphere that developed due to aerosol heating of the tropics. This in turn led to a strengthening of the polar vortex, which tends to modulate the planetary wave field in such a way that an anomalously positive Arctic Oscillation pattern is produced in the troposphere and at the surface, favouring warm conditions over the NH. During the summer, the model produced a more uniform cooling over the NH.     

Analyses and forecasts of a tornadic supercell outbreak using a 3DVAR system ensemble

As part of NOAA’s "Warn-On-Forecast" initiative, a convective-scale data assimilation and prediction system was developed using the WRF-ARW model and ARPS 3DVAR data assimilation technique. The system was then evaluated using retrospective short-range ensemble analyses and probabilistic forecasts of the tornadic supercell outbreak event that occurred on 24 May 2011 in Oklahoma, USA. A 36-member multi-physics ensemble system provided the initial and boundary conditions for a 3-km convective-scale ensemble system. Radial velocity and reflectivity observations from four WSR-88 Ds were assimilated into the ensemble using the ARPS 3DVAR technique. Five data assimilation and forecast experiments were conducted to evaluate the sensitivity of the system to data assimilation frequencies, in-cloud temperature adjustment schemes, and fixed- and mixed-microphysics ensembles. The results indicated that the experiment with 5-min assimilation frequency quickly built up the storm and produced a more accurate analysis compared with the 10-min assimilation frequency experiment. The predicted vertical vorticity from the moist-adiabatic in-cloud temperature adjustment scheme was larger in magnitude than that from the latent heat scheme. Cycled data assimilation yielded good forecasts, where the ensemble probability of high vertical vorticity matched reasonably well with the observed tornado damage path. Overall, the results of the study suggest that the 3DVAR analysis and forecast system can provide reasonable forecasts of tornadic supercell storms.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅱ： Numerical Simulation Experiments

Microphysics elements and vertical velocity retrieved were incorporated using the nudging method into the initial data assimilation of GRAPES （Global/Regional Assimilation and Prediction System） model. Simulation experiments indicated that nudging technique was effective in forcing the model forecast gradually consistent to the observations, yielding the thermodynamically and dynamically balanced analysis field. As viewed from the simulation results, water vapor is vital to precipitation, and it is a governing factor for the amount and duration of precipitation. The initial cloud water, rain water, and vertical velocity determine the strength distribution of convection and precipitation at the beginning time of forecast; the horizontal wind field steers the motion of the mesoscale weather system embedded in and impacts the position of precipitation zone to a large extent. The simulation experiments show that the influence of the initial retrieval data on prediction weakens with the increase of forecast time, and within the first hour of forecast, the retrieval data have an important impact on the evolution of the weather system, but its influence becomes trivial after the first three hours. Changing the nudging coefficient and the integral time-spacing of numerical model will bring some influences to the results. Herein only one radar reflectivity was used, the radar observations did not cover the whole model domain, and some empirical parameters were used in the retrieval method, therefore some differences still lie between simulation and observation to a certain extent, and further studies on several aspects are expected.     

HAZE-TO-FOG TRANSFORMATION DURING A LONG LASTING, LOW VISIBILITY EPISODE IN NANJING

Haze-to-fog transformation during a long lasting, low visibility episode was examined using the observations from a comprehensive field campaign conducted in Nanjing, China during 4-9 December 2013. In this episode, haze was transformed into fog and the fog lasted for dozens of hours. The impacts of meteorological factors such as wind, temperature (T) and relative humidity (RH) on haze, transition and fog during this episode were investigated. Results revealed significant differences between haze and fog days, due to their different formation mechanisms. Comparison was made for boundary-layer conditions during hazy days, haze-to-fog days and foggy days. Distributions of wind speed and wind direction as well as synoptic weather conditions around Nanjing had determinative impacts on the occurrences and characteristics of haze and fog. Weakened southerly wind in southern Nanjing resulted in high concentration of pollutants, and haze events occurred frequently during the study period. The wind speed was less than 1 m s-1 in the haze event, which resulted in a stable atmospheric condition and weak dispersion of the pollutants. The height of the temperature inversion was about 400 m during the period. The inversion intensity was weak and the temperature-difference was 4°C km-1 or less in haze, while the inversion was stronger, and temperature-difference was about 6°C km-1, approaching the inversion layer intensity in the fog event. Haze event is strongly influenced by ambient RH. RH values increased, which resulted in haze days evidently increased, suggesting that an increasing fraction of haze events be caused by hygroscopic growth of aerosols, rather than simply by high aerosol loading. When RH was above 90%, haze aerosols started to be transformed from haze to fog. This study calls for more efforts to control emissions to prevent haze events in the region.