Advancements to the Operational HWRF Modeling System at EMC

Similar documents
HMON (HNMMB): Development of a new Hurricane model for NWS/NCEP operations

Testing and Evaluation of GSI Hybrid Data Assimilation for Basin-scale HWRF: Lessons We Learned

Performance of the 2013 Operational HWRF

Advancements in Operations and Research on Hurricane Modeling and Ensemble Prediction System at EMC/NOAA

Tropical Cyclone Modeling and Data Assimilation. Jason Sippel NOAA AOML/HRD 2018 WMO Workshop at NHC

Retrospective and near real-time tests of GSIbased EnKF-Var hybrid data assimilation system for HWRF with airborne radar data

Transitioning Physics Advancements into the Operational Hurricane WRF Model

The EMC Mission.. In response to operational requirements:

Development and research of GSI based hybrid EnKF Var data assimilation for HWRF to improve hurricane prediction

Development of the Basin-scale HWRF Modeling System

COAMPS-TC 2015 Version, Performance, and Future Plans

NCEP Operational Hurricane Modeling System. HWRF Performance Verification in 2015

Motivation & Goal. We investigate a way to generate PDFs from a single deterministic run

NCEP Hurricane Modeling: HWRF 2016 Upgrade Plans. Sam Trahan Hurricane Team NOAA NCEP EMC and collaborators around the world

Impact of Assimilating Aircraft Reconnaissance Observations in Operational HWRF

2012 AHW Stream 1.5 Retrospective Results

HWRF Ocean: MPIPOM-TC

Developmental Testbed Center (DTC) Project for the Hurricane Forecast Improvement Program (HFIP)

Expansion of NCEP Operational Hurricane Weather Research and Forecast (HWRF) Model Forecast Guidance to all Global Tropical Cyclones

Operational Hurricane Modeling at NCEP/EMC

Discussion on HFIP RDITT Experiments. Proposal for extending the life of RDITT for one more year: Future Plans from Individual Groups

Overview of the NCEP Operational HWRF Modeling System

Assimilation of Airborne Doppler Radar Observations Using the Unified GSI based Hybrid Ensemble Variational Data Assimilation System for HWRF

Overview of HFIP FY10 activities and results

Numerical Weather Prediction: Data assimilation. Steven Cavallo

NCEP Operational Hurricane Modeling System. HWRF Performance Verification in 2016

Transition of Research to Operations

Performance and Verification of HWRF/HMON Ensemble Prediction System in 2017 Real time Parallel Experiment. Zhan Zhang, Weiguo Wang

PSU HFIP 2010 Summary: Performance of the ARW-EnKF Real-time Cloud-resolving TC Ensemble Analysis and Forecasting System.

WRF Model Simulated Proxy Datasets Used for GOES-R Research Activities

University of Miami/RSMAS

The role of testbeds in NOAA for transitioning NWP research to operations

Evaluation and Improvement of HWRF PBL Physics using Aircraft Observations

Hurricane Forecast Improvement Project (HFIP) Bob Gall HFIP Development Manager

ARW/EnKF performance for the 2009 Hurricane Season

An Overview of COAMPS-TC Development and Real-Time Tests

Case study validation of HWRF-HYCOM and HWRF-POM for Hurricane Isaac (2012)

DA/Initialization/Ensemble Development Team Milestones and Priorities

Track sensitivity to microphysics and radiation

Recent COAMPS-TC Development and Future Plans

BASIN-SCALE HWRF: Evaluation of 2017 Real-Time Forecasts

Coupled Ocean-Wave Model Team (Team 8) Report

Tropical cyclone simulations and predictions with GFDL s prototype global cloud resolving model

National Oceanic and Atmospheric Administration Hurricane Forecast Improvement Program Five-Year Strategic Plan

Evolution of the GFDL hurricane model in research and transition to NMC operations

Richard M. Yablonsky University of Rhode Island. WRF for Hurricanes Tutorial Boulder, CO 25 February 2010

HFIP Diagnostics Workshop Summary and Recommendations

Moving to a simpler NCEP production suite

HWRF Surface Layer Thermodynamics Evaluation. Eric W. Uhlhorn and Joseph J. Cione HFIP Hurricane Modeling Workshop September 2012

NOAA s Hurricane Forecast Improvement Project: Framework for Addressing the Weather Research Forecasting Innovation Act of 2017

Satellite Applications to Hurricane Intensity Forecasting

Experiments of Hurricane Initialization with Airborne Doppler Radar Data for the Advancedresearch Hurricane WRF (AHW) Model

Next Generation HWRF. HYCOM coupling. National Centers for Environmental Prediction. HWRF & HYCOM modeling team. Environmental Modeling Center

Hurricane Structure: Theory and Diagnosis

Developmental Testbed Center: Core Activities for HFIP

Recent Data Assimilation Activities at Environment Canada

2012 and changes to the Rapid Refresh and HRRR weather forecast models

Na#onal Hurricane Center Official Intensity Errors

The Developmental Testbed Center: Update on Data Assimilation System Testing and Community Support

HFIP Coupling Working Group

GSI Data Assimilation System Support and Testing Activities: 2013 Annual Update

Improvements to the NCEP Global and Regional Data Assimilation Systems

HWRF sensitivity to cumulus schemes

Objectives for meeting

Unidata Policy Meeting Key Program Status

Exploring the Use of Dynamical Weather and Climate Models for Risk Assessment

GFDL Hurricane Model Ensemble Performance During the 2012 Hurricane Season

Assimilation of GPS RO Data for Tropical Cyclone Prediction with HWRF Chunhua Zhou1, Hui Shao1, Bill Kuo1 and Ligia Bernardet2

High-Res Physics Tiger Team Report

The National Weather Service of the Future: Building a Weather-Ready Nation

AHW Ensemble Data Assimilation and Forecasting System

Assimilate W88D Doppler Vr for Humberto 05

Advanced Hurricane WRF (AHW) Physics

Assimilation of CrIS and ATMS Radiances in HWRF to Improve Hurricane/Typhoon Forecasts

Structure & Intensity Change: Future Directions IWTC VI, Topic 1

Masahiro Kazumori, Takashi Kadowaki Numerical Prediction Division Japan Meteorological Agency

Tropical Cyclone Mesoscale Data Assimilation

HWRF Ocean: The Princeton Ocean Model. HWRF Tutorial NCWCP, College Park, MD January 2018

Current and Future Experiments to Improve Assimilation of Surface Winds from Satellites in Global Models

A pseudo-ensemble hybrid data assimilation system for HWRF

Initialization of Tropical Cyclone Structure for Operational Application

NWS Operational Marine and Ocean Forecasting. Overview. Ming Ji. Ocean Prediction Center National Weather Service/NCEP. CIOSS/CoRP

Hurricane Structure: Theory and Application. John Cangialosi National Hurricane Center

Advanced diagnostics of tropical cyclone inner-core structure using aircraft observations

5.3 TESTING AND EVALUATION OF THE GSI DATA ASSIMILATION SYSTEM

NHC Ensemble/Probabilistic Guidance Products

ABSTRACT 3 RADIAL VELOCITY ASSIMILATION IN BJRUC 3.1 ASSIMILATION STRATEGY OF RADIAL

HFIP ENSEMBLE PLAN. Jun Du (EMC/NCEP), presenting on behalf of the HFIP Ensemble Team:

Fernando Prates. Evaluation Section. Slide 1

The Use of GPS Radio Occultation Data for Tropical Cyclone Prediction. Bill Kuo and Hui Liu UCAR

HFIP Annual Review Meeting November 5-7, 2018 Embassy Suites by Hilton Miami International Airport 3974 NW S River Dr, Miami, FL 33142

Unifying the NCEP Production Suite

The impact of assimilation of microwave radiance in HWRF on the forecast over the western Pacific Ocean

Have a better understanding of the Tropical Cyclone Products generated at ECMWF

Use of the GFDL Vortex Tracker

Ting Lei, Xuguang Wang University of Oklahoma, Norman, OK, USA. Wang and Lei, MWR, Daryl Kleist (NCEP): dual resolution 4DEnsVar

Ocean Model Impact Study Proposal for 2015

A proposal to establish a dedicated Hurricane Model Evaluation and Improvement Project. Joseph Cione Eric Uhlhorn NOAA/Hurricane Research Division

18A.2 PREDICTION OF ATLANTIC TROPICAL CYCLONES WITH THE ADVANCED HURRICANE WRF (AHW) MODEL

GSI DATA ASSIMILATION SYSTEM: COMMUNITY SUPPORT AND PRELIMINARY TESTING RESULTS

Transcription:

Advancements to the Operational HWRF Modeling System at EMC Vijay Tallapragada Hurricane Team Lead Environmental Modeling Center, NCEP/NOAA/NWS Camp Springs, MD 20746. WRF for Hurricanes Tutorial, Science Day, DTC, Boulder, CO. Tuesday April 26, 2011

Overview of the Operational HWRF HWRF Atmosphere Movable, two-way nested vortex following grid 9km inner domain and 27km outer domain), 42 vertical layers Advanced physics from GFDL/GFS Advanced vortex initialization with GSI/3DVAR HWRF Ocean Coupled to Princeton Ocean Model (POM) in the Atlantic Basin Feature based initialization of loop current and warm/cold core rings, cold wake specification during spin-up phase Operational HWRF products Numerical guidance on 6-hrly hurricane track and intensity for as many as five storms (both Atlantic and Eastern Pacific) High-resolution swaths (hourly, 10 th of a degree) for wind and precipitation along the projected storm path Simulated GOES synthetic satellite imagery (IR, VIS and WV) and radar reflectivity Four years into operations, since 2007

Design of movable nest All interpolations are done on a rotated lat-lon, E-grid with the reference lat-lon located at the centre of the parent domain. The nested domain can be freely moved anywhere within the grid points of the parent domain, the nested domain lat-lon lines will coincide with the lat-lon lines of the parent domain at integral parent-tonest ratio. Vortex following nest motion is achieved by tracking minimum 3

Why the HWRF? The HWRF is a next generation forecast system. It is an advancement over the GFDL in the following ways: In order to address intensity problem, we need higher resolution of hurricane core circulation to resolve more details of the storm structure. The hurricane circulation is driven by non-hydrostatic motions. HWRF is a non-hydrostatic model. GFDL is not. Also, the HWRF has its own analyses of the hurricane core. HWRF makes use of the operational data assimilation for the core initialization. This will allow the use of observations to help define the storm strength and structure. (GFDL does not have a vortex analyses and uses a bogusing technique for storm initialization)

Evolution of HWRF Initial implementation in 2007 hurricane season Model design and development of movable nested grid started in 2002 Initial HWRF workshop at NSF in 2004 28 different configurations tested individually (each with about 200 simulations) before initial implementation Extensive 3-season (2004-2006) pre-implementation testing of HWRF for all storms in the Atlantic and Eastern Pacific basins Vortex initialization upgrades in 2008 Address intensity bias for weaker systems, modifications to storm balance Infrastructure upgrade and transition to P6 in 2009 Capability enhancements to allow coupling to HyCOM and Wave Watch-III Script enhancements (identical scripts for NCO operations and EMC parallels) Physics and initialization upgrades in 2010 to improve the forecast skill. New baseline version with several bug fixes Modified surface physics formulation and use of Gravity Wave Drag parameterization Addition of satellite radiance data assimilation in the hurricane environment Focus on reducing intensity bias 5

HWRF GFDL Grid configuration 2-nests (coincident) 3-nests(not coincident) Nesting Force-feedback Interaction thru intra-nest fluxes Ocean coupling POM (Atlantic only) POM Convective parameterization SAS SAS Explicit condensation Ferrier Ferrier Boundary layer GFS non-local GFS non-local Surface layer GFDL..(Moon et. al.) GFDL..(Moon et. al.) Land surface model GFDL slab GFDL slab Dissipative heating Based on D-L Zhang Based on M-Y TKE 2.5 Radiation GFDL (cloud differences) GFDL 6

HWRF Operational Hurricane Forecast System NHC storm message Position domain Get OBS, Model Input initial & boundary conditions Ocean Initialization Initialize wake, loop currents & eddies Wrf si (used for topographical parameters) Wrf real:replace with interpolations from native model data storm analysis and data ingest 6hr 1st guess vortex relocation 3DVAR gsi for both nests Next cycle HWRF model Coupled with POM Synoptic fields for many variables Create file for track, intensity, etc 7

8

Hurricane Forecast System Requirements Document for T20 Requirements document co-written between EMC and NHC for transition to NCEP operations of a new candidate hurricane model or for annual or periodic upgrades to an existing operational hurricane model for transition to NCEP operations. Sets standards for: Performance testing case load, evaluation criteria System functionality initialization, coupling Computational performance 9

Research model forecasts of storm intensity for Hurricane Katrina showing benefit of high (1.6km) resolution (S. Chen) Intensity Forecast of Hurricane Katrina 0000 UTC 27 August 2005 Obs HWRF model High-res models (1.6km) GFDL model NHC Forecast HWRF Global models Katrina Landfall

08/22 18Z Cat. 2 Cat. 1 08/21 06Z Cat. 5 Cat. 5 08/20 18Z Cat. 4 Cat. 5 08/18 00Z Cat. 4 Cat. 3 08/17 12Z Cat. 2 Cat. 1 Tropical Storm 08/16 00Z Hurricane Dean 5 day forecasts of maximum winds starting from 8/19/06Z 90 kts Cat. 2 139 kts Cat. 5 08/22 18Z 08/20 18Z 08/19 06Z TS Cat 1 Cat 2 Cat 3 Cat 4 Cat 5

2010 HWRF operational configuration Based on WRF NMM Dynamic Core V2.0 Movable, two-way nested vortex following grid 9km inner domian (6 x 6 o domain) and 27km outer (75 x 75 domain) 42 vertical layers Advanced physics (tested in GFDL and GFS) Advanced vortex initialization that uses GSI 3D-var assimilation Coupled to Princeton Ocean Model (POM) with feature based initialization of loop current and warm/cold core rings and cold wake specification Upgrades/changes to the operational FY2009 HWRF configuration to improve the forecast skill: New baseline version with several bug fixes for radiation, land surface temperature and nonhydrostatic advection of vertical velocity Modified surface physics formulation based on observed profiles of drag coefficients for momentum and heat fluxes Gravity Wave Drag parameterization Addition of satellite radiance data assimilation in the hurricane environment Additional operational products Simulated GOES synthetic satellite imagery (IR, VIS and WV) and radar reflectivity MLD and OHC from POM and operational archival of ocean output Operational capability to run 5 th storm 12

Performance of the Operational HWRF for 2010 ATL hurricane season HWRF Atlantic Track Forecast Errors (less skillful at all times compared to GFS and GFDL) HWRF Atlantic Intensity Forecast Errors (less skillful at all times) Rapid growth of intensity errors Negative bias for HWRF Atlantic Intensity Forecasts 13

Performance of the Operational HWRF for 2010 EPAC hurricane season HWRF E-Pac Track Forecast Errors (better than GFS at all times, worse than GFDL except at 120 hr) HWRF E-Pac Intensity Forecast Errors (less skillful at all times) Rapid growth of intensity errors Negative bias for HWRF E-Pac Intensity Forecasts 14

Vortex Initialization Issues/Challenges HWRF produces too few strong storms Positive bias for weaker storms, negative for stronger storms (spinup/spindown) Insufficient vortex size correction Broad initial vertical structure 15

Issues with nest - topography interaction and feedback Issues with evolution of large-scale flow 16

HWRF Simulated-GOES Images

Band 3 Loop Simulated Band 3 Loop Real Convection and moisture transport in HWRF systematic dryness 18

Maximum wind speed (kts) Maximum Wind Speed (kts) Max Wind vs centeral Pressure at 6hr forecast Hurricane Alex and Earl 120 100 80 ctl: init ctl: +06hr diff: +06hr no-init:+06hr Linear (ctl: init) Linear (ctl: +06hr) Linear (diff: +06hr) Linear (no-init:+06hr) W-P Relationship 60 40 20 920 930 940 950 960 970 980 990 1000 1010 1020 Central Pressure (mb) Time series of max. wind speed (Celia: 2010062218) 85 80 75 70 9kts drop at the first 9min Rapid weakening of strong storms 65 60 55 50 84kts becomes 55kts around the first 2hours 45 40 1 10 19 28 37 46 55 64 73 82 91 100 109 118 127 136 Forecast Time (9min) 19

Model Upgrades (Atmosphere) Upgrade dynamical core to NMM community version V3.2 (EMC-DTC Collaboration) New GFS Deep Convection, Improved surface physics, and new FY2011 GSI/GFS IC/BC (EMC- GFDL Collaboration) Vortex initialization upgrades Improved storm size correction based on radius of 34 kt winds or ROCI and dynamical masswind consistency of the initial vortex (EMC- HRD collaboration) Modification of synthetic storm and its application in the initialization (vortex cycling) Upgrade HWRF GSI to V2.5 (community code) Ocean Upgrades HWRF Upgrades FY2011 (05/24/2011) Coupling to HYCOM Ocean Model (*Withdrawn from 2011 implementation plans, will be tested in parallel) Operational HWRF product enhancements Satellite angle corrections for simulated GOES WV and IR imagery, additional simulated microwave products New enhanced HWRF website for product display and navigation High-frequency output (3 hourly) and additional derived variables for diagnostics (EMC-NHC- CIRA Collaboration) (Supported by DTC) Improved storm size correction based on R34/ROCI

FY2011 HWRF Pre-Implementation Test Results Reduction of model spin-up/spin-down

IGOR 2010091300 HR20: PARENT HRDI: PARENT GFDL: PARENT HR20: NEST HRDI: NEST GFDL: NEST 22

IGOR 2010091300 10m wind speed HR20: NEST HRDI: NEST GFDL: NEST Size of the storm is much smaller with modified initialization and comparable to that of GFDL model

Ongoing Model development & collaborative efforts (FY2011 and beyond) Further advancements to the HWRF modeling system (EMC, HRD, URI) third nest capability, advanced diagnostic capability (diapost), idealized simulations Planned Demo during 2011 hurricane season (HFIP Stream 1.5) Multiple moving nests/improved model dynamics/ improved physics suitable for high resolution Fully coupled 3-way interactive oceanwave-atmosphere system Code management and community support (EMC, DTC) Setup HWRFV3.2+ and provide support to the community through DTC Development of R2O/O2R infrastructure and testing facility at DTC Advanced Data Assimilation (EMC, ESRL, OU & HRD) Hybrid EnKF/3DVAR DA Assimilation of real-time TDR data Assimilation of advanced satellite datasets 27 km 9 km 3 km NOAA Research with goals for transition to operations (27:09:03)

HWRF Physics (URI, GFDL, ESRL,HRD) Surface fluxes, sea spray and wave coupling Physics for high-resolution (convection, micro physics, PBL, LSM ) HWRF Diagnostics (HFIP, EMC, NHC, FSU, CIRA, HRD, UMBC/UMD) Identifying forecast errors from different components of model physics and dynamics Hurricane model diagnostics, evaluation and verification Develop a common and comprehensive diagnostics framework and tools to integrate model output with available observations for verification Enhanced real-time product display and navigation HWRF Ensembles Large Scale Flow Perturbations; Initial Storm Structure Perturbations; Physics-Based Perturbations High-Resolution HWRF and other parallels Real-time demo of triple nested (27/9/3) HWRF (HFIP Stream 1.5) Real-time demo of high-resolution 9:3 HWRF (HFIP Stream 2) Real-time demo of Doppler Radar DA experiments Real-time demo of NOAH LSM Coupled HWRF

Scales of Motion & Multi- Scale Models 26 2010 30 km 10 km 3 km 26

Why do we need higher resolution? Current operational resolution of HWRF (9km) cannot resolve critical hurricane features that can impact intensity and rainfall: Mesovortices ~10-20km wide in radial direction Hot towers ~5-15km wide in radial direction Inner rainbands ~10-30km wide Eye, eyewall of small storms etc. Recent results from HRD s high-resolution HWRF runs at 9/3 configuration have shown improvements over the operational HWRF runs (despite several differences, including lack of ocean coupling) Unified operational/research versions of HWRF modeling system (V3.2) enabled us to examine and evaluate the impacts of higher resolution, and assess the possibilities for transition to operations

Results from High-Resolution (9:3) HWRF Test operational HWRF with increased horizontal resolution (9km outer domain and 3km inner domain) Use operational HWRF physics, vortex initialization, vortex cycling and coupling to POM (no changes to ocean initialization or resolution) Run entire 2010 North Atlantic storms (>400 simulations) Use these results as a benchmark to test further advancements and to evaluate 27/9/3 configuration for operational purpose

Similar initial vortex. After Yucatan, vortex is temporarily destroyed. High-res recovers: 30km RMW low-res produces100-200km RMW Actual RMW is 9-27km. Example: Karl 13L 2010

Example: Alex 01L 2010 Similar situation as Karl. Correct RMW in 9:3 Much bigger storm in 27:9

Oper. HWRF High-Res HWRF High Resolution runs greatly improved track and intensity forecasts for Hurricane Alex Oper. HWRF High-Res HWRF

Example: Lisa 14L 2010 Low-res vortex too large again and misses real-world eyewall contraction to 18 km at ~48hrs

No significant gains in track forecasts from high-resolution runs 2010 Atlantic Seasonal Track and Intensity Statistics Systematic loss of intensity forecast skill compared to lowres and HWRF-X

Lessons learned so far.. Higher resolution alone is not enough. We need physics appropriate for high res. Started seeing better structure from high-res runs Vortex initialization requires several changes (we may be able to resolve RMW, eye and eye-wall appropriately). We should explore inner core DA with available obs. Operationalization is a big challenge. With optimal configuration, it is possible to fit triple nested domains within existing computational resources (work in progress) Providing high-resolution capability within the operational framework allows developers to test new and innovative methods to improve forecast skill.

Ongoing work Porting numerous fixes and upgrades from 2011 HWRF configuration Interoperability of several physics packages within WRF framework, idealized tests of NAM & NCAR phys WRF-NMM framework issues: NMM core assumes Ferrier scheme, microphysics arrays are zeroed out in many places. Working on fixing this Nest feedback issues, topography and land/sea mask related issues 27:9:3 simulation takes about two hours for a 5-day forecast using 128 processors. Bottlenecks we can fix: increase domain size to allow more processors base_state_parent wasted calculations in areas where nest is not present ~20% of runtime can be saved 3D boundary arrays NMM has 12 extra 3D arrays for nest boundary communication, while only 2D and 1D arrays are needed. ~5% of runtime gain removing these arrays. Multiprocessor coupler will allow faster interpolations and reduce overheads.

Planned 2012 Operational HWRF System Three atmospheric telescoping nested domains: 27km resolution 75x75 degree domain 9km resolution ~10x10 degree stormfollowing 3km resolution ~5x5 degree stormfollowing Include new nest motion algorithm and other dynamics improvements from HRD s developmental version of HWRF Coupled with POM/HyCOM ocean model. New coupler and modified HWRF vortex initialization for third nest Changes to HWRF physics appropriate for 3 km resolution

Coupler for Third Nest Extension of 2-way interactive NCEP coupler (dashed lines represent future development) 37

Initialization for Third Nest Extension of EMC/HRD 2011 HWRF operational initialization 38

Improved Vortex Initialization and Inner Core Data Assimilation Data assimilation and vortex initialization (EMC, ESRL, OU, AOML) Real-time test of the P3 TDR data flow from aircraft to NCO/TOC/AOC and assimilation using advanced GSI. Improved vortex initialization (model consistent 3-D balanced vortex) Ensemble data assimilation - hybrid EnKF (ESRL, OU, AOML) Planned Demo during 2011 hurricane season (HFIP Stream 2) HWRF Hybrid DA 39

Background HWRF vortex initialization/relocation Warm start: GFS environment + relocated and modified HWRF 6hr forecast vortex cold start: GFS environment + bogus vortex x (multiple time levels of background if FGAT) x t-1, x, x t+1, hybrid short range ensemble forecast x k en Assimilated obs rawinsonde, pibal, class sounding, profiler, dropsonde, AIRCFT, AIRCAR, GPSIPW, surface marine/land/splas h-level/mesonet satellite wind GSI (includes Hybrid DA Capability) Data selection, quality control, data thinning Iterative minimization 3DVAR cost function 2J(x) = x T B -1 x + (Hx - y) T R -1 (Hx - y) + J c Potential and future obs NOAA and USAF: flight level data, SFMR data new dropsonde cloud radiance data satellite radiance data: HIRS, AMSU-A, AMSU- B/MHS, AIRS, GOES sounder NOAA P3 TDR Pseudo-MSLP climatological isotropic B c flow-dependent anisotropic B a hybrid B=(β 1-1 B c + β 2-1 P en os) ensemble covariance P en = Σ x k en ( x k en ) T Analysis x a Post Balance EnKF Update ensemble members (x k en ) 126 hours forecast (3, 6, 9 hr forecast) Ensemble forecast 40

Planned FY2011 Real-Time Experiments on t-jet Experiment Physics Coupling Init Resources Stream 1.5* HWRF 27:9:3 Stream 2.0 HWRF-NOAH LSM Stream 2.0 HWRF-Hybrid DA + Stream 2.0 HWRF-HYCOM (Global) Stream 2.0 HWRF 27:9:3 HYCOM Stream 2.0 HWRF 27:9:3 with SC and PBL Stream 2.0 HWRF Choose from Table 1 NOAH LSM and 2011 Oper. Physics 2011 Oper. Physics 2011 Oper. Physics Modified Oper. Phys. GFS Shallow Convection and PBL Radiation & Microphysics POM POM POM HYCOM HYCOM POM/ HYCOM POM/ HYCOM Chose from Table 1 2011 Oper. Init w/ cycling 2011 Oper. Init w/ cycling 2011 Oper. Init w/ cycling 2011 Oper. Init w/ cycling 2011 Oper. Init w/ cycling 2011 Oper. Init w/ cycling +Tagged to real-time GFS Hybrid DA Ensemble Runs on CCS or tjet 172 cores per forecast for 4 hrs 15 GB archive space per forecast To be run on CCS 172 cores per forecast for 4 hrs 15 GB archive space per forecast 232 cores per forecast for 4 hrs 30 GB archive space per forecast Development in Progress To be run on CCS or Vapor Development in Progress To be run on CCS or Vapor Development in Progress To be run on CCS or Vapor 41

Analysis domain Outer domain: 80x80 degree, resolution 27 km Ghost domain: 20x20 degree, resolution 9 km Data assimilated * Quality controlled tail Doppler radar radial velocity data within [-3 hr, 3 hr] are assimilated * rawinsonde, pibal, class sounding, profiler, dropsonde, AIRCFT, AIRCAR, GPSIPW, surface data (marine/land/splashlevel/mesonet), satellite wind * satellite radiance data: AMSU-A, AMSU- B/MHS, HIRS, AIRS, GOES sounder * Pseudo-MSLP Assimilation of tail Doppler radar data GSI (3DVAR) * FGAT is used to represent obs at their appropriate time * Use flow-dependent anisotropic background error covariance [Riishøjgaard (1998), Purser 2003b, 2005]. TOMAS (2010 11 03 2356 Z - 04 0230 Z) 42

Post balance * large super gradient wind was found in HWRF s vortex initialization * Tangential wind from TDR data assimilation without additional constraint also shows large departure from gradient wind * Post balance is introduced to improve mass-wind balance for hurricanes * Modify temperature field based on gradient wind balance and hydrostatic balance and adjust surface pressure field corresponding to new temperature field without post balance with post balance 43

HWRF Hybrid DA with TDR Data TDR0: anisotropic B with post balance HYBG: beta1=0.5, B from global GSI Var/EnKF hybrid ensemble forecast (80 ensemble members), without post balance HYBH: beta1=0.5, B from hwrf ensemble forecast (48 ensemble members), without post balance The storm could not intensify as the analysis storm is quite shallow.

Need for a more comprehensive ocean model with DA capability GFDL-POM and HWRF-POM are proven systems. Why transition to HYCOM? Scientific / technical issues: Use operational RTOFS-Atlantic model for IC and BC: Realistic real-time ocean Full ocean model: Separate support community Solid base for physics based improvement. Management of NCEP production suite: No support (at NCEP) for further development of POM based system Cannot afford to maintain multiple systems 45

The HyHWRF Hurricane Prediction System GSI: Gridpoint Statistical Interpolation 46

HYCOM: one-way nested to RTOFS IC: RTOFS using data assimilation BC: RTOFS 6-day forecast 3-D Primitive governing equations, solved on a curvilinear mesh (8~20 km) and hybrid vertical coordinates (26 layers). Staggered Arakawa C-grid Numerical Integration: the leap-frog scheme, using 2 split time steps for fast barotropic and slow baroclinic mode. Nowcast: one-way Ocean alone 6-hr simulation, forced by GFS melt with parametric winds. 47

Model SST daily average GOES/AVHRR data POM-HWRF (R2) HYCOM-HWRF (R3) HR32 POM vs HYCOM IKE example POM on left, too warm under IKE. HYCOM on middle, too cool in general. Known issues in RTOFS, varies in general from season to season.. MMAB is working on addressing these issues HR32 (right): relatively warmer SST in the tropics, except more colder cold wake. 48

3-season HWRF-HYCOM testing Neutral impact on track and intensity errors Investigate the cause of negative bias/increased SST cooling and address SST drift issues Improved RTOFS DA and ocean model physics Extended diagnostics to evaluate the coupled system Need to develop a unified coupled system for both Atlantic and Eastern Pacific basins (Global HyCOM) to take full advantage Increased negative bias

Coupling to Wave-Watch III NOAA/NCEP in-house wave model, based on WAM. Operational global and (nested) regional model. Specialized Atlantic and Pacific hurricane wave models with blended winds from GFS and GFDL model. WAVEWATCH III will be coupled to HWRF Deep ocean model resolution dictated by GFS model Higher coastal model resolution dictated by model economy Highest model resolution in areas of special interest Hurricane nests moving with storm(s) like GFDL and HWRF

Advancing the HWRF System FY2011 & Beyond Resolution/ Infrastructure Physics DA/ Vortex Initialization Ocean 2011 2012 2013 2014* 2015* HWRFV3.2+ Triple nested HWRF (27/9/3 km) Deep Convection, Surface Physics Storm size correction, dynamic mass-wind consistenty HYCOM Coupling Increased vertical resolution, higher model top, upgrades to WRF infrastructure, NEMS/ESMF/NMM-B, community R2O efforts (HFIP) Horizontal diffusion, Shallow convection, Microphysics, Radiation, Surface Physics, Coupling to Waves and Land Surface, Physics for high-resolution Inner core DA (Doppler Radar, satellite), improved masswind balance Hybrid-EnKF DA, advanced vortex relocation procedure, improved GSI Improved ocean data assimilation, physics and resolution, unified coupled system for ATL & EPAC Waves One-way Wave Coupling Two-way wave coupling, multi-grid surf zone physics,. effects of sea spray Diagnostics and Product Development HWRF Ensembles, Coupling to Hydrological/ Surge/ Inundation models, diagnostics, product development Ongoing Work 2011 upgrades Planned developments *Potential Computer upgrade

Hurricane-Wave-Ocean-Surge- Inundation Coupled Models NCEP/Environmental Modeling Center Atmosphere- Ocean-Wave-Land HWRF SYSTEM NMM Upgrade V2.0 to Hurricane HWRF V3.2+ atmosphere in FY2011 Land and coastal waters Atmosphere/oceanic Boundary Layer wave spectra GFDL NOAH Slab LSM* Model winds air temp. WAVEWATCH III* Spectral wave model SST currents other fluxes runoff fluxes radiative fluxes FY2012 Princeton Upgrade Ocean to Model HYCOM (POM) Ocean for Model Atlantic Region wave fluxes High resolution Coastal, Bay & Estuarine hydrodynamic model* elevations currents 3D salinities temperatures Surge* Inundation* *Future developments

Next Generation Modeling Efforts with HWRF Next generation modeling efforts include basin scale regional domain with multiple moving nests. Downstream applications include serving WFOs for storm surge forecasts, landfall related flooding and inundation forecasts. NOAA has signed an MOU with India on application of HWRF modeling system for operational tropical cyclone prediction over Indian Seas (collaboration with AOML, U. Purdue) With Global HyCOM becoming operational at NCEP, it is possible to implement coupled HWRF-HYCOM system for all oceanic basins in the world. Work is in progress towards flexible vortex initialization procedure (collaboration with DTC).

HWRF Domain With Multiple Moving Nests Basin scale domain 7 days forecast SDA and cycling Regional ensembles/products Daily Tropical Outlook/genesis Computational Efficiency (27:9;about 2 h; 168 CPUs)

Coupling to Land Surface Model LSM Routing Model www.emc.ncep.noaa.gov/hwrf

Stream flow At 30.63N and 89.90W 1. The peak in WET run is much higher than in NLDAS, DRY run, NAM and the observed The peak time in DRY run is more comparable to the observation than in NLDAS. USGS: http://waterdata.usgs.gov

Sidr, 2007 Moving Nest Sidr, 2007 Intensity Forecast Sidr, 2007 Initial vortex Sidr, 2007 Track Forecast

Nargis, 2008 Track Forecast Moving Nest, Mala 2006

SUMMARY There has been lot of progress advancing the hurricane modeling capabilities at EMC, thanks to active collaboration between research and operations. Improving intensity/structure forecasts are orders of magnitude more difficult than was track forecasts. Requires substantial effort between research and operational hurricane communities This tutorial is one way to attract young generation to tackle these scientific challenges.

Real-time and pre-implementation T&E HWRF products: http://www.emc.ncep.noaa.gov/gc_wmb/vxt/index.html Acknowledgements: HWRF team at EMC EMC and HFIP Management Thanks for your attention Questions? Collaborations with NHC, DTC, HRD, GFDL, URI, CIRA and other HFIP partners