isu atmospheric component update – part i 3 rd rasm workshop 2012 fall

ISU Atmospheric Component Update – Part I

3rd RASM Workshop 2012 Fall

Justin GlisanIowa State University

2

Update

• Spectral nudging paper accepted by J. Clim.– Research funded through original DOE RACM grant– Co-authors: W.J. Gutowski, J.J Cassano, & M.E. Higgins

• Remaining two papers will be submitted to JGR• Postdoctoral SOMs research

J.Clim. Paper Background

• Long and short PAW simulations were run on the RACM domain

• A systematic, atmosphere-deep circulation bias formed within the northern Pacific storm track

• Various remedies tested, but with little success• Spectral or interior nudging was introduced

4

J. Clim. Paper Conclusions• Winter behavior more sensitive to nudging • Improve Cold Season Mean and Extreme Behavior– Stronger SN for precipitation– Weaker SN for surface temperatures

• Improve Warm Season Mean and Extreme Behavior– Weaker SN for precipitation– Stronger SN for surface temperatures

Optimal range for pan-Arctic simulations: 1/8th – 1/16th the WRF default

ISU Atmospheric Component Update – Part II

3rd RASM Workshop 2012 Fall

Brandon J. FiselIowa State University

Department of Geological and Atmospheric Sciences

Tucson, AZ 13-15 November

OVERVIEW

3rd RASM November 2012

Leap day in RASM Solved! WRF BC’s file name bug

ISU plans

RASM LEAP DAY

3rd RASM November 2012

Possible shift in WRF calendar? Compared lagged ERA-Interim (-48, -24,

+24 hrs) with RASM (bad leap, non-spectral, no Feb. 29). e.g., -48 hr lag -> ERA-Interim:1992-

Feb-28-00 UTC – RASM:1992-Mar-01-00 UTC.

Nothing suggestive of a WRF calendar problem (boundary discontinuities continue through simulation).

r31RBa2 (-48 hrs) r31RBa2 (-24 hrs)

r31RBa2 (+24 hrs)

RASM LEAP DAY

3rd RASM November 2012

LBC tendencies are computed on-the-fly in real.exe (share/module_bc -> stuff_bdytend subroutine).

Modifications made to met_em files. 2 methods:

1) blend (Seefeldt)

Method 1:met_em.d01.1992-02-28_00:00:00.nc -> met_em.d01.1992-02-28_00:00:00.ncmet_em.d01.1992-02-28_12:00:00.nc -> met_em.d01.1992-02-28_06:00:00.ncmet_em.d01.1992-02-29_00:00:00.nc -> met_em.d01.1992-02-28_12:00:00.nc

met_em.d01.1992-02-29_12:00:00.nc -> met_em.d01.1992-02-28_18:00:00.nc

met_em.d01.1992-03-01_00:00:00.nc -> met_em.d01.1992-02-29_00:00:00.nc

February 29 06, 12 and 18 UTC times left alone because those times are not important to RASM and will be removed from wrfbdy and wrffdda files after real.exe is completed.

RASM LEAP DAY

3rd RASM November 2012

LBC tendencies are computed on-the-fly in real.exe (share/module_bc -> stuff_bdytend subroutine).

Modifications made to met_em files. 2 methods:

1) blend (Seefeldt) 2) time-stagger

Method 1:met_em.d01.1992-02-28_00:00:00.nc -> met_em.d01.1992-02-28_00:00:00.ncmet_em.d01.1992-02-28_12:00:00.nc -> met_em.d01.1992-02-28_06:00:00.ncmet_em.d01.1992-02-29_00:00:00.nc -> met_em.d01.1992-02-28_12:00:00.nc

met_em.d01.1992-02-29_12:00:00.nc -> met_em.d01.1992-02-28_18:00:00.nc

met_em.d01.1992-03-01_00:00:00.nc -> met_em.d01.1992-02-29_00:00:00.nc

Method 2:met_em.d01.1992-02-28_12:00:00.nc -> met_em.d01.1992-02-28_12:00:00.ncmet_em.d01.1992-02-28_18:00:00.nc -> met_em.d01.1992-02-28_18:00:00.ncmet_em.d01.1992-03-01_00:00:00.nc -> met_em.d01.1992-02-29_00:00:00.nc

February 29 06, 12 and 18 UTC times left alone because those times are not important to RASM and will be removed from wrfbdy and wrffdda files after real.exe is completed.

After modification of met_em files, real.exe is run like usual. Following generation of wrfbdy and wrffdda files, February 29 is removed from files.

Scripts available on Garnet to make modifications to met_em and BC files.

r30RB1g r31RBa31 (method 1)

r31RBa32 (method 2) r31RBa31-r31RBa32

RASM LEAP DAY

3rd RASM November 2012

LBC tendencies are computed on-the-fly in real.exe (share/module_bc -> stuff_bdytend subroutine).

Modifications made to met_em files. 2 methods:

1) blend (Seefeldt) 2) time-stagger

both methods produced similar results; method 2 more simple

New files ready for use in RASM.

WRF BC NAME BUG

3rd RASM November 2012

Namelist setting allows user to specify a different name for BC files.

Bug – name overwritten to the default value (WRF v. 3.2 – 3.4.1). issue with new leap-day BC file

names WRF unable to find new Feb. files

Mediation_integrate.F

       IF ( .NOT. lbc_opened ) THEN         CALL construct_filename1 ( bdyname , 'wrfbdy' , grid%id , 2 )         WRITE(message,*)'Opening: ',TRIM(bdyname)         CALL wrf_debug(100,TRIM(message))         CALL open_r_dataset ( grid%lbc_fid, TRIM(bdyname) , grid , config_flags , "DATASET=BOUNDARY", ierr )          IF ( ierr .NE. 0 ) THEN            WRITE( message, * ) 'med_latbound_in: error opening ',TRIM(bdyname), ' for reading. IERR = ',ierr            CALL WRF_ERROR_FATAL( message )          ENDIF

       CALL construct_filename2a ( bdyname , grid%bdy_inname , grid%id , 2 , " " )       IF ( .NOT. lbc_opened ) THEN         WRITE(message,*)'Opening: ',TRIM(bdyname)         CALL open_r_dataset ( grid%lbc_fid, TRIM(bdyname) , grid , config_flags , "DATASET=BOUNDARY", ierr )          IF ( ierr .NE. 0 ) THEN            WRITE( message, * ) 'med_latbound_in: error opening ',TRIM(bdyname), ' for reading. IERR = ',ierr            CALL WRF_ERROR_FATAL( message )          ENDIF

WRF BC NAME BUG

3rd RASM November 2012

Namelist setting allows user to specify a different name for BC files.

Bug – name overwritten to the default value (WRF v. 3.2 – 3.4.1). issue with new leap-day BC file

names WRF unable to find new Feb. files

Mediation_integrate.F WRF folks will implement this fix in the

next version.

OVERVIEW

3rd RASM November 2012

Leap day in RASM ISU plans

Ongoing work Key questions Diagnostics Observational and GCM data Coupler/WRF output

ISU PLANS - ONGOING WORK

3rd RASM November 2012

Ongoing regimes and extremes work: Expanded the analysis to 17 years (1991 –

2007) using 6-member ensemble CORDEX simulations.

Key questions: Is there a similar breakdown in the

occurrence of seasonal regime behavior? Do we see more persistent regime behavior

becoming the “norm” during summer in recent years when sea ice is rapidly declining?

ISU PLANS - QUESTIONS

3rd RASM November 2012

ISU research questions: What effect does rapid loss in sea ice have

on the energy/moisture budgets? Also in the future?

Current and future spatial variations in energy/moisture budgets and relation to regime behavior change for same spatial areas? What effect is there on extremes in these regions?

ISU PLANS - QUESTIONS

3rd RASM November 2012

ISU research questions (continued): If sea ice continues to rapidly decline

thereby also delaying the re-freezing of sea ice, allowing more heat to be released by the ocean into the atmosphere, what effect does this have on future regime behavior?

How will regimes change in the future? Will the change in regime behavior be regionally dependent?

What effect would changing regime behavior have on extremes in the future Arctic? A more extreme Arctic?

ISU PLANS - DIAGNOSTICS

3rd RASM November 2012

Diagnostics: Monthly, seasonal and long-term (~20 years)

energy/moisture budget computations. Temporal and spatial variations on

energy/moisture budget terms on interannual and decadal time scales.

How do modeled (RASM) versus prescribed (stand-alone WRF) SST affect regime behavior?

Do RASM simulations have similar regime behavior over the same time period?

ISU PLANS - DIAGNOSTICS

3rd RASM November 2012

Diagnostics (continued): Sensitivity of RASM energy and moisture

budgets to model resolution. How does the NAO and AO affect regime

behavior? Influence of feedbacks in RASM versus stand-

alone WRF (cice+ocn+vic+datm). What is the interannual and decadal

variability of the net radiation from the ocean?

ISU PLANS - OBS DATA

3rd RASM November 2012

Observation data:ERA-Interim (BC’s)Global Precipitation Climatology Project

GEWEX Continental Scale International Project (GAPP) Surface Radiation Budget (SRB)

Pathfinder high resolution surface radiative fluxes

Canadian Meteorological Centre daily snow depth analysis data

Global Snow-Water-Equivalent Depth Coefficient-of-Variation Classification

ISU PLANS - OBS DATA

3rd RASM November 2012

Observation data (continued): Atmospheric Radiation Measurement (ARM)

Climate Research Facility AMSR-E/Aqua Daily 12.5 km Brightness

Temperature, Sea-Ice Concentration, & Snow Depth Polar Grids

ISU PLANS - GCM DATA

3rd RASM November 2012

GCM data: For future RASM simulations, what GCM BC’s

should be used? (e.g., ECHAM5 and HADCM3)

Sea ice declining faster than in GCM’s. Future sea-ice decline in RASM vs. GCM’s?

Image courtesy of:Stroeve et al., 2007. GRLdoi:10.1029/2007GL029703

Sea ice (brown) declining faster than GCM ensemble mean (solid black).

ISU PLANS - OTHER DATA

3rd RASM November 2012

Coupler and WRF data:Instantaneous vs. average output.

Need instantaneous values of state variables for budget studies. T, Q, U, V, ect.

SUMMARY

3rd RASM November 2012

WRF leap-day issue solved.WRF namelist bug submitted to wrfhelp.

ISU will evaluate energy/moisture budgets in RASM with focus on physical links between budgets and regimes/extremes.GCM data – what GCMs?

Need instantaneous values for energy/moisture budget studies.