Skip to content
Snippets Groups Projects

Migration to plugin template and add climpact functionality

Merged Mostafa Hadizadeh requested to merge migration into main
2 files
+ 267
2
Compare changes
  • Side-by-side
  • Inline
Files
2
+ 58
2
@@ -19,9 +19,10 @@ SPDX-License-Identifier: BSD-3-Clause
Barrier-free RCM data, representative for heat events, for urban impact models.
## Work assignment
## Description
This plugin extracts the development of the atmospheric temperature stratification during a specified heat/warm period from regional climate model data and tailors the data for the usage in urban impact models. For a selectable region of interest and period of time this plugin returns hourly temperature data at different layers/hights in a directly usable format (format specification). This plugin is based on the NUKLEUS ensemble of 3km horizontal grid.
Extension or composition of various existing plugins or program parts for the barrier-free transfer of meteorological data (representative of certain climatological heat events) into an urban climate model
Note that we are developing prototype here and have to make compromises for the first development stage.
## Input:
* predefined RCM data basis: NUKLEUS ensemble (9 member) → output is calculated for all 9 members
@@ -183,3 +184,58 @@ License management is handled with [``reuse``](https://reuse.readthedocs.io/).
If you have any questions on this, please have a look into the
[contributing guide][contributing] or contact the maintainers of
`heat2urbanimpact`.
# Rest of descriptions
## Input:
Input parameters for the plugin and help texts:
- **shape_file** (): "Select a geo reference file defining your region of interest, if None is selected (default), the whole region that is defined in the climate data will be taken. Note: You can either select a path on the HPC system, or a web url pointing to the file."
- **region** (selectfield): "Select a pre defined German municipality of interest. This selection has only effect if you don't chose a shape file."
- **split_by** (string): "If your selected geo reference file has multiple geometries and you whish to process each geometry separatly you can choose a split key that splits the geometries into different sub regions. The values for the split key will thereby used to distinguish the sub regions."
- **experiment** (selectfield): "Please select the model simulation of interest: historical or projection (global warming level +2K or +3K)."
- **event** (selectfield): "Please select the exceptionality of the heat/warm period by selecting percentiles of temperature probability distribution. Default is the 95th percentile."
- **length_of_event** (integer): "Please enter the number of days of the heat/warm period you would like to receive as output. An integer between 1 and 14 days would be plausible. Default is 3 days."
- **months_of_event** (string): "Please enter the month(s) when your heat/warm event is supposed to happen for your studies. E.g. enter for single month 'June' or for multiple months 'May:August'. Default is 'July'."
- **impact_model** (selectfield): "The format of the output will be tailored for the selected impact model (PALM or ENVI-met)."
Predefined parameters for the prototype:
- project: nukleus
- product: ceu-3
- models: all
- time_frequency: daily and hourly (issue of bias correction)
- variable: tas, tasmin, tasmax (temperature at different hights are not needed for statistical analysis of heat/warm days, but for the final output)
## Method --> still to be figured out!
### Data basis
It is still unclear if hourly bias corrected data of vertical air temperature profiles will be available. Therefore, we see three options to start the data anaylsis:
* option 0 (optimal) - bias-corrected data available: bias-corrected model data is used directly for picking out events and outputting data
* option 1: bias is (proven) not significant: Procedure as variant 0
* option 2*: bias-corrected data are available but only for daily data (mean, min, max): Difference: bias-corrected maximum - uncorrected maximum --> imprinting of this difference on the hourly values of the day under consideration or other method? Method applicable for all height levels? Method applicable for all variables?
* option 3 (most unfavorable): bias is unknown, is judged to be significant and no bias-corrected data: The raw model data are used to derive the climate change signal (e.g. 2K-World - reference period). The climate change signal is suitably imposed on measurement data. Challenge: Derivation of vertical profiles (only measurements at a height of 2 m are available for temperature).
*most plausible and feasible option at the moment
### "Filter out" heat/warm days and create times series of desired length
We are not looking for a specific heat/warm event or a concrete number of heat or summer days in the data. For the chosen experiment (historical, 2K, or 3K) we would like to know how the temperature stratification typically evolves when it is e.g. 2 days exceptionally hot/warm, or 5 days or 10 days. The exceptionality of the heat/warm period is defined with the input parameter _event_ (percentiles of temperature probability distribution) and the length of the heat period is defined with the input parameter _length_of_event_. To consider also the seasonality of heat/warm periods the user can select different months with the input parameter _months_of_event_.
To provide hourly time series is a must-have. We will select a heat/warm day based on bias corrected daily tmax data and then go on as described under option 2 (Data basis section) to imprint the bias correction to the hourly data and different hight levels.
**We still researching how to define realistic time series of heat events using a statistical data analysis or how to develop an artificial timeseries.**
## Final Plugin Output
As a result we get 9 timeseries (ensemble members) with a frequency of 1hr and a length of 1 to 14 days (selectable via _length_of_event_) of vertical temperature stratification (representative for selected heat/warm period via _event_) for the chosen historical, 2K, or 3K experiment which can be fed directly into an urban climate model (format is selectable via _impact_model_).
## Open issues:
* availability of hourly bias corrected data?
* availability of hourly topography corrected data?
* output options:
* offer plots of all 9 timeseries (+ reference, observation) to visualize ensemble bandwidth (+ difference to past)
* future wish list:
* timeseries of reference data and observation data for comparison of 2K/3K future with the past
* timeseries of additional variables: thermal stress (tas, wind speed, air humidity, short/long-wave irradiation surface), wind speed and direction (optimal: vertical profile or at least at AGS height geostrophic wind), shortwave (direct, diffuse) radiation at surface, air humidity (optimal: vertical profile), boundary conditon at ground surface (only initialisation): soil temperature and humidity
Loading