Soil Moisture Anomalies - ERA5_QM
FAIR Overall Score
Abstract
The Soil Moisture Anomaly is a drought indicator used to detect and monitor agricultural drought conditions, defined by a prolonged period of deficit in the availability of soil moisture to plants. The soil moisture anomalies provided as part of the Alpine Drought Observatory was derived from ERA5 Volumetric Soil Water Layers at different depths with Layer 1 (0-7cm), Layer 2 (7-28cm), Layer 3 (28–100cm), and Layer 4 (100-289cm). The input ERA5 soil moisture dataset was downscaled using a quantile mapping approach. Daily anomalies were calculated using a climatological mean and standard deviation of soil moisture from a smoothed time-series (running mean on a 10-day window) for a reference period of 1981–2020. The Soil Moisture Anomalies values range from -5 to +5, with negative values indicating drier than conditions while positive values indicate wetter than normal conditions, and -1 to +1 values indicates near-normal conditions. Datasets contains modified Copernicus Climate Change Service Information [1980–current year]; contains modified Copernicus Atmosphere Monitoring Service Information [1980-current year]. Neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus information or data it contains.
Keywords
Legal constraints
CC-BY-4.0
FAIR Overall Score: 77%
| Principle | Score | Earned | Level |
|---|---|---|---|
| Findable | 57% | 4 of 7 | moderate |
| Accessible | 100% | 7 of 7 | advanced |
| Interoperable | 67% | 4 of 6 | moderate |
| Reusable | 83% | 5 of 6 | moderate |
Evaluated by F-UJI web service: Anusuriya Devaraju, & Robert Huber. (2020). F-UJI - An Automated FAIR Data Assessment Tool. Zenodo. https://doi.org/10.5281/zenodo.6361400
Snippet code
##### ----Explore and download STAC data with Python ----- #####
from pystac_client import Client
## Read the example catalog
URL = 'https://stac.eurac.edu/'
catalog = Client.open(URL)
## List the Collections in the given Catalog
stac_collections = list(catalog.get_collections())
print(f"Number of collections: {len(stac_collections)}")
## Print collection IDs
print("Collections IDs:")
for collection in stac_collections:
print(f"- {collection.id}")
print("-------------------------")
## Retrieve a specific collection
collection = catalog.get_collection("ID_OfTheCollection")
## Search for items in the collection
collection_items = list(catalog.search(collections=['ID_OfTheCollection'], max_items=10).items())
print(collection_items)
## Retrieve a list of the first 10 items belonging to a specific collection
item = collection.get_item("NameOfTheItem")
#print(list(item.assets.items())[0:10])
##Print the band name and the href of a specific item retrieved from the previous list
print(item.assets["ID_OfTheCollection"].title)
print(item.assets["ID_OfTheCollection"].href)
## Print Item’s assets through the assets attribute, which is a dictionary
for asset_key in item.assets:
asset = item.assets[asset_key]
print("{}: {} ({})".format(asset_key, asset.href, asset.media_type))
## ------------- DONWLOAD COG files from a specific collection ------- ####
from pystac_client import Client
import requests
## Read the example catalog
URL = 'https://stac.eurac.edu/'
catalog = Client.open(URL)
## List the Collections in the given Catalog
stac_collections = list(catalog.get_collections())
# print(f"Number of collections: {len(stac_collections)}")
## Print collection IDs
print("Collections IDs:")
for collection in stac_collections:
print(f"- {collection.id}")
print("-------------------------")
## Retrieve a specific collection
collection = catalog.get_collection("ID_OfTheCollection")
## Search for items in the collection
collection_items = list(catalog.search(collections=['ID_OfTheCollection'], max_items=10).items())
print(collection_items)
url_basepath = "https://eurac-eo.s3-eu-west-1.amazonaws.com/"
collection_url = url_basepath + "ID_OfTheCollection" + "/"
print(collection_url)
for item_id in collection_items:
# item = item.assets[asset_key]
print(item_id.id)
cog_file = str(item_id.id + '.tif')
url = collection_url + cog_file
response = requests.get(url)
with open(cog_file, "wb") as f:
f.write(response.content)
Related docs
| # | Name | Description | Link | Date published | Category |
|---|---|---|---|---|---|
| 1 | MOOC Cubes and Clouds | Free Online Course teaching the concepts of data cubes, cloud platforms and open science in geospatial and EO. | Link | March 8, 2024 | OpenEO, STAC |
| 2 | STAC guidelines | Documentation to browse and download items of the STAC catalog | Link | Feb. 22, 2024 | STAC |