Towards a Convective Scale Reanalysis with a - - PowerPoint PPT Presentation

towards a convective scale reanalysis with a soil
SMART_READER_LITE
LIVE PREVIEW

Towards a Convective Scale Reanalysis with a - - PowerPoint PPT Presentation

Towards a Convective Scale Reanalysis with a Soil-Vegetation-Atmosphere- Transfer-Model International Symposium on Regional Reanalyses Clarissa Figura, Insa Thiele-Eich, Jan D. Keller, Wolfgang Kurtz, Clemens Simmer, Andreas Hense 19.07.2018


slide-1
SLIDE 1

International Symposium on Regional Reanalyses

Clarissa Figura, Insa Thiele-Eich, Jan D. Keller, Wolfgang Kurtz, Clemens Simmer, Andreas Hense

19.07.2018

Towards a Convective Scale Reanalysis with a Soil-Vegetation-Atmosphere- Transfer-Model

slide-2
SLIDE 2

Reanalyses offer spatially and temporally consistent data sets for global or regional grids and a certain vertical exent within the land surface and atmosphere. Applications:

  • Renewable Energy Sector
  • Regional climate analyses
  • Hydrology
  • Agricultural economics

Motivation

slide-3
SLIDE 3

Reanalyses offer spatially and temporally consistent data sets for global or regional grids and a certain vertical exent within the land surface and atmosphere. Applications:

  • Renewable Energy Sector
  • Regional climate analyses
  • Hydrology
  • Agricultural economics

Demand for regional reanalyses is growing!

Motivation

slide-4
SLIDE 4

Challenges 1

Misrepresentation of hydrology (soil moisture-evaporation- precipitation feedback) in atmospheric reanalyses, especially at catchment scale

  • Betts et al. 1998 and 2003: Investigations of water and

energy budget in ECMWFs reanalysis for Mississippi and subbasins → Coupled reanalysis approach

slide-5
SLIDE 5

Challenges 2

There exists a coupling between convection triggering and soil moisture (Cioni and Hohenegger, 2017)

  • idealized LES-Simulations show a strong coupling between

soil moisture and diurnal precipitation cycle  Convective scale setup

slide-6
SLIDE 6

Challenges 3

Models are not perfect due to parametrizations and the chaotic nature of the systems they represent (Lorenz, 1969).  Ensemble approach with perturbed realizations

slide-7
SLIDE 7

Experimental setup TerrSysMP

With kind permissions of P. Shrestha and M. Sulis

slide-8
SLIDE 8

KENDA (Kilometre Scale ENsemble Data Assimilation) is a Local Ensemble Transform Kalman Filter

Experimental setup

slide-9
SLIDE 9

Experimental setup

KENDA

slide-10
SLIDE 10

Experimental setup

KENDA

slide-11
SLIDE 11

Experimental setup

KENDA → Basic Cycling System (BaCy) Since March 2017 operational at Deutscher Wetterdienst COSMO Ensemble COSMO Ensemble LETKF LETKF

Observations Observations COSMO Ensemble COSMO Ensemble

slide-12
SLIDE 12

Experimental setup

KENDA with TerrSysMP ...work in progress LETKF LETKF

Observations Observations TerrSysMP Ensemble TerrSysMP Ensemble TerrSysMP Ensemble TerrSysMP Ensemble

slide-13
SLIDE 13

Experimental setup

EMVORADO (Efficient Modular VOlume RADar Operator) → Radar forward operator

QR QS … HEMVORADO Model

  • utput

= xb Simulation of radar beams = yb

slide-14
SLIDE 14

Experimental setup

TerrSysMP TerrSysMP TerrSysMP TerrSysMP LETKF LETKF

?

slide-15
SLIDE 15

Experimental setup

slide-16
SLIDE 16

Preliminary work

COSMO-REA6

  • COSMO 4.25/TERRA
  • ∆xy= 6 km
  • Atmospheric boundary: ERA-Interim
  • Data assimilation:

Nudging of conventional observations (e.g., buoys, radio soundings, aircraft) Soil Moisture / Snow / SST Analysis

slide-17
SLIDE 17

Preliminary work

Long-term deterministic TerrSysMP simulation (Mauro Sulis)

  • TerrSysMP fully coupled (COSMO 4.x, CLM 3.5)
  • ∆xy=1000 m/ 500 m
  • Atmospheric boundary: COSMO-DE-Analysis
  • Downscaling (no data assimilation)
slide-18
SLIDE 18

First results

Ensemble TerrSysMP downscaling

  • TerrSysMP fully coupled (COSMO 5.1, CLM 3.5)
  • ∆xy=1000 m/ 500 m
  • Atmospheric boundary: COSMO-DE-KENDA-Analysis
  • Time period: 16.05.-13.06.2014
slide-19
SLIDE 19

First results

Case study: Front

12 h accumulated rain [mm] for 21.05.2014 00 UTC - 21.05.2014 12 UTC

slide-20
SLIDE 20

First results

Case study: Convection

12 h accumulated rain [mm] for 22.05.2014 12 UTC - 23.05.2014 00 UTC

slide-21
SLIDE 21

First results

Feedback precipitation – soil moisture

IQR Stdev Precipitation (mm/24h) Soil moisture (pressure) 28.05.2014

slide-22
SLIDE 22

First results

Feedback precipitation – soil moisture

IQR Stdev Precipitation (mm/24h) Soil moisture (pressure) 26.05.2014

slide-23
SLIDE 23

First results

Contingency table

hit false alarm miss correct negatives Observed Model yes no yes no Observations: 68 DWD rain gauges stations in NRW domain

www.bremerhaven-wetter.de

slide-24
SLIDE 24

First results

Contingency table

Model BIAS ETS Log Odds COSMO-REA6 1.71

  • 0.03
  • 2.36

TSMP det 1.13

  • 0.04
  • 2.09

TSMP ens 1.04

  • 0.05
  • 1.78

Threshold: 0.1 mm/h

Model BIAS ETS Log Odds COSMO-REA6 2.95

  • 0.01

TSMP det 1.06

  • 0.03

TSMP ens 0.95

  • 0.03
  • 2.62

Threshold: 0.2 mm/h

16.-22.05.2014

slide-25
SLIDE 25

First results

Contingency table

Model BIAS ETS Log Odds COSMO-REA6 1.71

  • 0.03
  • 2.36

TSMP det 1.13

  • 0.04
  • 2.09

TSMP ens 1.04

  • 0.05
  • 1.78

Threshold: 0.1 mm/h

Model BIAS ETS Log Odds COSMO-REA6 2.95

  • 0.01

TSMP det 1.06

  • 0.03

TSMP ens 0.95

  • 0.03
  • 2.62

Threshold: 0.2 mm/h

16.-22.05.2014

slide-26
SLIDE 26

First results

Shortcomings in representation of precipitation: → Data assimilation (LETKF) is expected to enhance the results

slide-27
SLIDE 27

Conclusion

  • Qualitatively TerrSysMP with dynamic downscaling of

analyses is able to better reproduce small scale precipitation events in comparison to COSMO-REA6

  • Extended simulation time period necessary to evaluate the

quantitative accuracy of precipitation with verification scores

slide-28
SLIDE 28

Outlook

  • Extended time period with TerrSysMP-KENDA with

EMVORADO planned

  • Quantitative model comparison and verification in terms of

precipitation and soil moisture with independent observations

  • Evaluation of feedback of soil moisture on precipitation
  • Evaluate impact of initial conditions
slide-29
SLIDE 29

International Symposium on Regional Reanalyses

Clarissa Figura, Insa Thiele-Eich, Jan D. Keller, Wolfgang Kurtz, Clemens Simmer, Andreas Hense

19.07.2018

Towards a Convective Scale Reanalysis with a Soil-Vegetation-Atmosphere- Transfer-Model

slide-30
SLIDE 30

First results

Contingency table

Frequency Bias (BIAS):

  • BIAS=(hits+false alarms)/(hits+misses)
  • ratio of the frequency of modeled events to the frequency of
  • bserved events

→ BIAS<1: underforecast, BIAS>1: overforecast, perfect: 1 Equitable Threat Score (ETS):

  • ETS=(hits+hitsrandom)/(hits+misses+false alarms-hitsrandom)
  • fraction of observed and/or modeled events that were correctly

predicted, adjusted for hits associated with random chance → Range: -1/3 to 1, no skill:0, perfect: 1 Log Odds Ratio (LOR):

  • LOR=log((hits+correct negatives)/(misses+false alarms))
  • ratio of the odds of making a hit to the odds of making a false

alarm → Range: -∞ - +∞ no skill:0, perfect: +∞