ocada notes questions and discussion
Discussion goes here
Discussion goes here
Masayoshi ISHII, Hirotaka KAMAHORI, Hisayuki KUBOTA, Masumi ZAIKI, Ryo MIZUTA, Hiroaki KAWASE, Masaya NOSAKA, Hiromasa YOSHIMURA, Naga OSHIMA, Eiki SHINDO, Hiroshi KOYAMA, Masato MORI, Shoji HIRAHARA, Yukiko IMADA, Kohei YOSHIDA, Toru NOZAWA, Tetsuya TAKEMI, Takashi MAKI, Akio NISHIMURA, Global Historical Reanalysis with a 60-km AGCM and Surface Pressure Observations: OCADA, Journal of the Meteorological Society of Japan. Ser. II, 2024, Volume 102, Issue 2, Pages 209-240, Released on J-STAGE March 09, 2024, Advance online publication January 17, 2024, Online ISSN 2186-9057, Print ISSN 0026-1165, https://doi.org/10.2151/jmsj.2024-010, https://www.jstage.jst.go.jp/article/jmsj/102/2/102_2024-010/_article/-char/en,
Abstract:
A historical atmospheric reanalysis from 1850 to 2015 was performed using an atmospheric general circulation model assimilating surface pressure observations archived in international databases, with perturbed observational sea surface temperatures as a lower boundary condition. Posterior spread during data assimilation provides quantitative information on the uncertainty in the historical reanalysis. The reanalysis reproduces the evolution of the three-dimensional atmosphere close to those of the operational centers. Newly archived surface pressure observations greatly reduced the uncertainties in the present reanalysis over East Asia in the early 20th century. A scheme for assimilating tropical cyclone tracks and intensities was developed. The scheme was superior to the present several reanalyses in reproducing the intensity close to the observations and the positions. The reanalysis provides possible images of atmospheric circulations before reanalyses with full-scale observations become available, and opportunities for investigating extreme events that occurred before World War II. Incorporating dynamical downscaling with a regional model that includes detailed topography and sophisticated physics is an application of historical reanalysis to reveal the details of past extreme events. Some examples of past heavy rainfall events in Japan are shown using a downscaling experiment, together with dense rainfall observations over the Japanese islands.
A historical atmospheric reanalysis from 1850 to 2015 was performed using an atmospheric general circulation model assimilating surface pressure observations archived in international databases, with perturbed observational sea surface temperatures as a lower boundary condition. Posterior spread during data assimilation provides quantitative information on the uncertainty in the historical reanalysis. The reanalysis reproduces the evolution of the three-dimensional atmosphere close to those of the operational centers. Newly archived surface pressure observations greatly reduced the uncertainties in the present reanalysis over East Asia in the early 20th century. A scheme for assimilating tropical cyclone tracks and intensities was developed. The scheme was superior to the present several reanalyses in reproducing the intensity close to the observations and the positions. The reanalysis provides possible images of atmospheric circulations before reanalyses with full-scale observations become available, and opportunities for investigating extreme events that occurred before World War II. Incorporating dynamical downscaling with a regional model that includes detailed topography and sophisticated physics is an application of historical reanalysis to reveal the details of past extreme events. Some examples of past heavy rainfall events in Japan are shown using a downscaling experiment, together with dense rainfall observations over the Japanese islands.
An incomplete list can be found here:
https://www.cpc.ncep.noaa.gov/products/wesley/CORe/references.html
There are a number of variables of potential interest to researches investigating links between weather/climate and infectious diseases. The most common ones 2m temperature, 2m relative humidity, 2m specific humidity, and precipitation. There are other potential variables including 2m Tmax and Tmin. Most reanalyses make data available at monthly and daily (or sub daily time scales).
Here are ways to get these variables from different reanalysis is not an specific endorsement for this research. Most values will have an "error bar". Researchers should carefully read documentation on limitations and consider other sources as well.
20CRv3 1836-2015. Uses pressure as input.
ftp ftp.cdc.noaa.gov
cd to
2m temperature: /Datasets/20thC_ReanV3/Dailies/2mMO/air.2m.yyyy.nc (for previous 1981 use 2mSI).
/Datasets/20thC_ReanV3/Monthlies/2mSI-MO/air.2m.mon.mean.nc
2m specific humidity: /Datasets/20thC_ReanV3/Dailies/2mMO/shum.2m.yyyy.nc (for presvios1981 use 2mSI).
/Datasets/20thC_ReanV3/Monthlies/2mSI-MO/shum2m.mon.mean.nc
2m relative humidity: /Datasets/20thC_ReanV3/Dailies/2mMO/rhum2m.yyyy.nc (for previous 1981 use 2mSI).
/Datasets/20thC_ReanV3/Monthlies/2mSI-MO/rhum.2m.mon.mean.nc
precipitation rate: : /Datasets/20thC_ReanV3/Dailies/sfcMO/air.2m.yyyy.nc (for previous 1981 use sfcSI).
/Datasets/20thC_ReanV3/Monthlies/sfcSI-MO/air.2m.mon.mean.nc
JRA-55: 1958-present
Get a login
ncftp -u username -p password ds.data.jma.go.jp
cd /JRA-55/Hist/Monthly/
download files: go to directory/
2mT: anl_surf125: var name TMP_GDS0_HTGL_S113
rhum:anl_surf125: R_GDS0_ISBL_S113"
2mq anl_surf125: var name PFH_GDS0_HTGL_S113
prate fcst_phy2m125.TPRAT_GDS0_SFC_S130
Files are grib Convert to netCDF
ERA20-C: 1900-2010?
2m temperature: Month
2m specific humidity
2m relative humidity
precipitation rate:
NCEP Reanalysis 1948-2020.
Updated to near present and starts earlier than most. Lower resolution than some. Usees an older model.
ftp ftp.cdc.noaa.gov
cd to
2m temperature: /Datasets/ncep.reanalysis.dailyavgs/air.2m.yyyy.nc (day) /Datasets/ncep.reanalysis.derived/air.2m.mon.mean.nc (month)
2m specific humidity: /Datasets/ncep.reanalysis.dailyavgs/air.2m.yyyy.nc (day) /Datasets/ncep.reanalysis.derived/air.2m.mon.mean.nc (month)
2m relative humidity: /Datasets/ncep.reanalysis.dailyavgs/air.2m.yyyy.nc (day) /Datasets/ncep.reanalysis.derived/air.2m.mon.mean.nc (month)
precipitation rate:: /Datasets/ncep.reanalysis.dailyavgs/air.2m.yyyy.nc (day) /Datasets/ncep.reanalysis.derived/air.2m.mon.mean.nc (month)
ERA-5
Monthly
Not sure how to get daily.,
Hello Cathy,
Supersaturation with respect to ice can occur and values can be greater than 150%. However, the values shouldn't be negative. Given that pressure level RH is stored in spherical harmonics, it could be due to the spectral transformations or maybe a spectral truncation before the transformation. From where did you get the data and how did you process it (if at all)? If you tell me a date, time, level and location, I will check with the native grid.
Regards,
Paul
We checked Jan 1 0z 1979. Values had the range over all levels.
r_min = -6.01934528350829 ;
r_max = 162.923529684544 ;
We did get the files from NCAR but the values from Copernicus were the same. Grid was netCDF
dimensions:
longitude = 1440 ;
latitude = 721 ;
level = 37 ;
time = 1 ;
short r(time, level, latitude, longitude) ;
r:scale_factor = 0.00257798170338687 ;
r:add_offset = 78.450803209666 ;
etc.
from Copernicus. Let us know if you need further info.
Cathy
Also, we checked our specific humidity and values were also below 0 for pressure level values.
I know that can be the result of post processing. If it isn't in the native grid, should we truncate either rhum or shum at 0? If it is, should we truncate? Thanks.
Cathy
So should we truncate at <0 for RH and SH?
Cathy
Hello Cathy,
Values of relative and specific humidity should not be below zero, so yes, you should set them to zero. Relative humidity is archived as spherical harmonics, so even though the model should not have negative values, the spectral transformations can introduce them because of Gibbs oscillations. Specific humidity is archived in grid point space on a reduced Gaussian grid, but you are using data on a regular latitude/longitude grid, so the interpolation may have introduced the negative values. However, I would have thought the interpolation scheme would be designed not to do that, so I will check that out.
Best wishes,
Paul
ncep-notes-questions-and-discussion
What does the NCEP Reanalysis model use as its tropopause definition?
Current reference
The ERA5 Global Reanalysis Hersbach, H. et al. May 2020. QJRMS
please cite the article as doi: 10.1002/qj.3803.
Other References
How to cite
Dataset citable as: Copernicus Climate Change Service (C3S) (2017): ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate . Copernicus Climate Change Service Climate Data Store (CDS), date of access. https://cds.climate.copernicus.eu/cdsapp#!/home
MERRA-2 references
First pick the correct variable here:
https://disc.gsfc.nasa.gov/datasets?page=1&project=MERRA-2&keywords=%22MERRA-2%22
When you click on the correct variable, it will take you to a second webpage with tabs that you can click that include: (1) documentation papers you need to cite, and (2) the correct variable citation information. (And you apparently need both types of citation). For example, below is a link to the page for 3d assimilated, 3-hourly, model level data (i.e. the equivalent of what everyone uses for 6-hourly ERA-Interim variables like temperature).
https://disc.gsfc.nasa.gov/datasets/M2I3NVASM_V5.12.4/summary?keywords=%22MERRA-2%22
Dear all,
I wanted to know where to find wind speed data for a specified location in Kosovo?
Hope to hear from you soon.
With respect,
Bukurije Hoxha.
Hello,
I have two questions
1) what variables are used to calculate total reconstructed PM2.5(inst1_2d_bot_Nx: Bottom Layer Diagnostics) produced by MERRA-2 GMI?
2) When calculate the total reconstructed PM2.5, what variables are weighted? how much they are weighted?
Thank you in advance
I can't speak to MERRA-2 GMI. Can you provide a link to where you found this, and also if there is a documentation page?
some MERRA-2 PM2.5 data can be found here
https://giovanni.gsfc.nasa.gov/giovanni/#service=TmAvMp&starttime=&endtime=&variableFacets=dataProductTimeInterval%3Amonthly%3B&dataKeyword=PM%202.5
Valid ranges
Some pressure level RH values are not within 0-100 range and some as as far off as -200%. Should we be setting these values to 0-200 or do the negative values represent supersaturation or something like that?