Richard M. Yablonsky University of Rhode Island. WRF for Hurricanes Tutorial Boulder, CO 25 February 2010
|
|
- Maude Pitts
- 6 years ago
- Views:
Transcription
1 Richard M. Yablonsky University of Rhode Island WRF for Hurricanes Tutorial Boulder, CO 25 February
2 What is the Princeton Ocean Model? Three dimensional, primitive equation, numerical ocean model (commonly known as POM) Originally developed by Alan Blumberg and George Mellor in the late 1970 s Initially used for coastal ocean circulation applications (Blumberg and Mellor 1987) Open to the community during the 1990 s and 2000 s Many user generated changes incorporated into official code version housed at Princeton University 2
3 Developing HWRF Ocean (POM TC) Available POM code version transferred to University of Rhode Island (URI) in mid 1990 s POM code changes made at URI specifically to address ocean response to hurricane wind forcing This POM version coupled to GFDL hurricane model at URI Coupled GFDL/POM model operational at NCEP in 2001 Additional POM upgrades made at URI during 2000 s (e.g. initialization) and implemented in operational GFDL/POM Same version of POM coupled to operational HWRF in 2007 This POM version henceforth designated POM TC 3
4 Why Couple POM TC to HWRF? To create accurate SST field for input into the HWRF Evaporation (moisture flux) from sea surface provides heat energy to drive a hurricane Available energy decreases if storm c0re SST decreases Uncoupled hurricane models with static SST neglect SST cooling during integration high intensity bias One dimensional (vertical only) ocean models neglect upwelling, which can impact SST during integration (e.g. Yablonsky and Ginis 2009, MWR) 4
5 5
6 Physics of Storm Core SST Change 1) Vertical mixing/entrainment (Slide 7) 2) Upwelling (Slide 8) 3) Horizontal advection (Slide 9) 4) Heat flux to the atmosphere (Small by comparison) 6
7 1) Vertical mixing/entrainment Wind stress surface layer currents Current shear turbulence Turbulent mixing entrainment of cooler water POM-TC uses the Mellor-Yamada 2.5 turbulence closure submodel to parameterize vertical mixing Sea Warm surface sea temperature surface temperature decreases A T M O S P H E R E O C E A N Subsurface Cool subsurface temperature temperature increases This is a 1-D (vertical) process 7
8 2) Upwelling Cyclonic wind stress divergent surface currents Divergent currents upwelling Cyclonic hurricane Upwelling cooler water brought to surface vortex A T M O S P H E R E Warm sea surface temperature O C E A N Cool subsurface temperature This is a 3-D process 8
9 3) Horizontal advection Preexisting cold pool is located outside storm core Preexisting current direction is towards storm core Cyclonic hurricane Ocean currents advect cold pool under storm vortex core A T M O S P H E R E Warm sea surface temperature O C E A N Cool subsurface temperature This is a horizontal process 9
10 < < Prescribed propagation speed Cyclonic hurricane vortex < < A T M O S P H E R E < < Homogeneous initial SST < < O C E A N < < < Horizontally-homogeneous subsurface temperature < What is the impact of varying storm translation speed? 10
11 2.4 m s m s -1 11
12 3-D 1-D 2.4 m s 1 3-D 1-D 4.8 m s 1
13 POM TC Grid: United Region 47.5 N ~18-km grid spacing N 98.5 W W 13
14 POM TC Sigma VerJcal Coordinate 23 vertical sigma levels; free surface (η) Level placement scaled based on ocean bathymetry Largest vertical spacing occurs where ocean depth is 5500 m Location of 23 half-sigma levels when ocean depth is 5500 m: 5, 15, 25, 35, 45, 55, 65, 77.5, 92.5, 110, 135, 175, 250, 375, 550, 775, 1100, 1550, 2100, 2800, 3700, 4850, and 5500 m 14
15 Arakawa C Grid: External Mode Plan View Horizontal spatial differencing occurs on staggered Arakawa-C grid 2-D variables UA and VA are calculated at shifted location from η 15
16 Arakawa C Grid: Internal Mode Plan View Horizontal spatial differencing occurs on staggered Arakawa-C grid 3-D variables U and V are calculated at shifted location from T and Q T here represents variables T, S, and RHO Q here represents variables Km, Kh, Q2, and Q2l 16
17 VerJcal Grid: Internal Mode Elevation View Vertical spatial differencing also occurs on staggered grid 3-D variables W and Q are calculated at shifted depth from T and U T here represents variables T, S, and RHO Q here represents variables Km, Kh, Q2, and Q2l 17
18 Time Stepping POM TC has a split time step External (two dimensional) mode uses short time step: 22.5 seconds during pre coupled POM TC initialization 13.5 seconds during coupled POM TC integration Internal (three dimensional) mode uses long time step: 15 minutes during pre coupled POM TC initialization 9 minutes during coupled POM TC integration Horizontal time differencing is explicit Vertical time differencing is implicit 18
19 POM TC inijalizajon Prior to coupled model integration of HWRF/POM TC, POM TC must be initialized with a realistic, 3 D temperature (T) and salinity (S) field This T & S field must then be used to generate realistic ocean currents via geostrophic adjustment The spun up ocean must then incorporate the preexisting hurricane generated cold wake by applying the TC s wind stress using the NHC message file 19
20 Typical of Gulf of Mexico in Summer & Fall Typical of Caribbean in Summer & Fall 20
21 Approximate Locations of Oceanic Features During Hurricanes Katrina and Rita (2005) TX Rita LA Katrina MS AL GA Subsurface (75-m) ocean temperature during Katrina & Rita Gulf of Mexico warm core ring cold core ring Loop Current FL Warm Loop Current water and a warm core ring extend far into the Gulf of Mexico from the Caribbean Mexico Caribbean C Directly under Rita s & Katrina s track But how do we know the locations of (& how do we assimilate) these features in real-time? 21
22 How we modify GDEM T/S Climatology: Feature based modeling! Data-assimilated GDEM Altimetry Or define using real obs: e.g. AXBT 6 for LC and AXBT 13 (14) for WCR (CCR) Finally, assimilate SST and integrate ocean model for 48 hrs for geostrophic adjustment Start with GDEM T/S Look at altimetry/obs Define LC & ring positions Use Caribbean water along LC axis & in warm core ring center Make cold core ring center colder than env. Blend features w/ env. & sharpen fronts 22
23 Gustav : Ocean IniJalizaJon & Response (Next 6 Slides) 23
24 August GDEM T/S Climatology SST ~75-m Temperature Starting point is August GDEM T/S climatology August GDEM is then interpolated in time to start date by blending with September GDEM 24
25 Including Features & Sharpening SST ~75-m Temperature GDEM T/S climatology is modified using the feature-based model (see slide 19) This includes cross-frontal sharpening 25
26 00 hr Phase 1: GFS SST Assimilated SST ~75-m Temperature (Land/sea mask applied) At 00-hr phase 1, daily NCEP SST is assimilated into the upper ocean mixed layer T/S fields vertically-interpolated to POM σ-levels 26
27 48 hr Phase 1 / 00 hr Phase 2 SST and Surface Current Vectors ~75-m Temperature and Current Vectors (Land/sea mask applied) During 48-hr of phase 1 integration, SST is held constant while currents geostrophically adjust 48-hr phase 1 = 00-hour phase 2 27
28 72 hr Phase 2 / 00 hr Coupled SST and Surface Current Vectors ~75-m Temperature and Current Vectors (Land/sea mask applied) During 72-hr of phase 2 integration, cold wake is generated by applying NHC message file wind 72-hr phase 1 = 00-hour coupled HWRF/POM-TC 28
29 192 hr Phase 2 (like 120 hr Coupled) SST and Surface Current Vectors ~75-m Temperature and Current Vectors (Land/sea mask applied) During 120-hr of coupled HWRF/POM-TC run, cold wake is generated by HWRF wind + thermal forcing Cold wake is generated here by extending phase 2 29
30 HWRF/POM TC Coupling Wind speed (U a ) Temperature (T a ) Humidity (q a ) Hurricane Model Air Sea Interface Momentum flux (τ) Sensible heat flux (Q H ) Latent heat flux (Q E ) Momentum flux (τ) Temperature flux Shortwave radiation Ocean Model SST (T s ) 30
31 How to Run the POM TC Ocean IniJalizaJon: Technical Details (Next 13 Slides) 31
32 Scripts Scripts can be found in directory: ${HOME}/HWRF/ src/pomtc/ocean_scripts There are 5 scripts, 4 of which are kickit scripts Fifth script gfdl_pre_sortvit.sh is used internally by the fourth kickit script and should not be changed Four kickit scripts are designed to be run in order from kickit01 through kickit04 If only running default tutorial case , then no changes should be needed to the kickit scripts either Next slide explains how to run ocean initialization 32
33 Running the Ocean IniJalizaJon Go to the following directory: ${HOME}/HWRF/src/ pomtc/ocean_scripts Execute the four kickit scripts sequentially: [prompt]$./kickit01_sharpn.sh [prompt]$./kickit02_getsst.sh [prompt]$./kickit03_phase3.sh [prompt]$./kickit04_phase4.sh Next slide explains how to modify the kickit scripts if running a case other than default case
34 Modifying the kickit Scripts If running a case other than , edit the four kickit scripts before running them If running on Bluefire, the only two lines in each kickit script that need to be edited are: stormid=07l # e.g. SID = 07L start_date= # e.g. YYYYMMDDHH = If running on a machine other than Bluefire, these lines may also need to be edited, depending on the paths to your input datasets, source code, and output directories: data_d=/ptmp/hurrtutorial/datasets sorc_d=${home}/hwrf/src/pomtc work_d=/ptmp/${user}/hwrf/${stormid}/{start_date}/oceanprd 34
35 Procedures in kickit01_sharpn.sh Set the storm ID Set starting date for POM TC as YYYYMMDDHH Set which ocean climatology to use (default is GDEM) Set the ocean region (always use united with sharpening) Define the directories Create the work directory if it does not already exist Create the sharpening directory, overwriting old attempts Slice up the storm ID and starting date Use some pieces of starting date to select two climatology months Create parameter input_sharp depending on chosen climatology Create symbolic links for the input files Run the sharpening code Rename the sharpened climatology file in preparation for POM TC phase 1 End the sharpening script by returning to the ocean_script directory 35
36 Procedures in kickit02_getsst.sh Set the storm ID Set starting date for POM TC as YYYYMMDDHH Define the directories Create the work directory if it does not already exist Create the sst/mask/lonlat directory, overwriting old attempts Slice up the starting date, using HH to define the model cycle Create symbolic links for the GFS spectral input files Increase ulimit s to prevent a segmentation fault Run the getsst code Rename the sst/mask/lonlat files in preparation for POM TC phase 1 End the getsst script by returning to the ocean_script directory 36
37 Procedures in kickit03_phase3.sh Set the storm ID Set starting date for POM TC as YYYYMMDDHH Set the ocean region (always use united with sharpening) Define the directories Create the work directory if it does not already exist Create the phase3 directory, overwriting old attempts Slice up the starting date, using HH to define the model cycle Modify the phase3 parameter file by including the starting date Create symbolic links for all input files Run the POM TC code for phase 1 Rename the phase3 restart file in preparation for POM TC phase 2 End the phase3 script by returning to the ocean_script directory 37
38 Procedures in kickit04_phase4.sh Set the storm ID Set starting date for POM TC as YYYYMMDDHH Set the ocean region (always use united with sharpening) Define the directories Create the work directory if it does not already exist Create the phase4 directory, overwriting old attempts Slice up the starting date, using HH to define the model cycle Use programs gfdl_date2day and gfdl_day2date to manipulate the starting date, which is necessary if backing up three days to end phase 2 at the coupled starting date Run operational script gfdl_pre_sortvit.sh to extract track information for the specified storm from the yearly NHC hurricane message file (i.e. tcvitals); then append line of zeros Modify phase4 parameter file by including starting date, track, and phase3 restart files Create symbolic links for all input files Run the POM TC code for phase 2 Rename the phase4 restart file in preparation for coupled HWRF/POM TC run End the phase4 script by returning to the ocean_script directory 38
39 Executables There are 5 executable files associated with POM TC ocean initialization, which can be found in directory: $ {HOME}/HWRF/src/pomtc/ocean_exec The following slides list the function, input(s), output(s), and usage of each of these 5 executable files 39
40 gfdl_date2day.exe Function: Convert a date integer in YYMMDDHH format to Julian day in year YY Input(s): date2day.inp Output(s): fort.61 (date2day.dat) Usage: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_date2day.exe < date2day.inp 40
41 gfdl_day2date.exe Function: Convert year and Julian day input to date integer in YYMMDDHH format Input(s): day2date.inp Output(s): fort.61 (day2date.dat) Usage: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_day2date.exe < day2date.inp 41
42 gfdl_getsst.exe Function: Extract SST, land sea mask, and lon/lat data from the GFS spectral file Input(s): for11 (gfs.${start_date}.t${cyc}z.sfcanl) fort.11 (gfs.${start_date}.t${cyc}z.sfcanl) fort.12 (gfs.${start_date}.t${cyc}z.sanl) Output(s): fort.23 (lonlat.gfs) fort.74 (sst.gfs.dat) fort.77 (mask.gfs.dat) getsst.out Usage: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_getsst.exe >> getsst.out 42
43 gfdl_sharp_mcs_rf_l2m_rmy5.exe Function: Run the sharpening program, which takes the T/S climatology, horizontally interpolates it onto the POM TC grid for the United region domain, assimilates a land/sea mask and bathymetry, and employs the diagnostic, feature based modeling procedure Input(s): input_sharp fort.66 (gfdl_ocean_topo_and_mask.${region}) fort.8 (gfdl_gdem.${mm}.ascii) fort.90 (gfdl_gdem.${mmm2}.ascii) fort.24 (gfdl_ocean_readu.dat.${mm}) fort.82 (gfdl_ocean_spinup_gdem3.dat.${mm}) fort.50 (gfdl_ocean_spinup_gspath.${mm}) fort.55 (gfdl_ocean_spinup.bayuf) fort.65 (gfdl_ocean_spinup.fsgsuf) fort.75 (gfdl_ocean_spinup.sgyreuf) fort.91 (mmdd.dat) fort.31 (hwrf_gfdl_loop_current_rmy5.dat.${yyyymmdd}) fort.32 (hwrf_gfdl_loop_current_wc_ring_rmy5.dat.${yyyymmdd}) Output(s): fort.13 (gfdl_initdata.${region}.${mm}) sharpn.out Usage: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_sharp_mcs_rf_l2m_rmy5.exe < input_sharp > sharpn.out 43
44 gfdl_ocean_united.exe Function: Run POM TC ocean phase 1 or phase 2 (also known historically as ocean phase 3 and phase 4, respectively, as in the model code) Input(s): fort.10 (parameters.inp) fort.15 (nullfile if phase 1; track if phase 2) fort.21 (sst.gfs.dat) fort.22 (mask.gfs.dat) fort.23 (lonlat.gfs) fort.13 (gfdl_initdata.${region}.${mm}) fort.66 (gfdl_ocean_topo_and_mask.${region}) fort.14 (not used if phase 1; RST.phase3.${region} if phase 2) Output(s): RST.yymmddhh (RST.phase3.${region} if phase 1; RST.final if phase 2) phase3.out if phase 1; phase4.out if phase 2 Usage: Phase 1: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_ocean_${region}.exe > phase3.out Phase 2: ${HOME}/HWRF/src/pomtc/ocean_exec/gfdl_ocean_${region}.exe > phase4.out 44
45 Key References Gopalakrishnan, S., and Coauthors, 2010: Hurricane Weather and Research and Forecasting (HWRF) Model Scientific Documentation. L. Bernardet, Ed., 75 pp. Mellor, G. L., 2004: User s guide for a threedimensional, primitive equation, numerical ocean model (June 2004 version). Prog. in Atmos. And Ocean. Sci., Princeton University, 56 pp. Yablonsky, R. M., and I. Ginis, 2008: Improving the ocean initialization of coupled hurricane ocean models using feature based data assimilation. Mon. Wea. Rev., 136,
HWRF Ocean: MPIPOM-TC
HWRF v3.7a Tutorial Nanjing, China, December 2, 2015 HWRF Ocean: MPIPOM-TC Ligia Bernardet NOAA SRL Global Systems Division, Boulder CO University of Colorado CIRS, Boulder CO Acknowledgement Richard Yablonsky
More informationHWRF Ocean: The Princeton Ocean Model. HWRF Tutorial NCWCP, College Park, MD January 2018
HWRF Ocean: The Princeton Ocean Model Isaac Ginis Graduate School of Oceanography University of Rhode Island HWRF Tutorial NCWCP, College Park, MD 23-25 January 2018 1 1 Why Couple a 3-D Ocean Model to
More informationAtmosphere-Ocean Interaction in Tropical Cyclones
Atmosphere-Ocean Interaction in Tropical Cyclones Isaac Ginis University of Rhode Island Collaborators: T. Hara, Y.Fan, I-J Moon, R. Yablonsky. ECMWF, November 10-12, 12, 2008 Air-Sea Interaction in Tropical
More information6C.4 USING AXBTS TO IMPROVE THE PERFORMANCE OF COUPLED HURRICANE-OCEAN MODELS
6C.4 USING AXBTS TO IMPROVE THE PERFORMANCE OF COUPLED HURRICANE-OCEAN MODELS Richard M. Yablonsky* and Isaac Ginis Graduate School of Oceanography, University of Rhode Island, Narragansett, Rhode Island
More information12A.2 NUMERICAL SIMULATIONS OF THE HURRICANE INTENSITY RESPONSE TO A WARM OCEAN EDDY
12A.2 NUMERICAL SIMULATIONS OF THE HURRICANE INTENSITY RESPONSE TO A WARM OCEAN EDDY Richard M. Yablonsky* and Isaac Ginis University of Rhode Island, Narragansett, Rhode Island 1. INTRODUCTION Hurricanes
More information6A.4 IMPACT OF A WARM OCEAN EDDY S CIRCULATION ON HURRICANE-INDUCED SEA SURFACE COOLING WITH IMPLICATIONS FOR HURRICANE INTENSITY
6A.4 IMPACT OF A WARM OCEAN EDDY S CIRCULATION ON HURRICANE-INDUCED SEA SURFACE COOLING WITH IMPLICATIONS FOR HURRICANE INTENSITY Richard M. Yablonsky* and Isaac Ginis University of Rhode Island, Narragansett,
More informationEvolution of the GFDL hurricane model in research and transition to NMC operations
Evolution of the GFDL hurricane model in research and transition to NMC operations Robert E. Tuleya (Saic, EMC & CCPO/ODU) Morris Bender (GFDL) Isaac Ginis (URI) Tim Marchok (GFDL) URI/GSO Y. Kurihara
More informationLimitation of One-Dimensional Ocean Models for Coupled Hurricane Ocean Model Forecasts
4410 M O N T H L Y W E A T H E R R E V I E W VOLUME 137 Limitation of One-Dimensional Ocean Models for Coupled Hurricane Ocean Model Forecasts RICHARD M. YABLONSKY AND ISAAC GINIS Graduate School of Oceanography,
More informationImpact of a Warm Ocean Eddy s Circulation on Hurricane-Induced Sea Surface Cooling with Implications for Hurricane Intensity. Richard M.
Impact of a Warm Ocean Eddy s Circulation on Hurricane-Induced Sea Surface Cooling with Implications for Hurricane Intensity Richard M. Yablonsky Isaac Ginis Graduate School of Oceanography University
More informationCoupled Ocean-Wave Model Team (Team 8) Report
Coupled Ocean-Wave Model Team (Team 8) Report George Halliwell (co-lead, NOAA/AOML/PhOD) Hendrik Tolman (co-lead, NOAA/NCEP) Isaac Ginis (URI) Chris Fairall (NOAA/ESRL) Shaowu Bao (NOAA/ESRL) Jian-Wen
More information18A.2 PREDICTION OF ATLANTIC TROPICAL CYCLONES WITH THE ADVANCED HURRICANE WRF (AHW) MODEL
18A.2 PREDICTION OF ATLANTIC TROPICAL CYCLONES WITH THE ADVANCED HURRICANE WRF (AHW) MODEL Jimy Dudhia *, James Done, Wei Wang, Yongsheng Chen, Qingnong Xiao, Christopher Davis, Greg Holland, Richard Rotunno,
More informationP1M.4 COUPLED ATMOSPHERE, LAND-SURFACE, HYDROLOGY, OCEAN-WAVE, AND OCEAN-CURRENT MODELS FOR MESOSCALE WATER AND ENERGY CIRCULATIONS
P1M.4 COUPLED ATMOSPHERE, LAND-SURFACE, HYDROLOGY, OCEAN-WAVE, AND OCEAN-CURRENT MODELS FOR MESOSCALE WATER AND ENERGY CIRCULATIONS Haruyasu NAGAI *, Takuya KOBAYASHI, Katsunori TSUDUKI, and Kyeongok KIM
More informationCase study validation of HWRF-HYCOM and HWRF-POM for Hurricane Isaac (2012)
Case study validation of HWRF-HYCOM and HWRF-POM for Hurricane Isaac (2012) Pat Fitzpatrick and Yee Lau, Mississippi State University Hyun-Sook Kim, Marine Modeling and Analysis Branch, NOAA/NWS/NCEP/EMC
More information5. General Circulation Models
5. General Circulation Models I. 3-D Climate Models (General Circulation Models) To include the full three-dimensional aspect of climate, including the calculation of the dynamical transports, requires
More informationMotivation & Goal. We investigate a way to generate PDFs from a single deterministic run
Motivation & Goal Numerical weather prediction is limited by errors in initial conditions, model imperfections, and nonlinearity. Ensembles of an NWP model provide forecast probability density functions
More informationA New Ocean Mixed-Layer Model Coupled into WRF
ATMOSPHERIC AND OCEANIC SCIENCE LETTERS, 2012, VOL. 5, NO. 3, 170 175 A New Ocean Mixed-Layer Model Coupled into WRF WANG Zi-Qian 1,2 and DUAN An-Min 1 1 The State Key Laboratory of Numerical Modeling
More informationDevelopment of Ocean and Coastal Prediction Systems
Development of Ocean and Coastal Prediction Systems Tal Ezer Program in Atmospheric and Oceanic Sciences P.O.Box CN710, Sayre Hall Princeton University Princeton, NJ 08544-0710 phone: (609) 258-1318 fax:
More informationMODELING EAST ASIAN SUMMER MONSOON IN 1998 WITH A COUPLED REGIONAL AIR-SEA MODEL. Xuejuan Ren
INYS 2008 MODELING EAST ASIAN SUMMER MONSOON IN 1998 WITH A COUPLED REGIONAL AIR-SEA MODEL Xuejuan Ren Yaocun Zhang, Yongfu Qian, Huijuan Lin Department of Atmospheric Sciences, Nanjing University, Nanjing,
More informationHWRF Surface Layer Thermodynamics Evaluation. Eric W. Uhlhorn and Joseph J. Cione HFIP Hurricane Modeling Workshop September 2012
HWRF Surface Layer Thermodynamics Evaluation Eric W. Uhlhorn and Joseph J. Cione HFIP Hurricane Modeling Workshop 17-18 September 2012 Special thanks. HRD HWRF modeling team Gopal, Xuejin Zhang, Thiago
More informationOcean Data Assimilation and Initialization Procedure for the Coupled GFDL/URI Hurricane Prediction System
1918 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 22 Ocean Data Assimilation and Initialization Procedure for the Coupled GFDL/URI Hurricane Prediction System ALEKSANDR FALKOVICH NOAA/NCEP/EMC,
More informationPRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Event Response
PRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK June 2014 - RMS Event Response 2014 SEASON OUTLOOK The 2013 North Atlantic hurricane season saw the fewest hurricanes in the Atlantic Basin
More informationImproving Surface Flux Parameterizations in the NRL Coupled Ocean/Atmosphere Mesoscale Prediction System
Improving Surface Flux Parameterizations in the NRL Coupled Ocean/Atmosphere Mesoscale Prediction System LONG-TERM GOAL Shouping Wang Naval Research Laboratory Monterey, CA 93943 Phone: (831) 656-4719
More informationAdvanced Hurricane WRF (AHW) Physics
Advanced Hurricane WRF (AHW) Physics Jimy Dudhia MMM Division, NCAR 1D Ocean Mixed-Layer Model 1d model based on Pollard, Rhines and Thompson (1973) was added for hurricane forecasts Purpose is to represent
More informationThis supplementary material file describes (Section 1) the ocean model used to provide the
P a g e 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 The influence of ocean on Typhoon Nuri (2008) Supplementary Material (SM) J. Sun 1 and L.-Y. Oey *2,3 1: Center for Earth System
More informationImpact of a Warm Ocean Eddy s Circulation on Hurricane-Induced Sea Surface Cooling with Implications for Hurricane Intensity
MARCH 2013 Y A B L O N S K Y A N D G I N I S 997 Impact of a Warm Ocean Eddy s Circulation on Hurricane-Induced Sea Surface Cooling with Implications for Hurricane Intensity RICHARD M. YABLONSKY AND ISAAC
More informationMomentum Flux Budget Across Air-sea Interface under Uniform and Tropical Cyclone Winds
Momentum Flux Budget Across Air-sea Interface under Uniform and Tropical Cyclone Winds Yalin Fan 1, Isaac Ginis 2, Tetsu Hara 2 1. AOS, Princeton University / GFDL, NOAA, Princeton, New Jersey 2. Graduate
More informationTesting and Evaluation of GSI Hybrid Data Assimilation for Basin-scale HWRF: Lessons We Learned
4th NOAA Testbeds & Proving Ground Workshop, College Park, MD, April 2-4, 2013 Testing and Evaluation of GSI Hybrid Data Assimilation for Basin-scale HWRF: Lessons We Learned Hui Shao1, Chunhua Zhou1,
More informationThe Operational GFDL Coupled Hurricane Ocean Prediction System and a Summary of Its Performance
DECEMBER 2007 BENDER ET AL. 3965 The Operational GFDL Coupled Hurricane Ocean Prediction System and a Summary of Its Performance MORRIS A. BENDER NOAA/GFDL, Princeton, New Jersey ISAAC GINIS Graduate School
More informationNumerical Weather Prediction: Data assimilation. Steven Cavallo
Numerical Weather Prediction: Data assimilation Steven Cavallo Data assimilation (DA) is the process estimating the true state of a system given observations of the system and a background estimate. Observations
More informationNext Generation HWRF. HYCOM coupling. National Centers for Environmental Prediction. HWRF & HYCOM modeling team. Environmental Modeling Center
HWRF tutorial 1 Next Generation HWRF HYCOM coupling Hyun-Sook Kim HWRF & HYCOM modeling team 2018 HWRF Tutorial Environmental Modeling Center National Centers for Environmental Prediction US National Weather
More informationHurricane Intensity: Governing Factors and Forecasting Challenges. EAS 470 Final Paper Allison Wing
Hurricane Intensity: Governing Factors and Forecasting Challenges EAS 470 Final Paper Allison Wing Tropical cyclones are undoubtedly among the mostly deadly and destructive natural phenomena found on Earth
More informationApplications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations
Applications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations Miyazawa, Yasumasa (JAMSTEC) Collaboration with Princeton University AICS Data
More informationHFIP Coupling Working Group
HFIP Coupling Working Group Progress report with contributions of Tolman (NCEP), Chen (NRL), Ginis (URI) Hendrik L. Tolman Chief, Marine Modeling and Analysis Branch NOAA / NWS / NCEP / EMC Hendrik.Tolman@NOAA.gov
More information2013 ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Cat Response
2013 ATLANTIC HURRICANE SEASON OUTLOOK June 2013 - RMS Cat Response Season Outlook At the start of the 2013 Atlantic hurricane season, which officially runs from June 1 to November 30, seasonal forecasts
More informationOcean Model Impact Study Proposal for 2015
Ocean Model Impact Study Proposal for 2015 OMITT-1 Background: HWRF w/ 3D POM-TC Yablonsky et al. (2010 IHC) confirmed POM tended to under-cool in response to prescribed wind stress based on observed TC
More informationThe Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science
The Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science Outline Basic Dynamical Equations Numerical Methods Initialization
More informationCOUPLED TROPICAL CYCLONE-OCEAN MODELING FOR CAPABILITIES. Graduate School of Oceanography. Phone: Fax:
COUPLED TROPICAL CYCLONE-OCEAN MODELING FOR TRANSITION TO OPERATIONAL PREDICTIVE CAPABILITIES Isaac Ginis and Lewis M. Rothstein Graduate School of Oceanography University of Rhode Island, Narragansett,
More informationNear-surface weather prediction and surface data assimilation: challenges, development, and potential data needs
Near-surface weather prediction and surface data assimilation: challenges, development, and potential data needs Zhaoxia Pu Department of Atmospheric Sciences University of Utah, Salt Lake City, Utah,
More informationThe HIGHTSI ice model and plans in SURFEX
Air Snow T in Ice with snow cover T sfc x T snow Q si F si h s h s The HIGHTSI ice model and plans in SURFEX Bin Cheng and Laura Rontu Water Ice T ice h i Finnish Meteorological Institute, FI-11 Helsinki,
More informationImproving Air-Sea Coupling Parameterizations in High-Wind Regimes
Improving Air-Sea Coupling Parameterizations in High-Wind Regimes PI: Dr. Shuyi S. Chen Co-PI: Dr. Mark A. Donelan Rosenstiel School of Marine and Atmospheric Science, University of Miami 4600 Rickenbacker
More informationToward Accurate Coastal Ocean Modeling
Toward Accurate Coastal Ocean Modeling Peter C. Chu Naval Postgraduate School Monterey, CA 93943, USA Email: pcchu@nps.edu http://www.oc.nps.navy.mil/~chu International Council for Sciences, Scientific
More informationWRF Modeling System Overview
WRF Modeling System Overview Louisa Nance National Center for Atmospheric Research (NCAR) Developmental Testbed Center (DTC) 27 February 2007 1 Outline What is WRF? WRF Modeling System WRF Software Design
More informationCold air outbreak over the Kuroshio Extension Region
Cold air outbreak over the Kuroshio Extension Region Jensen, T. G. 1, T. Campbell 1, T. A. Smith 1, R. J. Small 2 and R. Allard 1 1 Naval Research Laboratory, 2 Jacobs Engineering NRL, Code 7320, Stennis
More informationWater Stratification under Wave Influence in the Gulf of Thailand
Water Stratification under Wave Influence in the Gulf of Thailand Pongdanai Pithayamaythakul and Pramot Sojisuporn Department of Marine Science, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
More informationThe Ocean-Atmosphere System II: Oceanic Heat Budget
The Ocean-Atmosphere System II: Oceanic Heat Budget C. Chen General Physical Oceanography MAR 555 School for Marine Sciences and Technology Umass-Dartmouth MAR 555 Lecture 2: The Oceanic Heat Budget Q
More informationJuly Forecast Update for North Atlantic Hurricane Activity in 2018
July Forecast Update for North Atlantic Hurricane Activity in 2018 Issued: 5 th July 2018 by Professor Mark Saunders and Dr Adam Lea Dept. of Space and Climate Physics, UCL (University College London),
More informationTHE INFLUENCE OF HIGHLY RESOLVED SEA SURFACE TEMPERATURES ON METEOROLOGICAL SIMULATIONS OFF THE SOUTHEAST US COAST
THE INFLUENCE OF HIGHLY RESOLVED SEA SURFACE TEMPERATURES ON METEOROLOGICAL SIMULATIONS OFF THE SOUTHEAST US COAST Peter Childs, Sethu Raman, and Ryan Boyles State Climate Office of North Carolina and
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season : August 12, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met Service) is
More informationComments by William M. Gray (Colorado State University) on the recently published paper in Science by Webster, et al
Comments by William M. Gray (Colorado State University) on the recently published paper in Science by Webster, et al., titled Changes in tropical cyclone number, duration, and intensity in a warming environment
More informationAssimilation Impact of Physical Data on the California Coastal Ocean Circulation and Biogeochemistry
Assimilation Impact of Physical Data on the California Coastal Ocean Circulation and Biogeochemistry Yi Chao, Remote Sensing Solutions (RSS)/UCLA; John D. Farrara, RSS; Fei Chai, University of Maine; Hongchun
More informationWeather Research and Forecasting Model. Melissa Goering Glen Sampson ATMO 595E November 18, 2004
Weather Research and Forecasting Model Melissa Goering Glen Sampson ATMO 595E November 18, 2004 Outline What does WRF model do? WRF Standard Initialization WRF Dynamics Conservation Equations Grid staggering
More informationImproved Fields of Satellite-Derived Ocean Surface Turbulent Fluxes of Energy and Moisture
Improved Fields of Satellite-Derived Ocean Surface Turbulent Fluxes of Energy and Moisture First year report on NASA grant NNX09AJ49G PI: Mark A. Bourassa Co-Is: Carol Anne Clayson, Shawn Smith, and Gary
More informationKevin E Trenberth NCAR
Attribution of Recent Increases in Atlantic Hurricane Activity Kevin E Trenberth NCAR Issues for detection and attribution of changes in hurricanes What has happened? How good is the observational record?
More informationRecent Developments in the Navy Coastal Ocean Model and its application as the ocean component in regional coupled forecast models
Recent Developments in the Navy Coastal Ocean Model and its application as the ocean component in regional coupled forecast models Tommy Jensen, Paul Martin, Clark Rowley, Tim Campbell, Richard Allard,
More informationEnergy transport and transfer in the wake of a tropical cyclone
Energy transport and transfer in the wake of a tropical cyclone Claudia Pasquero Department of Earth and Environmental Sciences Università degli Studi di Milano - Bicocca Agostino Meroni, Francesco Ragone
More informationExploring the Use of Dynamical Weather and Climate Models for Risk Assessment
Exploring the Use of Dynamical Weather and Climate Models for Risk Assessment James Done Willis Research Network Fellow National Center for Atmospheric Research Boulder CO, US Leverages resources in the
More informationENSO Outlook by JMA. Hiroyuki Sugimoto. El Niño Monitoring and Prediction Group Climate Prediction Division Japan Meteorological Agency
ENSO Outlook by JMA Hiroyuki Sugimoto El Niño Monitoring and Prediction Group Climate Prediction Division Outline 1. ENSO impacts on the climate 2. Current Conditions 3. Prediction by JMA/MRI-CGCM 4. Summary
More informationEnhancing predictability of the Loop Current variability using Gulf of Mexico Hycom
Enhancing predictability of the Loop Current variability using Gulf of Mexico Hycom Matthieu Le Hénaff (1) Villy Kourafalou (1) Ashwanth Srinivasan (1) Collaborators: O. M. Smedstad (2), P. Hogan (2),
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season Forecast : June 11, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met Service)
More informationLarge-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling Eric D. Skyllingstad
More informationObjectives for meeting
Objectives for meeting 1) Summarize planned experiments 2) Discuss resource availability Aircraft Instrumentation Expendables 3) Assign working groups to complete each experiment plan Flight planning and
More informationImpact of frontal eddy dynamics on the Loop Current variability during free and data assimilative HYCOM simulations
Impact of frontal eddy dynamics on the Loop Current variability during free and data assimilative HYCOM simulations Matthieu Le Hénaff (1) Villy H. Kourafalou (1) Ashwanth Srinivasan (1) George R. Halliwell
More informationNumerical Simulations and Observations of Surface Wave Fields under an Extreme Tropical Cyclone
SEPTEMBER 2009 F A N E T A L. 2097 Numerical Simulations and Observations of Surface Wave Fields under an Extreme Tropical Cyclone YALIN FAN, ISAAC GINIS, AND TETSU HARA Graduate School of Oceanography,
More informationAn Introduction to Coupled Models of the Atmosphere Ocean System
An Introduction to Coupled Models of the Atmosphere Ocean System Jonathon S. Wright jswright@tsinghua.edu.cn Atmosphere Ocean Coupling 1. Important to climate on a wide range of time scales Diurnal to
More informationLink between Hurricanes and Climate Change: SST
Link between Hurricanes and Climate Change: SST FRIDAY: PROJECT 2 Proposals due! Bring laptops on weds (at least one per group)! Discuss assessment on Weds. Roger Francois will talk on Friday about the
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season Forecast Issued: May 10, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met
More informationNHC Ensemble/Probabilistic Guidance Products
NHC Ensemble/Probabilistic Guidance Products Michael Brennan NOAA/NWS/NCEP/NHC Mark DeMaria NESDIS/STAR HFIP Ensemble Product Development Workshop 21 April 2010 Boulder, CO 1 Current Ensemble/Probability
More informationPre-Season Forecast for North Atlantic Hurricane Activity in 2018
Pre-Season Forecast for North Atlantic Hurricane Activity in 2018 Issued: 30 th May 2018 by Professor Mark Saunders and Dr Adam Lea Dept. of Space and Climate Physics, UCL (University College London),
More informationLarge-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling Eric D. Skyllingstad
More informationUse of the GFDL Vortex Tracker
Use of the GFDL Vortex Tracker Tim Marchok NOAA / Geophysical Fluid Dynamics Laboratory WRF Tutorial for Hurricanes January 24, 2018 Outline History & description of the GFDL vortex tracker Inputs & Outputs
More informationSIMULATION OF ARCTIC STORMS 7B.3. Zhenxia Long 1, Will Perrie 1, 2 and Lujun Zhang 2
7B.3 SIMULATION OF ARCTIC STORMS Zhenxia Long 1, Will Perrie 1, 2 and Lujun Zhang 2 1 Fisheries & Oceans Canada, Bedford Institute of Oceanography, Dartmouth NS, Canada 2 Department of Engineering Math,
More information1. Introduction. In following sections, a more detailed description of the methodology is provided, along with an overview of initial results.
7B.2 MODEL SIMULATED CHANGES IN TC INTENSITY DUE TO GLOBAL WARMING Kevin A. Hill*, Gary M. Lackmann, and A. Aiyyer North Carolina State University, Raleigh, North Carolina 1. Introduction The impact of
More informationRelated Improvements. A DFS Application. Mark A. Bourassa
Related Improvements in Surface Turbulent Heat Fluxes A DFS Application Center for Ocean-Atmospheric Prediction Studies & Department of Earth, Ocean and Atmospheric Sciences, The Florida State University,
More informationThe role of testbeds in NOAA for transitioning NWP research to operations
ECMWF Workshop on Operational Systems November 18, 2013 The role of testbeds in NOAA for transitioning NWP research to operations Ligia Bernardet 1* and Zoltan Toth 1 1 NOAA ESRL Global Systems Division,
More informationAnalysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Analysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes Carol Anne Clayson Dept.
More informationAssimilation of GPS RO Data for Tropical Cyclone Prediction with HWRF Chunhua Zhou1, Hui Shao1, Bill Kuo1 and Ligia Bernardet2
Assimilation of GPS RO Data for Tropical Cyclone Prediction with HWRF Chunhua Zhou1, Hui Shao1, Bill Kuo1 and Ligia Bernardet2 1 National Center for Atmospheric Research (NCAR) 2 NOAA/ESRL/GSD Sixth FORMOSAT-3/COSMIC
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationJuly Forecast Update for Atlantic Hurricane Activity in 2017
July Forecast Update for Atlantic Hurricane Activity in 2017 Issued: 4 th July 2017 by Professor Mark Saunders and Dr Adam Lea Dept. of Space and Climate Physics, UCL (University College London), UK Forecast
More informationATOC 5051 INTRODUCTION TO PHYSICAL OCEANOGRAPHY. Lecture 19. Learning objectives: develop a physical understanding of ocean thermodynamic processes
ATOC 5051 INTRODUCTION TO PHYSICAL OCEANOGRAPHY Lecture 19 Learning objectives: develop a physical understanding of ocean thermodynamic processes 1. Ocean surface heat fluxes; 2. Mixed layer temperature
More information6.1 CASE STUDY: SLOSH USING GRIDDED WIND FIELDS FOR HURRICANES IRENE-2011 AND SANDY-2012
6.1 CASE STUDY: SLOSH USING GRIDDED WIND FIELDS FOR HURRICANES IRENE-2011 AND SANDY-2012 Dongming Yang 1* and Arthur Taylor 2 1. Ace Info Solutions, Reston, Virginia 2. NOAA / NWS / Office of Science and
More informationApplying satellite-derived ocean measurements for tropical cyclone intensity studies and forecasts. Gustavo Jorge Goni
Applying satellite-derived ocean measurements for tropical cyclone intensity studies and forecasts Gustavo Jorge Goni National Oceanic and Atmospheric Administration Atlantic Oceanographic and Meteorological
More informationJuly Forecast Update for Atlantic Hurricane Activity in 2016
July Forecast Update for Atlantic Hurricane Activity in 2016 Issued: 5 th July 2016 by Professor Mark Saunders and Dr Adam Lea Dept. of Space and Climate Physics, UCL (University College London), UK Forecast
More informationOutput Module (last updated 8/27/2015)
Output Module (last updated 8/27/2015) Goals: Provide a common output format and content from all modeling groups for Tier 2 files. Provide a data file that is compatible with the GFDL Vortex Tracker program
More informationCoastal Antarctic polynyas: A coupled process requiring high model resolution in the ocean and atmosphere
Coastal Antarctic polynyas: A coupled process requiring high model resolution in the ocean and atmosphere Mike Dinniman and John Klinck Center for Coastal Physical Oceanography Old Dominion University
More informationCoupled Ocean Circulation and Wind-Wave Models With Data Assimilation Using Altimeter Statistics
Coupled Ocean Circulation and Wind-Wave Models With Data Assimilation Using Altimeter Statistics Abstract ZHANG HONG, SANNASIRAJ, S.A., MD. MONIRUL ISLAM AND CHOO HENG KEK Tropical Marine Science Institute,
More informationImportance of air-sea interaction on the coupled typhoon-wave-ocean modeling
Importance of air-sea interaction on the coupled typhoon-wave-ocean modeling Collaborators: I. Ginis (GSO/URI) T. Hara (GSO/URI) B. Thomas (GSO/URI) H. Tolman (NCEP/NOAA) IL-JU MOON ( 文一柱 ) Cheju National
More informationJohn Steffen and Mark A. Bourassa
John Steffen and Mark A. Bourassa Funding by NASA Climate Data Records and NASA Ocean Vector Winds Science Team Florida State University Changes in surface winds due to SST gradients are poorly modeled
More informationObservations and Modeling of SST Influence on Surface Winds
Observations and Modeling of SST Influence on Surface Winds Dudley B. Chelton and Qingtao Song College of Oceanic and Atmospheric Sciences Oregon State University, Corvallis, OR 97331-5503 chelton@coas.oregonstate.edu,
More informationRobert Weaver, Donald Slinn 1
1 1 Robert Weaver, Donald Slinn 1 Department of Civil and Coastal Engineering, University of Florida, Gainesville, Florida Supported by the US Office of Naval Research AGU Fall Meeting 2002 Poster OS72A-0342
More informationA Reformulation of the Logistic Growth Equation Model (LGEM) for Ensemble and Extended Range Intensity Prediction
A Reformulation of the Logistic Growth Equation Model (LGEM) for Ensemble and Extended Range Intensity Prediction Mark DeMaria NOAA/NESDIS, Fort Collins, CO Andrea Schumacher and Kate Musgrave CIRA/CSU,
More informationApril Forecast Update for Atlantic Hurricane Activity in 2018
April Forecast Update for Atlantic Hurricane Activity in 2018 Issued: 5 th April 2018 by Professor Mark Saunders and Dr Adam Lea Dept. of Space and Climate Physics, UCL (University College London), UK
More informationOcean Model Development for COAMPS
Ocean Model Development for COAMPS Paul Martin Naval Research Laboratory Stennis Space Center, MS 39529 phone: (228) 688-5447 fax: (228) 688-4759 email: martin@nrlssc.navy.mil Award #: N0001498WX30328
More informationExtratropical transition of North Atlantic tropical cyclones in variable-resolution CAM5
Extratropical transition of North Atlantic tropical cyclones in variable-resolution CAM5 Diana Thatcher, Christiane Jablonowski University of Michigan Colin Zarzycki National Center for Atmospheric Research
More informationNumerical simulations of tropical cyclone ocean interaction with a high resolution coupled model
University of Rhode Island DigitalCommons@URI Graduate School of Oceanography Faculty Publications Graduate School of Oceanography 1993 Numerical simulations of tropical cyclone ocean interaction with
More informationAtmosphere Ocean Coupled Dynamics of Cyclones in the Midlatitudes
2432 MONTHLY WEATHER REVIEW VOLUME 132 Atmosphere Ocean Coupled Dynamics of Cyclones in the Midlatitudes XUEJUAN REN Department of Atmospheric Sciences, Nanjing University, Nanjing, China WILLIAM PERRIE
More informationIntroduction. One way to monitor the status of this energy is though monitoring sea surface
Sears 1 Understanding Hurricane Intensity Using Sea Surface Height and Temperature Information John Sears (Plymouth State University) Robbie Hood (NASA-MSFC) Frank LaFontaine (Raytheon) Abstract Warmer
More informationDYNAMICAL DOWNSCALING OF COUPLED MODEL HISTORICAL RUNS
FINAL REPORT FOR PROJECT 1.5.4 DYNAMICAL DOWNSCALING OF COUPLED MODEL HISTORICAL RUNS PRINCIPAL INVESTIGATOR: DR. JOHN MCGREGOR, CSIRO Marine and Atmospheric Research, John.McGregor@csiro.au, Tel: 03 9239
More informationAdvances in Coastal Inundation Simulation Using Unstructured-Grid Coastal Ocean Models
Advances in Coastal Inundation Simulation Using Unstructured-Grid Coastal Ocean Models Bob Beardsley (WHOI) Changsheng Chen (UMass-Dartmouth) Bob Weisberg (U. South Florida) Joannes Westerink (U. Notre
More informationImproving Ocean Model Initialization for Coupled Tropical Cyclone Forecast Models Using GODAE Nowcasts
2576 M O N T H L Y W E A T H E R R E V I E W VOLUME 136 Improving Ocean Model Initialization for Coupled Tropical Cyclone Forecast Models Using GODAE Nowcasts G. R. HALLIWELL JR. AND L. K. SHAY MPO/RSMAS,
More informationAn upper to mid level circulation (Low) in the NW Caribbean is
WEEKLY WEATHER OUTLOOK BELIZE, CENTRAL AMERICA PERIOD: Monday-Monday, Oct. 1 October 8, 2012 DATE ISSUED: Monday, Oct 1, 2012 6:00 am RFrutos EcoSolutions & Services SYNOPSIS: An upper to mid level circulation
More information