Atmosphere

JRA-55: Notes, Questions, and Discussion

Created by mapr8844 on - Updated on 01/29/2020 12:54

JRA-55: Notes, Questions, and Discussion

Quick Overview and useful links

 

General Reanalysis Notes, Questions, and Discussion

01/29/2020: Some tropical cyclones occurring over the Northeast Pacific and the North Atlantic from 1959 to 1987 were erroneously represented as anti-cyclonic vortices in the Japanese 55-year Reanalysis (JRA-55) dataset. A detailed report is available at the JRA-55 website:



# Issue with tropical cyclone analysis in JRA-55

https://jra.kishou.go.jp/JRA-55/index_en.html#quality


Note: the following comments have been added to the body of this page due to them not migrating during our recent site upgrade:

 

JRA-55: 2016 Data

Submitted by Johanna Yepes (not verified) on Fri, 05/27/2016 - 08:30. 

Hi, Does the 2016 hourly data available? Thank you

Re: JRA-55: 2016 Data

Submitted by Kazutoshi Onogi (not verified) on Wed, 06/01/2016 - 01:53. 

Unfortunately hourly data is not available in JRA-55. Please visit http://jra.kishou.go.jp/JRA-55/index_en.html and read the documents "JRA-55 Product Users' Handbook" at the bottom of the page.

Re: JRA-55: 2016 Data

Submitted by Chiaki Kobayashi (not verified) on Mon, 05/30/2016 - 04:00. 

JRA data are available a couple of days behind the real time from JMA's web site (http://jra.kishou.go.jp/JRA-55/index_en.html) except for 3-dimensional average diagnostic fields and model level forecast fields. Thanks, Chiaki KOBAYASHI Climate Research Department Meteorological Research Institute Tsukuba, Ibaraki, 305-0052, Japan.

Re: JRA-55: 2016 Data

Submitted by gilbert.p.compo on Fri, 05/27/2016 - 13:32. 

It appears that the data are not yet at NCAR in http://rda.ucar.edu/datasets/ds628.0/ . I recommend contacting the Data Specialist mentioned on the page directly.

Gilbert P. Compo

University of Colorado/CIRES

NOAA/ESRL/Physical Sciences Division

http://www.esrl.noaa.gov/psd/people/gilbert.p.compo

JRA-55 assmilation of horizontal wind from land stations

Submitted by gilbert.p.compo on Thu, 05/19/2016 - 13:25. 

Does JRA-55C or JRA-55 assimilate horizontal winds from land surface stations? Thanks.

Gilbert P. Compo

University of Colorado/CIRES

NOAA/ESRL/Physical Sciences Division

http://www.esrl.noaa.gov/psd/people/gilbert.p.compo

Re: JRA-55 assmilation of horizontal wind from land stations

Submitted by Chiaki Kobayashi (not verified) on Fri, 05/20/2016 - 07:47. 

Dear Gil Horizontal winds from land surface stations are assimilated only in the screen-level analysis component, not in the atmospheric analysis component. Screen-level analysis fields are not used as the initial conditions for forecasts, therefore have no impact on subsequent assimilation cycles. Best regards, Chiaki -- Chiaki KOBAYASHI Climate Research Department Meteorological Research Institute Tsukuba, Ibaraki, 305-0052, Japan.

JRA-55C assimilation of horizontal wind from aircraft

Submitted by gilbert.p.compo on Thu, 05/19/2016 - 13:24. 

Does JRA-55C assimilate horizontal wind observations from aircraft? I do not think so based on Kobayashi et al. 2014 but wanted to double check. Thanks.

Gilbert P. Compo

University of Colorado/CIRES

NOAA/ESRL/Physical Sciences Division

http://www.esrl.noaa.gov/psd/people/gilbert.p.compo

Re: JRA-55C assimilation of horizontal wind from aircraft

Submitted by Chiaki Kobayashi (not verified) on Fri, 05/20/2016 - 07:50. 

Dear Gil, No, The observation from air crafts are not assimilated in JRA-55C. I think you are understanding correctly. Best regards, Chiaki -- Chiaki KOBAYASHI Climate Research Department Meteorological Research Institute Tsukuba, Ibaraki, 305-0052, Japan.

nenghan wan (not verified)

Mon, 06/06/2022 - 22:40

Hi, Sir,

From previous comments, I konw Brightness temperature is the SST, but I'd like to know if the SST over Arctic and Antarctica are the water temperature or sea ice temperature, the trend is so different when I compare the monthly time series of SST for ERA5 and JRA55 over Arctic and Antarctica, as I understand the SST in EAR5 is for water temperature. Would you mind provide some references? Thank you for your help, I really appreciate.

best,

nenghan

Do you mean references for JRA55?

https://climate.mri-jma.go.jp/pub/ocean/JRA55-do/docs/v1_5-manual/User_…

As I understand it, the COBE2-SST is the model SST. COBE2 references are

  • Hirahara, S., Ishii, M., and Y. Fukuda, 2014: Centennial-scale sea surface temperature analysis and its uncertainty. J of Climate27, 57-75.
  • Folland, C. K. and D. E. Parker, 1995: Correction of instrumental biases in historical sea surface temperature data. Q. J. R. Meteorol. Soc., 121, 319-367.
  • Ishii, M., A. Shouji, S. Sugimoto, and T. Matsumoto, 2005: Objective Analyses of Sea-Surface Temperature and Marine Meteorological Variables for the 20th Century using ICOADS and the Kobe Collection. Int. J. Climatol., 25, 865-879.
  • Japan Meteorological Agency, 2006: Characteristics of Global Sea Surface Temperature Analysis Data (COBE-SST) for Climate Use. Monthly Report on Climate System Separated Volume, 12, 116pp.

 

Cathy Smith

Dear Aaron,

Here are some papers that intercompare upper tropospheric winds field from different reanalyses, though the latest reanalyses such as ERA5 are not covered.

 

Long et al., 2017: Climatology and interannual variability of dynamic variables in multiple reanalyses evaluated by the SPARC Reanalysis Intercomparison Project (S-RIP). Atmos. Chem. Phys., 17, 14593–14629, https://doi.org/10.5194/acp-17-14593-2017.

Manney et al., 2017: Reanalysis comparisons of upper tropospheric–lower stratospheric jets and multiple tropopauses. Atmos. Chem. Phys., 17, 11541–11566, https://doi.org/10.5194/acp-17-11541-2017.

Manney, and Hegglin, 2018: Seasonal and regional variations of long-term changes in upper-tropospheric jets from reanalyses. J. Climate, 31, 423-448, https://doi.org/10.1175/JCLI-D-17-0303.1.

 

I hope this helps.

 

Regards,

Shinya

 

Dear Aaron,

To my knowedge, no detailed reanalysis intercomparison of upper-tropospheric/lower-stratospheric winds similar to what Long et al. and manney et al. have done has been published yet that includes ERA5, unfortunately. The report from the SPARC Reanalysis Intercomparison Project (S-RIP), which should be published this fall and whose Chap 3 provided the basis for the Long et al. paper, may include some evaluation of ERA5 winds (but this is still to be confirmed). In addition, some insight into the skill of ERA5 upper-level winds can be found in Hersbach et al. (2020) and - with focus on the QBO - in Simmons et al. (2020).

Hersbach, H., et al., 2020: The ERA5 global reanalysis. Q.J.R. Meteorol Soc., 146, 1999-2049, https://doi.org/10.1002/qj.3803.

Simmons, A., et al., 2020: Global stratospheric temperature bias and other stratospheric aspects of ERA5 and ERA5.1. ECMWF Technical Memorandum number 859, https://doi.org/10.21957/rcxqfmg0.

I hope this helps.

Kind regards,
Julien
(ERA5 Team)

Jih-Wang Aaron Wang (not verified)

Thu, 05/02/2019 - 11:03

Dear All:

Sorry, I made a mistake in the previous post.

I don't understand why JRA55's resolution is ~55km. If it uses T319 spectral model, isn't the horizontal resolution=[earth circumference]/[320*2]=40000km/640=62.5km.

Thanks,

Aaron

Dear Aaron,

I'm sorry for the delay in replying to you.

Your calculation is basically correct for the grid intervals in the tropics and for those in the north-south direction. The value of ~55km represents the grid intervals in the west-east direction in the mid-latitudes since they are narrower in higher latitudes.

Regards,

Shinya 

Cathy.Smith@noaa.gov

Mon, 02/25/2019 - 11:00

Brightness temperature is actually SST. This SST comes from the COBE daily data, but is not exactly the same because it includes a parameterized diurnal cycle. So it differs from the COBE data in terms of the temporal variability as well as the resolution.

This brightness temperature is computed from surface upward longwave radiation assuming that the surface is a black body. As long as this assumption is valid, the brightness temperature can be used as skin temperature as well.

David Carvalho (not verified)

Tue, 03/27/2018 - 15:33

Hi, I need to know what variables (T, UV winds, surface pressure, etc.) JRA55 assimilates from land met stations and moored buoy arrays such as TAO, TRITON and RAMA. I saw on this forum that it doesn't assimilate winds from land met stations, if I understood correctly.

Also, does JRA55 assimilates any variables measured by moored buoys offshore the Spanish Atlantic coast, operated by the Spanish Agency Puertos del Estado?

Any help would be much appreciated, thanks!

Cathy.Smith@noaa.gov

Mon, 10/24/2016 - 10:09

Users should be aware that for 2mT, 2mq, 10m u and v and some others, that they should use the "forecast" variable and not the "analysis" variable. The WRIT pages now have both.

----------------------------------------------------------------------------------------------------------

Table 5 and section 3.2a "screen level analysis" of

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. Meteorol. Soc. Japan, 93, 5-48, doi: 10.2151/jmsj.2015-001

describes the details of this "separate" analysis:

The analysis of screen-level variables (2-m temperatures, 2-m relative humidities and 10-m winds) is performed separately from the atmospheric analysis component. These variables are analyzed with a univariate 2-dimensional optimal interpolation (2D-OI). In JRA-25, departures were estimated by comparing observations with the first guess at the analysis time. In JRA-55, observations are compared with the first guess at the actual observation time and the departures are applied at the analysis time (the so-called first guess at the appropriate time (FGAT) approach). Temperature and wind observations from islands are not used because they are not necessarily representative at the grid scale of JRA-55. Determining whether an observation is from an island is based on the 0.25-degree resolution land cover data; consequently, observations from the coast are also excluded. Screen-level analysis fields are not used as initial conditions for forecasts.

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

ERA-20C Notes, Questions, and Discussion

Created by mapr8844 on - Updated on 01/22/2019 17:46

ERA-20C: Notes, Questions, and Discussion

Quick overview and useful links

Direct dataset questions to ECMWF User Support https://climate.copernicus.eu/help-support

General Reanalysis Notes, Questions, and Discussion

Greetings

I read a recent reply to Yiyi's question (10/2017 comment #5401) regarding ERA-20C by Paul Berrisford, "For most applications you would use the monthly means of daily means. You would only use the synoptic means if you wanted to restrict your attention to a particular point in, or portion of, the diurnal cycle."

Can a month of monthly means (forecasts) and a month of synoptic means (forecasts) can be used interchangeably? Are there applications where one might be preferred? I assume that the comparison must be MODA forecast to SYN forecast, and AN to AN, if the forecasts contain only model output and the analysis contains assimilation?

-thanks in advance

Good moring!

We are working with 20CR and ERA20C on Shclimate dynamics. We observe some important differences prior to 1950 between both reanalysis. We would like to know when and where ship log data  is included in ERA20C between 1905 and 1950 south of 30S.

Thanks!

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

ERA-Interim: Notes, Questions, and Discussion

Created by mapr8844 on - Updated on 09/01/2020 19:08

 

 

 

Direct dataset questions to ECMWF User Support https://climate.copernicus.eu/help-support

Note: the following three comments have been added to this page due to them not migrating during our upgrade:

 

Re: ERA-Interim: Notes, Questions, and Discussion

Submitted by Natalie T (not verified) on Tue, 07/12/2016 - 13:14. 

Dear all, Is acceleration due to gravity already accounted for in the vertical integral of northward/eastward water vapour flux? I didn't find mention of anything of the sort in the documentation, but didn't want to repeat the action in my analysis if it has already been done. Thank you in advance.

Re: ERA-Interim: Notes, Questions, and Discussion

Submitted by Namaoui (not verified) on Tue, 07/12/2016 - 01:58. 

Hello, Iam Beginner in Era interim , I want to know how download a frid of Era interim ? Regards

ERA-Interim: Questions

Submitted by Johanna Yepes (not verified) on Fri, 05/13/2016 - 09:33. 

Dear Era-Interim staff, I would like to ask about the units of the runoff data. I downloaded the monthly means of runoff and I can see that the units are meters but I don't know the time units for that, those are meters per day or meters per month? Thank you.

YEONWOO CHOI (not verified)

Thu, 10/22/2020 - 14:39

Dear, Sir/Madam

I have question about the surface downward UV radiation data (i.e., UV-B).

I found that the range of UV values extracted from ERA-INT and ERA5 is different, even though their units are the same (W/m2).

For example, in Jan 1979, global UV-B derived from ERA5 ranges from 0 to 45, but that derived from ERA-INT ranges from 20 to 130.

What is the exact difference between ERA-INT and ERA5 data on the value of the UVB?

 

Best

YW 

 

 

DearYW,

thank you for your interest in ERA5 and ERA-Interim.

The difference between ERA5 and ERA-Interim is that for fluxes (and other accumulated quantities) is that for ERA5 accumulations are since the previous post-processing step (hourly for ERA5 and 3-hourly for its uncertainty estimate) while it is since from the start of the originating forecast for ERA-Interim.

More info can be found at: 

https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation#ERA5:datadocumentation-Meanratesandaccumulations

SUCHISMITA CHOUDHURY

Wed, 04/22/2020 - 09:18

How do we process ERA-Interim data in Matlab and obtain the plot for any of the parameters. Have downloaded the NetCDF file and was cable of reading the file but unable to plot the data. I am a Research Student working on Trend analysis of Total column Water Vapor.

Hello,

ERA-Interim provides the parameter "Total column water vapour" (code 137.128). In addition, there is the parameter "Total column water" (code 136.128). Both of these parameters can be downloaded from the webpage https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/ at sub-daily temporal resolution. They are also available as monthly means, see the menu option "Monthly means of daily means" on the left hand side of the above page.

Regards,

Paul

lais schmalfuss (not verified)

Fri, 12/07/2018 - 06:33

Hello,

Somebody could help me with this?

I downloaded data from the Interim ERA and forgot to rename it.
Now I know the data years, but I do not know which months each data represents.

Thank's

Laís

Dear Lais,

If you retrieved GRIB data, the date information (or century, year, month, day), time and length of the forecast, for forecast data, is contained within the GRIB header of each message. If you retrieved netCDF data, one of the coordinates of the file, is time since the reference date/time.

Regards,

Paul

Cathy.Smith@noaa.gov

Thu, 07/12/2018 - 15:21

I used to get the last date of the monthly ERA-Interim reanalysis by going to a ERAInterim monthly webpage and noting the last month listed. Now I don't see a last month. How can I find the last month available? Is there a page that has it? Tx.

Dear Cathy,

The page you are looking for is the page for daily data:

https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/

The data are updated on a monthly basis, so if the daily data is available for a month, then the monthly mean data will be available too.

Regards,

Paul

I see the last date again. But I have had no luck getting the (monthly) data via a python script for the last 2 months or so. The code was working and suddenly just hangs. I suspect it's an overloaded server. I was able to get a file via the web interface though it took a while. That's not a great option for me but I could try to make it work. 

Hello,

We have looked at this and we can download data. Could you give us more details about what you are trying to obtain and how?

The easiest way to do this might be to create an issue on the C3S website: go to https://climate.copernicus.eu/help-and-support, click "login to the C3S Enquiry Portal", then "Any question" and fill in the form.

When the problem is resolved, we can post a summary here.

Thanks,

Paul

Kyle Itterly (not verified)

Fri, 05/25/2018 - 12:01

Hello,

I have a question and a quick note to those looking at the diurnal cycle of accumulated forecast parameters in general.

Let's use total precipitation as the example. I seem to be getting a few negative values for 3-hourly precipitation rate when subtracting the 3-hourly accumulated steps and dividing by the interval. It turns out that it's relatively rare and the non-negative values appear to match other results and to be of sufficient accuracy. Applying this same methodology to other accumulated parameters e.g. radiative fluxes appears to be fine. My theory is that bilinearly interpolating the accumulated precipitation (before subtraction of steps) to 1 degree by 1 degree (done by the WebAPI server) is the core of my issue. So, what would be the proper resolution and order of computation to ensure the 3-hourly precipitation rate is being calculated as intended?

My note to others looking at the diurnal cycle of ERA forecast parameters is to consider using the WebAPI so you can concatenate later forecast steps 6-18 hours from 00 and 12 UTC instead of using the first few steps 3-12, which can be affected by model spinup errors. Additionally, using forecast steps 6-18 hours provides a full diurnal cycle of model CAPE without erroneous values of zero in forecast step 3. This took me a while to figure out, and it turns out that the diurnal cycle of some parameters (e.g. OLR) appears to be less biased using the later steps.

Thanks in advance for any insight,

Kyle

David Carvalho (not verified)

Tue, 03/27/2018 - 15:31

Hi, I need to know what variables (T, UV winds, surface pressure, etc.) ERA-Interim assimilates from land met stations and moored buoy arrays such as TAO, TRITON and RAMA.

Also, does ERA-Interim assimilates any variables measured by moored buoys offshore the Spanish Atlantic coast, operated by the Spanish Agency Puertos del Estado?

Any help would be much appreciated, thanks!

Paul Berrisford (not verified)

Mon, 09/04/2017 - 11:37

Dear Andrew,

From the data, you could calculate the geopotential height on the model levels and interpolate the winds to 100gpm. This would give you the most accurate answer. For calculating geopotential, see equation 2.21 in Part III "Dynamics and numerical procedures" of the ERA-Interim model documentation:

https://www.ecmwf.int/sites/default/files/elibrary/2007/9220-part-iii-dynamics-and-numerical-procedures.pdf

To a first approximation the near surface model levels in ERA-Interim are at a constant height above the surface. However, that approximate height is proportional to temperature (in K), so the standard atmosphere gives slightly different results to those assuming a uniform scale height, see:

https://www.ecmwf.int/en/forecasts/documentation-and-support/60-model-levels

With the standard atmosphere, 100m is located between levels 58 and 57, though slightly closer to 57. So you could assume the standard atmosphere and interpolate between these two levels.

If you're happy with the discrepancy due to temperature, you could just use the winds at level 57, following the result from a 7km scale height.

Regards,

Paul

 

Hi Paul,

Thanks for your reply, and sorry it's been awhile. I ended up doing some sensitivity tests for different model levels around 100m on the annual time scale (yearly means) for the geographic area of Texas. From what I can tell, there is no discernible difference between the different methods (using model level 57 versus interpolation of closest model levels). 

Thank you for your help.

Best,

Andrew

Andrew Kumler (not verified)

Tue, 08/22/2017 - 16:42

Hello,

I am currently doing research using ERA-Interim on hub-height winds (100 meters). ERA-Interim outputs data at model levels, one of these being model level 57 (1000 mb, or 100 m agl). This assumes a surface pressure of 1013.25 mb and scale height of the atmosphere of 7 km. Would this be an appropriate level to use to assess 100 meter winds over a large geographic area, or should alternative levels be used? The wind power law can be used to calculate 100 meter winds using winds at a lower heights as well. Or should I use geopotential heights to account for changes in pressure at 100 m agl? I assume changes in geopotential at 100 m would be fairly negligible over long periods of time.

Thanks,

Andrew

Dear Cathy,

as Paul said, there is no automatic way of checking the latest available date. The ERA-Interim release cycle is described here: https://software.ecmwf.int/wiki/display/CKB/Why+can%27t+I+download+ERA-Interim+data+for+recent+months

Regards

Karl

Cathy.Smith@noaa.gov

Wed, 04/12/2017 - 08:52

How can I automatically find the last date of the ERA-inerim? If I request a date past the last date, the ecmwf perl module  fails with many errors. I can look at the webpage but then I can't update automatically.

 

Dear Cathy,

There is no automatic way of checking the latest available date. In order to find this date, you need to look at the ERA-Interim data webpage:

http://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/

Regards,

Paul

Lu Sun (not verified)

Thu, 11/17/2016 - 12:05

Dear Sir/Madam,

I am wondering how the CAPE in ERA-Interim data be calculated? Because I want to compared the results with other reanalysis data but most of them need us to calculate CAPE individually. Some important information such as: what level should be treated us the beginning point (1000mb? 1017mb? or the surface pressure?) Thank you for your replying.

 

Dear Lu Sun,

There is a description of the definition of CAPE in the ERA-Interim model documentation, Part IV: Physical processes:

http://www.ecmwf.int/sites/default/files/elibrary/2007/9221-part-iv-physical-processes.pdf

In particular, look at equation 5.61 in Section 5.11 Diagnostics for post processing: CAPE, which is in Chapter 5 Convection.

The calculation is done from model (hybrid sigma/pressure) level data.

Regards,

Paul

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

CFSR: Notes, Questions, and Discussion

Created by mapr8844 on - Updated on 05/31/2017 13:43

CFSR: Notes, Questions, and Discussion

Quick overview and useful links

 

General Reanalysis Notes, Questions, and Discussion


Notice: the following two comments have been placed on the body of this page as they didn't migrate during the site upgrade:

1.) 

Re: CFSR: Notes, Questions, and Discussion

Submitted by Anonymous (not verified) on Thu, 05/12/2016 - 18:09. 

What is the Ocean data input for CFSR after 1997 which seems the data volume are more than ICOADS 2.5? cited by SURANJANA SAHA et al., 2010: All of the observed marine data from 1948 through 1997 have been supplied by the COADS datasets.

2.) 

Re: CFSR: Notes, Questions, and Discussion

Submitted by gilbert.p.compo on Mon, 05/16/2016 - 11:26. 

Courtesy of David Behringer: ---------- Forwarded message ---------- From: David Behringer - NOAA Affiliate Date: Mon, May 16, 2016 at 9:37 AM Subject: Re: reanalyses.org question about ocean input data in CFSR Hi Gil, After 1997 we pulled in data from the NODC data base. Beginning in 2000 we merged profiles from the GODAE and NODC databases. We found the merged set was larger than either alone. We've also done our own QC. Dave

Gilbert P. Compo

University of Colorado/CIRES

NOAA/ESRL/Physical Sciences Division

 

Lilia Caballero (not verified)

Wed, 01/27/2021 - 00:00

Hi. I would like to know which type of cartographic projection does CFSR use?

Wesley Ebisuzaki (not verified)

Thu, 01/28/2021 - 09:12

In reply to by Lilia Caballero (not verified)

   The CFSR uses a spectral model of a spherical Earth for the atmospheric fields.  The

post-processor converts the spectral fields to the (global) gridded fields for the users.

The global fields can be converted to various projections like Mercator, polar stereographic,

and Lambert conical by using programs like copygb, copygb2, and wgrib2. 

   The land-surface and ocean models are grid models and use different grids.

David (not verified)

Tue, 03/27/2018 - 15:36

Hi, I need to know what variables (T, UV winds, surface pressure, etc.) CFSR assimilates from land met stations and moored buoy arrays such as TAO, TRITON and RAMA.

Also, does CFSR assimilates any variables measured by moored buoys offshore the Spanish Atlantic coast, operated by the Spanish Agency Puertos del Estado?

Any help would be much appreciated, thanks!

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

Arctic System Reanalysis: Notes, Questions, and Discussion

Created by mapr8844 on - Updated on 05/31/2017 13:42

Arctic System Reanalysis: Notes, Questions, and Discussion
 

Quick overview and useful links

 

General Reanalysis Notes, Questions, and Discussion

George Victor (not verified)

Fri, 03/12/2021 - 09:48

Hi,

I'm George Victor, a student at Uppsala University, Sweden. I'm interested in using the arctic system reanalysis (15km) dataset for wave model simulations.

For my study, I expect to have a regular rectangular grid data. However, since ASR data is not based on regular rectangular projection, I'd like to know what projection is it based on. I'd also be grateful if you could suggest me some tools to convert the  data to rectrangular  projection. Thank you

The information below deals with the projection question. Regret I cannot assist with the tool query. If George Victor contacts us directly and provides more information, we may be able to direct him to the necessary software.

Dave Bromwich, Ohio State University. bromwich.1@osu.edu

The WPS parameter for ASRv2 that uses the Advanced Research WRF (ARW):

    The projection is polar stereographic.

    The true latitude is 60N.

    The domain center latitude is 90N and standard longitude is -175E.

 REF_LAT= 90N : A real value specifying the latitude part of a (latitude, longitude) location whose (i,j) location in the simulation domain is known. For ARW, ref_lat gives the latitude of the center-point of the coarse domain by default (i.e., when ref_x and ref_y are not specified).

REF_LON=-175E : A real value specifying the longitude part of a (latitude, longitude) location whose (i, j) location in the simulation domain is known. For ARW, ref_lon gives the longitude of the center-point of the coarse domain by default (i.e., when ref_x and ref_y are not specified). West longitudes are negative, and the value of ref_lon should be in the range [-180, 180]. 

STAND_LON=-175E : A real value specifying, for ARW, the longitude that is parallel with the y-axis in the Lambert conformal and polar stereographic projections. For the regular latitude-longitude projection, this value gives the rotation about the earth's geographic poles. 

The WRF Preprocessing System (WPS) document:

https://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.0/users_guide_chap3.html

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

Notes, Questions and Discussion

Created by esrl_admin on - Updated on 03/29/2022 08:53

Users have requested workspaces and comment areas for general reanalysis questions and for asking specific questions or posting comments about specific reanalysis datasets.

The links below are good places to post notes, ask questions, comment, or start a discussion about each dataset.

Or, leave a comment at the bottom of this page. 

Links to Notes, Questions, and Discussions for each reanalysis dataset:

Reanalyses.org Overview Comparison Table

 

Reanalyses.org Overview of each reanalysis

 

SPARC-Reanalysis Intercomparison Project (S-RIP) Comparison Paper and its Tables 

  • the paper includes several tables comparing assimilating models, resolution, major physical parameterizations, data input, and boundary conditions

Bennet Foli (not verified)

Sat, 01/02/2021 - 11:01

Hello!

I would like to know if PIRATA buoy data are assimilated in NCEP and ERA5 reanalysis winds and if so which specific buoys (locations) and which parameters from the buoys are utilized. Can I also have some references to these.

Thank you.

Thank you for the reply.

The following are the station IDs and corresponding locations

GL_TS_MO_13010; % Lat 0N, Lon 0E
GL_TS_MO_15002; % Lat 0N, Lon 10W
GL_TS_MO_31006; % Lat 4N, Lon 23W
GL_TS_MO_13002; % Lat 20N, Lon 23W
 

Thank you.

Bennet

I tried to estimate the 2m horizontal wind using the log wind profile with 10 m horizontal wind, and found it is quite different to the 2m horizontal wind provided in MERRA2 output. So, I am wondering how did you calculate the horizontal wind at 2m? Thanks.

Posted from Arlindo da Silva and Andrea Molod, both of NASA

 

From Dr. da Sliva-

From the glossary: U10M, U2M, V10M, V2M: The eastward and northward wind component at 10 m, and 2m in the surface layer, in m s -1 . These values are defined above displacement height (DISPH) and are dependent on the surface layer parameterization being used.

 

From Dr. Molod

hi there,

MERRA-2 used what we call the Helfand (Helfand and Schubert, 1995) Monin-Obukhov surface layer parameterization. The charnock coefficients (relating stress to speed over the ocean) were modified substantially from what is in the manuscript....  The 2-m and 10-m winds are “picked off” the non-dimensional vertical wind profiles that the local stability dictates. It will not use a neutral (log wind) profile.

Our surface layer thickness in MERRA-2 is approximately 60m, and so in M2 the 2-m and 10-m winds we output are essentially interpolated along the correct profile between 0 at the sfc and the model lowest level wind at 60 m.

If you are interested in repeating this calculation, it would be possible to approximate the profile using the Cu, u* and Ri outputs from MERRA-2 and getting the Psi values (non-dimensional vertical wind profile) from those.

I can send you the algebra or you can find it in the Helfand and Schubert paper in the appendix....

Let me know if there is any more information you would like.

andrea

Hi,

Recently I'd like to use MERRA2 data to do some analysis. I read some papers and knew that most output variables were calculated by

parameterizations in the GEOS-5 AGCM model it used and assimilated observations. 

May I ask how they made the output variables physically consistent? there are so many parameterizations that seem not related together.

 

Best,

Orrin

Hi Orrin,

The parameterizations in any GCM or atmospheric model generally depend on each other. There are prognostic state fields: T, Q, U, V (sometime liquid and ice water too) and say aerosols and chemical constituents. clouds are formed depending on the states, radiation is computed from the states and the clouds, and each then feed back into the states.  Land and ocean surfaces connect with the atmosphere, exchanging information through fluxes.  This is how the models are designed.  The code can be complicated in places. Best is to have an advisor who has experience with models and work with one. (without knowing your background)

Mike

Richard Yablonsky (not verified)

Fri, 05/12/2017 - 12:23

Hello... we are attempting to simulate the coastal flooding in Great Britain from the 1953 North Sea flood, and we are having trouble finding an atmospheric reanalysis product that provides reasonably accurate surface wind and pressure forcing for our simulations.  Our time period of interest is January 30, 1953 to February 2, 1953.  Any advice on which reanalysis product(s) would be most appropriate for our purposes would be much appreciated.

Dear Richard,

Which datasets have you already tried in this context? 

There are only 5 possibilities: 

20CR v2

20CR v2c

ERA-20C

CERA-20C

NCEP-NCAR 

Links to all of these are available at http://reanalyses.org/atmosphere/overview-current-atmospheric-reanalyses. 

For the ensemble-based 20CR v2 and v2c, you should use individual members, though using the ensemble mean for 1953 Great Britain should be alright.  CERA-20C also has an ensemble, and individual members would probably be better for your purposes. 

best wishes,

gil

20CR2c has a decent qualitative representation of the storm - https://vimeo.com/165142838 - I've not looked at the quantitative accuracy. I think David Hein, at the Met Office has has a regional model run covering the UK downscaled from one of the 20CR ensemble members. That might also be worth a look. 

Cathy.Smith@noaa.gov

Fri, 04/07/2017 - 09:57

The WRIT page is a good start.

For specific countries, The World Bank has data by country but not to 2015. You might contact them and see if they have it available.

http://data.worldbank.org/data-catalog/cckp_historical_data

They use CRU data.

You could also look at the URL below  (they have different years/months available). NCEI has many datasets and could have it lesehwre. You can look/ask them.

https://www.ncdc.noaa.gov/sotc/global/201513


good luck!

Cathy

Tracy Carty (not verified)

Thu, 04/06/2017 - 08:26

Hi

I work for a leading development NGO.  I'm writing a paper on the East Africa drought and I'm trying to get the following information, or close to it based on what's available.  I'm not up to speed with the data sets or how to use them and hence this message is a call to see if anyone would be able to lend me a hand.  I'm hoping this is a quick job for someone who has the data to hand and knows what they are doing??

[Mean annual temperatures increased from 1960-2016 by X in Kenya, X in Ethiopia and X in Somalia][would also be interested for 2015 and first few months of 2017 if available]

FYI - I've also emailed merra-questions@lists.nasa.gov to see if anyone can help with this.

Best wishes,

Tracy

 

Dear Tracy,

If you go to the Web-based reanalysis intercomparison tools https://reanalyses.org/atmosphere/writ and select the WRIT Timeseries button, you can make these calculations yourself for regions that encompass your countries of interest. The JRA-55 Screen dataset is probably the best for your purposes, as it includes actual near-surface thermometer values and is fairly up to date. The JRA-55 itself does not directly include the near surface temperature observations (it does include the upper-air, which start near the surface). 

Enter your latitude and longitude box approximating each country and you will get a time series that includes statistics about how much the temperature has changed.

You can also use the WRIT Maps and zoom in on Africa, and show the difference between the current period and an earlier period.

Hope that this is helpful.

Best wishes,

gil

 

Ashwin (not verified)

Thu, 09/01/2016 - 02:31

Does anyone have any code samples that are newer than this one - http://www.esrl.noaa.gov/psd/data/gridded/readgeneral.F ?

 

I want to be able to read a reanalysis 2 file using fortran and netcdf .

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

steps for data downloads

Created by Cathy.Smith@noaa.gov on - Updated on 03/07/2022 18:48

How to download files: Step-by-Step

These methods may not  be the only way! See the source for more information. Other sites may also have the data (NCAR, for example) and those methods would be different than those below. 

MERRA-2 (Using ftp subsetter)

  1. ​Go to page http://disc.sci.gsfc.nasa.gov/daac-bin/FTPSubset2.pl
  2. From the list, find type of data you are looking for. 3D, daily, monthly, etc. Note assimilation is term used for what many researchers know as analysis.
  3. Optionally select a region to subset.
  4. Select​ a start and end month. You cannot get all years so select a small subset to start with.
  5. Select a variable(s)
  6. Optional- select a remapping type and grid to specify grid size.  Some match other datasets. Some may not be obvious and you may need to try them and see what you get.
  7. Select an output format. netcdf4-classic can be read by many applications.
  8. Select "start search". You will get a list of files. You can save this list and run wget --content-disposition -i wget_n3LWiFa9
  9. You can also download a file and look at it to see if it's what you need.
  10. To get all files, you can run a script to get each file with a wget -content-disposition command and fill in the dates you need. Or you can use an editor to generate the list, changing the filename dates.
  11. Use name/password: You need to use a name and password to get the data. wget --user username --password user password
  12. Check to see that the files are downloaded completely.

A few notes:

  • There are different streams so you have to change those depending on the year.
  • You can change the output file names to whatever you want. 
  • You must close the window of the subsetted data search page in order to perform a new search!

JRA-55:

This out out of date. NCAR has the data there. 

  1. Register for dataset. You will get a user name and passwordhttp://jra.kishou.go.jp/JRA-55/index_en.html#usage
  2. Get the data via ftp with your name/password.  See http://jra.kishou.go.jp/JRA-55/index_en.html#download
  3. From the command line, ncftpls -u name -p password ftp://ds.data.jma.go.jp/JRA-55/Hist/Monthly/anl_p125 for example. Or use ncftpget
  4. Data is divided into directories: http://jra.kishou.go.jp/JRA-55/document/JRA-55_handbook_LL125_en.pdf

ERA-Interim:

  1. Go to this webpage: http://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/
  2. Login before getting data. You will need to apply for a login. Safari may have issues with ERA pages on MACS; Firefox will work.
  3. On the left column of the webpage above, select the type of analysis you need (time scale, data type).
  4. Select a month, variable. View MARS request. You will see a python script sample. That script will list the parameters you need to download the data files from a script. The parameters can be used by python or Perl scripts or similar languages.
  5. Download required modules depending on what language you are using. See this webpage https://software.ecmwf.int/wiki/display/WEBAPI/Access+ECMWF+Public+Datasets
  6. Edit your script to download the complete set of dates/variables you need. Include the required modules in your script. 

20CRV2c:

  1. ftp to ftp.cdc.noaa.gov
  2. cd /Datasets/20CRV2c/
  3. get files.  ncftp or similar command line client is preferable as the whole file is more likely to download than using a browser.  Data is divided into directories by time scales and type.
  4. For the browser link, use https://downloads.psl.noaa.gov/Datasets/20thC_ReanV2c/

20CRV3:

  1. ftp to ftp.cdc.noaa.gov
  2. cd /Datasets/20CRV3/
  3. get files.  ncftp or similar command line client is preferable as the whole file is more likely to download than using a browser.  Data is divided into directories by time scales and type.
  4. For the browser link, use https://downloads.psl.noaa.gov/Datasets/20thC_ReanV3/

 

 

Jeannine St-Jacques (not verified)

Fri, 06/23/2017 - 11:04

Hello

 

  I'm trying to find some web-based reanalysis tools that I used to use. I need to look at average 500 mb GPH for a few months and do some correlation plots. There used to be a webtool here that I used but there has been a big reorganization.

 

Thanks

Jeannine

Dear Jeannine,

What tools are you not able to find? Perhaps what you are looking for is linked at http://reanalyses.org/atmosphere/how-obtainplotanalyze-data ?

Best wishes,

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

20CRv2c: Notes, Questions, and Comments

Created by Cathy.Smith@noaa.gov on - Updated on 07/18/2016 10:13

20CRv2c: Notes, Questions, and Comments

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

ERA-20C references

Created by paul.poli@shom.fr on - Updated on 10/28/2016 14:21

The following references describe system components:

  • Initial, ensemble production, from which background errors were derived:

    Poli, P, Hersbach, H, Tan, DGH, Dee, DP, Thepaut, J-J, Simmons, A, Peubey, C, Laloyaux, P, Komori, T, Berrisford, P, Dragani, R, Trémolet, Y, Hólm, EV, Bonavita, M, Isaksen, L, Fisher, M. (2013) The data assimilation system and initial performance evaluation of the ECMWF pilot reanalysis of the 20th-century assimilating surface observations only (ERA-20C). ERA Report Series 14
  • Observation feedback:

    Poli, P, Hersbach, H, Berrisford, P, Dee, DP. (2015) The observation feedback archive for the ICOADS and ISPD data sets. ERA Report Series 18

 

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

Reanalysis Model Output: Similarity of format and attributes?

Created by Cathy.Smith@noaa.gov on - Updated on 07/18/2016 10:13

Background: Model output can differ in many ways including

  • variables
  • attributes
  • units
  • levels
  • file divisions for time-series (single time step in a file,  single day, single month...)
  • output grid resolution
  • file names
  • time format in netcDF (unlimited vs not unlimited)

Is it reasonable to suggest that any of these be standardized so that comparisons are easier. So that tools that analyze the data are easier to use?

gerald.potter

Wed, 05/27/2015 - 00:04

The NASA Climate Model Data Service (CDS) and the NASA Center for Climate Simulation (NCCS) has initiated a new project is called the Collaborative REAnalysis Technical Environment – Intercomparison Project (CREATE-IP). The primary object of this project is to convert the existing reanalyses to a format consistent with the CMIP5 directives which includes standardized variable names, attributes, file names and time format. The data is distributed though the Earth System Grid Federation (ESGF) in netCDF CF compliant format at: https://www.earthsystemcog.org/projects/create-ip/. Each variable is stored separately and contains all time values. The horizontal resolution of each reanalysis is unchanged from the way it was produced but guidance is provided to assist with modification and all available model levels are included. Currently, the monthly variables provided are those with a one-to-one correspondence to the CMIP5 archive. High frequency output is now being processed to provide 1 hour, 3 hour and 6 hour data for a select set of 2D and 3D variables. In addition, plans are underway to publish the O-A, O-F terms and selected observations using the same formatting.

Add new comment

CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.