Topo: New topography file (with no "leaking" into ocean")

[cecilehannay]

There was a lot of relevant info in the run request using the new topography file (with no "leaking" into ocean"). We are saving this info into a new discussion.

Description of the run: #117

We have produced a new topography file where we use a Laplacian smoother (instead of the default distance weighted smoother) and turn off smoothing where the land fraction is 0 (i.e. no smoothing over ocean points; partial ocean points are still smoothed). We are curious if not smoothing over ocean points has an impact in coupled model setup since the temperatures passed to the ocean, where previously the topography was "leaking into the ocean", could be biased cold (assuming temperature decreasing with height on average).

The Figure below shows energy spectra for different topographies:

A couple of observations regarding the energy spectra:

• The energy spectra is the same when using a 3km intermediate cubed-sphere field (ncube3000, orange line) compared to using a 16km intermediate cubed-sphere grid (ncube540, red line) when computing energy spectra on 1 degree finite-volume grid. That means, we can assess smoothing using the ncube540 intermediate cubed-sphere grid which is much more computationally efficient to run the smoother on compared to the ncube3000 intermediate cubed-sphere grid.
• Interpolation matters for computing energy spectra: bi-linear interpolation from 1 degree spectral-element grid to 1 degree finite-volume grid produces damping (compare blue and red lines). Please note that the two grid "see" the exact same smoothed topography on the intermediate cubed-sphere grid. Hence when comparing different levels of smoothing always use the same target grid.
• The energy spectra with new and old smoother are very similar (compare red and green lines).
• Not smoothing over ocean does add some energy at many scales (compared cyan and red).
• Black line is unsmoothed topography (for reference)

The next two plots show Greenland area and South America topographies (all mapped to fv 1 degree grid):

• ncube540 raw: raw topography on ~16km intermediate cubed-sphere (ncube540)
• CESM2: Current default topography
• nc540: Laplacian smoothed topography (nu=28E7) - old CESM1.5 smoothing level when the SE topography was smoothed using SE (and not in the topography software as done now; which is needed for the ridge scheme variables)
• nc540;Co11: Laplacian smoothed topography (nu=20E7) that matches C011 smoothing (default CESM2) on the ncube540 grid (for the ncube3000 grid that is equivalent to C060 smoothing radius)
• nc540; no-ocn: Laplacian smoothed topography (nu=20E7) that matches C011 smoothing (default CESM2) but no smoothing over ocean.

The cross section plots (lower right) are though 30S and 65N.

The following observations are made:

• Smoothing over peaks seems to be roughly the same for all smoothing options -> not smoothing over ocean does not seem to affect peaks a couple of grid-points inland.
• One can visually see no smoothing over ocean.

Global plot from @adamrher showing difference between new PHIS minus the current namelists_defaults file:

As expected lower topography near coastlines. However, another difference is the use of bedmachine topography dataset so that added changes over Greenland.

Some mathematical details on the smoothing operator:

Topo file

bnd_topo = "/glade/p/cgd/amp/pel/topo/files/ne30pg3_gmted2010_bedmachine_nc3000_Laplace0100_20220516.nc"

Contact

@PeterHjortLauritzen
@JulioTBacmeister
@adamrher