You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 4 Next »

1. Ensemble version

TO BE UPDATED



Ensemble identifier codeJMA GEPS2203JMA GEPS2103JMA GEPS2003JMA GEPS1701GSM1403C
Short DescriptionGlobal ensemble system that simulates initial uncertainties using the Local Ensemble Transform Kalman Filter (LETKF), the singular vectors and lagged averaging forecasts, model uncertainties due to physical parameterizations using a stochastically perturbed physics tendencies (SPPT) scheme, and uncertainties on surface boundary conditions using sea surface temperature (SST) perturbation. Although it is uncoupled model, SST is relaxed to the ensemble-mean SST based on prediction using coupled model after lead time of 12 days. Ensembles are based on 50 members run once a week with 2 start dates (Tuesday and Wednesday at 12Z) up to day 34.Global ensemble system that simulates initial uncertainties using the Local Ensemble Transform Kalman Filter (LETKF), the singular vectors and lagged averaging forecasts, model uncertainties due to physical parameterizations using a stochastically perturbed physics tendencies (SPPT) scheme, and uncertainties on surface boundary conditions using sea surface temperature (SST) perturbation. Although it is uncoupled model, SST is relaxed to the ensemble-mean SST based on prediction using coupled model after lead time of 12 days. Ensembles are based on 50 members run once a week with 2 start dates (Tuesday and Wednesday at 12Z) up to day 34.Global ensemble system that simulates initial uncertainties using the Local Ensemble Transform Kalman Filter (LETKF), the singular vectors and lagged averaging forecasts, model uncertainties due to physical parameterizations using a stochastically perturbed physics tendencies (SPPT) scheme, and uncertainties on surface boundary conditions using sea surface temperature (SST) perturbation. Although it is uncoupled model, SST is relaxed to the ensemble-mean SST based on prediction using coupled model after lead time of 12 days. Ensembles are based on 50 members run once a week (Tuesday and Wednesday at 00Z and 12Z) up to day 34.Global ensemble system that simulates initial uncertainties using the Local Ensemble Transform Kalman Filter (LETKF), the singular vectors and lagged averaging forecasts, model uncertainties due to physical parameterizations using a stochastically perturbed physics tendencies (SPPT) scheme, and uncertainties on surface boundary conditions using sea surface temperature (SST) perturbation. Although it is uncoupled model, SST is relaxed to the ensemble-mean SST based on prediction using coupled model after lead time of 12 days. Ensembles are based on 50 members run once a week (Tuesday and Wednesday at 00Z and 12Z) up to day 34.Global ensemble system that simulates initial uncertainties using the bred vectors and lagged averaging forecasts and model uncertainties due to physical parameterizations using a stochastic scheme. Ensembles are based on 50 members, run once a week (Tuesday, Wednesday at 12Z) up to day 34.
Research or operationalOperationalOperationalOperationalOperationalOperational
Data time of first forecast run30/03/202130/03/202124/03/202022/03/201705/03/2014

2. Configuration of the EPS






Is the model coupled to an ocean model?NoNoNoNoNo
If yes, please describe ocean model briefly including frequency of coupling and any ensemble perturbation appliedN/AN/AN/AN/AN/A
Is the model coupled to a sea Ice model?NoNoNoNoNo
If yes, please describe sea-ice model briefly including any ensemble perturbation appliedN/AN/AN/AN/AN/A
Is the model coupled to a wave model?NoNoNoNoNo
If yes, please describe wave model briefly including any ensemble perturbation appliedN/AN/AN/AN/AN/A
Ocean modelN/AN/AN/AN/AN/A
Horizontal resolution of the atmospheric modelTL479 (about 40 km) up to 18 days, TL319 (about 55 km) after 18 days.TL479 (about 40 km) up to 18 days, TL319 (about 55 km) after 18 days.TL479 (about 40 km) up to 18 days, TL319 (about 55 km) after 18 days.TL479 (about 40 km) up to 18 days, TL319 (about 55 km) after 18 days.TL319 (about 55 km)
Number of model levels12812812812860
Top of model0.01 hPa0.01 hPa0.01 hPa0.01 hPa
Type of model levelshybrid (sigma-p) coordinatehybrid (sigma-p) coordinatehybrid (sigma-p) coordinatehybrid (sigma-p) coordinatesigma
Forecast length34 days (816 hours), but archived up to 32.5 days (780 hours)34 days (816 hours), but archived up to 32.5 days (780 hours)34 days (816 hours), but archived up to 32.5 days (780 hours)34 days (816 hours), but archived up to 32.5 days (780 hours)34 days (816 hours)
Run Frequencyonce a week (combination of Tuesday and Wednesday at 12Z)once a week (combination of Tuesday and Wednesday at 12Z)once a week (combination of Tuesday and Wednesday at 00Z and 12Z)once a week (combination of Tuesday and Wednesday at 00Z and 12Z)once a week (combination of Tuesday and Wednesday at 00Z and 12Z)
Is there an unperturbed control forecast included?YesYesYesYesYes
Number of perturbed ensemble members48 (totally 2 controls from 2 initial dates)48 (totally 2 controls from 2 initial dates)46 (4 controls from each initial date), but archived as 49 (1 control from each initial date)46 (4 controls from each initial date), but archived as 49 (1 control from each initial date)48 (2 controls from each initial date)
Integration time step12 minutes up to 18 days and 20 minutes after 18 days12 minutes up to 18 days and 20 minutes after 18 days12 minutes up to 18 days and 20 minutes after 18 days12 minutes up to 18 days and 20 minutes after 18 days20 minutes

3. Initial conditions and perturbations






Data assimilation method for control analysishybrid 4D Var-LETKFhybrid 4D Var-LETKFhybrid 4D Var-LETKF4D Var4D Var
Resolution of model used to generate Control AnalysisTL959L128TL959L128TL959L100TL959L100TL959L100
Ensemble initial perturbation strategyLETKF + singular vectors (initial SV; Northern Hemisphere, Tropics and Southern Hemisphere) + Lagged Average ForecastingLETKF + singular vectors (initial SV; Northern Hemisphere, Tropics and Southern Hemisphere) + Lagged Average ForecastingLETKF + singular vectors (initial SV; Northern Hemisphere, Tropics and Southern Hemisphere) + Lagged Average ForecastingLETKF + singular vectors (initial SV; Northern Hemisphere, Tropics and Southern Hemisphere) + Lagged Average ForecastingBred vectors (extratropics (NH) plus tropics) + Lagged Average Forecasting
Horizontal and vertical resolution of perturbationsTL319L128 (LETKF), T63L40 (SV)TL319L128 (LETKF), T63L40 (SV)TL319L100 (LETKF), T63L40 (SV)TL319L100 (LETKF), T63L40 (SV)TL319L60
Perturbations in +/- pairsYes (SV), No (LETKF), No (SST)Yes (SV), No (LETKF), No (SST)Yes (SV), No (LETKF), No (SST)Yes (SV), No (LETKF), No (SST)Yes
Initialization of land surface




3.1 What is the land surface model (LSM) and version used in the forecast model, and what are the current/relevant references for the model? Are there any significant changes/deviations in the operational version of the LSM from the documentation of the LSM?

LSM is the JMA SIB. The JMA SIB is based on the Simple Biosphere (SiB) developed by Sellers et al.(1986) and implemented by Sato et al.(1989a) and Sato et al.(1989b) for the land surface process in forecast model, and overall specifications are comprehensively updated and refined schemes are introduced. Significant changes from SiB are;

  • replacement the force-restore method for soil temperature prediction with a multilayer soil heat and water flux model and separate layers for snow

  • introduction of an equation of heat conduction with seven soil levels

  • consideration of the release or absorption of latent heat from phase change for soil temperature prediction

  • introduction of a new snow model with up to four layers with consideration for thermal diffusion, increased density from snow compaction and reduction of albedo due to snow aging

The Simple Biosphere (SiB) developed by Sellers et al.(1986), Sato et al.(1989a)

LSM is the JMA SIB. The JMA SIB is based on the Simple Biosphere (SiB) developed by Sellers et al.(1986) and implemented by Sato et al.(1989a) and Sato et al.(1989b) for the land surface process in forecast model, and overall specifications are comprehensively updated and refined schemes are introduced. Significant changes from SiB are;

  • replacement the force-restore method for soil temperature prediction with a multilayer soil heat and water flux model and separate layers for snow

  • introduction of an equation of heat conduction with seven soil levels

  • consideration of the release or absorption of latent heat from phase change for soil temperature prediction

  • introduction of a new snow model with up to four layers with consideration for thermal diffusion, increased density from snow compaction and reduction of albedo due to snow aging

The Simple Biosphere (SiB) developed by Sellers et al.(1986), Sato et al.(1989a)

LSM is the JMA SIB. The JMA SIB is based on the Simple Biosphere (SiB) developed by Sellers et al.(1986) and implemented by Sato et al.(1989a) and Sato et al.(1989b) for the land surface process in forecast model, and overall specifications are comprehensively updated and refined schemes are introduced. Significant changes from SiB are;

  • replacement the force-restore method for soil temperature prediction with a multilayer soil heat and water flux model and separate layers for snow

  • introduction of an equation of heat conduction with seven soil levels

  • consideration of the release or absorption of latent heat from phase change for soil temperature prediction

  • introduction of a new snow model with up to four layers with consideration for thermal diffusion, increased density from snow compaction and reduction of albedo due to snow aging

The Simple Biosphere (SiB) developed by Sellers et al.(1986), Sato et al.(1989a)

LSM is the JMA SIB. The JMA SIB is based on the Simple Biosphere (SiB) developed by Sellers et al.(1986) and implemented by Sato et al.(1989a) and Sato et al.(1989b) for the land surface process in forecast model, and overall specifications are comprehensively updated and refined schemes are introduced. Significant changes from SiB are;

  • replacement the force-restore method for soil temperature prediction with a multilayer soil heat and water flux model and separate layers for snow

  • introduction of an equation of heat conduction with seven soil levels

  • consideration of the release or absorption of latent heat from phase change for soil temperature prediction

  • introduction of a new snow model with up to four layers with consideration for thermal diffusion, increased density from snow compaction and reduction of albedo due to snow aging

The Simple Biosphere (SiB) developed by Sellers et al.(1986), Sato et al.(1989a)

The Simple Biosphere (SiB) developed by Sellers et al.(1986), Sato et al.(1989a) and Sato et al.(1989b) has been implemented for the land surface process in forecast model.
3.2 How is soil moisture initialized in the forecasts? (climatology / realistic / other)? If “realistic”, does the soil moisture come from an analysis using the same LSM as is coupled to the GCM for forecasts, or another source? Please describe the process of soil moisture initialization. Is there horizontal and/or vertical interpolation of initialization data onto the forecast model grid? If so, please give original data resolution(s).  Does the LSM differentiate between liquid and ice content of the soil? If so, how are each initialized? If all model soil layers are not initialized in the same way or from the same source, please describe.Soil moisture is initialized realistically by JMA’s soil moisture analysis based on Simplified Extended Kalman Filter (SEKF). Soil moisture and other initialization data are interpolated from the resolution of the deterministic system (TL959) to the resolution of the EPS forecast model grid (TL479). No vertical interpolation is applied. The LSM defines two soil moisture state variables, liquid and ice. The amounts of liquid and ice are determined by soil temperature and the ratio of the two soil moisture state variables cycled from the model forecast in the JMA Global Analysis (GA). To prevent unrealistic soil moisture drifting after long-term integration, initial soil moisture deeper than the 4th layer is set to the climatological value. See Ochi (2020) for details.Soil moisture is initialized realistically by JMA’s soil moisture analysis based on Simplified Extended Kalman Filter (SEKF). Soil moisture and other initialization data are interpolated from the resolution of the deterministic system (TL959) to the resolution of the EPS forecast model grid (TL479). No vertical interpolation is applied. The LSM defines two soil moisture state variables, liquid and ice. The amounts of liquid and ice are determined by soil temperature and the ratio of the two soil moisture state variables cycled from the model forecast in the JMA Global Analysis (GA). To prevent unrealistic soil moisture drifting after long-term integration, initial soil moisture deeper than the 4th layer is set to the climatological value. See Ochi (2020) for details.Soil moisture is initialized realistically by JMA’s soil moisture analysis based on Simplified Extended Kalman Filter (SEKF). Soil moisture and other initialization data are interpolated from the resolution of the deterministic system (TL959) to the resolution of the EPS forecast model grid (TL479). No vertical interpolation is applied. The LSM defines two soil moisture state variables, liquid and ice. The amounts of liquid and ice are determined by soil temperature and the ratio of the two soil moisture state variables cycled from the model forecast in the JMA Global Analysis (GA). To prevent unrealistic soil moisture drifting after long-term integration, initial soil moisture deeper than the 4th layer is set to the climatological value. See Ochi (2020) for details.Soil moisture is initialized realistically by JMA’s soil moisture analysis based on Simplified Extended Kalman Filter (SEKF). Soil moisture and other initialization data are interpolated from the resolution of the deterministic system (TL959) to the resolution of the EPS forecast model grid (TL479). No vertical interpolation is applied. The LSM defines two soil moisture state variables, liquid and ice. The amounts of liquid and ice are determined by soil temperature and the ratio of the two soil moisture state variables cycled from the model forecast in the JMA Global Analysis (GA). To prevent unrealistic soil moisture drifting after long-term integration, initial soil moisture deeper than the 4th layer is set to the climatological value. See Ochi (2020) for details.Initial soil moisture data (which consists of three layers) is produced by the offline simulations of the land surface model. The land processes in this simulations are similar to the one set in forecast model, and no horizontal and vertical interpolations are introduced in this analysis. Soil ice is handled as soil water, and no soil ice is used as initial condition specifically.
3.3 How is snow initialized in the forecasts? (climatology / realistic / other)? Is there horizontal and/or vertical interpolation of data onto the forecast model grid? If so, please give original data resolution(s) Are snow mass, snow depth or both initialized? What about snow age, albedo, or other snow properties?The global snow depth is analyzed using SYNOP snow depth data on the day. The forecasted snow depth corrected by satellite-estimated snow area is employed as first guess of the snow depth analysis. A two-dimensional optimum interpolation (OI) is employed for the analysis method. The analyzed snow depth is converted to the EPS forecast model grid (TL479). In addition, snow depth on Japan land grids is replaced with interpolated latest observation value of AMeDAS (Automated Meteorological Data Acquisition System) snow depth as well as Japanese SYNOP. The analyzed snow depth is converted to snow mass using snow density cycled from the GA model forecast and the snow mass is used as an initial condition for the LSM. Snow age is cycled from the GA model forecast and albedo is calculated with the snow age.The global snow depth is analyzed using SYNOP snow depth data on the day. The forecasted snow depth corrected by satellite-estimated snow area is employed as first guess of the snow depth analysis. A two-dimensional optimum interpolation (OI) is employed for the analysis method. The analyzed snow depth is converted to the EPS forecast model grid (TL479). In addition, snow depth on Japan land grids is replaced with interpolated latest observation value of AMeDAS (Automated Meteorological Data Acquisition System) snow depth as well as Japanese SYNOP. The analyzed snow depth is converted to snow mass using snow density cycled from the GA model forecast and the snow mass is used as an initial condition for the LSM. Snow age is cycled from the GA model forecast and albedo is calculated with the snow age.The global snow depth with 1.0 degree latitude/longitude resolution is analyzed using SYNOP snow depth data on the day. A two-dimensional optimum interpolation (OI) is employed for the analysis method. The analyzed snow depth is interpolated to the EPS forecast model grid (TL479). In addition, snow depth on Japan land grids are replaced with interpolated latest observation value of AMeDAS (Automated Meteorological Data Acquisition System) snow depth as well as Japanese SYNOP. The analyzed snow depth is converted to snow mass using snow density cycled from the GA model forecast and the snow mass is used as an initial condition for the LSM. Snow age is cycled from the GA model forecast and albedo is calculated with the snow age.The global snow depth with 1.0 degree latitude/longitude resolution is analyzed using SYNOP snow depth data on the day. A two-dimensional optimum interpolation (OI) is employed for the analysis method. The analyzed snow depth is interpolated to the EPS forecast model grid (TL479). In addition, snow depth on Japan land grids are replaced with interpolated latest observation value of AMeDAS (Automated Meteorological Data Acquisition System) snow depth as well as Japanese SYNOP. The analyzed snow depth is converted to snow mass using snow density cycled from the GA model forecast and the snow mass is used as an initial condition for the LSM. Snow age is cycled from the GA model forecast and albedo is calculated with the snow age.Initial snow data is also produced by using the offline simulation of the land surface model, and no horizontal interpolation is introduced. In the offline simulation, snow depth data is updated once a day (00UTC) by the two-dimensional Optimal Interpolation using the SYNOP snow depth data. The first guess is calculated by snow depth data from the offline simulation and snow cover data estimated by satellite observation. Forecast model calculates the snow water equivalent, so snow depth data is converted to the snow water equivalent. Snow density is set as a function related to the snow water equivalent (Verseghy 1991). Snow albedo is set as a function of wavelength and snow temperature. The age of snow is not considered.
3.4 How is soil temperature initialized in the forecasts? (climatology / realistic / other) Is the soil temperature initialized consistently with soil moisture (frozen soil water where soil temperature ≤0°C) and snow cover (top layer soil temperature ≤0°C under snow)? Is there horizontal and/or vertical interpolation of data onto the forecast model grid? If so, please give original data resolution(s)  If all model soil layers are not initialized in the same way or from the same source, please describe.Soil temperature cycled from the GA model forecast with resolution of TL959 is used as an initial condition for the LSM. It is interpolated to the EPS forecast model grid. The liquid and ice content of the soil initialized in the forecasts are determined by the soil temperature cycled from the GA model forecastSoil temperature cycled from the GA model forecast with resolution of TL959 is used as an initial condition for the LSM. It is interpolated to the EPS forecast model grid. The liquid and ice content of the soil initialized in the forecasts are determined by the soil temperature cycled from the GA model forecastSoil temperature cycled from the GA model forecast with resolution of TL959 is used as an initial condition for the LSM. It is interpolated to the EPS forecast model grid. The liquid and ice content of the soil initialized in the forecasts are determined by the soil temperature cycled from the GA model forecastSoil temperature cycled from the GA model forecast with resolution of TL959 is used as an initial condition for the LSM. It is interpolated to the EPS forecast model grid. The liquid and ice content of the soil initialized in the forecasts are determined by the soil temperature cycled from the GA model forecastSoil temperature is also initialized by the offline simulations of the land surface model. No horizontal and vertical interpolation is implemented. Note that soil layers are three for soil moisture, while it is only one layer for soil temperature. Snow cover is updated at each 00UTC based on the snow depth analysis. At the same time, soil temperature for all grids where snow exists is set as less than 0 deg. Once the soil temperature becomes less than 0deg, soil water changes soil ice (No consideration about freeze latent heat). No soil water scatters or moves in the freezing soil.
3.5 How are time-varying vegetation properties represented in the LSM? Is phenology predicted by the LSM? If so, how is it initialized? If not, what is the source of vegetation parameters used by the LSM? Which time-varying vegetation parameters are specified (e.g., LAI, greenness, vegetation cover fraction) and how (e.g., near-real-time satellite observations? Mean annual cycle climatology? Monthly, weekly or other interval?)There is no consideration about time-varying of vegetation properties. The climatology for each vegetation based on Dorman and Sellers (1989) is set as outer parameters (Some parameters are monthly data, such as LAI, the ratio of green leaves, the ratio of vegetation)There is no consideration about time-varying of vegetation properties. The climatology for each vegetation based on Dorman and Sellers (1989) is set as outer parameters (Some parameters are monthly data, such as LAI, the ratio of green leaves, the ratio of vegetation)There is no consideration about time-varying of vegetation properties. The climatology for each vegetation based on Dorman and Sellers (1989) is set as outer parameters (Some parameters are monthly data, such as LAI, the ratio of green leaves, the ratio of vegetation)There is no consideration about time-varying of vegetation properties. The climatology for each vegetation based on Dorman and Sellers (1989) is set as outer parameters (Some parameters are monthly data, such as LAI, the ratio of green leaves, the ratio of vegetation)There is no consideration about time-varying of vegetation properties. The climatology for each vegetation based on Dorman and Sellers (1989) is set as outer parameters (Some parameters are monthly data, such as LAI, the ratio of green leaves, the ratio of vegetation)
3.6 What is the source of soil properties (texture, porosity, conductivity, etc.) used by the LSM?

Soil thermal and hydraulic properties are calculated using soil texture from the HWSD (Harmonized World Soil Database).

Soil thermal and hydraulic properties are calculated using soil texture from the HWSD (Harmonized World Soil Database).

Soil thermal and hydraulic properties are calculated using soil texture from the HWSD (Harmonized World Soil Database).

Soil thermal and hydraulic properties are calculated using soil texture from the HWSD (Harmonized World Soil Database).

The source of soil properties is different depending on the property. For example, soil porosity is set as outer parameters for each type of vegetation, while soil heat conductivity is set as a function related to the porosity and soil moisture in the first soil layer. In addition, some parameters are not considered such as difference of soil texture.
3.7 If the initialization of the LSM for re-forecasts deviates from the procedure for forecasts, please describe the differences.Land surface values are estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55 (Kobayashi et al. 2015)Land surface values are estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55 (Kobayashi et al. 2015)Land surface values are estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55 (Kobayashi et al. 2015)Land surface values are estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55 (Kobayashi et al. 2015)There is no difference between re-forecast and operational forecast about the procedure for initialization of land surface. For the re-forecast, land analysis data of JRA-55 is utilized as the land initial data and it is derived from the offline system forced by JRA-55 atmospheric field. This system is similar to the operational system, but note that the atmospheric forcing for operational offline system is given from the operational Global Analysis.

4. Model uncertainties perturbations






Is model physics perturbed?  If yes, briefly describe methodsStochastically perturbed physics tendencies (SPPT) schemeStochastically perturbed physics tendencies (SPPT) schemeStochastically perturbed physics tendencies (SPPT) schemeStochastically perturbed physics tendencies (SPPT) schemeStochastic physics
Do all ensemble members use exactly the same model version?SameSameSameSameSame
Is model dynamics perturbed?NoNoNoNoNo
Are the above model perturbations applied to the control forecast?NoNoNoNoNo

5. Surface boundary perturbations






Perturbations to sea surface temperature?YesYesYesYesNo
Perturbation to soil moisture?NoNoNoNoNo
Perturbation to surface stress or roughness?NoNoNoNoNo
Any other surface perturbation?NoNoNoNoNo
Are the above surface perturbations applied to the Control forecast?NoNoNoNoNo
Additional commentsNoneNoneNoneNoneNone

6. Other details of the models






Description of model gridLinear gridLinear gridLinear gridLinear gridLinear grid
List of model levels in appropriate coordinatesSee appendixSee appendixSee appendixSee appendixhttp://jra.kishou.go.jp/JRA-55/document/JRA-55_handbook_TL319_v2_en.pdf
What kind of large scale dynamics is used?Spectral semi-lagrangianSpectral semi-lagrangianSpectral semi-lagrangianSpectral semi-lagrangianSpectral semi-lagrangian
What kind of boundary layer parameterization is used?Mellor and Yamada level 2 and diffusive coefficients based on Han and Pan (2011) for stable BL, JMA (2019)Mellor and Yamada level 2 and diffusive coefficients based on Han and Pan (2011) for stable BL, JMA (2019)Mellor and Yamada level 2 and diffusive coefficients based on Han and Pan (2011) for stable BL, JMA (2019)Mellor and Yamada level 2 and diffusive coefficients based on Han and Pan (2011) for stable BL, Yonehara et al. (2014)Mellor and Yamada level 2.5
What kind of convective parameterization is used?Arakawa and Schubert (1974), JMA (2019Arakawa and Schubert (1974), JMA (2019Arakawa and Schubert (1974), JMA (2019)Arakawa and Schubert (JMA 2013), Yonehara et al. (2014), Yonehara et al. (2017)Arakawa and Schubert (JMA 2013)
What kind of large-scale precipitation scheme is used?Sundqvist (1978), JMA (2019)Sundqvist (1978), JMA (2019)Sundqvist (1978), JMA (2019)Sundqvist (1978), Yonehara et al. (2017)Sundqvist (1978)
What cloud scheme is used?Smith (1990), Kawai and Inoue (2006), JMA (2019)Smith (1990), Kawai and Inoue (2006), JMA (2019)Smith (1990), Kawai and Inoue (2006), JMA (2019)Smith (1990), Kawai and Inoue (2006), Yonehara et al. (2017)Smith (1990), Kawai and Inoue (2006), Yonehara et al. (2017)
What kind of land-surface scheme is used?JMA-SIB, JMA (2019), Yonehara et al. (2020)JMA-SIB, JMA (2019), Yonehara et al. (2020)JMA-SIB, JMA (2019), Yonehara et al. (2020)JMA-SIB, Yonehara et al. (2017)SiB (Sato et al. 1989)
How is radiation parametrized?
  • Longwave radiation: JMA (2019)
  • Shortwave radiation: JMA (2019)
  • Longwave radiation: JMA (2019)
  • Shortwave radiation: JMA (2019)
  • Longwave radiation: JMA (2019)
  • Shortwave radiation: JMA (2019)
  • Longwave radiation: Yabu (2013)
  • Shortwave radiation: JMA (2013), Nagasawa (2012)
Outline of the operational numerical weather prediction at the Japan Meteorological Agency

7. Re-forecast configuration






Number of years covered40 years (1981-2020)40 years (1981-2020)30 years (1981-2010)32 years (1981-2012)30 years (1981-2010)
Produced on the fly or fix re-forecasts?Fixed re-forecasts in advanceFixed re-forecasts in advanceFixed re-forecasts in advanceFixed re-forecasts in advanceFixed re-forecasts in advance
FrequencyThe re-forecasts consist of a 13-member ensemble starting the 15th and the last dates of calendar months.The re-forecasts consist of a 13-member ensemble starting the 15th and the last dates of calendar months.The re-forecasts consist of a 13-member ensemble starting the 15th and the last dates of calendar months.The re-forecasts consists of a 5 member ensemble starting the 10th, 20th, the last dates of calendar months.The re-forecasts consists of a 5 member ensemble starting the 10th, 20th, the last dates of calendar months.
Ensemble size13 members13 members13 members5 members5 members
Initial conditionsJRA-55 (TL319L60) + land surface values estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55JRA-55 (TL319L60) + land surface values estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55JRA-55 (TL319L60) + land surface values estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55JRA-55 (TL319L60) + land surface values estimated with the land-surface model in the Global EPS using atmospheric forcing from JRA-55JRA-55 (TL319L60) + JRA-55 land analysis (TL319)
Is the model physics and resolution the same as for the real-time forecastsYesYesYesYesYes
If not, what are the differencesN/AN/AN/AN/AN/A
Is the ensemble generation the same as for real-time forecasts?NoNoNoNoYes, except for lagged average forecasting.
If not, what are the differencesLETKF perturbations are not used and singular vectors (initial SV + evolved SV) are usedLETKF perturbations are not used and singular vectors (initial SV + evolved SV) are usedLETKF perturbations are not used and singular vectors (initial SV + evolved SV) are usedLETKF perturbations are not used and singular vectors (initial SV + evolved SV) are usedN/A
Other relevant information

The JMA re-forecasts dataset is a "fixed" dataset which means that the re-forecasts are produced once from a "frozen" version of the model and are used for a number of years to calibrate real-time forecast. The JMA re-forecasts consist of a 13-member ensemble running twice a month from 1981 to 2020. The start dates correspond to 1st and 16st of each month at 00Z minus 12 hours (28 February instead of 29 February). Here is the complete list of re-forecast start dates:

15/31 January - 15/28 February - 15/31 March - 15/30 April - 15/31 May - 15/30 June - 15/31 July - 15/31 August - 15/30 September - 15/31 October - 15/30 November and 15/31 December 1981-2020

The S2S database contains the complete JMA re-forecast dataset.


The JMA re-forecasts are archived in the S2S database with 2 date attributes:

  • hdate which corresponds to the actual starting date of the re-forecast
  • date which corresponds to the ModelVersionDate. Since the JMA re-forecasts are "fixed" re-forecasts this ModelVersiondate is the same for all the re-forecasts and equal to 20210331. This variable will change when a new version of the JMA model for extended-range forecasts will be implemented.

The JMA re-forecasts dataset is a "fixed" dataset which means that the re-forecasts are produced once from a "frozen" version of the model and are used for a number of years to calibrate real-time forecast. The JMA re-forecasts consist of a 13-member ensemble running twice a month from 1981 to 2020. The start dates correspond to 1st and 16st of each month at 00Z minus 12 hours (28 February instead of 29 February). Here is the complete list of re-forecast start dates:

15/31 January - 15/28 February - 15/31 March - 15/30 April - 15/31 May - 15/30 June - 15/31 July - 15/31 August - 15/30 September - 15/31 October - 15/30 November and 15/31 December 1981-2020

The S2S database contains the complete JMA re-forecast dataset.


The JMA re-forecasts are archived in the S2S database with 2 date attributes:

  • hdate which corresponds to the actual starting date of the re-forecast
  • date which corresponds to the ModelVersionDate. Since the JMA re-forecasts are "fixed" re-forecasts this ModelVersiondate is the same for all the re-forecasts and equal to 20210331. This variable will change when a new version of the JMA model for extended-range forecasts will be implemented.

The JMA re-forecasts dataset is a "fixed" dataset which means that the re-forecasts are produced once from a "frozen" version of the model and are used for a number of years to calibrate real-time forecast. The JMA re-forecasts consist of a 13-member ensemble running twice a month from 1981 to 2020. The start dates correspond to 1st and 16st of each month at 00Z minus 12 hours (28 February instead of 29 February). Here is the complete list of re-forecast start dates:

15/31 January - 15/28 February - 15/31 March - 15/30 April - 15/31 May - 15/30 June - 15/31 July - 15/31 August - 15/30 September - 15/31 October - 15/30 November and 15/31 December 1981-2010

The S2S database contains the complete JMA re-forecast dataset.


The JMA re-forecasts are archived in the S2S database with 2 date attributes:

  • hdate which corresponds to the actual starting date of the re-forecast
  • date which corresponds to the ModelVersionDate. Since the JMA re-forecasts are "fixed" re-forecasts this ModelVersiondate is the same for all the re-forecasts and equal to 20200331. This variable will change when a new version of the JMA model for extended-range forecasts will be implemented.

The JMA real-time ensemble forecasts include 4 start dates (Tuesdays at 00Z and 12Z + Wednesdays at 00Z + 12Z). The ensemble size is 13 members for Tuesdays 12Z and Wednesdays 00 and 12Z and 11 members for Tuesdays 00Z. For user convenience, ensemble data from 4 start dates are archived as a single 50-member ensemble starting on Wednesdays at 12Z. List of initialize date and index of archived member in final 50 member ensemble are as follows;

  • Wednesday 12 Z : 0,1, …, 12
  • Wednesday 00 Z : 13,14, …, 25
  • Tuesday 12 Z : 26, 27, …, 38
  • Tuesday 00 Z : 39, 40, …, 49

Note that member 0, 13, 26 and 39 are control forecast from each initial date, however, they are encoded as pf for user convenience.

The JMA re-forecasts dataset is a "fixed" dataset which means that the re-forecasts are produced once from a "frozen" version of the model and are used for a number of years to calibrate real-time forecast. The JMA re-forecasts consist of a 5-member ensemble running three times a month from 1981 to 2012. The start dates correspond to 1st / 11th and 21st of each month at 00Z minus 12 hours (28 February instead of 29 February). here is the complete list of re-forecast start dates:

10/20/31 January - 10/20/28 February - 10/20/31 March - 10/20/30 April - 10/20/31 May - 10/20/30 June - 10/20/31 July - 10/20/31 August - 10/20/30 September - 10/20/31 October - 10/20/30 November and 10/20/31 December 1981-2012

The S2S database contains the complete JMA re-forecast dataset.

The JMA real-time ensemble forecasts include 4 start dates (Tuesdays at 00Z and 12Z + Wednesdays at 00Z + 12Z). The ensemble size is 13 members for Tuesdays 12Z and Wednesdays 00 and 12Z and 11 members for Tuesdays 00Z. For user convenience, ensemble data from 4 start dates are archived as a single 50-member ensemble starting on Wednesdays at 12Z. List of initialize date and index of archived member in final 50 member ensemble are as follows;

  • Wednesday 12 Z : 0,1, …, 12

  • Wednesday 00 Z : 13,14, …, 25

  • Tuesday 12 Z : 26, 27, …, 38

  • Tuesday 00 Z : 39, 40, …, 49

Note that member 0, 13, 26 and 39 are control forecast from each initial date, however, they are encoded as pf for user convenience.

The JMA re-forecasts are archived in the S2S database with 2 date attributes:

  • hdate which corresponds to the actual starting date of the re-forecast

  • date which correspond to the ModelVersionDate. Since the JMA re-forecasts are "fixed" re-forecasts this ModelVersiondate is the same for all the re-forecasts and equal to 20170131. This variable will change when a new version of the JMA model for extended-range forecasts will be implemented.

The JMA re-forecasts dataset is a "fixed" dataset which means that the re-forecasts are produced once from a "frozen" version of the model and are used for a number of years to calibrate real-time forecast. The JMA re-forecasts consist of a   5-member ensemble running three times a month from 1981 to 2010. The start dates correspond to 1st / 11th and 21st of each month at 00Z minus 12 hours (28 February instead of 29 February). here is the complete list of re-forecast start dates:

10/20/31 January - 10/20/28 February - 10/20/31 March - 10/20/30 April - 10/20/31 May - 10/20/30 June - 10/20/31 July - 10/20/31 August - 10/20/30 September - 10/20/31 October - 10/20/30 November and 10/20/31 December 1981-2010

The S2S database contains the complete JMA re-forecast dataset.

As for the other models, JMA re-forecasts are archived in the S2S database with 2 date attributes:

  • hdate which corresponds to the actual starting date of the re-forecast
  • date which correspond tot he ModelVersionDate.Since the JMA re-forecasts are "fixed" re-forecasts this ModelVersiondate is the same for all the re-forecasts and equal to  20140304. This variable will change when a new version of the JMA model for extended-range forecasts will be implemented.

8. References

  • Arakawa, A. and W. H. Schubert, 1974: Interaction of a Cumulus Cloud Ensemble with the Large-Scale Environment, Part I. J. Atmos. Sci., 31, 674-701.
  • Dorman, J. L. and P. J. Sellers, 1989: A global climatology of albedo, roughness length and stomatal resistance for atmospheric general circulation models as represented by the simple biosphere model (SiB). J. Appl.Meteor., 28, 833–855.
  • Han, J. and H-L, Pan. 2011: Revision of convection and vertical diffusion schemes in the NCEP global forecast system. Weather and Forecasting, 26.4, 520-533.
  • Japan Meteorological Agency, 2019: Outline of the operational numerical weather prediction at the Japan Meteorological Agency.
  • Kobayashi, S., Y. Ota, Y. Harada, A. Ebita, M. Moriya,H. Onoda, K. Onogi, H. Kamahori, C. Kobayashi, H. Endo, K. Miyaoka, and K. Takahashi, 2015: The JRA-55 reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan. 93, 5-48.
  • Lott, F. and M.J. Miller, 1997: A new subgrid-scale orographic drag parameterization: Its formulation and testing. Quart. J. Roy. Meteor. Soc., 123, 101–127.
  • Ochi, K., 2020: Preliminary results of soil moisture data assimilation into JMA Global Analysis. CAS/JSC WGNE Res. Activ. Atmos. Oceanic Modell., 50, 1.15-1.16.
  • Sato, N., P. J. Sellers, D. A. Randall, E. K. Schneider, J. Shukla, J. L. Kinter III, Y-T Hou, and E. Albertazzi, 1989a: Effects of implementing the simple biosphere model (SiB) in a general circulation model. J. Atmos. Sci., 46, 2757–2782.
  • Sato, N., P. J. Sellers, D. A. Randall, E. K. Schneider, J. Shukla, J. L. Kinter III, Y-T Hou, and E. Albertazzi, 1989b: Implementing the simple biosphere model (SiB) in a general circulation model: Methodologies and results. NASA contractor Rep. 185509, NASA. 76pp.
  • Scinocca, J. F. 2003: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models. J. Atmos. Sci., 60(4), 667-682.
  • Sellers, P. J., Y. Mintz, Y. C. Sud, and A. Dalcher, 1986: A simple biosphere model (SiB) for use within general circulation models. J. Atmos. Sci., 43, 505–531.
  • Simmons, A. J. and D. M. Burridge, 1981: An Energy and Angular-Momentum Conserving Vertical Finite-Difference Scheme and Hybrid Vertical Coordinates. Mon. Wea. Rev., 109, 758-766.
  • Takakura, T., and T. Komori, 2020: Two-tiered sea surface temperature approach implemented to JMA’s Global Ensemble Prediction System, CAS/JSC WGNE Res. Activ. Atmos. Oceanic Modell., 50, 6.15-6.16.
  • Yonehara, H., C. Matsukawa, T. Nabetani, T. Kanehama, T. Tokuhiro, K. Yamada, R. Nagasawa, Y. Adachi, and R. Sekiguchi, 2020: Upgrade of JMA’s Operational Global Model. CAS/JSC WGNE Res. Activ. Atmos. Oceanic Modell., 50, 6.19-6.20.

Appendix. Hybrid coordinates

Model level fields are produced for 128 hybrid levels. Each hybrid level is defined with half-levels  as the boundary;

where  is the surface pressure. Coefficients A and B are given in Table A for k = 0, 1, 2, …, 128. The following equation by Simmons and Burridge (1981) gives full-level pressure;

where C=1 and k=1, 2, …, 127. The full-level pressure for the uppermost level (k=128) is given by

Table A gives half-level and full-level pressures with a surface pressure of 1000hPa.


Table A. Model level from 1 to 128.

k

A[Pa]

B

Ph [Pa]

Pf [Pa]

1

0.000000000000

1.000000000000

100000.000000000000

99906.149328801500

2

0.367279053361

0.998119607566

99812.328035645300

99707.078615143900

3

1.651295433209

0.996002149193

99601.866214743700

99480.841827595100

4

4.263935162251

0.993556025633

99359.866498482900

99219.297574726400

5

8.818642865415

0.990699763605

99078.795003366000

98915.476159705600

6

16.153564274720

0.987360935873

98752.247151551600

98563.497767488800

7

27.352655050102

0.983475161368

98374.868791829000

98158.495693304200

8

43.760653705889

0.978985208099

97942.281463645000

97696.543075913600

9

66.989186167242

0.973840213666

97451.010552792600

97174.582116438100

10

98.912437237234

0.967995031094

96898.415546672800

96590.355173546100

11

141.651809373079

0.961409701475

96282.621956901400

95942.337440485500

12

197.549799934382

0.954049049521

95602.454752054800

95229.671150837900

13

269.133969108446

0.945882393768

94857.373345948900

94452.101429317500

14

359.072347557138

0.936883359805

94047.408328036600

93609.914021442900

15

470.121954299976

0.927029782579

93173.100212176600

92703.875201491000

16

605.072274485724

0.916303682526

92235.440527052200

91735.174189560300

17

766.685600539283

0.904691299836

91235.815584124600

90705.368407497400

18

957.636088881686

0.892183171556

90175.953244477800

89616.331877189700

19

1180.449250079540

0.878774237209

89057.872970937400

88470.207023681400

20

1437.443395530380

0.864463960067

87883.839402250300

87269.360090425000

21

1730.674330567250

0.849256452964

86656.319626988400

86016.340312364600

22

2061.884332303790

0.833160599383

85377.944270597100

84713.842927258300

23

2432.456198180200

0.816190162450

84051.472443213600

83364.676041107800

24

2843.372912386710

0.798363876191

82679.760531518600

81971.731299931800

25

3295.183263110750

0.779705514941

81265.734757235400

80537.958263339100

26

3787.973561424270

0.760243938082

79812.367369670600

79066.342322563600

27

4321.345466897210

0.740013108242

78322.656291106200

77559.885961130400

28

4894.399817081290

0.719052081756

76799.607992665500

76021.593118835800

29

5505.726286743880

0.697404970581

75246.223344854400

74454.456390393700

30

6153.398665018050

0.675120874964

73665.486161433700

72861.446768227600

31

6834.975529818620

0.652253786076

72060.354137414400

71245.505624873700

32

7547.506113219540

0.628862457583

70433.751871486600

69609.538623254500

33

8287.541182649680

0.605010244769

68788.565659596100

67956.411242983100

34

9051.148804124690

0.580764909450

67127.639749119800

66288.945616008800

35

9833.934898763650

0.556198388528

65453.773751564300

64609.918376429900

36

10631.068546365100

0.531386523789

63769.720925242600

62922.059244645800

37

11437.312024228700

0.506408750331

62078.187057363800

61228.050086197400

38

12247.055590877300

0.481347741058

60381.829696701600

59530.524208817900

39

13054.357028916700

0.456289004840

58683.257512887600

57832.065687358200

40

13852.985946051900

0.431320436397

56985.029585744600

56135.208534621000

41

14636.472796443700

0.406531816609

55289.654457338000

54442.435566027300

42

15398.162525485300

0.382014262845

53599.588809964400

52756.176837013500

43

16131.272660216100

0.357859630043

51917.235664546200

51078.807563275800

44

16828.955566715800

0.334159864585

50244.942025257400

49412.645465099300

45

17484.364477778600

0.311006314494

48584.995927132800

47759.947507769900

46

18090.722762945100

0.288489001104

46939.622873350700

46122.906038724200

47

18641.395773489400

0.266695859038

45310.981677293800

44503.644350057200

48

19129.964452997900

0.245711952979

43701.159750914300

42904.211719761400

49

19550.299766124400

0.225618681387

42112.167904844900

41326.578007900800

50

19896.636870744900

0.206492978760

40545.934746732300

39772.627903007100

51

20163.980086975300

0.188403206951

39004.300782028000

38244.154930225900

52

20351.944344566100

0.171370679911

37489.012335655400

36742.855344864500

53

20463.342581962100

0.155383728414

36001.715423361600

35270.322043458200

54

20501.374046673500

0.140425756613

34543.949708016400

33828.038628992400

55

20469.538812267800

0.126476038663

33117.142678580100

32417.373767254300

56

20371.603412767300

0.113510007752

31722.604187988600

31039.575968092000

57

20211.564215127300

0.101499572658

30361.521480941200

29695.768918140300

58

19993.608937901800

0.090413458959

29034.954833767600

28386.947481304800

59

19722.076787294600

0.080217571292

27743.833916544400

27113.974469751500

60

19401.417732755300

0.070875372401

26488.954972838800

25877.578272274200

61

19036.151481018200

0.062348274144

25270.978895384400

24678.351408969000

62

18630.826728478400

0.054596035287

24090.430257184200

23516.750061661900

63

18189.981276200600

0.047577160614

22947.697337580200

22393.094609378300

64

17718.103579492800

0.041249295828

21843.033162334100

21307.571177633000

65

17219.596275363200

0.035569612810

20776.557556335600

20260.234190198700

66

16698.742187492900

0.030495180003

19748.260187780400

19251.009892823500

67

16159.673251338300

0.025983313127

18758.004564073600

18279.700800639100

68

15606.342733896800

0.021991901889

17805.532922819900

17345.991005115600

69

15042.501046123400

0.018479709003

16890.471946455900

16449.452262833600

70

14471.675363861100

0.015406638528

16012.339216706500

15589.550777349100

71

13897.153188406700

0.012733971269

15170.550315345700

14765.654577045400

72

13321.969893417000

0.010424565774

14364.426470857500

13977.041386529500

73

12748.900223548600

0.008443024230

13593.202646582800

13222.906886473900

74

12180.453634514400

0.006755823303

12856.035964764400

12502.373257022700

75

11618.873296265500

0.005331410662

12152.014362450400

11814.497902710800

76

11066.138522444300

0.004140268568

11480.165379269900

11158.282262021800

77

10523.970341313200

0.003154946421

10839.464983428800

10532.680612040700

78

9993.839886724720

0.002350064638

10228.846350531500

9936.608787727700

79

9476.979262550330

0.001702292571

9647.208519692700

9368.952745786510

80

8974.394520051310

0.001190303424

9093.424862452510

8828.576914590710

81

8486.880384221110

0.000794709276

8566.351311844540

8314.332283627920

82

8015.036371079950

0.000497979401

8064.834311166120

7825.064198118570

83

7559.283951848330

0.000284345023

7587.718454119770

7359.619836325270

84

7119.884440234330

0.000139693594

7133.853799600960

6916.855358210030

85

6696.957303929350

0.000051455511

6702.102855046850

6495.642724032220

86

6290.498628904240

0.000008486026

6291.347231482130

6094.876189775000

87

5900.493980737070

0.000000000000

5900.493980737070

5713.478492448980

88

5528.481630297230

0.000000000000

5528.481630297230

5350.406741907360

89

5174.285933388900

0.000000000000

5174.285933388900

5004.658036279840

90

4836.925350720490

0.000000000000

4836.925350720490

4675.274815057110

91

4515.466275302810

0.000000000000

4515.466275302810

4361.349956723390

92

4209.028002503950

0.000000000000

4209.028002503950

4062.031616233200

93

3916.787433562190

0.000000000000

3916.787433562190

3776.527781221010

94

3637.983481907540

0.000000000000

3637.983481907540

3504.110504425720

95

3371.921127743180

0.000000000000

3371.921127743180

3244.119743480640

96

3117.975037623500

0.000000000000

3117.975037623500

2995.966708364220

97

2875.592632860050

0.000000000000

2875.592632860050

2759.136582338770

98

2644.296454685800

0.000000000000

2644.296454685800

2533.190445678590

99

2423.685637127690

0.000000000000

2423.685637127690

2317.766195254820

100

2213.436263278170

0.000000000000

2213.436263278170

2112.578220408930

101

2013.300350929100

0.000000000000

2013.300350929100

1917.415570860470

102

1823.103194212330

0.000000000000

1823.103194212330

1732.138341042090

103

1642.738784867620

0.000000000000

1642.738784867620

1556.672003523160

104

1472.163056683520

0.000000000000

1472.163056683520

1390.999458895190

105

1311.384746465990

0.000000000000

1311.384746465990

1235.150637357940

106

1160.453751007800

0.000000000000

1160.453751007800

1089.189593919690

107

1019.446986748220

0.000000000000

1019.446986748220

953.199187899828

108

888.451928822532

0.000000000000

888.451928822532

827.263627663412

109

767.548215932542

0.000000000000

767.548215932542

711.449386955591

110

656.787947376753

0.000000000000

656.787947376753

605.785246379858

111

556.175551249932

0.000000000000

556.175551249932

510.242460795943

112

465.648342355864

0.000000000000

465.648342355864

424.716270813441

113

385.059081069088

0.000000000000

385.059081069088

349.010127626377

114

314.161951254557

0.000000000000

314.161951254557

282.824045691573

115

252.603356889854

0.000000000000

252.603356889854

225.748400577188

116

199.918760089339

0.000000000000

199.918760089339

177.264223682977

117

155.536429527503

0.000000000000

155.536429527503

136.750604725015

118

118.788443103480

0.000000000000

118.788443103480

103.499217296971

119

88.928627808076

0.000000000000

88.928627808076

76.735284754437

120

65.156391730376

0.000000000000

65.156391730376

55.643586362110

121

46.644705017873

0.000000000000

46.644705017873

39.397472794456

122

32.569931284291

0.000000000000

32.569931284291

27.188427315847

123

22.140906946455

0.000000000000

22.140906946455

18.253569151567

124

14.624692644256

0.000000000000

14.624692644256

11.898702444718

125

9.366805063801

0.000000000000

9.366805063801

7.515060935547

126

5.804438422468

0.000000000000

5.804438422468

4.588710053849

127

3.472094398676

0.000000000000

3.472094398676

2.702510420768

128

2.000000000000

0.000000000000

2.000000000000

1.000000000000

129

0.000000000000

0.000000000000

0.000000000000


 

  • No labels