At continental mid-latitude, soil moisture (SM) is a key component of the climate systems and land surface initialization is crucial for subseasonal-to-seasonal (S2S) predictions. We introduce a new stochastic global land reanalysis system called the Norwegian Climate Prediction Model Land (NorCPM-Land), which will be used to initialize the land component of the Norwegian Climate Prediction Model (NorCPM). We assimilate the blended SM from the European Space Agency’s Climate Change Initiative (ESA CCI) into a 30-member offline simulation of the land surface Community Land Model (CLM). Fluxes are provided by 30-member historical simulations of the full coupled NorCPM. The Ensemble Kalman Filter (EnKF) updates daily the soil column from the SM data using the cumulative density function matching method.The NorCPM-Land is currently produced for 40 years from 1980 to 2019. Assimilation significantly improves the land surface state variability and reduces error by 10.5% when validated using independent SM observations and by reanalysis estimates from ERA5-Land. It also yields an improvement of land surface energy, runoff and net primary production. We demonstrate that adjusting the underlying soil moisture considerably enhances the ability to simulate land surface state dynamics.
Atmospheric gravity waves (GWs) span a broad range of length scales. As a result, the un-resolved and under-resolved GWs have to be represented using a sub-grid scale (SGS) parameterization in general circulation models (GCMs). In recent years, machine learning (ML) techniques have emerged as novel methods for SGS modeling of climate processes. In the widely-used approach of supervised (offline) learning, the true representation of the SGS terms have to be properly extracted from high-fidelity data (e.g., GW-resolving simulations). However, this is a non-trivial task, and the quality of the ML-based parameterization significantly hinges on the quality of these SGS terms. Here, we compare three methods to extract 3D GW fluxes and the resulting drag (GWD) from high-resolution simulations: Helmholtz decomposition, and spatial filtering to compute the Reynolds stress and the full SGS stress. In addition to previous studies that focused only on vertical fluxes by GWs, we also quantify the SGS GWD due to lateral momentum fluxes. We build and utilize a library of tropical high-resolution ($\Delta x =3~km$) simulations using weather research and forecasting model (WRF). Results show that the SGS lateral momentum fluxes could have a significant contribution to the total GWD. Moreover, when estimating GWD due to lateral effects, interactions between the SGS and the resolved large-scale flow need to be considered. The sensitivity of the results to different filter type and length scale (dependent on GCM resolution) is also explored to inform the scale-awareness in the development of data-driven parameterizations.
The Tropical cyclone (TC) forecast skill of the eight global medium-range forecast models which are participating in the DIMOSIC (DIfferent Models, Same Initial Conditions) project is investigated in this study. Each model was used to generate 10-day forecasts from the same initial conditions provided by the European Centre for Medium-Range Weather Forecasts. There are a total of 123 initial dates spanning in one year from June 2018 to June 2019 with a 3-day interval. The TC track and intensity forecasts are evaluated against the best track dataset. TC-related precipitation and tropical cyclogenesis forecasts are also compared to explore the differences and similarities of TC forecasts across the models. This comparison of TC forecasts allows model developers in different centers to benchmark their model against other models, with the impact of the initial condition quality removed. The verifications reveal that most models show slow-moving and right-of-track biases in their TC track forecasts. Also, a common dry bias in TC-related precipitation indicates a general deficiency in TC intensity and convection in the models which should be related to insufficient model resolution. These findings provide important references for future model developments.
Simulating whole atmosphere dynamics, chemistry, and physics is computationally expensive. It can require high vertical resolution throughout the middle and upper atmosphere, as well as a comprehensive chemistry and aerosol scheme coupled to radiation physics. An unintentional outcome of the development of one of the most sophisticated and hence computationally expensive model configurations is that it often excludes a broad community of users with limited computational resources. Here, we analyze two configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6)) with simplified “middle atmosphere” chemistry at nominal 1 and 2 degree horizontal resolutions. Using observations, a reanalysis, and direct model comparisons, we find that these configurations generally reproduce the climate, variability, and climate sensitivity of the 1 degree nominal horizontal resolution configuration with comprehensive chemistry. While the background stratospheric aerosol optical depth is elevated in the middle atmosphere configurations as compared to the comprehensive chemistry configuration, it is comparable between all configurations during volcanic eruptions. For any purposes other than those needing an accurate representation of tropospheric organic chemistry and secondary organic aerosols, these simplified chemistry configurations deliver reliable simulations of the whole atmosphere that require 35% to 86% fewer computational resources at nominal 1 and 2 degree horizontal resolution, respectively.
The development of any atmosphere or ocean model warrants a suite of test cases to verify its spatial and temporal discretizations, order of accuracy, stability, reproducability, portability, scalability, etc. In this paper, we present a suite of shallow water test cases designed to verify the barotropic solver of atmosphere and ocean models. These include the non-dispersive coastal Kelvin wave; the dispersive inertia-gravity wave; the dispersive planetary and topographic Rossby waves; the barotropic tide; and a non-linear manufactured solution. These test cases check the implementation of the linear pressure gradient term; the linear constant or variable-coefficient Coriolis and bathymetry terms; and the non-linear advection terms. Simulation results are presented for a variety of time-stepping methods as well as two spatial discretizations: a mimetic finite volume method based on the TRiSK scheme, and a high-order discontinuous Galerkin spectral element method. We explain the strategies that need to be adopted for specifying initial and non-periodic boundary conditions on hexagonal meshes. Convergence studies of every test case are conducted with refinement in both space and time, only in space, and only in time. The convergence slopes match the expected theoretical predictions.