Value-added Impact of Geostationary Hyperspectral Infrared Sounders on Local Severe Storm Forecasts—via a Quick Regional OSSE

Zhenglong LI,Jun LI∗,Pei WANG,Agnes LIM,Jinlong LI,Timothy J.SCHMIT,Robert ATLAS,Sid-Ahmed BOUKABARA,and Ross N.HOFFMAN

1Cooperative Institute for Meteorological Satellite Studies,University of Wisconsin-Madison,1225 West Dayton Street,Madison,WI 53706,USA

2Center for Satellite Applications and Research,NOAA,Madison,WI 53706,USA

3Atlantic Oceanographic and Meteorological Laboratory,OAR/NOAA,Miami,FL 33149,USA

4Center for Satellite Applications and Research,NOAA,College Park,MD 20740,USA

ABSTRACT Accurate atmospheric temperature and moisture information with high temporal/spatial resolutions are two of the key parameters needed in regional numerical weather prediction(NWP)models to reliably predict high-impact weather events such as local severe storms(LSSs).High spectral resolution or hyperspectral infrared(HIR)sounders from geostationary orbit(GEO)provide an unprecedented source of near time-continuous,three-dimensional information on the dynamic and thermodynamic atmospheric fields—an important benefit for nowcasting and NWP-based forecasting.In order to demonstrate the value of GEO HIR sounder radiances on LSS forecasts,a quick regional OSSE(Observing System Simulation Experiment)framework has been developed,including high-resolution nature run generation,synthetic observation simulation and validation,and impact study on LSS forecasts.Results show that,on top of the existing LEO(low earth orbit)sounders,a GEO HIR sounder may provide value-added impact[a reduction of 3.56%in normalized root-mean-square difference(RMSD)]on LSS forecasts due to large spatial coverage and high temporal resolution,even though the data are assimilated every 6 h with a thinning of 60 km.Additionally,more frequent assimilations and smaller thinning distances allow more observations to be assimilated,and may further increase the positive impact from a GEO HIR sounder.On the other hand,with denser and more frequent observations assimilated,it becomes more difficult to handle the spatial error correlation in observations and gravity waves due to the limitations of current assimilation and forecast systems(such as a static background error covariance).The peak reduction of 4.6%in normalized RMSD is found when observations are assimilated every 3 h with a thinning distance of 30 km.

Key words:OSSE,hyperspectral,sounder,high-impact weather

1.Introduction

High spectral,or hyperspectral,resolution instruments like the advanced infrared(IR)sounders(Menzel et al.,2018)onboard the low earth orbit(LEO)satellites,such as the Atmospheric InfraRed Sounder(AIRS)(Chahine et al.,2006),Infrared Atmospheric Sounding Interferometer(IASI)(Hilton et al.,2012),and Crosstrack Infrared Sounder(CrIS)(Gambacorta et al.,2014),are able to provide global atmospheric temperature and moisture pro file information with high vertical resolution for numerical weather prediction(NWP)models.According to reports from a number of global NWP centers,hyperspectral IR sounders have the largest positive impact,from a single instrument,to the assimilation and prediction results of current operational assimilation forecast systems(Cardinali,2009;Joo et al.,2013;Mallick et al.,2017).In combination with advanced microwave sounder data,advanced IR sounder data have been successfully assimilated with significant positive impacts at many operational NWP centers.These include the European Centre for Medium-Range Weather Forecasts(McNally et al.,2006;Collard and McNally,2009;McNally et al.,2014),the National Centers for Environmental Prediction(Le Marshall et al.,2005,2006a,2006b),M´et´eo-France(Aulign´e et al.,2003),and the United Kingdom Met Office(Hilton et al.,2009).

Monitoring moisture that has large spatial and temporal variability needs more frequent observations.Such high temporal resolution(1 h or better)atmospheric observation can be provided by geostationary satellites,such as the current geostationary operational environmental satellite(GOES)Sounder(Menzel et al.,1998)and advanced imagers.The imagers include the Advanced Baseline Imager(Schmit et al.,2005,2017)onboard GOES-16 and-S,-T and-U;the Advanced Himawari Imager onboard the Japanese Himawari-8/-9 satellites(Bessho et al.,2016;Ma et al.,2017);and the Advanced Geosynchronous Radiation Imager onboard the Chinese FengYun-4(FY-4)series(Yang et al.,2017).Those sensors usually lack vertical resolution due to broader spectral response functions(SRFs)and limited channels.On the other hand,high spectral resolution sensors have a much greater vertical resolving power for temperature,moisture,winds and trace gases than low spectral resolution sensors.

A hyperspectral IR sounder onboard a satellite at geostationary orbit can provide four-dimensional atmospheric temperature,moisture,and wind pro files.Schmit et al.(2009)outlined the advantages of a geostationary satellite–based advanced IR sounder on nowcasting,forecasting and NWP applications over a polar orbit satellite–based sounder by providing rapid refresh measurements for severe weather events.An advanced IR sounder called InfraRed Sounder has been planned for the EUMETSAT Meteosat Third Generation sounding mission to be launched in 2022.The Geostationary Interferometric InfraRed Sounder onboard the Chinese FY-4 series(Yang et al.,2017)is the first geostationary satellite–based high spectral resolution IR sounder.Launched on 11 December 2016,FY-4A is the first satellite of the FY-4 series,though it is experimental.Although the Hyperspectral Environmental Suite was removed from the GOES-R series,NOAA continues to have valid justification for measurements from an advanced IR sounder in geostationary orbit(Schmit et al.,2009).In order to understand the unique value of a geostationary advanced IR sounder for high-impact weather(HIW),such as local severe storms(LSSs),a quick regional Observing System Simulation Experiment(rOSSEs)framework has been developed.This includes generating the Nature Run(NR),simulating LEO orbits,simulating and validating the synthetic radiances,impact study of data assimilation,and NWP forecasting of HIW events.

An OSSE,a type of sensitivity experiment,can answer questions related to the value of new or proposed observing systems and how these and other data are assimilated(Atlas,1997;Atlas et al.,2015;Hoffman and Atlas,2016).In an OSSE,the NR is very important as a proxy of reality.It should be a free run of a forecast model without observation data assimilation,have realistic phenomenology and variability compared with nature,and should be as independent as possible from the data assimilation system model.With NR data,the new and existing observing systems data can be simulated for assimilation.For a more comprehensive description and explanation of OSSEs,see Hoffman and Atlas(2016).OSSEs are classified into two categories depending on the NR and the forecast model.A global OSSE uses a global NR and global NWP forecast model with a focus on large-scale weather systems.A quick rOSSE uses a high-resolution regional NR(or a global NR if the resolution requirement is met or embedded in the global NR,if possible)and a regional NWP model with a focus on mesoscale or small-scale HIW events,such as LSSs and tropical cyclones.

A global OSSE with a geostationary advanced IR sounder has been developed by NOAA’s Joint Center for Satellite Data Assimilation(Boukabara et al.,2016a,2018),with a focus on global-scale weather systems.The quick rOSSE method has been widely used to understand the value-added impact of future observations,such as the NASA Cyclone Global Navigation Satellite System(Zhang et al.,2017),and geostationary hyperspectral IR sounder(Wang et al.,2013;Jones et al.,2017).In this study,a quick rOSSE is carried out to study the impact of a geostationary advanced IR sounder on HIW forecasts compared with existing LEO sounders,including advanced IR sounders such as AIRS,IASI and CrIS.

In this work,the orbit simulator and radiance simulation/validation are introduced in section 2.The assimilation experiments and impact study are outlined in section 3,followed by a discussion in section 4.Finally,a summary and outlook are given in section 5.

2.Simulation and validation

In rOSSE-based studies,the first step is the simulation,including the NR,the orbits,and the synthetic observations.This section provides details on how the simulations are carried out and how they are validated.

2.1.NR

In an OSSE study,the NR represents the real atmospheric state,which requires coverage over a large domain and high spatial and temporal resolutions to provide optimal atmospheric information for both the simulation and assimilation.For an LSS OSSE,a very high resolution NR is needed to capture the mesoscale structures in the storms.For this study,the Weather Research and Forecasting-Nonhydrostatic Mesoscale Model(WRF-NMM),version 3.6.1,is used as the regional NWP model to generate the high-resolution NR(HRNR),with a spatial resolution of 2 km.The domain covers almost the whole of the continental United States and part of the eastern Pacific Ocean,with grid points of 3200(lon)×1320(lat)(see Fig.1 for the NR domain coverage,along with the imagery of storm relative helicity at 0600 UTC 27 May 2006).The model is set up with 51 vertical layers from the surface to a model top at 10 hPa.The schemes used include the Eta(Ferrier)microphysics scheme,the GFDL longwave and shortwave schemes,the Eta similarity surface layer scheme,the Noah land surface model,and the MYJ planetary boundary layer scheme.No cumulus parameterization is needed due to the high resolution(2 km)of the HRNR(Skamarock and Weisman,2009).

Fig.1.Model domain for the NR(outer black lines)and experiments(inner black lines).The values shown are storm relative helicity at 0600 UTC 27 May 2006.Note the pattern of storm relative helicity is similar to the accumulated precipitation shown in the bottom panel in Fig.2.

To avoid the impact of the real atmosphere,the initial and boundary data are from the global OSSE six-hourly Global Forecast System(GFS)analysis(Lim et al.,2017).The global OSSE uses the Goddard Earth Observing System Model Version 5 NR(G5NR,Gelaro et al.,2015).Most existing observation systems are assimilated in the global OSSE,such as the conventional radiosonde observations(RAOBs),surface observations,aircraft,pro filers,radiances from AMSU-A onboard the MetOp and NOAA satellites,Advanced Technology Microwave Sounder(ATMS),AIRS,IASI,CrIS,and the GPS radio occultation(GPSRO)bending angle.The global OSSE output files have a 1°×1°resolution,with 27 vertical levels from the surface to 10 hPa.The original resolutions are downscaled for the HRNR by the WRF model.

The HRNR initializes at 1800 UTC 25 May 2006,and forecasts for 54 h,ending at 0000 UTC 28 May 2006.The storm runs from 0000–1200 UTC 27 May.The HRNR is 30 h ahead of the storm event,to fully prepare the data simulation and assimilation experiments.Figure 2 shows the 3-h precipitation from the HRNR at 0600 UTC 27 May,along with that from the G5NR.The storm precipitation is mainly in southeastern Nebraska,western Iowa,central Kansas and southward to southern Oklahoma and Texas.The HRNR clearly catches the storm precipitation pattern and location as reasonably well as G5NR,but the HRNR has much finer spatial resolution and shows much better spatial structures of the LSS than the G5NR.Due to the high spatial resolution,the precipitation intensity of the HRNR is also stronger than that of the G5NR.Note that the pattern of the accumulated precipitation shown in the upper panel of Fig.2 is similar to the storm relative helicity shown in Fig.1.

2.2.LEO orbit simulator

A satellite orbit simulator was developed at the University of Wisconsin-Madison to simulate LEO satellite orbital information based on track information(Frederick W.Nagle,personal communication).For given starting and ending times,the simulator calculates the orbital information,such as geolocation,satellite geometry(both zenith and azimuth angle),solar geometry,time,etc.For CrIS and IASI,the field of view(FOV)numbers are also provided.Such information is needed for simulating LEO radiance observations.The simulator has excellent agreement with real observations,as shown in Fig.3 for CrIS and IASI/MetOp-A for half of the scan line of a randomly selected granule.Besides,the satellite orbit simulator has been expanded to include many other instruments,such as AIRS,IASI/MetOp-B,AMSU-A from the NOAA satellites,Aqua and MetOp-A/B,ATMS,etc.The orbit simulator can also simulate orbits for future sensors as well,like the microwave sensors for TROPICS(Time-Resolved Observations of Precipitation structure and storm Intensity with a Constellation of Smallsats)mission(https://tropics.ll.mit.edu/CMS/tropics/Mission-Overview).

2.3. Radiative transfer model

Fig.2.Three-hour precipitation from HRNR(upper)and G5NR(lower)from 0300–0600 UTC 27 May 2006.Note the pattern of 3-h accumulated precipitation is similar to the storm relative helicity at 0600 UTC shown in Fig.1.

Fig.3.Validation of simulated(a)CrIS and(b)IASI/MetOp-A orbits using observations.The red dots represent the simulated locations of CrIS and IASI/MetOp-A FOVs,and blue dots represent observations.The CrIS time is 2354:10 UTC 23 September 2014,and the IASI/MetOp-A time is 0023:58 UTC 24 September 2014.Only half of the scan lines are shown due to symmetry.

An efficient hyperspectral IR radiative transfer model(HIRTM)has been developed to simulate synthetic radiance observations in both clear and cloudy skies(Li et al.,2017)for AIRS,CrIS and IASI.The HIRTM uses the Stand-alone AIRS Radiative Transfer Algorithm(SARTA;Hannon et al.,1996;Strow et al.,2003)for atmospheric transmittance calculation.The cloud absorption and scattering are calculated using the cloudy model developed by Wei et al.(2004).Figure 4 shows a clear-sky RTM comparison between SARTA and the Community RTM(CRTM)for CrIS.The University of Wisconsin-Madison(UW-Madison)SeeBor database[15 704 global pro files(Seemann et al.,2003,2008)]are used as inputs in the radiance calculations.The statistical results show that the CRTM[version 2.1.3,with optical depth in pressure space coefficients(Chen et al.,2010,2012)]agrees with SARTA reasonably well,especially for the window re-gion,without any obvious angle dependency.The bias shows some spectral variations,especially in the longwave CO2region(＞ 13 µm).However,most of the bias is within ±0.3 K.The water vapor channels have a slightly larger bias,but in general are within±0.5 K.The standard deviation(STD)shows very good agreement between SARTA and CRTM.There is no strong spectral dependency shown,except on some of the absorption lines(water vapor,ozone and CO2).These small but substantial differences between SARTA and CRTM are important for OSSE studies.In an OSSE,it is generally preferred that the same RTM is not used in the radiance simulation and assimilation—the equivalent of assuming a perfect RTM,which is not true in reality.The HIRTM requires four parameters of clouds in order to perform the cloudy radiative transfer calculation:the effective cloud particle diameter;cloud optical thickness;effective cloud-top pressure;and cloud phase(Wei et al.,2004).All four parameters can be calculated from the hydrometer pro files in the NR(Li et al.,2017).In addition to using the HIRTM for AIRS,CrIS and two IASIs,CRTM V2.1.3,with optical depth in absorber space(ODAS)coefficients,are used to simulate the radiances for ATMS and six AMSU-As from MetOp-A/B and the NOAA satellites.In total,11 existing LEO instruments are simulated to represent the existing LEO capability(Table 1).

Fig.4.Clear-sky RTM comparison(bias in upper panel;STD in lower panel)between SARTA and CRTM with different satellite viewing angle zones(1.9°–17.1°,20.9°–36.5° and 40.6°–57.8°)for CrIS.The UW SeeBor database containing 15 704 global pro files are used in the simulation.

Table 1.Simulated LEO satellite radiances assimilated in the experiments.

2.4. Synthetic observation simulation

Key parameters from the NR used for the synthetic radiance simulation are listed in Table 2.Additional ancillary data include the climatological ozone pro file and the UW High Spectral Resolution IR emissivity database.All pro files are linearly interpolated from 51 WRF model levels to 101 RTM levels.Levels above the model top are filled with climatological temperature and water vapor pro files.An overall shift of the climatological pro files is performed to ensure smooth change around the model top.Use of climatological pro files in those upper levels has very minor impacts on the radiative transfer calculations.

Table 2.Key parameters from the WRF-NMM forecast used in the radiative transfer calculation.

Once the NR and orbits are ready,the synthetic radiance can be simulated via the following steps.For a given geolocation,the FOV size is calculated based on the viewing angle and the provided nadir FOV size (the spatial resolution).With the FOV size known,all NR pixels falling into that FOV will be used to calculate the mean thermodynamic and hydrometer pro files as well as surface parameters,which are used as inputs for the RTM.Random noise is added based on the instrument specification of the Noise Equivalent delta Radiance for IR or Noise Equivalent delta Temperature for microwave.All LEO radiances are converted to Binary Universal Form for the Representation of meteorological data(BUFR)format.

Simulating a future GEO hyperspectral IR sounder follows a more straightforward two-step process.In the first step,the radiance simulation is performed for every grid point in the NR.The satellite and solar geometries for each grid point can be calculated with the known satellite longitude and observation time.In the second step,a linear spatial interpolation is used to convert the simulated radiances to satellite observation locations.This step actually decreases the spatial resolution from the NR to the satellite sensor.All hyperspectral IR sounder radiances,from both GEO and LEO,are simulated using the HIRTM.All others are simulated using CRTM ODAS coefficients.

For synthetic Global Telecommunication System(GTS)observations,only the RAOB pro files of temperature,moisture,and U/V(zonal/meridional)winds are simulated using the nearest pro file from the NR at real RAOB locations extracted from real GTS data with random Gaussian noise added.The T/Q(temperature/moisture)random noise is generatedusing eigenvectors derived from the UW-Madison See-Bor database(Seemann et al.,2008).A similar technique has been used before to generate simulated GFS forecast pro files from RAOB pro files(Jin et al.,2008;Li et al.,2009).The overall STDs of error for temperature and moisture pro files are 1K and10%in relative humidity—the same as those used in the Community Gridpoint Statistical Interpolation(GSI)for Prep BUFR data(Keyser,2018).The dry bias is not considered in the moisture simulation.Random noise for wind is added in a similar way,but with eigenvectors derived from the GFS forecast.The advantage of eigenvectors instead of random noise at each level is that the former includes vertical correlation,which is more realistic than the latter.The synthetic RAOB data are generated every 6 h,with the majority of the data from 0000 and 1200 UTC.

2.5.Verification of simulated radiances

Validating the simulated synthetic radiances follows two steps:validating the RTM and validating the simulated radiances.The first step ensures the radiance simulation is carried out properly by the RTM,in both clear and cloudy skies.Usually,the validation is performed by comparing with real observations or simulations from other RTMs,like the lineby-line RTM.Both the CRTM(Chen et al.,2010,2012)and the HIRTM(Li et al.,2017)have been validated before.The second step ensures the simulated radiances are reasonable.This means that both the NR and the RTM have to be reasonable.

In this study,the future GEO hyperspectral IR sounder is assumed to be IASI at the GOES-16 on-orbit check-out location at 89.5°W,with a spatial resolution of 4 km and temporal resolution of 1 h.The validation of the simulated synthetic radiances is difficult because the NR has nothing to do with reality;there are no real observations that are close to the HRNR that can be used for direct validation.Besides,there are no geostationary satellite radiance observations from 89.5°W at the HRNR time.However,one can still use real observations to indirectly verify the simulated GEO IASI radiances.

To demonstrate that,another set of GEO IASI radiances are simulated by placing IASI at the location of GOES-12(75°W).The synthetic GEO IASI radiances,which are simulated from the HRNR using the HIRTM,are further convolved with the GOES-12 Imager SRF to generate the simulated GOES-12 Imager radiances.Instead of a direct comparison between the simulation and the observation(the two are totally different),the focus is on whether the temporal variation of the simulated radiances is similar to reality.If the temporal variation of the HRNR is not realistic,or if the simulated radiances are not realistic,it is unlikely the temporal variation of the simulated radiances will be in a pattern similar to the real observations.Figure 5 shows that the temporal variations of the simulations for GOES-12 Imager“water vapor”band 3(6.48 µm)and the“longwave window”band 4(10.71µm)are highly similar to those of the observations.The temporal variation is de fined as Tbt2−Tbt1,where t2−t1 is 60 minutes for the GOES-12 Imager.Comparisons of band 3(6.48µm)and band 6(13.31µm),and band 4(10.71µm)and band 6(13.31 µm),also show similar agreement.The similarity of the temporal variations between the simulation and the observation indicates that:(1)the two simulated GOES-12 Imager channels are consistent with each other(both are reasonably simulated);(2)the simulated radiancescapture the real temporal variation reasonably well;and(3)the temporal variation of the thermodynamic information from the HRNR is realistic.With this agreement between the simulation and observation for the GOES-12 Imager,it is reasonable to believe the GEO IASI synthetic radiances from 89.5°W are reasonably simulated.

Fig.5.Temporal variation of GOES-12 Imager bands 3 and 4 from the(a)simulation and(b)observations.The similarity between the simulation and observation indicates channels 3 and 4 are reasonably simulated from the NR.

2.6.Synthetic sounding retrievals

Similar to Jones et al.(2017),synthetic atmospheric sounding retrievals are generated using a simple linear regression technique in clear sky.While assimilating radiances and assimilating sounding retrievals for the hyperspectral IR sounders are theoretically equivalent(Migliorini,2012),the latter does have some advantages at the current stage.For example,in the retrieval,more channels could be used,while the radiance assimilation is limited by channel correlation and computational resources.In particular,for future GEO hyperspectral IR sounders, which have finer spatial resolution and higher temporal resolution,there will be significantly more observations than for existing LEO sounders.The retrieval greatly reduces the computational burden,whilst at the same time retaining most of the information.The retrieved atmospheric temperature and moisture soundings,together with the simulated RAOBs,are converted to PrepBUFR.

3.Assimilation experiments and impact study

For assessing the impact on LSS forecasts,all simulated synthetic observations are converted to either BUFR or Prep-BUFR.The BUFR data include simulated LEO radiances from IASI on MetOp-A/B,CrIS,AIRS,ATMS,and AMSUA on the MetOp-A/B and NOAA satellites(Table 1).The PrepBUFR data include GEO IASI sounding retrievals and RAOBs.While a calibration is preferred for a rigorous OSSE,it is noted that the rOSSE system used in this study is not calibrated.The calibration requires carefully tuning the simulated observation errors so that they have similar impact as in the real world.The short lifetime of the HRNR(less than three days)makes it difficult to perform a statistically sound calibration.For that reason,many rOSSE studies do not perform a calibration(Cintineo et al.,2016;Jones et al.,2017;Zhang et al.,2017),especially when the two respective models for the NR and forecast impact are different.

3.1.Experimental design

The GSI three-dimensional variational system,V3.3,is used as the data assimilation model,which is supported by the Development Testbed Center for community research and study based on the NOAA operational data assimilation system.It can assimilate various kinds of observations,including conventional data,satellite radiances,radar reflectivity,etc.The GEO IASI soundings are assimilated as radiosonde pro files in this study in the GSI model.Considering this study focuses on severe weather forecasting in a regional model,the background error covariance from the North American Mesoscale Forecast System(NAM)is adopted in the experiments.The enhanced bias correction scheme(Zhu et al.,2014)is applied for the assimilation of satellite radiances,and the bias coefficients are updated for every cycling run.The same bias correction scheme is used for all experiments.

To avoid the identical twin problem in OSSEs,the Advanced Research version of the WRF model(WRF-ARW),V3.6.1,is chosen as the regional NWP forecast model for the LSS simulation.The global OSSE analysis data are used to initialize the forecast domain and provide the boundary conditions for the rOSSE.No synthetic observations from the global OSSE are used in this study.Compared with WRF-NMM,used for NR generation,the WRF-ARW has a coarser horizontal spatial resolution(9 km),smaller coverage(450×280 grid points)(Fig.1),and different schemes are used,including the New Thompson scheme for microphysics,Rapid Radiative Transfer Model-Global(RRTMG)for the longwave radiation scheme,RRTMG for the shortwave radiation scheme,the Yonsei University scheme for the planetary boundary layer,and the Kain–Fritsch scheme for cumulus parameterization.

Figure 6 is a flow chart of the experiments for LSS simulation.The model starts at 0600 UTC 26 May,with a 6-h spin-up,and the cycled data assimilation is performed from 1200 UTC 26 May to 0000 UTC 27 May,followed by an 18-h forecast to 1800 UTC 27 May.The 6-h cycling is chosen to be consistent with the NAM and GFS.For the control(CNTRL)experiments,the data assimilated include simulated radiosonde observations(temperature,moisture,and U/V winds)and simulated satellite radiance observations from the sensors listed in Table 1.This represents the existing capability.The experiment with GEO IASI soundings(GEO,for short),as additional observations,represents the capability from one future geostationary hyperspectral IR sounder.The temperature and moisture soundings from GEO IASI are assimilated only from the surface to 200 hPa,with a relaxed observational error covariance[four times the original rootmean-square difference(RMSD)calculated from comparison with the NR].Due to the high density of the GEO IASI data(4-km spatial resolution),the pro files are thinned every 60 km at first.Impacts from different combinations of the assimilation cycle frequency and thinning distance are further discussed in a later section.

3.2.Impact on LSS forecasts

The focus of this study is to show the impact of a GEO hyperspectral IR sounder over existing LEO sounders for LSS forecasts.Compared with LEO,GEO increases the spatial coverage for the regional model.Despite the global coverage from LEO,GEO is capable of providing more frequent full-disk coverage for an LSS study.Figure 7 shows the data coverage of GEO IASI temperature at 500 hPa at 1200 UTC 26 May(with a 60-km thinning).The GEO IASI data are available for the whole model domain,and only data under clear skies are assimilated in the GSI model.The assimilated GEO IASI soundings directly affect thermodynamic and dynamic fields at the analysis time,although the percentage of data assimilated is only 0.125%(see Table 3).The temperature field and the differences between the HRNR and experiments at 1200 UTC 26 May are shown in Fig.8.Compared with the HRNR,the temperature fields of the CNTRL are warmer in the western part of the model domain.The temperature RMSD is 0.55 K,and the STD is 0.54 K.With the assimilation of the GEO IASI soundings,the warm bias at thewestern end of the model domain is smaller,and the temperature RMSD is reduced to 0.43 K,and the STD to 0.39 K.The RMSD is reduced by 22.16%,and the STD reduced by 27.72%.A similar reduction in RMSD can be seen at the two other assimilation times.

Table 3.Reduction of the final normalized RMSD in percentage terms(%, first number)by different experiments(thinning distance and refresh rate)compared to the CNTRL run,and the percentage(%,second number)of the assimilated observations compared with all available clear-sky GEO IASI soundings within the model domain for all cycle steps.

Fig.6.Flow chart of the experimental design.

Fig.7.Data coverage of the retrieval temperature(units:K)at 500 hPa at 1200 UTC 26 May.

Fig.8.Temperature at 500 hPa of(a)NR,(b)CNTRL and(c)GEO IASI at 1200 UTC 26 May,and the differences between the NR and experiments for temperature at 500 hPa:(d)NR−CNTRL;(e)NR−GEO;(f)GEO−CNTRL(units:K).

To further assess the impact of GEO IASI soundings on the analyzed temperature field,the temperature RMSD profile between the HRNR and experiments is calculated from 200 hPa to 1000 hPa at 1200 and 1800 UTC 26 May,and 0000 UTC 27 May(Fig.9).The blue dashed line represents CNTRL and the red line GEO IASI.The temperature fields in the whole model domain are used for verification of the analysis fields.At all three analysis times,GEO has a smaller RMSD than CNTRL for the whole vertical pro file from 200 to 1000 hPa.The largest RMSD differences are always found at around 200 hPa,which indicates the temperature adjustments with assimilating GEO IR soundings are larger near the top of the troposphere than at lower levels.For the RMSD of GEO IASI,an average reduction of 0.15 K throughout the pro files is found over CNTRL.Based on the comparison with the HRNR,assimilating GEO IASI soundings brings the temperature fields in the troposphere closer to the HRNR,and reduces the difference between CNTRL and the HRNR.

Instead of showing further results from individual parameters,an overall evaluation strategy is carried out.The purpose is to use one single parameter to characterize the overall impact for the LSS simulations in the whole domain,which includes most of the important analysis/forecast parameters.For that purpose,the following parameters are selected:the temperature,relative humidity,and U/V winds at four standard atmospheric levels(250 hPa,500 hPa,700 hPa and 850 hPa)to represent the atmospheric thermodynamic fields of the whole domain.In addition,convective available potential energy(CAPE),convective inhibition(CIN),and helicity are used to reflect the environmental instabilities,which are often used by forecasters to make predictions,as well as precipitation,which is an important forecast parameter for LSSs.

The RMSDs are calculated for each parameter for the whole model domain using the HRNR as a reference,except for precipitation;the results are shown in Table 4.Due to the high temporal and spatial variability of precipitation,the equitable threat score(ETS)is often used to characterize both the precipitation amount and pattern(Hamill,1999).To be consistent with the RMSD for which lower values indicate better results, 1-ETS, which is the difference between one and ETS scores of 0.1 mm as the threshold,is used.Notice the temperature,relative humidity and winds are the mean value of the four standard levels.In Table 4,the smaller values between the two experiments are in bold,and the percentage changes are also listed in the last column.For all parameters except CAPE,GEO IASI shows better results than CNTRL.For example,for temperature,assimilating GEO IASI soundings reduces the RMSD by 7.15%,relative humidity by 4.86%,U-wind by 5.9%,V-wind by 7%,CIN by 1.59%,etc.However,it is still impossible to combine different parameters together due to the different units for the parameters.Therefore,a normalization process is carried out to ensure the sum of the square equals 1 for each parameter,including 1-ETS.And the final nominalized RMSD is calculated using a weighted average:

•Thermodynamic parameters(T/Q/U/V),50%

•CAPE,10%

•CIN,10%

•Helicity,10%

•Precipitation 1-ETS,20%

The selected instabilities(10%)and precipitation 1-ETS(20%)are given relatively large weights because of their importance in LSS simulation,especially the precipitation.T/Q/U/V are averaged together with a 50%weight.The calculation is performed at the last analysis time(0000 UTC 27 May)and every 6-h forecast(0600,1200,and 1800 UTC 27 May).The precipitation ETSs are calculated for every 6-h accumulated rainfall at 0600,1200 and 1800 UTC 27 May.

Fig.9.Temperature RMSD between NR and CNTRL(blue),and between NR and GEO(red)from 200 hPa to 1000 hPa at(a)1200 UTC 26 May 26,(b)1800 UTC 26 May,and(c)0000 UTC 27 May.

Table 4.RMSD of the parameters between the NR and experiments.Bold values are the smaller ones between the two experiments.Bold values for percentage change indicate improvement by GEO.

Using the above method,the final normalized RMSD of CNTRL and GEO IASI are calculated.Assimilating GEO IASI soundings is effective in reducing the normalized RMSD from 0.72 to 0.69,which reduces the total normalized RMSD by around 3.56%.The reduction is mainly from the thermodynamic information,which indicates assimilating GEO IASI soundings may benefit the temperature,relative humidity and winds substantially.Overall,the GEO IASI soundings are adding value to the LSS simulation for both the analysis and forecast compared with CNTRL.

4.Analysis and discussion

Despite a spatial resolution of 4 km and temporal resolution of 1 h,the GEO IASI soundings are assimilated with a thinning of 60 km and a 6-h cycling.This means that about 0.1%of the clear-sky observations are assimilated;the majority of the observations are not used.This is not the best way to take advantage of the high temporal and spatial resolutions from GEO.Thinning and relatively large cycling hours can help reduce the impacts from the error correlation(Liu and Rabier,2002,2003)in observations,the gravity wave,and initial imbalances in the model(Lynch and Huang,1992;Benjamin et al.,2004).However,it is important to explore how current assimilation systems and forecast models can take advantage of the high temporal and spatial resolutions from a future GEO hyperspectral IR sounder.

To study the impact of different thinning distances and cycling hours,more experiments are carried out with different thinning distances:15,30,60,and 120 km.For each experiment,the cycling hours of both 3 and 6 h are tested.In total,there are eight experiments in addition to the CNTRL run.Table 3 shows the percentage of observations assimilated,as compared with all available clear-sky GEO IASI soundings within the model domain for all cycled steps.With the same thinning distance,the cycling hour has little impact on the percentage of observations assimilated.However,for a fixed cycling hour,the thinning distance has significant impacts;as the thinning distance halves,the percentage of assimilated observations increases at a rate slightly less than four due to clouds.To quantify the impact of the increased assimilated observations,a final normalized RMSD(not shown)is calculated for each experiment and the CNTRL.The reductions of RMSD in percentage terms for all experiments,compared with the CNTRL,are shown in Table 3.Overall,all experiments show a positive impact from assimilating GEO IASI soundings(the final normalized RMSD is always smaller than the CNTRL),regardless of the thinning distance and cycling hours.When the assimilation cycle is every 6 h,a smaller thinning distance leads to better analysis/forecast results.This is reasonable because there is no error correlation added in the simulated radiance observations.Some error correlations might exist in the GEO IASI sounding retrievals from the retrieval algorithm.However,with a cycle of 6 h,WRF-ARW is able to manage it.When the assimilation cycle is every 3 h,a smaller thinning distance appears to help from 120 km down to 30 km.With a thinning distance of 30 km,the GEO IASI soundings have the largest impact of 4.60%over the CNTRL.This means a reduction in the normalized RMSD of the analysis and forecast by 4.60%from the GEO IASI soundings.Compared with a cycle of 6 h and a thinning of 60 km,there is a slightly more than 1%RMSD reduction.A further smaller thinning distance appears to reduce the positive impact from GEO IASI soundings.A couple of possible reasons may explain this result.One is probably due to the gravity wave and spatial error correlations,which become more difficult for WRF-ARW to handle within 3 h.To further shorten the cycling hour and the thinning distance may require some numerical techniques,such as a digital filter(Huang et al.,2007;Peckham et al.,2016),to remove the exhibited noise from the imbalance in the model forecasts,which is beyond the scope of this study and therefore not discussed any further.Another is the WRF-ARW model resolution(9 km)used in the study.Thinning GEO IASI at 15 km is very close to the model resolution.Combining these dense observations with the high frequency cycle assimilation(here,3 h)could make the model dynamical adjustment difficult from the data in fluence.There should be an optimal balance among the data density,data assimilation cycle and model resolution.In our cases,30-km thinning with 3-h cycling assimilation is the best,and 15-km thinning with 6-h cycling assimilation is the second best.In general,with an appropriate model resolution,the higher spatial resolution GEO IASI soundings and the higher assimilation cycle may increase the forecast skill.

OSSEs can provide quantitative information on observing system impacts for new as well as current observing systems.It can also be used for data assimilation system diagnosis and methodology improvement.Information from OSSEs can lead to better planning and decisions.Quick rOSSEs have the advantage of being cheap and fast to perform,and can sometimes be used to answer questions relating to a particular storm or to demonstrate potential.However,quick rOSSEs are limited in utility since they typically cannot be used to establish statistically significant quantitative results.Both global full OSSEs and quick rOSSEs are helpful for planning and decision making,while the planning and decision making cannot fully rely on OSSE results for three reasons.First,OSSEs have limitations:an NR might not necessarily reflect nature in reality,especially for mesoscale and small-scale weather systems.For example,an NR might have difficulty in simulating the storm intensity due to a limitation of the NWP model used to create the NRs.In addition,the current assimilation system used for impact studies might not be advanced enough to account for the full benefits of future advanced observing systems.Third,OSSEs mainly target NWP applications,while observations from an instrument may have much wider applications.For example,highdensity observations can be used to provide warning information in the pre-convection environment(Ferretti and Faccani,2005,Li et al.,2011),and can also be used for environment and climate monitoring.

5.Summary and future work

A quick rOSSE framework has been developed,including NR generation,orbit simulator,synthetic observation simulation and validation,and impact study.The application to a case study of an LSS using GEO IASI observations has been carried out.GEO IASI has larger spatial coverage and higher temporal resolution,compared with LEO advanced IR sounders,for regional NWP models.The impact study shows that GEO IASI may add substantial positive value to the analysis and forecasting of the thermodynamic parametersof T/Q/U/V.The evaluation based on the combined normalizedRMSD from different parameters indicates an overall positive impact on the LSS analysis and forecast for the case demonstrated.The reduction in the normalized RMSD in the analysis and forecast is 3.56%.In addition,experiments with different thinning distances and assimilation cycling hours reveal that more frequent assimilation and a shorter thinning distance may further increase the positive impact from GEO IASI with the current GSI/WRF-ARW system.However,denser and more frequent assimilated observations also poses a significant burden on the assimilation and forecast systems,especially in handling the spatial error correlation in observations and gravity waves with a static background error covariance matrix,which could hinder the positive impact from GEO IASI.The peak reduction of 4.60%in the normalized RMSD is found with a thinning distance of 30 km and a cycle of 3 h.Any further reduction in thinning distance and cycling hour should also consider the model resolution used.

It should be noted that this is a single case study;the results do not provide statistical significance.Ideally,more cases should be used,but identifying typical LSS cases from a global NR is not easy due to the fact that global NRs are limited in simulating strong LSSs.Nevertheless,the case study demonstrates the potential value-added impact from GEO hyperspectral IR sounder data on LSSs.

To improve background error representativeness,future work should use a hybrid/ensemble data assimilation system with a flow-dependent background error covariance.Improving sounding retrievals using the one-dimensional variational method(Li et al.,2000),expanding soundings into cloudy skies(Li et al.,2009,2016;Wang et al.,2014,2015,2017),and including geostationary imager radiances(Honda et al.,2018a,2018b)for assimilation will also be conducted.Comparisonbetween direct sounding radiance assimilation and retrieval assimilation should also be included,with an appropriate experimental design.Another important future aim is to assimilate atmospheric motion vectors(AMVs)derived from water vapor soundings;the AMVs from continuous water vapor measurements are an important feature of GEO advanced IR sounders,which could further benefit LSS forecasts.

Acknowledgements.This work is partly supported by the NESDIS OPPA OSSE program(Grant No.NA15NES4320001).The views,opinions,and findings contained in this report are those of the authors and should not be construed as an official NOAA or U.S.government position,policy,or decision.Thanks to the Joint Center for Satellite Data Assimilation for providing the“S4”supercomputer[Supercomputer for Satellite Simulations and Data AssimilationBoukabara et al.(2016b)],physically located at SSEC,University of Wisconsin-Madison,as the main computational resource for this research study.