560 likes | 579 Views
Explore the impact of topographic parameters on surface radiation with correction factors and case studies. Contact Matteo Buzzi for insights and summaries.
E N D
Grid scale parameterization of topographic effects on radiation Matteo Buzzi, Mathias W. Rotach and Marco Arpagaus MeteoSwiss
Overview • Introduction • The topographic problem • General aspects • Topographic parameters • Parameterization scheme • Corrections factors • Corrected surface radiation components • Case studies • Summary, conclusions and outlook Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
The problem: shortwave radiation . Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
The problem: thermal downward radiation 4 esT eff 4 esT 4 esT Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 slope slope
Introduction • NWP mesoscale models resolution tends to the order of magnitude of 1 km: topographic effects become important • Topographic shadowing (horizon, sun elevation) • Slope angle • Slope aspect • Sky view (horizon) • aLMo radiation scheme (Ritter and Geleyn,1992): surface radiation is computed on horizontal surfaces, topographic effects are not considered • Idea: try a simple approach with correction factors for surface radiation components Müller and Scherrer (2005) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Overview • The topographic problem • Introduction • Topographic parameters • Parameterization scheme • Corrections factors • Corrected surface radiation components • Case studies • Summary, conclusions and outlook Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topographic parameters: slope aspect (7 km) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topographic parameters: slope aspect (2 km) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topographic parameters: slope angle (7 km) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topographic parameters: slope angle (2 km) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Overview • The topographic problem • Introduction • Topographic parameters • Parameterization scheme • Correction factors • Corrected surface radiation components • Case studies • Summary, conclusions and outlook Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Sky view factor: 7 km Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Sky view factor: 2 km Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
2004-12-11 08 UTC (7 km) Correction factor for direct radiation radiation on horizontal surface radiation on inclined and oriented surface Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
2004-12-11 08 UTC (2 km) Correction factor for direct radiation SE slope increased radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 Shadowing NW slope
2004-12-11 12 UTC (7 km) Correction factor for direct radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
2004-12-11 12 UTC (2 km) Correction factor for direct radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
2004-12-11 15 UTC (7 km) Correction factor for direct radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
2004-12-11 15 UTC (2 km) Correction factor for direct radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Direct solar radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 Müller and Scherrer (2005)
Diffuse downwelling solar radiation Reduced sky view Reflections by adjacent areas Assumptions: • Diffuse downwelling radiation considered isotropic • Adjacent areas have the same upwelling diffuse radiation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 Müller and Scherrer (2005)
Thermal downwelling radiation Reduced sky view Radiation from adjacent areas Assumptions: • fac = ratio of the upwelling thermal radiation of a grid-point and its adjacent points (at higher elevation): considers a horizontal length scale (7 or 4 km), a mean slope angle (sky view) and a mean vertical temperature lapse rate (6.5°C/1000m) Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 Modified from Müller and Scherrer (2005)
Overview • Introduction • Topographic parameters • Parameterization scheme • Correction factors • Corrected surface radiation components • Case studies • 2004-12-11: aLMo 3.15, 7 km and 2 km • 2003-06-22: aLMo 3.15, 7 km and 2 km • Summary, conclusions and outlook Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2004-12-11 00 UTC +09 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2004-12-11 00 UTC +12 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2004-12-11, grid-point Correction factor NET solar radiation Surface temperature Sensible heat flux at the surface Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 corrected control
Results: 7 km, 2004-12-11, grid-point Correction factor NET solar radiation Surface temperature Snow amount Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 corrected control
Results: 2 km, 2004-12-11 00 UTC +09 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2004-12-11 00 UTC +12 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2004-12-11 00 UTC +14 Total cloud cover difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2004-12-11 00 UTC +12 10 m wind vectors difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 corrected control
Results: 7 km, 2003-06-22 00 UTC +08 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2003-06-22 00 UTC +12 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2003-06-22 00 UTC +08 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2003-06-22 00 UTC +12 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2003-06-22 00 UTC +12 Total cloud cover difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2003-06-22 00 UTC +14 Total cloud cover difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Overview • Introduction • Topographic parameters • Parameterization scheme • Correction factors • Corrected surface radiation components • Case studies • Summary, conclusions and outlook Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Summary and conclusions • The Müller and Scherrer (2005) scheme for topographic effects on radiation has been implemented into the aLMo • Some sensitivity case studies have been carried out • The impact of the topographic effects (shadowing, slope angle, slope aspect and sky view) is substantial at high resolution • Some significant indirect impacts (feedbacks) even at 7km: related to snow melt, stability (turbulence) and low clouds • Higher impact during the winter time due to sun elevation and snow conditions • Significant secondary effects on clouds cover at 2km • Surface thermal changes at high resolution have also secondary effects on thermal circulation Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Outlook • Next step: systematic validation with observations • Extension of the domain to the entire aLMo operational domain (7 km and 2 km) and start test suite • Compare implementation results with the sub-grid-scale version of the Müller and Scherrer (2005) scheme (differences in preprocessor) • Indirect effects (feedbacks) on clouds and stability (turbulence) have to be further investigated • How to solve the problem of clouds at high resolution? Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Thank you very much for your attention! matteo.buzzi@meteoswiss.ch
Implementation into the LM • Preprocessor: prior to model integration • Computes corrections factors (time dipendent) + sky view factor (time indipendent) and save results in binary file • LM reads the correction factors (parallelized) • Corrections take place in the fesft.incf radiation subroutine • In spectral loop for solar radiation (direct and diffuse) • After the spectral loop for thermal radiation: the two components (up and down) have to be recomputed from the thermal radiation balance, corrected and added to a new radiation balance Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topography: 7 km Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Topography: 2 km Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
The problem: clouds at high resolution Model: no cloud cover Reality: cloud cover Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 2 km, 2003-06-22 00 UTC +12 10 m wind vectors difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005 corrected control
Results: 2 km, 2004-12-11 00 UTC +12 Total cloud cover difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2003-06-22 00 UTC +14 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2004-12-11 00 UTC +15 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005
Results: 7 km, 2003-08-07 00 UTC +15 Surface temperature difference Matteo.buzzi@MeteoSwiss.ch COSMO Meeting, WG 3, 20.09.2005