We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR (Seamless System for Prediction and EArth System Research) is built from component models recently developed at GFDL ‐ the AM4 atmosphere model, MOM6 ocean code, LM4 land model and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0o (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1o to 0.25o. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM 4 models, but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time‐mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM 4 models.
Lu, Lv, Shaoqing Zhang, Stephen G Yeager, Gokhan Danabasoglu, P Chang, Lixin Wu, Xiaopei Lin, Anthony Rosati, and Feiyu Lu, September 2020: Impact of Coherent Ocean Stratification on AMOC Reconstruction by Coupled Data Assimilation with a Biased Model. Journal of Climate, 33(17), DOI:10.1175/JCLI-D-19-0735.1. Abstract
The Atlantic meridional overturning circulation (AMOC) is of great importance in Earth’s climate system, and reconstructing its structure and variability by combining observations with a coupled model is a key step in understanding historical and future states of AMOC. However, models always have systematic errors called bias owing to imperfect numerical representation of the real world. Model bias and the sparse nature of ocean observations, particularly in deep oceans, make it difficult to generate a complete historical picture of AMOC structure and variability. Here, two coupled models that are biased with respect to each other are used to design “twin” experiments to systematically study the influence of model bias on AMOC reconstruction. One model is used to produce the “observations” that sample the “true” solution of the AMOC to be reconstructed, while the other model is used to incorporate the “observations” to reconstruct the “truth” through coupled data assimilation (CDA). The degree to which the “truth” is recovered by a CDA scheme assesses the critical role of coherent (both upper- and deep-ocean incorporate enough observations to mitigate stratification instability) ocean stratification on AMOC reconstruction. Results show that balancing restoration of climatology and assimilation of observations is vital to better reconstruct AMOC structure and variability, given that most ocean observations are only available in the upper 2000 m. The gained results serve as a guideline in ocean-state estimation with a balance of deep restoring and upper data constraint for climate prediction initialization, especially for decadal predictions.
The next‐generation seasonal prediction system is built as part of the Seamless System for Prediction and EArth System Research (SPEAR) at the Geophysical Fluid Dynamics Laboratory (GFDL) of NOAA. SPEAR is an effort to develop a seamless system for prediction and research across timescales. The ensemble‐based ocean data assimilation (ODA) system is updated for Modular Ocean Model version 6 (MOM6), the ocean component of SPEAR. Ocean initial conditions for seasonal predictions, as well as an ocean state estimation, are produced by the MOM6 ODA system in coupled SPEAR models. Initial conditions of the atmosphere, land, and sea ice components for seasonal predictions are constructed through additional nudging experiments in the same coupled SPEAR models. A bias correction scheme called Ocean Tendency Adjustment (OTA) is applied to coupled model seasonal predictions to reduce model drift. OTA applies the climatological temperature and salinity increments obtained from ocean data assimilation as 3‐dimensional tendency terms to the MOM6 ocean component of the coupled SPEAR models. Based on preliminary retrospective seasonal forecasts, we demonstrate that OTA reduces model drift—especially SST forecast drift—in coupled model predictions and improves seasonal prediction skill for applications such as El Niño‐Southern Oscillation (ENSO).
Sun, Jingzhe, Zhengyu Liu, Feiyu Lu, Weimin Zhang, and Shaoqing Zhang, June 2020: Strongly Coupled Data Assimilation Using Leading Averaged Coupled Covariance (LACC). Part III: Assimilation of Real World Reanalysis. Monthly Weather Review, 148(6), DOI:10.1175/MWR-D-19-0304.1. Abstract
Recent studies proposed LACC (leading averaged coupled covariance) as an effective strongly coupled data assimilation (SCDA) method to improve the coupled state estimation over weakly coupled data assimilation (WCDA) in a coupled general circulation model (CGCM). This SCDA method, however, has been previously evaluated only in the perfect model scenario. Here, as a further step towards evaluating LACC for real world data assimilation, LACC is evaluated for the assimilation of reanalysis data in a CGCM. Several criterions are used to evaluate LACC against the benchmark WCDA. It is shown that despite significant model bias, LACC can improve the coupled state estimation over WCDA. Compared to WCDA, LACC increases the globally averaged anomaly correlation coefficients (ACCs) of sea surface temperature (SST) by 0.036 and atmosphere temperature at the bottom level (Ts) by 0.058. However, there also exist regions where WCDA outperforms LACC. Although the reduction in the anomaly root-mean-square error (RMSE) is not as consistently clear as the increase in ACC, LACC can largely correct the biased model climatology.
Zhang, Shaoqing, Zhengyu Liu, X-F Zhang, Xinrong Wu, G Han, Y Zhao, X Yu, C Liu, Y Liu, S Wu, and Feiyu Lu, et al., June 2020: Coupled data assimilation and parameter estimation in coupled ocean–atmosphere models: a review. Climate Dynamics, 54(11-12), DOI:10.1007/s00382-020-05275-6. Abstract
Recent studies have started to explore coupled data assimilation (CDA) in coupled ocean–atmosphere models because of the great potential of CDA to improve climate analysis and seamless weather–climate prediction on weekly-to-decadal time scales in advanced high-resolution coupled models. In this review article, we briefly introduce the concept of CDA before outlining its potential for producing balanced and coherent weather–climate reanalysis and minimizing initial coupling shocks. We then describe approaches to the implementation of CDA and review progress in the development of various CDA methods, notably weakly and strongly coupled data assimilation. We introduce the method of coupled model parameter estimation (PE) within the CDA framework and summarize recent progress. After summarizing the current status of the research and applications of CDA-PE, we discuss the challenges and opportunities in high-resolution CDA-PE and nonlinear CDA-PE methods. Finally, potential solutions are laid out.
Liu, Zhengyu, C He, and Feiyu Lu, August 2018: Local and Remote Responses of Atmospheric and Oceanic Heat Transports to Climate Forcing: Compensation versus Collaboration. Journal of Climate, 31(16), DOI:10.1175/JCLI-D-17-0675.1. Abstract
We present a theoretical study on local and remote responses of atmosphere and ocean meridional heat transports (AHT and OHT, respectively) to climate forcing in a coupled energy balance model. We show that, in general, a surface heat flux forces opposite AHT and OHT responses in the so-called compensation response, while a net heat flux into the coupled system forces AHT and OHT responses of the same direction in the so-called collaboration response. Furthermore, unless the oceanic thermohaline circulation is significantly changed, a remote climate response far away from the forcing region tends to be dominated by the collaboration response, because of the effective propagation of a coupled ocean–atmosphere energy transport mode of collaboration structure. The relevance of our theory to previous CGCM experiments is also discussed. Our theoretical result provides a guideline for understanding of the response of heat transports and the associated climate changes.
Lu, Feiyu, and Zhengyu Liu, November 2018: Assessing Extratropical Influence on Observed El Nino-Southern Oscillation Events Using Regional Coupled Data Assimilation. Journal of Climate, 31(21), DOI:10.1175/JCLI-D-17-0849.1. Abstract
The extratropical influence on the observed events of El Niño–Southern Oscillation (ENSO) variability from 1948 to 2015 is assessed by constraining the extratropical atmospheric variability in a coupled general circulation model (CGCM) using the regional coupled data assimilation (RCDA) method. The ensemble-mean ENSO response to extratropical atmospheric forcing, which is systematically and quantitatively studied through a series of RCDA experiments, indicates robust extratropical influence on some observed ENSO events. Furthermore, an event-by-event quantitative analysis shows significant differences of the extratropical influence among the observed ENSO events, both in its own strength and in its relation to tropical precursors such as the equatorial Pacific heat content anomaly. This study provides the first dynamic quantitative assessment of the extratropical influence on observed ENSO variability on an event-by-event basis.
Lu, Feiyu, Zhengyu Liu, Y Liu, Shaoqing Zhang, and R Jacob, May 2017: Understanding the control of extratropical atmospheric variability on ENSO using a coupled data assimilation approach. Climate Dynamics, 48(9), DOI:10.1007/s00382-016-3256-7. Abstract
The control of extratropical atmospheric variability on ENSO variability is studied in a coupled general circulation model (CGCM) utilizing an ensemble-based coupled data assimilation (CDA) method in the perfect-model framework. Assimilation is limited to the desired model components (e.g. atmosphere) and spatial areas (e.g. the extratropics) to study the ensemble-mean model response (e.g. tropical response to “observed” extratropical atmospheric variability). The CDA provides continuously “corrected” extratropical atmospheric forcing and boundary conditions for the tropics and the use of ensemble optimizes the observational forcing signal over internal variability in the model component or region without assimilation. The experiments demonstrate significant control of extratropical atmospheric forcing on ENSO variability in the CGCM. When atmospheric “observations” are assimilated only poleward of 20° in both hemispheres, most ENSO events in the “observation” are reproduced and the error of the Nino3.4 index is reduced by over 40 % compared to the ensemble control experiment that does not assimilate any observations. Further experiments with the assimilation in each hemisphere show that the forced ENSO variability is contributed roughly equally and independently by the Southern and Northern Hemisphere extratropical atmosphere. Further analyses of the ENSO events in the southern hemisphere forcing experiment reveal robust precursors in both the extratropical atmosphere over southeastern Pacific and equatorial Pacific thermocline, consistent with previous studies of the South Pacific Meridional Mode and the discharge-recharge paradigm, respectively. However, composite analyses based on each precursor show that neither precursor alone is sufficient to trigger ENSO onset by itself and therefore neither alone could serve as a reliable predictor. Additional experiments with northern hemisphere forcing, ocean assimilation or different latitudes are also performed.
This paper tests the idea of substituting the atmospheric observations with atmospheric reanalysis when setting up a coupled data assimilation system. The paper focuses on the quantification of the effects on the oceanic analysis resulted from this substitution and designs four different assimilation schemes for such a substitution. A coupled Lorenz96 system is constructed and an ensemble Kalman filter is adopted. The atmospheric reanalysis and oceanic observations are assimilated into the system and the analysis quality is compared to a benchmark experiment where both atmospheric and oceanic observations are assimilated. Four schemes are designed for assimilating the reanalysis and they differ in the generation of the perturbed observation ensemble and the representation of the error covariance matrix. The results show that when the reanalysis is assimilated directly as independent observations, the root-mean-square error increase of oceanic analysis relative to the benchmark is less than 16% in the perfect model framework; in the biased model case, the increase is less than 22%. This result is robust with sufficient ensemble size and reasonable atmospheric observation quality (e.g., frequency, noisiness, and density). If the observation is overly noisy, infrequent, sparse, or the ensemble size is insufficiently small, the analysis deterioration caused by the substitution is less severe since the analysis quality of the benchmark also deteriorates significantly due to worse observations and undersampling. The results from different assimilation schemes highlight the importance of two factors: accurate representation of the error covariance of the reanalysis and the temporal coherence along each ensemble member, which are crucial for the analysis quality of the substitution experiment.
Lu, Feiyu, Zhengyu Liu, Shaoqing Zhang, and Y Liu, September 2015: Strongly Coupled Data Assimilation Using Leading Averaged Coupled Covariance (LACC). Part I: Simple Model Study. Monthly Weather Review, 143(9), DOI:10.1175/MWR-D-14-00322.1. Abstract
This paper studies a new Leading Averaged Coupled Covariance (LACC) method for the strongly coupled data assimilation (SCDA). The SCDA not only uses the coupled model to generate the forecast and assimilates observations into multiple model components like the weakly coupled version (WCDA), but also applies cross update using the coupled covariance between variables from different model components. The cross update could potentially improve the balance and quality of the analysis, but its implementation has remained a great challenge in practice due to different timescales between model components.
In a typical extra-tropical coupled system, the ocean-atmosphere correlation shows a strong asymmetry with the maximum correlation occurring when the atmosphere leads the ocean by about the decorrelation time of the atmosphere. The LACC method utilizes such asymmetric structure by using the leading forecasts and observations of the fast atmospheric variable for cross update, therefore increasing the coupled correlation and enhancing the signal-to-noise ratio in calculating the coupled covariance. Here it is applied to a simple coupled model with the Ensemble Kalman Filter (EnKF). With the LACC method, the SCDA reduces the analysis error of the oceanic variable by over 20% compared to the WCDA and 10% compared to the SCDA using simultaneous coupled covariance. The advantage of the LACC method is more notable when the system contains larger errors, such as in the cases with smaller ensemble size, bigger timescale difference or model biases.
Lu, Feiyu, Zhengyu Liu, Shaoqing Zhang, Y Liu, and R Jacob, November 2015: Strongly Coupled Data Assimilation Using Leading Averaged Coupled Covariance (LACC). Part II: CGCM experiments. Monthly Weather Review, 143(11), DOI:10.1175/MWR-D-15-0088.1. Abstract
This paper uses a fully coupled general circulation model (CGCM) to study the Leading Averaged Coupled Covariance (LACC) method in a strongly coupled data assimilation (SCDA) system. Our previous study in a simple coupled climate model (Lu et al. 2015) has shown that, by calculating the coupled covariance using the leading averaged atmospheric states, the LACC method enhances the signal-to-noise ratio and improves the analysis quality of the slow model component compared to both the traditional weakly coupled data assimilation without cross-component adjustments (WCDA) and the regular SCDA using the simultaneous coupled covariance (SimCC).
Here in Part II, we test the LACC method with a CGCM in a perfect-model framework. By adding the observational adjustments from the low-level atmosphere temperature to the sea surface temperature (SST), the SCDA using LACC significantly reduces the SST error compared to WCDA over the globe; it also improves from the SCDA using SimCC, which performs better than the WCDA only in the deep tropics. The improvement in SST analysis is a result of the enhanced signal-to-noise ratio in the LACC method, especially in the extra-tropical regions. The improved SST analysis also benefits the subsurface ocean temperature and low-level atmosphere temperature analyses through dynamic and statistical processes.