Search

The ability of two downscaling experiments to correctly simulate thermodynamic conditions over the Iberian Peninsula (IP) is compared in this paper. To do so, three parameters used to evaluate the unstable conditions in the atmosphere are evaluated: the total totals index (TT), convective available potential energy (CAPE), and convective inhibition (CIN). The Weather and Research Forecasting (WRF) model is used for the simulations. The N experiment is driven by ERA-Interim's initial and boundary conditions. The D experiment has the same configuration as N, but the 3DVAR data assimilation step is additionally run at 00:00, 06:00, 12:00, and 18:00 UTC. Eight radiosondes are available over the IP, and the vertical temperature and moisture profiles from the radiosondes provided by the University of Wyoming and the Integrated Global Radiosonde Archive (IGRA) were used to calculate three parameters commonly used to represent atmospheric instability by our own methodology using the R package aiRthermo. According to the validation, the correlation, standard deviation (SD), and root mean squared error (RMSE) obtained by the D experiment for all the variables at most of the stations are better than those for N. The different methods produce small discrepancies between the values for TT, but these are larger for CAPE and CIN due to the dependency of these quantities on the initial conditions assumed for the calculation of a lifted air parcel. Similar results arise from the seasonal analysis concerning both WRF experiments: N tends to overestimate or underestimate (depending on the parameter) the variability of the reference values of the parameters, but D is able to capture it in most of the seasons. In general, D is able to produce more reliable results due to the more realistic values of dew point temperature and virtual temperature profiles over the IP. The heterogeneity of the studied variables is highlighted in the mean maps over the IP. According to those for D, the unstable air masses are found along the entire Atlantic coast during winter, but in summer they are located particularly over the Mediterranean coast. The convective inhibition is more extended towards inland at 00:00 UTC in those areas. However, high values are also observed near the southeastern corner of the IP (near Murcia) at 12:00 UTC. Finally, no linear relationship between TT, CAPE, or CIN was found, and consequently, CAPE and CIN should be preferred for the study of the instability of the atmosphere as more atmospheric layers are employed during their calculation than for the TT index. [ABSTRACT FROM AUTHOR]

A new diagram is proposed for the verification of vector quantities generated by multiple models against a set of observations. It has been designed with the objective, as in the Taylor diagram, of providing a visual diagnostic tool which allows an easy comparison of simulations by multiple models against a reference dataset. However, the Sailor diagram extends this ability to two-dimensional quantities such as currents, wind, horizontal fluxes of water vapour and other geophysical variables by adding features which allow us to evaluate directional properties of the data as well. The diagram is based on the analysis of the two-dimensional structure of the mean squared error matrix between model and observations. This matrix is separated in a part corresponding to the bias and the relative rotation of the two orthogonal directions (empirical orthogonal functions; EOFs) which best describe the vector data. Since there is no truncation of the retained EOFs, these orthogonal directions explain the total variability of the original dataset. We test the performance of this new diagram to identify the differences amongst the reference dataset and a series of model outputs by using some synthetic datasets and real-world examples with time series of variables such as wind, current and vertically integrated moisture transport. An alternative setup for spatially varying time-fixed fields is shown in the last examples, in which the spatial average of surface wind in the Northern and Southern Hemisphere according to different reanalyses and realizations from ensembles of CMIP5 models are compared. The Sailor diagrams presented here show that it is a tool which helps in identifying errors due to the bias or the orientation of the simulated vector time series or fields. The R implementation of the diagram presented together with this paper allows us also to easily retrieve the individual diagnostics of the different components of the mean squared error and additional diagnostics which can be presented in tabular form. [ABSTRACT FROM AUTHOR]

The ability of two downscaling experiments to correctly simulate the instability conditions that can trigger thunderstorms over the Iberian Peninsula is compared in this paper. To do so, three instability indices are evaluated: TT index, CAPE and CIN. The WRF model is used for the simulations. The N experiment is driven by ERA-Interim's initial and boundary conditions; The D experiment has the same configuration as N, but the 3DVAR data assimilation step is additionally run at 00, 06, 12 and 18 UTC. Eight radiosondes are available over the IP, and the values for these indices calculated from the University of Wyoming were chosen as reference in the validation of both simulations. Additionally, measured variables at different pressure levels from the radiosondes provided byWyoming were used to calculate the three instability indices by our own methodology using the R package aiRthermo. According to the validation, the correlation, SD and RMSE obtained by the experiment D for all the indices in most of the stations are better than those for N. The different methodologies produce small discrepancies between the values for TT, but these are larger for CAPE and CIN due to the dependency of these indices on the initial conditions assumed for the calculation of an air parcel's vertical evolution. Similar results arise from the seasonal analysis concerning both WRF experiments: N tends to overestimate or underestimate (depending on the index) the variability of the reference values, but D is able to capture it in most of the seasons. The heterogeneity of the indices is highlighted in the mean maps over the Iberian Peninsula. According to those from D, the ingredients for the development of convective precipitation during winter are found along the entire Atlantic coast, but in summer they are located particularly in the Mediterranean coast. The chances of developing thunderstorms in those areas at 12 UTC is much higher than at 00 UTC; The convective inhibition is more extended towards inland at 00 UTC in those areas, which prevents storms from developing. However, high values are observed near Murcia also at 12 UTC. [ABSTRACT FROM AUTHOR]

A new diagram is proposed for the verification of vector quantities generated by multiple models against a set of observations. It has been designed with the idea of extending the Taylor diagram to two dimensional quantities such as currents, wind, or horizontal fluxes of water vapour, salinity, energy and other geophysical variables. The diagram is based on the analysis of the two-dimensional structure of the mean squared error matrix between model and observations. This matrix is separated in a part corresponding to the bias and the relative rotation of the empirical orthogonal functions of the data. We test the performance of this new diagram to identify the differences amongst the reference dataset and the different model outputs by using examples with wind, current, vertically integrated moisture transport and wave energy flux time series. An alternative setup is shown in the last examples, where the spatial average of surface wind in the Northern and Southern Hemispheres according to different reanalyses and realizations of CMIP5 models are compared. The examples of use of the Sailor diagram presented show that it is a tool which helps in identifying errors due to the bias or the orientation of the simulated vector time series or fields. The R implementation of the diagram presented together with this paper allows also to easily retrieve the individual diagnostics of the different components of the mean squared error. [ABSTRACT FROM AUTHOR]