Gemclim 3.3.3.1 nudging : Différence entre versions

Un article de Informaticiens département des sciences de la Terre et l'atmosphère
Aller à: navigation, charcher
m (References:)
m (References:)
 
(16 révisions intermédiaires par le même utilisateur non affichées)
Ligne 2: Ligne 2:
  
 
Spectral nudging is a way to force (nudge) the larger scales of an RCM towards the ones of its driving GCM resp. RCM.<br>  
 
Spectral nudging is a way to force (nudge) the larger scales of an RCM towards the ones of its driving GCM resp. RCM.<br>  
 +
<br>
  
The information below is for '''CRCM5/GEMCLIM v_3.3.3''' and '''v_3.3.3.1''' only!<br>Older version of CRCM5/GEMCLIM do not use the same spectral nudging, but a veeeery costly one which is controlled differently and might or might not work. - Sorry.
+
=== Spectral nudging parameters in GEM:<br>  ===
 
 
=== References:<br>  ===
 
 
 
*Riette, S. (2000): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2000_Riette_Pilotage_Spectral.pdf '''Technique de Pilotage Spectral du Modele Regional Canadien de Climat''']
 
  
*Šeparović L, de Elía R, Laprise R (2012) Impact of spectral nudging and domain size in studies of RCM response to parameter modification. Clim Dyn 38(7-8): 1325-1343 [http://link.springer.com/article/10.1007%2Fs00382-011-1072-7 http://link.springer.com/article/10.1007%2Fs00382-011-1072-7]<br>'''Exerpt:'''<br>"Nested limited-area Regional Climate Models (RCMs) are models that dynamically downscale global General Circulation Model (GCM) simulations or objective analyses to high-resolution computational grids, using a high-resolution representation of the surface forcing and model dynamics. RCMs require the information on some prognostic variables as their lateral boundary conditions (LBC). The choices of integration domains and nesting techniques are free parameters of RCMs. The optimal integration domain depends on the particular situation, although there are some general recommendations that can facilitate user’s judgment (e.g., Laprise et al. 2008). For example, Leduc and Laprise (2008) showed that the use of a too small domain could result in the simulations being deficient in fine-scale variance. It has been also noted that in large continental-scale domains RCM large-scale variables can considerably drift from the driving fields, which can then result in appearance of large spurious gradients in the vicinity of the outflow boundaries. Spectral nudging (SN; von Storch et al. 2000; Biner et al. 2000) has been employed to ensure that the model solution remains close to the large-scale components of the driving fields over the entire domain. However, the use of spectral nudging remains an open issue. …"
+
&nbsp;&nbsp;&nbsp;[[spn GEM3.3.3|GEM3.3.3.1 (CRCM5)]]<br>
 
+
&nbsp;&nbsp;&nbsp;[[spn GEM4.8.lts12|GEM4.8.lts12]]
*von Storch H, Langenberg H, Feser F (2000): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2000_vonStorch_Spectral_Nudging.pdf '''A spectral nudging technique for dynamical downscaling purposes.'''] Mon Weather Rev 128: 3664–3673
 
*Biner S, Caya D, Laprise R, Spacek L (2000) Nesting of RCMs by imposing large scales. Research activities in Atmospheric and Oceanic Modelling, WMO/TD 987, Report 30: 7.3-7.4
 
  
 
<br>
 
<br>
  
== Grid size restrictions<br>  ==
+
=== References:<br> ===
 
 
To be able to use spectral nudging in CRCM5, the '''core''' '''grid needs to match the FFT criteria in x- and y-direction'''. You can use 'findfft' to get the possible values to set 'Grd_ni' and 'Grd_nj' in your gemclim_settings.nml. (You need to set the model environment first to have access to 'find'fft'.)<br>For example:  
 
 
 
&nbsp;&nbsp;&nbsp; findfft -gnimax 500 -cfl ''Pil_maxcfl''
 
 
 
where the key '-cfl' needs to be set to value of 'Pil_maxcfl' from your gemclim_settings.nml!!!!!<br> Your 'Grd_ni' and 'Grd_nj' will then have to be set to one of the values findfft is listing.<br> If your grid is larger than 500 points, just set '-gnimax' to a number large enough.
 
 
 
== Activate spectral nudging in CRCM5  ==
 
  
To activate spectral nudging in CRCM5 one has to set several parameters in the file 'gemclim_settings.nml' in the namelist 'gem_cfgs'.<br>
+
*'''Xu''' Z. and Z.-L. Yang (2015): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2015_Xu_Spectral_Nudging_Bias_Correction.pdf '''A new dynamical downscaling approach with GCM bias corrections and spectral nudging.'''] J. Geophys. Res. Atmos., 120, 3063–3084
  
=== Fields to nudge<br>  ===
+
*'''Berg''' P., R. Döscher, and T. Koenigk (2013): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2013_Berg_RCAO_nudging.pdf '''Impacts of using spectral nudging on regional climate model RCA4 simulations of the Arctic'''] Geosci. Model Dev., 6, 849–859, 2013
  
In CRCM5 one has the option of nudging u-wind, v-wind, and/or temperature.<br>Tests have shown that it is sufficient to nudge u-wind and v-wind and that nudging the temperature as well does not change much the results. Therefore I&nbsp;suggest to set:
+
*'''Wang''' J., and V.R. Kotamarthi (2013): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2013_Wang_Kotamarthi_Nudging.pdf '''Assessment of Dynamical Downscaling in Near-Surface Fields with Different Spectral Nudging Approaches Using the Nested Regional Climate Model (NRCM)'''] JAMC Vol 52, 1576-1591
  
&nbsp; &nbsp; '''Lam_spn_nudging_S = 'UV'&nbsp;&nbsp;&nbsp;&nbsp; ,'''  
+
*'''Bowden''' J.H., C.G. Nolte, and T.L. Otte (2012): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2013_Bowden_WRF_Spectral_nudging.pdf '''Simulating the impact of the large-scale circulation on the 2-m temperature and precipitation climatology.'''] Clim Dyn (2013) 40:1903–1920
  
'''One needs to set this parameter to switch the spectral nudging on!'''
+
*'''Šeparović''' L., R. de Elía, and R. Laprise (2012) [http://link.springer.com/article/10.1007%2Fs00382-011-1072-7 '''Impact of spectral nudging and domain size in studies of RCM response to parameter modification.'''] Clim Dyn 38(7-8): 1325-1343 <br>'''Exerpt:'''<br>"Nested limited-area Regional Climate Models (RCMs) are models that dynamically downscale global General Circulation Model (GCM) simulations or objective analyses to high-resolution computational grids, using a high-resolution representation of the surface forcing and model dynamics. RCMs require the information on some prognostic variables as their lateral boundary conditions (LBC). The choices of integration domains and nesting techniques are free parameters of RCMs. The optimal integration domain depends on the particular situation, although there are some general recommendations that can facilitate user’s judgment (e.g., Laprise et al. 2008). For example, Leduc and Laprise (2008) showed that the use of a too small domain could result in the simulations being deficient in fine-scale variance. It has been also noted that in large continental-scale domains RCM large-scale variables can considerably drift from the driving fields, which can then result in appearance of large spurious gradients in the vicinity of the outflow boundaries. Spectral nudging (SN; von Storch et al. 2000; Biner et al. 2000) has been employed to ensure that the model solution remains close to the large-scale components of the driving fields over the entire domain. However, the use of spectral nudging remains an open issue. …"
  
== Control the strength of the spectral nudging  ==
+
*'''Mohino''' E., S. Janicot, H. Douville, L.Z.X. Li (2011): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2012_Mohino_Nudging.pdf '''Impact of the Indian part of the summer MJO on West Africa using nudged climate simulations.'''] Clim Dyn (2012) 38:2319–2334
  
There are 3 aspects controlling the strength of the spectral nudging which will all get applied at each time step, one after the other:<br>
+
*'''Alexandru''' A. (2009): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2009_Adelina_SN.pdf '''Sensitivity Study of Regional Climate Model Simulations to Large-Scale Nudging Parameters.'''] Mon Weather Rev 137: 1666-1686
  
#a filter in spectral space (horizontal)<br>
+
*'''Omrani''' H., P. Drobinski, and T. Dubos (2009): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2011_Drobinski_QJRMS_spectral_nudging.pdf '''Spectral nudging in regional climate modeling: how strong should we nudge?'''] Q.J.R. Meteorol. Soc. 00: 1–6 (2009)
#the shape of the vertical profile<br>
 
#the intensity of nudging applied at each time step<br>
 
  
=== 1) Half response wavelength<br>  ===
+
*'''Radu''' R., M. Déqué, and S. Somot (2008): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2008_Radu_ALADIN_Nudging.pdf '''Spectral nudging in a spectral regional climate model'''] Tellus (2008), 60A, 898–910
  
This is the wavelength at which equal weight is given to the driving data and the RCM solution. For longer wavelengths more weight is given to the driving data, for shorter wavelengths more weight is given to the RCM. A recommanded value is about 10 or 15 times the grid point size of the driving data. The half response wavelength needs to be given in kilometers. <br>
+
*'''Riette''' S., and D. Caya (2002): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2002_Riette_Caya.pdf '''Sensitivity of short simulations to the various parameters in the new CRCM spectral nudging. Research Activities in Atmospheric and Oceanic Modelling'''] WMO/TD-1105, Rep. 32, H. Ritchie, Ed., WMO, 7.39–7.40
  
For example set:<br>
+
*'''Riette''' S. (2000): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2000_Riette_Pilotage_Spectral.pdf '''Technique de Pilotage Spectral du Modele Regional Canadien de Climat''']
  
'''&nbsp; &nbsp; Lam_spn_half_respon = 1000.&nbsp; ,'''<br>
+
*'''Biner''' S., D. Caya, R. Laprise, L. Spacek (2000) Nesting of RCMs by imposing large scales. Research activities in Atmospheric and Oceanic Modelling, WMO/TD 987, Report 30: 7.3-7.4
  
The default is 1177.
+
*'''von Storch''' H., H. Langenberg, and F. Feser (2000): [http://pascal.sca.uqam.ca/~winger/Documents/Spectral_Nudging/2000_vonStorch_Spectral_Nudging.pdf '''A spectral nudging technique for dynamical downscaling purposes.'''] Mon Weather Rev 128: 3664–3673
 
 
http://pascal.sca.uqam.ca/~winger/Images/spn_half_respn.png
 
 
 
=== 2) Controlling the vertical profile<br>  ===
 
 
 
There are 3 parameters to control the nudging profile:
 
 
 
&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev&nbsp;&nbsp;&nbsp;&nbsp;: above this level (in hybrid values) the profile is set to 1, to full strength nudging (default 0., model top).<br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp;: below this level (in hybrid values) the profile is set to 0, no nudging (default 1., surface).<br>&nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S&nbsp;: Between these two hybrid levels the profile goes from 0 to 1. Either linearly (''''LINEAR'''') (default) or with a cos2 (''''COS2'''') function. <br>
 
 
 
For example set:
 
 
 
'''&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev &nbsp; = 0.&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp; &nbsp; , <br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; = 0.5 &nbsp; &nbsp;&nbsp; &nbsp;&nbsp; ,<br> &nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S = 'LINEAR' ,'''
 
 
 
<br> http://pascal.sca.uqam.ca/~winger/Images/spn_profile.png
 
 
 
=== 3) Controlling the strength  ===
 
 
 
The spectral nudging gets applied at each time step. The difference between RCM and driving data gets calculated (for different wave lengths) and with that an increment, by which the model results would need to get nudged to match the driving data. But this increment is usually not fully applied, but only a fraction of it. Since the idea is to "nudge" the large scales of the model towards the ones of the driving data and not to impose them. <br><br>To control how strongly the model is nudged, one does not say "I want to apply n% of this increment", but one rather says "I want this increment to be applied over h hours". So if, for example, one wants to apply only 5% of the increment and the model time step is 30 minutes, the time needed to apply this increment fully (if applied at each time step) would be 10 hours. So in this case one would set:<br>
 
 
 
&nbsp;&nbsp;&nbsp; '''Lam_spn_relax_hours = 10. ,'''
 
 
 
10 is also the default value.
 
  
 
<br>
 
<br>
 
== Modifications to your executable  ==
 
 
Unfortunately, the spectral nudging in GEMCLIM v_3.3.3 and v_3.3.3.1 is not exactly doing what it is supposed to do.<br>You therefore need a patch.
 
 
If you are running '''v_3.3.3''', please talk to me (winger.katja@uqam.ca).
 
 
If you are running '''v_3.3.3.1''', please follow the instructions below:<br>
 
 
#Set the model environment (with '''3331'''&nbsp;&nbsp; resp.&nbsp;&nbsp; . s.ssmuse.dot gemclim_3.3.3.1-dev )<br>
 
#Go into the directory in which you create your executables
 
#Get a copy of the corrected function:<br>&nbsp;&nbsp;&nbsp; '''cp ${gemclim}/Sources/Patches/spn_calfiltre.ftn .'''
 
#Recreate the Makefile etc:'''<br>&nbsp;&nbsp; r.make_exp'''
 
#Compile the corrected function:<br>&nbsp;&nbsp;&nbsp; '''make spn_calfiltre.o'''
 
#Link and load the model executable:<br>&nbsp;&nbsp;&nbsp; '''make gemclimdm'''<br>
 
 
Then launch your simulation as usual.<br>
 

Version actuelle datée du 12 de février 2019 à 19:18

Spectral nudging

Spectral nudging is a way to force (nudge) the larger scales of an RCM towards the ones of its driving GCM resp. RCM.

Spectral nudging parameters in GEM:

   GEM3.3.3.1 (CRCM5)
   GEM4.8.lts12


References:

  • Šeparović L., R. de Elía, and R. Laprise (2012) Impact of spectral nudging and domain size in studies of RCM response to parameter modification. Clim Dyn 38(7-8): 1325-1343
    Exerpt:
    "Nested limited-area Regional Climate Models (RCMs) are models that dynamically downscale global General Circulation Model (GCM) simulations or objective analyses to high-resolution computational grids, using a high-resolution representation of the surface forcing and model dynamics. RCMs require the information on some prognostic variables as their lateral boundary conditions (LBC). The choices of integration domains and nesting techniques are free parameters of RCMs. The optimal integration domain depends on the particular situation, although there are some general recommendations that can facilitate user’s judgment (e.g., Laprise et al. 2008). For example, Leduc and Laprise (2008) showed that the use of a too small domain could result in the simulations being deficient in fine-scale variance. It has been also noted that in large continental-scale domains RCM large-scale variables can considerably drift from the driving fields, which can then result in appearance of large spurious gradients in the vicinity of the outflow boundaries. Spectral nudging (SN; von Storch et al. 2000; Biner et al. 2000) has been employed to ensure that the model solution remains close to the large-scale components of the driving fields over the entire domain. However, the use of spectral nudging remains an open issue. …"
  • Biner S., D. Caya, R. Laprise, L. Spacek (2000) Nesting of RCMs by imposing large scales. Research activities in Atmospheric and Oceanic Modelling, WMO/TD 987, Report 30: 7.3-7.4