Adding four-dimensional data assimilation by analysis nudging to the Model for Prediction Across Scales – Atmosphere (version 4.0)

The Model for Prediction Across Scales – Atmosphere (MPAS-A) has been modified to allow four-dimensional data assimilation (FDDA) by the nudging of temperature, humidity, and wind toward target values predefined on the MPAS-A computational mesh. The addition of nudging allows MPAS-A to be used as a global-scale meteorological driver for retrospective air quality modeling. The technique of “analysis nudging” developed for the Penn State/National Center for Atmospheric Research (NCAR) Mesoscale Model, and later applied in the Weather Research and Forecasting model, is implemented in MPAS-A with adaptations for its polygonal Voronoi mesh. Reference fields generated from 1°×1° National Centers for Environmental Prediction (NCEP) FNL (Final) Operational Global Analysis data were used to constrain MPAS-A simulations on a 92–25km variable-resolution mesh with refinement centered over the contiguous United States. Test simulations were conducted for January and July 2013 with and without FDDA, and compared to reference fields and near-surface meteorological observations. The results demonstrate that MPAS-A with analysis nudging has high fidelity to the reference data while still maintaining conservation of mass as in the unmodified model. The results also show that application of FDDA constrains model errors relative to 2m temperature, 2m water vapor mixing ratio, and 10m wind speed such that they continue to be at or below the magnitudes found at the start of each test period.


Introduction
Combining data at various times in a dynamical model to provide time continuity and dynamic balance among the model fields was first suggested by Charney et al. (1969) and has become known as four-dimensional data assimilation (FDDA). Nudging, also known as Newtonian relaxation, involves the use of special terms in forecast equations to "nudge" an atmospheric model simulation toward observations or some reference state. It was originally employed for dynamic initialization (Anthes, 1974;Hoke and Anthes, 1976). Many approaches to data assimilation have been developed for dynamical atmospheric models including simple interpolation and more complex variational and stochastic methods. However, these have usually been intended to initialize prognostic models. Nudging was tested as a means for improving diagnostic simulations by Stauffer and Seaman (1990) in the Penn State/National Center for Atmospheric Research (NCAR) Mesoscale Model -version 4 (MM4) (Anthes et al., 1987). In that study, nudging was applied in two ways. The model solutions were nudged toward either gridded analyses ("analysis nudging") or individual observations ("obs nudging"). After MM4, both forms of nudging were implemented in its successor MM5 (Grell et al., 1995) and subsequently in the Weather Research and Forecasting (WRF) model (Skamarock and Klemp, 2008).
In addition to basic variables like temperature, humidity, wind, cloud cover, and precipitation, meteorological simulations guided by nudging provide factors critical to air quality modeling that are not easily observed, such as stability, turbulence, mixing height, etc. Nudging applied in MM4, MM5, and WRF has been used at the U.S. Environmental Protection Agency (U.S. EPA) for almost three decades to support air quality modeling, first with the Regional Acid Deposition Model (RADM) (Chang et al., 1987) and continuing to the present day with the Community Multiscale Air Quality (CMAQ) model (Byun and Schere, 2006;Appel et al., 2017).
It has long been recognized that air quality at any particular location can be affected by pollution sources on local to global scales (NRC, 1998(NRC, , 2010. Air quality models are often applied with relatively coarse horizontal resolution on hemispheric and global scales to provide boundary information for nested, higher-resolution regional models (Bullock Jr. et al., 2008;Jacobson and Ginnebaugh, 2010;Schere et al., 2012;Mathur et al., 2014). In various applications, this nested modeling strategy has created unrealistic simulations at the lateral boundaries of internal model domains due to discontinuities in horizontal and vertical resolution and/or differing modeling assumptions between separate models used at each scale (Warner et al., 1997;Bullock Jr. et al., 2009;Tudor and Termonia, 2010;Mathur et al., 2017).
To address the need for a global-to-local air quality modeling system that can avoid boundary problems associated with model domain nesting, this work adapts the Model for Prediction Across Scales -Atmosphere (MPAS-A) (Skamarock et al., 2012) for use as the meteorological component of a future coupled meteorological-chemical modeling system. MPAS-A, which features a global computational mesh based on a centroidal Voronoi tessellation (Du et al., 1999), offers gradual mesh refinement rather than discrete nesting to a focal region. For retrospective air quality modeling, an FDDA approach based on analysis nudging has been developed and tested in MPAS-A as described below.

Experimental design, implementation, and testing
FDDA by way of analysis nudging, similar to that described in Stauffer and Seaman (1990), has been added to MPAS-A. Unlike MM4, MM5, and WRF, which are limited-area models with rectangular computational grids, MPAS-A has a polygonal computational mesh as illustrated in Fig. 1. Nonetheless, once the required "target" fields (i.e., reference data for nudging) are defined to match the MPAS-A prognostic variable array, analysis nudging in MPAS-A is similar to its ancestral implementations.

Description of analysis nudging in MPAS-A
Analysis nudging is applied in MPAS-A by the addition of a nudging tendency term to the normal predictive equation. The nudging tendency for prognostic variable α is calculated as where G α is a nudging inverse timescale or "nudging coefficient", W PBL and W layer are special binary terms (1 or 0), and α target is the target or reference value for α. Trusted reference fields are typically only available at certain times each day and temporal interpolation is required to provide target values at each model time step. It may be advantageous to avoid nudging in the planetary boundary layer (PBL) so as not to disrupt simulated diurnal processes. In this case, W PBL can be set equal to 1 in layers above the simulated PBL top and set equal to 0 in layers below or containing the PBL top. Otherwise, W PBL can be set equal to 1 in all layers. W layer is a similar binary term to allow the exclusion of nudging near the surface based simply on vertical layer number. G α , W PBL , and W layer are all defined independently for each of the nudged variables.
Analysis nudging has been applied for potential temperature (Θ), water vapor mixing ratio (q v ), and wind. Treating wind involves extra complications because of the way it is represented in the MPAS-A mesh. As illustrated by Fig. 1, scalar prognostic variables including Θ and q v are defined at the cell centers. However, the prognostic variable for wind in MPAS-A is the component perpendicular to the cell faces (U). To nudge wind, meridional and zonal decompositions at the cell centers are used. These model variables URe-constructZonal and UReconstructMeridional already exist to treat the influence of parameterized convection and PBL processes on the wind field. While the wind component across cell edges (U) could be nudged directly, this method would require 50% more comparisons between prognostic and target values since there are 3 times as many cell edges as there are cells. Nudging tendencies for UReconstructZonal and UReconstructMeridional are translated to cell edges in the same manner as the tendencies for PBL and convection processes.

Creating target fields
The MPAS-A modeling system already provides model initialization software, namely the executable program init_atmosphere_model. For this study, initialization fields were created at each time where nudging target fields were desired using 1°×1° National Centers for Environmental Prediction (NCEP) FNL (Final) Operational Global Analysis data (ds083.2) (NCEP/NWS/NOAA/U.S. Department of Commerce, 2000). Target fields could be based on other analytical methods such as three-and four-dimensional variational assimilation (3D-Var, 4D-Var) or an ensemble Kalman filter (EnKF). However, the NCEP FNL data are already produced using techniques similar to 4D-Var. A full description of the NCEP FNL data product is available from the "Documentation" tab at https://doi.org/10.5065/ D6M043C6. The spatial resolution of the 1°×1° NCEP FNL data approximates that of the coarse portion of the MPAS-A mesh used in this study. The nudged variables Θ, q v , UReconstructZonal, and UReconstructMeridional were extracted from each initialization file, renamed th_ fdda_new, qv_fdda_new, u_fdda_new, and v_ fdda_new, respectively, and used to compile the necessary FDDA input files. The modified MPAS-A reads in new FDDA targets every 6h when NCEP FNL data are available, specifically at 00:00, 06:00, 12:00, and 18:00UTC. Target values at intervening times during the simulation are computed using linear time interpolation.
The FDDA target variable names contain "new" to indicate that, for the time increment at which they are read, the values represent the target value at the end of the upcoming 6h FDDA time interval. Unlike WRF which reads in "old" and "new" targets for each FDDA interval, the modified MPAS-A reads only "new" values. The "new" values from the previous time interval are recycled to be used as the "old" values, thus reducing the FDDA target file size by half. At simulation start time, initial values for Θ, q v , UReconstructZonal, and UReconstructMeridional are used to set th_fdda_old, qv_ fdda_ old, u_fdda_old, and v_fdda_old, respectively. Thus, simulation start time must be at 00:00, 06:00, 12:00, or 18:00UTC in order to maintain the 6h FDDA data interval.
Scripts have been written to automate the process of running init_atmosphere_model for each FDDA time, extracting the four nudged variables, and composing the FDDA target input file. They perform variable extraction and FDDA input file composition using NetCDF Operators (NCO) software available at http://nco.sourceforge.net/ (last access: 9 July 2018).
As mentioned before, model initialization and FDDA inputs were produced from 1°×1° NCEP FNL data using the init_atmosphere_model software included in the MPAS-A version 4.0 public distribution. For this study, init_atmosphere_model was slightly modified to allow finer vertical resolution near the surface. Air quality models typically require fine vertical resolution in the PBL in order to better simulate pollutant emissions which are commonly near the surface. To produce sufficiently thin layers near the surface, the unmodified init_atmosphere_model required an unreasonable number of layers due to the 1.5-power function used to define layer boundary heights. Modified Fortran code described in Appendix A was developed such that only 50 layers were required with the model top specified at 30km. Using the modified code and given a surface elevation at sea level, layer thickness is 18m at the bottom, 232m at 1.5km elevation, 1000m at 12.5km elevation, and 1729m at the top. This layer structure was used to produce model initialization and FDDA target data files for all MPAS-A model simulations described below.
The modified init_atmosphere_model was also used to produce update fields for sea surface temperature and sea ice at 6h intervals throughout each test period. For this purpose, the new layer generation function had no bearing, but a problem was discovered in the original MPAS-A model code where sea ice was being analyzed over land areas. This problem was solved with additional code modifications described in Appendix B.
Once the required initialization, surface update, and FDDA target fields were in place, MPAS-A simulations were performed with the atmosphere_model program.  Table 3 shows the new namelist options added as part of the FDDA implementation and the values used for testing in this study. These options are similar to those used in the WRF model for FDDA application. WRF contains them within a special fdda subset of namelist options. For now, they have been added to the physics namelist input variable list for MPAS-A. The primary option to invoke FDDA is config_fdda_scheme. As a default, config_fdda_scheme=off, and FDDA is not invoked in the modified MPAS-A model. If FDDA is invoked with a value of analysis, then the other options in Table 3 become applicable. The modified MPAS-A code also includes a second option for FDDA called scaled which allows the user to adjust nudging strength based on MPAS cell size. This option is still under development and has not been investigated as a part of this study.
FDDA can be selectively applied or omitted for each meteorological variable (Θ, q v , U). If applied, the nudging strength is controlled by a variable-specific nudging coefficient. FDDA for each variable can be applied throughout all vertical layers, or only above a particular layer number. In many previous applications of WRF, it was common for FDDA to only be applied above the PBL so as not to disrupt the diurnal evolution of the PBL with data from a linear time interpolation Bowden et al., 2012Bowden et al., , 2013Bullock Jr. et al., 2014). Table 3 shows namelist options provided to avoid nudging in the PBL, avoid nudging below a specified layer number, or both. It is also worth noting here that the default nudging coefficients implemented in MPAS-A are equal to 3.0×10 −4 s −1 for all variables, just like in WRF. The theoretical reasoning comes from Stauffer and Seaman (1990) where they equate this nudging timescale to that of meteorological phenomena at the meso-α spatial scale. Unlike for temperature and wind, nudging water vapor concentration perturbs atmospheric mass in the simulation. Previous studies using WRF have chosen to employ smaller nudging coefficients for q v versus other variables Bowden et al., 2012Bowden et al., , 2013Bullock Jr. et al., 2014). Results discussed later in this work show some benefit from doing so. For this study, the nudging coefficient for q v was 1 order of magnitude smaller than that for the other variables, except for a special test where the value was kept equal.

Results
To evaluate FDDA in MPAS-A, test simulations for January and July of 2013 were performed with both the standard version of the model and the modified model using FDDA by analysis nudging. The modified MPAS-A was also applied with FDDA turned off to verify agreement with the results obtained from the standard model. Model results from the standard and modified versions were first compared to the FDDA target fields. Obviously, nudging strongly toward the target fields should produce good agreement with those fields. The intent of these first comparisons was to verify that using nudging coefficients for temperature, humidity, and wind similar to those used in WRF would constrain MPAS-A simulations in a reasonable manner. To further test the capabilities of FDDA in MPAS-A, simulated surface-level data for temperature, humidity, and wind speed from both the standard and modified MPAS-A models were then compared to observational data from the Meteorological Assimilation Data Ingest System (MADIS) (https://madis.noaa.gov, last access: 9 July 2018). Finally, total dry air, total water vapor, and total atmospheric mass calculations were performed to test for any corruption of mass conservation by the implementation of FDDA in MPAS-A. Differences are calculated as FDDA target field minus MPAS-A with FDDA to reflect the polarity of the nudging force. The center and right columns show the same information for layers 28 and 45, respectively. Layer 1 extends from 0 to ~18m above the surface where the surface is at mean sea level (m.s.l.), and to ~15m above the surface over the highest resolved terrain. The vertical span of layer 28 varies from 5002-5551 to 9449-9784ma.m.s.l. depending on the resolved terrain height which varies from −82 to 5425m. So layer 28 represents a 330-550m thick layer somewhere in the middle troposphere. The span of layer 45 varies from 20622-22024 to 20682-22048ma.m.s.l. Layer 45 varies only slightly due to the MPAS-A hybrid vertical coordinate system having shifted almost completely to a height coordinate at that altitude in the lower stratosphere.

Comparisons to FDDA target fields
By 10 days into the simulation, the simulation without FDDA already shows significant potential temperature differences from the FDDA target fields in all three layers shown in Fig. 2. These differences are especially noticeable in layer 45 where an apparent stratospheric warming event is stronger in the "no FDDA" simulation and is longitudinally displaced about 120° from the location in the reanalysis-based target field. The unconstrained simulation also shows a high-latitude cold pool that is not in the target field and much colder stratospheric temperatures over the tropics. There are also some interesting differences in layer 28 around the high terrain of the Himalayas where the "no FDDA" simulation resulted in much warmer temperatures than the target values. The simulation with FDDA matches the target fields for Θ almost perfectly for layers 28 and 45. Near the surface, maps for layer 1 show difference from the Θ target field in both MPAS-A simulations, mostly in Arctic regions where simulated temperatures are generally colder. But the simulation with FDDA shows surface temperatures closer to the target values even though FDDA was not applied in the PBL.
To demonstrate that FDDA continues to constrain MPAS-A simulations through longer time periods, Fig. 3 shows the same information as Fig. 2, but this time the Θ fields are for 00:00UTC, 31 January 2013, 30 days into the simulation. The deviation of the "no FDDA" simulation from the target fields is larger than at day 10, but the results from the simulation with FDDA continue to follow the target fields closely for layers 28 and 45. However, layer 1 continues to be too cold across the Arctic with FDDA applied above the PBL. The simulation without FDDA is too cold in some parts of the Arctic and too warm in others, and significant deviations from the target field are apparent in many locations around the globe. A quick investigation of observed surface temperatures at Utqiaġvik (formerly Barrow), Alaska, found that the simulated surface temperatures are about 10K too cold at that location in the FDDA simulation. This points to the fact that FDDA applied only above the PBL can keep simulated temperatures at the surface from being too high due to the effect of convective mixing, but it cannot prevent them from being too cold. Work is ongoing to remedy both warm and cold biases in surface temperature through the use of other land surface and PBL models which nudge soil temperature and moisture towards known conditions.
To investigate temporal variability of the nudging effects on Θ, the simulated and target values for January 2013 were plotted for layer 28 at Research Triangle Park, NC (35.93° N, 78.96° W), along with the value of the nudging term. The results in Fig. 4 show that the simulated values tracked the FDDA target values quite closely but do deviate at times. Simultaneously, the nudging term increases in magnitude to counteract these deviations. A significant perturbation is evident at around 00:00UTC on 18 January. To investigate this occurrence and to illustrate the spatial variability of the nudging, spatial plots of the Θ nudging term were made for layer 28 at 00:00UTC on 18 January, one with a focus on North America and the second showing the entire global domain. Figure 5 shows the spatial extent of the perturbation and the spatial patterns of positive and negative forcing at that time.
These patterns match the spatial scale of meteorological features where model error might grow over time were it not for the corrective effects of analysis nudging.
The array of maps in Fig. 6 shows water vapor mixing ratio at 00:00UTC, 31 January 2013, in the same arrangement as for potential temperature in Fig. 3. Even with weak q v nudging, the simulation with FDDA matches the target fields well for all three layers. Without FDDA, the simulated pattern of water vapor deviates significantly from the target in all layers. In layer 45, the "no FDDA" case shows higher q v values all across the tropics than exist in the target field or in the simulation with FDDA where water vapor is practically absent. These results suggest that even weak nudging of water vapor can mitigate what appears to be artificial vertical diffusion of tropospheric water vapor into the stratosphere. Figure 7 shows layer 1 fields at 00:00UTC, 31 January 2013, for potential temperature and water vapor mixing ratio, but this time focused on the contiguous United States (CONUS). At this point in time, a strong cold front stretching from western Pennsylvania, down the Appalachian Mountains, and into the Gulf of Mexico was advancing from the west. It is important to note once again that analysis nudging was only applied above the PBL. Nonetheless, not only does the simulation with FDDA place the front in the correct location, but the simulated front shows sharper detail for Θ and q v than in the target fields. This is especially true for q v where it appears that weaker nudging above the PBL allows simulated details at the surface to be better conserved. Finer detail is evident, not only along the cold front but also in other locations where it appears variations in terrain and land cover type may be important. It stands to reason that if nudging had been applied within the PBL, these simulated fine details at the surface would have been blurred somewhat by blending with the less resolved target fields. Of course, finer detail does not in itself indicate better accuracy.
To address that issue, simulation results were also compared to observational data as described later in Sect. 3.2.
Vertical cross-section plots of water vapor mixing ratio show that FDDA can constrain undesirable model behavior even when the nudging strength is quite weak. Figure 8 shows vertical cross sections along longitude 80° E from 55° N to 55° S for 00:00UTC, 31 January 2013 (day 30). This particular location was chosen to investigate the effect of the Himalayas on MPAS-A simulation results. The height of the Himalayas (~5km) as resolved by the 92km mesh size is shown by the white area at the bottom of each plot. Figure 8a shows the standard MPAS simulation result, Fig. 8b shows the FDDA target field, and Fig. 8c shows the results from MPAS-A with FDDA using moisture nudging at 1/10 the strength used for the other variables. As previously seen in Fig. 6, the unconstrained simulation shows signs of upward transport and/or diffusion of water vapor into the lower stratosphere with water vapor mixing ratios over 2 orders of magnitude higher than the target values. The simulation with FDDA almost completely eliminates this deviation from the target, even with the weak nudging strength.
Wind velocities and flow patterns from MPAS-A simulations and FDDA target fields were investigated with streamline plots. Figure 9 shows global streamline analyses for layer 28 (~500-300hPa), again for 00:00UTC, 31 January 2013. As might be expected from a 30-day forecast, the flow field from the simulation without FDDA (Fig. 9a) differs significantly from the FDDA target flow field (Fig. 9b). However, the simulation with FDDA (Fig. 9c) follows the FDDA target data almost perfectly. To show the ability of FDDA to maintain finer-scale fidelity, similar streamline analyses focused on the CONUS are shown in Fig. 10. Again, the simulation with FDDA is almost identical to the target flow field. Streamline analysis for layer 1 focused on the southeastern US. Figure 11 shows some noticeable differences between the simulation with FDDA and the target field. These differences are not surprising given that no nudging was applied in the PBL. Also, the FDDA target fields above the PBL were derived from 1° FNL reanalysis data, while the simulation cell size is ~25km in this region. Terrain effects on wind flow direction and speed appear to be more significant in the simulation than in the FDDA target field.
Similar comparisons of MPAS-A simulations to FDDA target data were made for the July 2013 test period. All of these comparisons showed essentially the same results as were found for January 2013. While weather systems and patterns were generally more quiescent, at least in the Northern Hemisphere, simulations unconstrained by FDDA still deviated significantly from target fields after a few days, while those constrained by FDDA maintained their fidelity relative to the target data.

Comparisons to observational data
MPAS-A simulation results were compared to observations of 2m temperature, 2m humidity, and 10m wind speed. To assure data quality, only aviation routine weather reports (METAR) and surface airways observation (SAO) reports from the MADIS repository were used. This comparison was made using the Atmospheric Model Evaluation Tool (AMET) described in Appel et al. (2011). AMET was configured to calculate daily evaluation statistics for the entire global domain and for a subdomain confined within 25 to 50° N latitude and 67 to 125° W longitude, basically covering the CONUS where the horizontal mesh size was 25km. Daily statistics were calculated for both the January and July 2013 test periods. Figure 12 shows the time series of 2m temperature root mean squared error (RMSE) for January 2013. The top graph shows results for the entire global domain, while the bottom graph shows results for the CONUS subdomain. Three MPAS-A simulations were analyzed, the standard model denoted as no FDDA, FDDA applied using relatively weak q v nudging denoted as FDDA, and FDDA applied using equal nudging strength for all variables denoted as FDDA (equal). The no FDDA cases show error increasing right from the start, both globally and over the CONUS. For the global domain, the FDDA and FDDA (equal) cases both show RMSE actually decreasing somewhat over the first 10 days but generally holding steady throughout the month. Globally, model performance with weak q v nudging is slightly better than with equal nudging strength for Θ, q v , and U. In the CONUS analysis, 2m temperature RMSE for the FDDA case decreases more significantly over the first 10 days than in the global analysis and remains below the starting values throughout the remainder of the month. Once again, relatively weak q v nudging improves model error statistics to some degree. Figure 13 shows the same information as Fig. 12 Even though FDDA nudging was not applied within the PBL for any variable, the results above show model errors near the surface were constrained quite well, except where simulated surface temperatures were too cold. Also, near the surface is where finer horizontal resolution of the model relative to the FDDA target data source has its greatest effect.
Further study is anticipated to better identify optimal FDDA nudging strengths for Θ, q v , and U in MPAS-A, and to better understand the vertical levels of the atmosphere where nudging should be applied.

Mass conservation tests
In addition to the comparisons described above, the ability of MPAS-A to conserve simulated atmospheric mass was also tested. For each month-long test period, all model simulations reported total atmospheric mass and total water vapor mass at each 150s simulation time step. This was accomplished with minor additions to the Fortran code in the time integration module (./src/core_atmosphere/dynamics/mpas_atm_time_ integration.F ). The modified module is included in the model code repository (https://doi.org/10.5281/ zenodo.1101204).
As shown in Fig. 18, all simulations including those with FDDA conserved total moist air within the global model domain to within 5 parts in 100000 of their starting values. Total water vapor mass (Fig. 19) varied more significantly in time in each simulation, and this is to be expected due to evaporation and precipitation processes. There is a diurnal signal evident in the water vapor mass total from all simulations, most likely due to longitudinal variations in evaporation and precipitation potential under solar radiation caused by the geographic distribution of continents and oceans. For January 2013, the no FDDA simulation lost over 5% of its initial quantity of water vapor. This could be indicative of too much simulated precipitation, too little simulated evaporation, or both. Also, the model initialization based on the NCEP FNL analysis could be too moist. The simulations with FDDA all tended to maintain more total water vapor relative to the standard model. In general, the unconstrained simulations tended to lose water vapor at the start of the simulation and come to an equilibrium point significantly lower than the simulations with FDDA. The FDDA (equal) simulations tended to quickly establish and then maintain the most total water vapor. Obviously, the FNL analysis indicated a moister atmosphere than the unconstrained MPAS-A simulations could maintain.
It is interesting to note that the FDDA and FDDA (equal) cases had almost identical trends in total moist air mass, with any differences in total water vapor almost perfectly canceled out by differences in dry air (Fig. 20). Also, the January 2013 no FDDA simulation gained total dry air mass to about 13 parts in 100000 of the initial value. This is the same simulation that lost a significant fraction of its initialized water vapor, once again showing an opposite conservation response between dry air and water vapor.
Overall, the results in Figs. 18-20 demonstrate that the addition of FDDA does not degrade mass conservation relative to the standard MPAS-A. Conservation of dry air mass is most important if MPAS-A is to be used as the meteorological driver for air quality modeling. These results show that using FDDA could actually offer some improvement in that regard.

Summary and conclusions
The U.S. EPA is working to make MPAS-A suitable for use as the meteorological component of an integrated meteorology and air quality modeling system for global-to-finescale applications. The ability to constrain simulated meteorology to resemble historical reanalysis fields at comparable spatial scales is crucial to making this integrated modeling system a practical diagnostic tool for air quality research. FDDA applied through analysis nudging has been used for decades to provide this constraint in other models such as MM4, MM5, and WRF. The results shown here demonstrate that it also works quite well in MPAS-A. Comparison of MPAS-A simulations of January and July 2013 with and without FDDA demonstrates that unconstrained simulations deviate significantly from historical conditions in only a few days, while those constrained through analysis nudging follow historical conditions well in most situations. Due to surface-layer decoupling, analysis nudging applied only above the PBL was not able to constrain the development of excessively cold surface temperatures in Arctic areas during the January 2013 simulation period. However, this can be addressed with the use of land surface models that also employ FDDA. Further study is already underway at the U.S. EPA to determine the best strength with which to nudge temperature, humidity, and wind in MPAS-A and the levels of the atmosphere to best apply that nudging.
The target fields toward which MPAS-A state variables are nudged could come from a number of sources. Historical meteorological reanalysis products have previously been used for this purpose in regional and hemispheric modeling with WRF, and the results here suggest they can also be used with MPAS-A on the global scale. This study applied MPAS-A with a variable 92-25km mesh with the refined region centered on North America. Target fields used here were based on the 1°×1° NCEP FNL reanalysis product. As such, there was not a great disparity in horizontal resolution between the simulations and the target fields where the MPAS-A mesh size was 92km. However, where the mesh was more refined, MPAS-A was capable of delivering additional horizontal detail, and the results shown here indicate weaker nudging may produce superior results, at least when nudging water vapor mixing ratio. Finally, adding FDDA did not disrupt the ability of MPAS-A simulations to conserve mass, and this is an important point when considering its use for air quality modeling.
Code and data availability.
The MPAS-A model software used in this project is a subset of the complete Model for Prediction Across Scales (MPAS) developed by Los Alamos National Security, LLC (LANS) and the University Corporation for Atmospheric Research (UCAR), and distributed under a three-clause BSD license allowing distribution of original and derivative works under conditions that have been satisfied here. The full text of this BSD license can be found at http://mpas-dev.github.io/files/documents/MPAS-DevelopersGuide.pdf (last access: 9 July 2018 The init_atmosphere core previously only activated packages to read in static fields when config_init_case==7 (and either config_vertical_grid or config_met_interp were true). However, for config_init_case==8, we need the landmask field for interpolating, e.g., sea ice. Since we may potentially need other static fields as well when creating surface update files, the simplest solution seems to be to simply read all static fields when config_ init_case==8 by activating the vertical_stage_in package for this case.
The basis for these changes is the MPAS-        Simulation results from MPAS-A using FDDA and FDDA target fields for potential temperature (K) and water vapor mixing ratio (kgkg −1 ) in layer 1 at 00:00UTC on 31                config_radtlw_interval Temporal interval between calls to the parameterizations of longwave radiation, format "yyyy-mm-dd_hh:mm:ss".