Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

The Design Document page provides a description of the algorithms, implementation and planned testing including unit, verification, validation and performance testing. Please read  Step 1.3 Performance Expectations that explains feature documentation requirements from the performance group point of view. 

Design Document

Expand
titleClick here for instructions to fill up the table below ......

The first table in Design Document gives overview of this document, from this info the Design Documents Overview page is automatically created.

In the overview table below 4.Equ means Equations and Algorithms, 5.Ver means Verification, 6.Perf - Performance, 7. Val - Validation

  • Equations: Document the equations that are being solved and describe algorithms
  • Verification Plans: Define tests that will be run to show that implementation is correct and robust. Involve unit tests to cover range of inputs as well as benchmarks.
  • Performance expectations: Explain the expected performance impact from this development
  • Validation Plans: Document what process-based, stand-alone component, and coupled model runs will be performed, and with what metrics will be used to assess validity

Use the symbols below (copy and paste) to indicate if the section is in progress or done or not started.

...

Page Properties
idFeature_PR


Info

Overview table for the owner and an approver of this feature

1.Description

SNICAR-AD Shortwave Radiation for Sea Ice and Land

2.Owner

Cheng Dang (Unlicensed)
3.Created 
4.Equ(tick)
5.Ver(warning)(tick)
6.Perf(warning)(tick)
7.Val(warning)(tick)
8.Approver
9.Approved Date
V2.0Pending



...

Expand
titleClick here for Table of Contents ...


Panel

Table of Contents

Table of Contents





Title: SNICAR-AD Shortwave Radiation for Sea Ice and Land

Requirements and Design

E3SM Cryosphere Group

Date:  


Summary


The goal of this task is to improve remove the discrepancy of shortwave radiative transfer simulation of snow-covered sea ice by adopting the updated radiative for snowpacks on sea ice and on land in E3SM. We implement the radiative transfer scheme SNICAR-AD in the sea-ice component MPAS-seaice . A parallel modification has also been made to and the land component ELM , together they provide consistent for E3SM to predict consistent and improved shortwave radiative transfer computation for properties of snow-covered surfaces.


Requirements

...

ELM: modify the land shortwave code to include a subroutine that adopts the new radiative transfer algorithm.

MPAS-seaice: modify the sea-ice shortwave radiative transfer code to include a new subroutine SNICAR_AD that

  1. provides new input files for snow and ice optics adopted by ELM;
  2. expands the shortwave radiative calculation from three to five subdivided solar bands adopted by ELM;
  3. removes empirical adjustment for radiative-effective snow grain size;
  4. adopts a new algorithm to correct for the bias of near-infrared albedo and absorption when solar zenith angles are larger than 75 degrees.


Date last modified: Jun 1st, 2019 


Algorithmic Formulations

Design solution: described in the following article

Dang, C., Zender, C. S., and Flanner, M. G.: Inter-comparison and improvement of 2-stream shortwave radiative transfer models for a unified treatment of cryospheric surfaces in ESMs, The Cryosphere Discuss., https://doi.org/10.5194/tc-2019-22, in review, 2019.

Date last modified: July 1st, 2019
Contributors: Cheng Dang (Unlicensed) Charlie Zender Mark Flanner




Design and Implementation

Implementation: 

In the sea-ice and land shortwave code (ice_shortwace.F90, we add a new subroutine SnowSnicarMod.F90), we implemented SNICAR_AD to compute the shortwave radiative properties of snow-covered sea ice, which that can be turned on by setting config_use_snicar_ad = true in the namelist.Source /off by configuration toggles in the model namelists.

MPAS-source code: https://github.com/chengdang/MPAS-Model/tree/seaice/snowice_sw_snicar. The code modification is merged to the sea-ice developing branch, with default config_use_snicar_ad = false.

The pseudocode in ice_shortwave.F90

if snow-covered sea ice:
    if (config_use_snicar_ad):
          call snicar_ad
    else 
          call default shortwave algorithm
    end
end
subroutine snicar_ad

Date last modified: July 1st, 2019

Contributors: Cheng Dang (Unlicensed)




Planned Verification and Unit Testing 

Verification and Unit Testing:

  • The offline radiative transfer model has been tested against a multi-stream benchmark model. (Dang, C., Zender, C. S., and Flanner, M. G.: Inter-comparison and improvement of 2-stream shortwave radiative transfer models for a unified treatment of cryospheric surfaces in ESMs, The Cryosphere Discuss., https://doi.org/10.5194/tc-2019-22, in review, 2019.)


Date last modified: July 1st, 2019

Contributors: Cheng Dang (Unlicensed)

Results of the




Planned Validation Testing 

Validation Testing: 

Date last modified: July 1st, 2019

Contributors: Cheng Dang (Unlicensed)

Planned Validation Testing 

Validation Testing: 

  • Results of the modified E3SM model are tested against observation provided by mpas-analysis and e3sm diagnostics.
  • Results of the modified E3SM model will be tested against observations/reanalysis/model simulations of snow depth, sea-ice concentration/fraction, dust and black carbon concentration (Fully-coupled historical simulations are in progress).
    • BC/Dust in snow: 
      • The Arctic: Doherty et al., 2010
      • China: Wang et al., 2013
      • North America: Doherty et al., 2014
      • Data summary: Dang et al., 2017
    • Sea-ice concentration/extent: PIOMAS
    • Snow depth: to be added.


Date last modified: July 1st, 2019
Contributors: Cheng Dang (Unlicensed)



Planned Performance Testing 

Performance Testing: 

The new implementation is not expected to have a significant impact on the computational cost/performance of E3SM.

Date last modified: July 1st, 2019

Contributors: Cheng Dang (Unlicensed)