# Atmospheric variables below model terrain

**URL:** <https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225>\
**Category:** C3S - Datasets and Usage\
**Tags:** c3s, era5, global\_reanalysis\
**Created:** [30 September 2019 15:32 UTC](https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225 "2019-09-30T15:32:21Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![Michela](https://forum.ecmwf.int/letter_avatar_proxy/v4/letter/m/898d66/32.png) [@Michela](https://forum.ecmwf.int/u/Michela)\
**Post date:** [30 September 2019 15:32 UTC](https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225/1 "2019-09-30T15:32:21Z")

</div>

ERA5 reanalysis data on pressure levels at lower levels, eg at 1000 or 850 hPa, can show atmospheric variables that are below the model terrain. This happens in regions such as the Himalayas or Andes and is due to the interpolation of data from model to pressure levels&nbsp; (the interpolation documentations is at [http://www.umr-cnrm.fr/gmapdoc/spip.php?article157](http://www.umr-cnrm.fr/gmapdoc/spip.php?article157)).  
 The data which is below the surface can be masked out using the surface pressure or the [surface geopotential](https://confluence.ecmwf.int/display/CKB/ERA5%3A+surface+elevation+and+orography).  
 For more information about ERA5 levels see the [ERA5 documentation](https://confluence.ecmwf.int/pages/viewpage.action?pageId=76414402).

---

<div class="post-metadata">

**Author:** ![Henrique\_Barbosa](https://forum.ecmwf.int/letter_avatar_proxy/v4/letter/h/839c29/32.png) [@Henrique\_Barbosa](https://forum.ecmwf.int/u/Henrique_Barbosa)\
**Post date:** [2 March 2021 12:04 UTC](https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225/2 "2021-03-02T12:04:00Z")

</div>

The link you sent for the ERA5 documentation says the orography can be downloaded from the era5-single-levels dataset.&nbsp;

I'm confused, because that dataset varies with time. I've tried, and I can indeed download "hourly fields" of geopotential at the surface (aka topography). Shouldn't the orography be an invariant field, like in ERA-Interim for instance?

[https://apps.ecmwf.int/datasets/data/interim-full-invariant/](https://apps.ecmwf.int/datasets/data/interim-full-invariant/)

Or does the topography in ERA5 really changes?&nbsp;

I would be a good idea to clarify this point in the documentation.

---

<div class="post-metadata">

**Author:** ![Michela](https://forum.ecmwf.int/letter_avatar_proxy/v4/letter/m/898d66/32.png) [@Michela](https://forum.ecmwf.int/u/Michela)\
**Post date:** [3 March 2021 16:03 UTC](https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225/3 "2021-03-03T16:03:42Z")

</div>

Hi,

the surface geopotential (orography in the CDS) is listed in&nbsp; **Table 1: surface and single level parameters: invariants** at the following link:

[ERA5: data documentation - Copernicus Knowledge Base - ECMWF Confluence Wiki](https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation#ERA5:datadocumentation-Parameterlistings)&nbsp;

Please also have a look at the description for orography available on the overview page&nbsp;[ERA5 hourly data on single levels from 1979 to present (copernicus.eu)](https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview)&nbsp;where it states that the parameter does not vary in time.

Thanks

Michela

---

<div class="post-metadata">

**Author:** ![Felix\_Bracht](https://forum.ecmwf.int/letter_avatar_proxy/v4/letter/f/ec9cab/32.png) [@Felix\_Bracht](https://forum.ecmwf.int/u/Felix_Bracht)\
**Post date:** [13 October 2022 11:22 UTC](https://forum.ecmwf.int/t/atmospheric-variables-below-model-terrain/1225/4 "2022-10-13T11:22:01Z")

</div>

Hi [Michela Giusti](https://confluence.ecmwf.int/display/~usmg),

I am working with **_ERA5 reanalysis data on pressure levels_** and try to mask out&nbsp;atmospheric variables that are below the model terrain using the invariant geopotential from the era5-single-levels dataset.

However, I have realized that the geopotential in **_ERA5 reanalysis data on pressure levels&nbsp;_** varies over time.

Hence, I observe different values for the geopotential at different hours (and dates) for the same location (coordinates).

For instance, I observe the following values in my data:

| level | time | geopotential | longitude | latitude |
| 1 | 01.01.2020 10:00 | 458472 | -100 | 40 |
| 1 | 01.01.2020 22:00 | 460877 | -100 | 40 |
| 1 | 01.11.2020 10:00 | 458277 | -100 | 40 |
| 1 | 01.11.2020 22:00 | 460718 | -100 | 40 |
| ... | ... | ... | ... | ... |

Shouldn't the geopotential be time-invariant within the same pressure level and coordinates?

How can I mask out atmospheric variables below the model terrain based on changing geopotentials in the ERA5 reanalysis data on pressure levels data?
