dart tutorial sec’on 15: dart experiments: control and design · 1. can just ‘restart’ from...

9
The Na’onal Center for Atmospheric Research is sponsored by the Na’onal Science Founda’on. Any opinions, findings and conclusions or recommenda’ons expressed in this publica’on are those of the author(s) and do not necessarily reflect the views of the Na’onal Science Founda’on. ©UCAR DART Tutorial Sec’on 15: DART Experiments: Control and Design

Upload: others

Post on 23-Sep-2020

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

TheNa'onalCenterforAtmosphericResearchissponsoredbytheNa'onalScienceFounda'on.Anyopinions,findingsandconclusionsorrecommenda'onsexpressedinthispublica'onarethoseoftheauthor(s)anddonotnecessarilyreflecttheviewsoftheNa'onalScienceFounda'on.

©UCAR

DARTTutorialSec'on15:DARTExperiments:ControlandDesign

Page 2: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

DARTobserva'onsequence(obs_sequence)files

Containa!me-orderedlistofobserva'ondefini'ons:1.Typeofobserva'on(radiosondetemperature,radarreflec'vity),2.Loca'onofobserva'on,3.Timeofobserva'on,4.Observa'onerrorvariance,5.Addi'onalstuffforcomplexobserva'ontypes,6.Qualitycontrolvalues.Whendrivingfilterassimila'on,alsocontainobservedvalues.

(Anobserva'onsequencecanhave0ormorevaluesassociatedwitheachobserva'on;sec'on17).

DARTTutorialSec'on15:Slide2

Page 3: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

DARTexperimentsaredrivenbyinputobserva'onsequences

model_moddefinesashortest_!me_between_assimila!ons,Δt.Ini'al'meofensemble(frominputfileornamelist)isti.

Δt/2 Δt/2

Δt/2 Δt/2

Δt/2 Δt/2

ModelAdvance

ModelAdvance

ti

ti+Δt

ti+2Δt

Allobserva'onsinthis'meintervalupdateprioratti.

Filtercon'nuesun'lallobserva'onsinsequencehavebeenused.Firstobserva'oncanNOTbeearlierthanfirst‘window’.Ifa‘window’hasnoobserva'ons,noassimila'onisdone.

DARTTutorialSec'on15:Slide3

Page 4: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

DARTexperimentsaredrivenbyinputobserva'onsequences

Thiscapabilityistoolimited:Wouldliketoallow:1.Dynamicallyadjus'ngmodelΔtformodelswiththiscapability;2.WindowwidthsthatcouldbelessthanΔt;

Δt/2 Δt/2

Δt/2 Δt/2

Δt/2 Δt/2

ModelAdvance

ModelAdvance

ti

ti+Δt

ti+2Δt

Allobserva'onsinthis'meintervalupdateprioratti.

3.Timeinterpola'onforwardoperators.

DARTTutorialSec'on15:Slide4

Page 5: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

Dealingwithmodelswithmul'-level'medifferencing

t

t

Example:Leapfrog

1.Canjust‘restart’fromsingle'me-leveladereachassimila'on.

Thiscanleadtonumericalinstabilityif‘restarts’toofrequent.Limitofdenseobserva'onsin'me,becomesforwarddifferencing..

DARTTutorialSec'on15:Slide5

Page 6: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

Dealingwithmodelswithmul'-level'medifferencing

t

Example:Leapfrog

2.Canexpanddefini'onofmodelstatetoincludemul'ple'mes.Statevectorincludes'mestiandti+1forobserva'onswith

'mesinthisinterval.Thiscanimproveperformance.Alsopermitseasy'meinterpola'onforwardoperators.Worksincurrentimplementa'onifmodel_modinterfaceismodified.Mostmodelsdonotimplementthiscapability.

ti ti+1

DARTTutorialSec'on15:Slide6

Page 7: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

Experimenttypes

1.Realdatafilteringassimila'ons:observa'onsfrominstruments.2.ObservingSystemSimula'onExperiments(OSSEs):

Observa'onsaresynthe'c.Modelintegra'onsubs'tutesfortruth.Forwardoperatorforeachobserva'onappliedtotruthstate.RandomsamplefromN(0,σobs)addedin.σobsfromobserva'onsequencefile.

3.ObservingSystemExperiments(OSEs):

Userealobserva'ons,butwithholdsomewithpurpose.

DARTTutorialSec'on15:Slide7

Page 8: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

Experimenttypes

4.MixedOSEs/OSSEs:Addsynthe'cobserva'onstorealobserva'ons.‘Truth’forsynthe'ccomesfrommodelintegra'onfromlastassimilatedstatees'mate.

5.Observa'ontarge'ng:

GivenOSEorOSSE,Addobserva'onsinfuturetoimprovefutureperformance,Alreadydoneopera'onallyforweatherpredic'on.Whereshouldaplaneflytogetmostvaluableobserva'ons?

6.Smoothing:notcurrentlysupportedinDARTManhagan

Useobserva'onsinfuturetoimprovestatees'mate,Canbeturnedonwith&smoother_nml:num_lags.

DARTTutorialSec'on15:Slide8

Page 9: DART Tutorial Sec’on 15: DART Experiments: Control and Design · 1. Can just ‘restart’ from single ’me-level aer each assimilaon. This can lead to numerical instability if

1.   FilteringForaOneVariableSystem2.   TheDARTDirectoryTree3.   DARTRun>meControlandDocumenta>on4.   Howshouldobserva>onsofastatevariableimpactanunobservedstatevariable?

Mul>variateassimila>on.5.   ComprehensiveFilteringTheory:Non-Iden>tyObserva>onsandtheJointPhaseSpace6.   OtherUpdatesforAnObservedVariable7.   SomeAddi>onalLow-OrderModels8.   DealingwithSamplingError9.   MoreonDealingwithError;Infla>on10.   RegressionandNonlinearEffects11.   Crea>ngDARTExecutables12.   Adap>veInfla>on13.   HierarchicalGroupFiltersandLocaliza>on14.   QualityControl15.   DARTExperiments:ControlandDesign16.   Diagnos>cOutput17.   Crea>ngObserva>onSequences18.   LostinPhaseSpace:TheChallengeofNotKnowingtheTruth19.   DART-CompliantModelsandMakingModelsCompliant20.   ModelParameterEs>ma>on21.   Observa>onTypesandObservingSystemDesign22.   ParallelAlgorithmImplementa>on23.  Loca'onmoduledesign(notavailable)24.  Fixedlagsmoother(notavailable)25.   Asimple1Dadvec>onmodel:TracerDataAssimila>on

DARTTutorialIndextoSec'ons

DARTTutorialSec'on15:Slide9