Recommended for Plenary | Recommended for Group Session | ACME or Collab. | Speaker/Affiliation | Title | Abstract | |
|---|---|---|---|---|---|---|
| 1 | collab. | Matthew Long
Matthew Long <mclong@ucar.edu> | MARBL talk: Ocean biogeochemistry in the Earth system modeling framework: applications and approach | Marine biogeochemical cycles play a fundamental role regulating climate, most directly through impacts on atmospheric carbon dioxide. This has motivated the development of ocean biogeochemistry modules as components of Earth system models, primarily for the purpose of simulating changes in ocean sinks for carbon dioxide under future emissions scenarios. These models, however, can be applied to a range of other interesting problems related to marine ecology and biogeochemistry. In this talk, I provide a brief overview of some research questions in ocean biogeochemistry and ecology that I find interesting. I then describe our efforts to develop the Marine Biogeochemistry Library (MARBL), which is a modular implementation of ocean biogeochemistry that aims to be flexible and capable of operating within different physical frameworks---thereby enabling research across a broad array of questions. | ||
| 2 | collab. | Michal A. Kopera (UCSC), Wieslaw Maslowski (NPS), Francis X. Giraldo (NPS)
Michal A. Kopera <makopera@ucsc.edu> | A new ice sheet / ocean interaction model for Greenland Fjords using discontinuous Galerkin method | One of the key outstanding challenges in modeling of climate change and sea-level rise is the ice-sheet/ocean interaction in narrow, elongated and often geometrically complicated fjords around Greenland. The goal of the Fjord-DG (FDG) project is to build a separate, highresolution module for use in Earth System Models (ESMs) to realistically represent the fjord bathymetry and coastlines and the fine-scale processes occurring within the fjord and at the ice shelf interface, using discontinuous Galerkin (DG) methods. FDG is currently at the first stage of development. We used NUMA (Non-hydrostatic Unified model of the Atmosphere) framework to develop the incompressible Navier-Stokes equation (INSE) solver, which will be used as a dynamical core in the fjord ocean model. We will present some preliminary results of idealized INSE test cases, and discuss further avenues of the project progress. The key features of the FDG module will be high-order accuracy, geometrical flexibility and nonconforming adaptive mesh refinement to resolve the processes occurring near the ice-sheet/ ocean interface without introducing prohibiting computational cost. The non-hydrostatic model will account for the stationary ice-shelf with sub-shelf ocean interaction, basal melting and subglacial meltwater influx, with boundary conditions at the surface to account for floating sea ice. The boundary conditions will be provided to the model via CPL7 coupler to emulate the integration with ESM. FDG will be tested initially on Sermilik Fjord using real bathymetry, boundary and initial conditions, and evaluated against observations and other model results for this fjord. The overarching goal of the project is to be able to resolve the ice-sheet/ocean interactions around the entire Greenland’s coast and two-way couple with climate models like ACME. | ||
| 3 | collab. | Michael Prather, Juno Hso, Alex Viedenbaum, Alex Nicolau; UC Irvine | Solar-J - RRTMG comparisons An examination of systematic biases in ACME solar heating rates - comparison with multi-stream Solar-J | As part of the development of a combined photolysis-heating module for DOE's ACME, we have taken the better-resolved photolysis code (Cloud-J, < 778 nm) and merged it with the RRTMG bins longward of 778 nm. For clear-sky calculations, the results are similar as expected, but for cirrus and stratus decks the 2-stream models in RRTMG produce significant biases that vary as a function of solar zenith angle and cloud optical depth when compared with the 8-stream Solar-J model. This is expected, but the size is significant, and may affect the mean meteorology. We present comparisons of the different options for solar heating codes used in CESM and ACME with those from the newly developed Solar-J. | ||
| 4 | collab. | Xianglei Huang (the University of Michigan) | Incorporate realistic spectral emissivity of surfaces into the CESM and the influence on simulated radiation budget, mean climate, and climate changes. | The actual surface emissivity is spectrally dependent, for both oceans and lands. Like many state-of-the-art GCMs, the atmospheric radiation scheme in the CESM still assumes the surface being blackbody. A few recent studies showed that the surface spectral emissivity could have impact on the simulated radiation budget and climate, which motivated this study. First, we developed a global surface spectral emissivity dataset suitable for the GCM and NWP models. The dataset is based on the first-principle calculation for the far-IR. It is anchored on MODIS surface emissivity retrievals for the mid-IR and validated against IASI retrieved surface spectral emissivity. Using LBLRTM, we also carry out a benchmark study to understand the errors in the RRTMG_LW scheme when surface spectral emissivity is included.
We implement this dataset into the CESM and modify code to ensure consistency between land module and atmosphere module of the CESM. Comparing to the standard CESM, the slab-ocean run shows a change of 1.9 Wm-2 for the global-mean TOA net radiative flux (downward positive). The atmospheric column net radiativing cooling changes by -0.6 Wm-2. The changes have distinct spatial patterns. The results from 20-year fully coupled CESM run show similar but generally smaller changes than the slab-ocean run. The impact on simulated precipitation and surface temperature will be also discussed. | ||
| 5 | collab. | Michael S. Pritchard, Hossein Parishani, Mathew C. Wyant, Marat Khairoutdinov, Balwinder Singh, Christopher S. Bretherton Mike Pritchard [mailto:mspritch@uci.edu] | Towards low cloud-permitting cloud superparameterization. | Results are shown from a prototype “ultra-parameterized” version of the Community Atmosphere Model. In this approach, O(10k) embedded cloud resolving models (CRMs) are used within a 2-degree global climate model (i.e. a Multiscale Modeling Framework, or MMF) but with radically refined interior resolution relative to standard superparameterization — with CRM horizontal resolution approaching 250 m and vertical resolution approaching 20 m in the marine boundary layer inversion zone. Philosophically, the goal is to explicitly capture the outer scales of turbulence involved in marine boundary layer (MBL) dynamics towards a more robust treatment of global low cloud feedback in climate models. To make the approach computationally feasible, an algorithm for accelerating the mean state evolution of CRMs is implemented, as well as software engineering for enhanced parallel scalability and GPU co-processing. Pilot hindcast tests reveal a more satisfying representation of MBL vertical decoupling dynamics and surface fluxes but also an unintended consequences of chronically low liquid water concentration linked to overentrainment symptoms, inconsistent with offline LES benchmark studies. Fixing an error in how surface fluxes are transmitted to CRMs in the MMF helps recover some missing coastal cloud fraction and using CRM grids with an extreme aspect ratio (~200) helps recover much of the missing LWP in association with modified turbulent kinetic energy statistics. Overall this highlights emerging issues and potential tuning strategies at a previously unencountered frontier within the grey zone of quasi-resolved turbulence for multi-scale climate modeling. | ||
| 6 | collab. | Ethan Butler [mailto:eebutler@umn.edu] | Estimating Global Maps of Trait Distributions | Plant traits have been measured at diverse sites, but there are still wide areas with minimal or non-existent measurements. This spatial limitation necessitates some method of interpolation to construct continuous trait surfaces for use in independent modeling or as input to global land surface models. Here we present two distinct methodologies to leverage a sub-set of the global database of plant traits, TRY, to create continuous maps of trait distributions for leaf nitrogen, specific leaf area, and mass based maximum photosynthesis rate. First, a categorical method, leveraging classification of the species observed in TRY and satellite estimates of their plant functional type abundances - analogous to how traits are currently assigned to PFTs in land surface models. Second, a spatial statistical method which additionally estimates how the distribution of a trait changes in accord with both climate and soil covariates. These methods produce global maps of full trait distributions at a resolution commensurate with many other observational datasets and Earth System Models and may be used as input or cross validation for such models. | ||
| 7 | collab. | François Primea and Keith Moore, UCI | Developing fast implicit spin-‐up capabilities for CESM to study the dynamics o the ocean’s nitrogen and carbon cycles | We are investigating the impacts of non-‐Redfield plankton stoichiometry on the marine carbon cycle using CESM. We tested two methods to generate non-‐Redfield patterns (th Redfield N/P ratio of 16/1 is used i th standard CESM):
We are evaluating the different formulations by comparing the simulations to field observations o [PO3 an N = [NO3]-‐16[PO4]. The model evaluation is still ongoing, but is showing promising results, which we will present along with comparisons to recent inverse model results. We will also present an update on the development of the fast implicit biogeochemical tracer solver and our plan for using the solver to evaluate the impacts of non-‐Redfield stoichiometry on the oxygen-‐, phosphorous-‐, nitrogen-‐, and carbon-‐cycle dynamics of the ocean. | ||
| 8 | collab. | Chaopeng Shen, William J. Riley and John M. Melack Civil and Environmental Engineering, Pennsylvania State University | Scale-aware hydrologic and biogeochemical modeling for the Amazon and the world: model enhancement, multi-scale strategies and dataset generation | Modeling hydrology and biogeochemistry faces several prominent challenges including insufficient subsurface flow representation, hydrologic and coupled biogeochemical scaling issue and lack of high quality verification datasets. Our project addresses these issues with a rich set of activities which also target generic improvements in hydrologic modeling in future Earth System modeling, especially in data-poor regions. We have made progress in (a) model improvement: we implemented a Global Process-based Adaptive Watershed Simulator with Community Land Model, (PAWS+CLM4.5) package with Fortran-native reflection capabilities to enable detailed, process-based hydrologic modeling in the world and convenient modeling coupling; (b) multi-scale methods: we investigated several scaling methods including moment matching, reduced-order modeling and soil-moisture-fractal-based scaling, and are testing a perturbed prototype scaling method which seems promising; (c) channel-land interactions: we created data extraction algorithm for physically-based channel modeling, and examined the influence of channel processes on simulated water and carbon fluxes and states; (d) Amazon water budgets: we modeled hydrology in several basins in the Amazon, which highlighted the importance of groundwater flow in streamflow and evapotranspiration (ET) seasonality; and (e) production of “hydrologic-model-free” verification data using satellite products, including a GRACE-assisted-Budyko-curve global ET product and a GRACE-based streamflow change estimate (we focus on GRACE as it is not sensitive to the heavy canopy in the Amazon, unlike MODIS and SMAP). These activities have resulted in a list of publications and will serve to broaden the capabilities of the Earth System models. | ||
| 9 | ACME | Bryce Harrop Phil Rasch (pnl.gov), Po-Lun Ma Atmosphere | Evaluating monsoon circulations in ACME v1 ne30 experiments | The monsoon circulations are the dominant mode of seasonal variability in the tropics. Marked by its strong annual cycle in precipitation, the global monsoon system brings water to approximately half of the world’s population. Despite the importance of the monsoon in the current climate, several shortcomings still exist in our understanding as well as our ability to accurately model monsoon behavior. As part of the development of the Accelerated Climate Modeling for Energy model version 1 (ACME v1), it is important to evaluate the model’s ability to represent the monsoon (both globally and regionally). To evaluate the monsoon circulation within ACME v1, we compare with an older version of the model (ACME v0), another commonly used climate model (the Community Atmosphere Model), and a number of observational datasets. Several important metrics have been identified from the literature as important markers of the monsoon circulation and these are used to quantify the evaluation of the monsoon. Globally, the monsoon simulated by ACME v1 compares well with observations, but regional biases remain. | ||
| 10 | ACME | Land | Hypsometric analysis improves topography-based subgrid structures for the ACME Land Model | Topography exerts a major control on land surface processes through its influence on atmospheric forcing, soil and vegetation properties, network topology and drainage area. Land surface spatial structure that captures spatial heterogeneity influenced by topography is expected to improve representation of land surface processes in land surface models. For example, land surface modeling using subbasins instead of regular grids as computational units has demonstrated improved scalability of simulated runoff and streamflow processes. Two methods (Global and Local) are applied to derive new land surface spatial structures by further dividing subbasins into subgrid units based on topographic properties to take advantage of the emergent patterns and scaling properties of atmospheric, hydrologic, and vegetation processes in land surface models. The Global method utilizes the elevation classification scheme employed in Leung and Ghan (1995; 1998) combined with classifications of topographic slope and aspect to discretize each subbasin into multiple subgrid units. While, in the Local method, each subbasin is divided into multiple subgrid units using elevation classes derived based on hypsometric characteristics combined with classes of topographic aspect. In this study, the relative merits of using hypsometric characteristics in deriving topography-based subgrid structures are evaluated over the topographically contrasting regions of the Northwestern United States. Results highlight the relative advantages of the Local method over the Global method in capturing topographic heterogeneity and spatial patterns of atmospheric forcing and land cover. | ||
| 11 | ACME | Khachik Sargsyan, Daniel Ricciuto Land | Quantifying the Impacts of Parametric Uncertainty on Biogeochemistry in the ACME Land Model | A surrogate construction is a routine approach for highly expensive models enabling studies that otherwise require an infeasible number of model evaluations. Polynomial chaos machinery is a convenient tool for representing uncertain parametric inputs, propagating them through a model of interest, as well as for surrogate construction. However, large number of input parameters pose significant challenges. We develop advanced UQ and machine learning methods for construction of high-dimensional model surrogates. In particular, this poster will highlight Weighted Iterative Bayesian Compressed Sensing algorithm that enables efficient surrogate construction and uncertainty decomposition. We applied the technique to the ACME Land Model, focusing on biogeochemistry, for several output QoIs at nearly 100 FLUXNET sites covering multiple plant functional types and climates, varying 65 input parameters over broad ranges of possible values. | ||
| 12 | ACME | Jinyun Tang, William Riley (Unlicensed) Land | Huge divergence in land-atmosphere carbon exchange resulting from ambiguous numerical coupling between carbon and nitrogen dynamics | With the V0 ACME land model, we analyzed the land-atmosphere carbon exchange as simulated with carbon and nitrogen dynamics coupled using three legitimate but different numerical implementations of nitrogen limitation: (1) mineral nitrogen based limitation (MNL), (2) net nitrogen uptake based limitation (NUL), and (3) proportional nitrogen flux based limitation (PNL). For the 1850-2000 period, the three approaches resulted in very similar global distributions of carbon and nitrogen stocks, and many almost overlapping mass and energy fluxes. However, a strong divergence occurred in the simulated land-atmosphere carbon exchange for the 2001-2300 period under the RCP4.5 atmospheric CO2 forcing. Quantitatively, this divergence is as large as that of the CMIP5 models by 2100 and is about ~1000 Pg C by year 2300. Detailed analysis indicates that this divergence resulted from (1) the MNL, NUL, and PNL schemes predict progressively weaker nitrogen limitation, so that the PNL scheme leads to higher nitrogen loss through aerobic and anaerobic denitrification and surface and subsurface hydrological transport and (2) the usually high carbon to nitrogen ratio. Therefore, considering that the ratios of carbon to nutrients (other than nitrogen) are often high, we expect small inconsistencies in imposing nutrient limitation would likely lead to large divergence in predicted ecosystem carbon stocks. Since the MNL scheme was used in both V0 and V1 versions of the ACME land module, we recommend next versions to use more robust numerical coupling (e.g., PNL scheme) between carbon and nutrients. A paper on this work has been submitted for review. | ||
| 13 | ACME | Hui Wan, Kai Zhang, Phil Rasch (pnl.gov), Balwinder Singh, Xingyuan Chen, Jim Edwards | Non-bit-for-bit solution reproducibility: a new test based on time step convergence | |||
| 14 | ACME | |||||
| 15 | ACME | |||||
| 16 | ACME |