Engineering topics
Publications and source records attributed to Liu, Zhengyu.
We have successfully developed the first working version of online ensemble coupled data assimilation (ECDA) on E3SM for both version 1 and version 2. Our schemes have also been tested successfully for assimilating SST data in both synthetic observations and real world observations, laying a foundation for future development.
The long-term response of the Atlantic meridional overturning circulation (AMOC) to anthropogenic forcing has been difficult to detect from the short direct measurements available due to strong interdecadal variability. Here, we present observational and modeling evidence for a likely accelerated weakening of the AMOC since the 1980s under the combined forcing of anthropogenic greenhouse gases and aerosols. This likely accelerated AMOC weakening signal can be detected in the AMOC fingerprint of salinity pileup remotely in the South Atlantic, but not in the classic warming hole fingerprint locally in the North Atlantic, because the latter is contaminated by the “noise” of interdecadal variability. Our optimal salinity fingerprint retains much of the signal of the long-term AMOC trend response to anthropogenic forcing, while dynamically filtering out shorter climate variability. Given the ongoing anthropogenic forcing, our study indicates a potential further acceleration of AMOC weakening with associated climate impacts in the coming decades.
The role of ocean forcing on Atlantic multidecadal variability (AMV) is assessed from the (downward) heat flux–SST relation in the framework of a new stochastic climate theory forced by red noise ocean forcing. Previous studies suggested that atmospheric forcing drives SST variability from monthly to interannual time scales, with a positive heat flux–SST correlation, while heat flux induced by ocean processes can drive SST variability at decadal and longer time scales, with a negative heat flux–SST correlation. Here, first, we develop a theory to show how the sign of heat flux–SST correlation is affected by atmospheric and oceanic forcing with time scale. In particular, a red noise ocean forcing is necessary for the sign reversal of heat flux–SST correlation. Furthermore, this sign reversal can be detected equivalently in three approaches: the low-pass correlation at lag zero, the unfiltered correlation at long (heat flux) lead, and the real part of the heat flux–SST coherence. Second, we develop a new scheme in combination with the theory to assess the magnitude and time scale of the red noise ocean forcing for AMV in the GFDL SPEAR model (Seamless System for Prediction and Earth System Research) and observations. In both the model and observations, the ocean forcing on AMV is in general comparable with the atmospheric forcing, with a 90% probability greater than the atmospheric forcing in observations. In contrast to the white noise atmospheric forcing, the ocean forcing has a persistence time comparable or longer than a year, much longer than the SST persistence of ~3 months. This slow ocean forcing is associated implicitly with slow subsurface ocean dynamics.
Numerical stability is of critical importance in general circulation models because it affects the design of algorithms, time to solution, and computational costs associated with the simulations, which are very expensive in practice. In this paper we extend the stability analysis for ocean-atmosphere coupling proposed in [Zhang et al., J. Sci. Comput. 84, 44(2020)] to a more realistic model that includes horizontal advection. We analyze various time-stepping strategies. We find that advection has a stabilizing effect in scenarios common to climate models when bulk interface condition and explicit flux coupling are used. We also show that our method can be used to study the stability impact of advection for other interface conditions such as Dirichlet-Neumann conditions.
At the Last Glacial Maximum (LGM), records suggest drier conditions in the northwest United States and wetter conditions in the southwest United States relative to present-day as well as widespread changes in the isotopic composition of water. However, the mechanisms responsible for these changes remain ambiguous. Here, we explore differences in western United States hydroclimate between the LGM and preindustrial with a water isotope tracer enabled Earth System Model. We then use proxy forward models to compare simulated and recorded δ 18 O in speleothems. We find that the pattern of hydroclimate response in the western United States at the LGM relates to a combination of 1) increased frequency and southward shifted wintertime extratropical cyclones in the North Pacific, 2) greater rainout of moisture as it moves over the continent, and 3) reduced evaporation in the cooler LGM climate. The simulated lower δ 18 O of precipitation at the LGM relates predominantly to an increase in cool season moisture removal efficiency, with a secondary contribution from relatively more cool season precipitation. Both surface temperatures and North American ice sheets contribute to these hydroclimate changes at the LGM. Comparisons between δ 18 O from proxy forward models and speleothem records in the western United States show general agreement at the LGM, with increasing depletion moving towards the continental interior. Finally, this study highlights the similarities and differences between hydrologic and δ 18 O changes at the LGM and emphasizes the utility of model-proxy comparison for interpretation.
Speleothem records from the South Asian summer monsoon (SASM) region display variability in the ratio of O-18 and O-16 (δ O-18) in calcium carbonate at orbital frequencies. The dominant mode of variability in many of these records reflects cycles of precession. There are several potential explanations for why SASM speleothem records show a strong precession signal, including changes in temperature, precipitation, and circulation. Here we use an Earth system model with water isotope tracers and water-tagging capability to deconstruct the precession signal found in SASM speleothem records. Our results show that cycles of precession-eccentricity produce changes in SASM intensity that correlate with local temperature, precipitation, and δ O-18. However, neither the amount effect nor temperature differences are responsible for the majority of the SASM δ O-18 variability. Instead, changes in the relative moisture contributions from different source regions drive much of the SASM δ O-18 signal, with more nearby moisture sources during Northern Hemisphere summer at aphelion and more distant moisture sources during Northern Hemisphere summer at perihelion. Further, we find that evaporation amplifies the δ O-18 signal of soil water relative to that of precipitation, providing a better match with the SASM speleothem records. This work helps explain a significant portion of the long-term variability found in SASM speleothem records.
Studying the El Nino Southern Oscillation (ENSO) in the past can help us better understand its dynamics and improve its future projections. However, both paleoclimate reconstructions and model simulations of ENSO strength at the Last Glacial Maximum (LGM; 21 ka B.P.) have led to contradicting results. Here we perform model simulations using the recently developed water isotope-enabled Community Earth System Model (iCESM). For the first time, model-simulated oxygen isotopes are directly compared with those from ENSO reconstructions using the individual foraminifera analysis (IFA). We find that the LGM ENSO is most likely weaker comparing with the preindustrial. The iCESM suggests that total variance of the IFA records may only reflect changes in the annual cycle instead of ENSO variability as previously assumed. Furthermore, the interpretation of subsurface IFA records can be substantially complicated by the habitat depth of thermocline-dwelling foraminifera and their vertical migration with a temporally varying thermocline.