Chen, X., Zhang, Y., Zhang, M., Feng, Y., Wu, Z., Qiao, F., et al. (2013). Intercomparison between observed and simulated variability in global ocean heat content using empirical mode decomposition, part I: modulated annual cycle.
Clim Dyn, 41(11-12), 2797–2815.
Deng, J., Wu, Z., Zhang, M., Huang, N. E., Wang, S., & Qiao, F. (2018). Using Holo-Hilbert spectral analysis to quantify the modulation of Dansgaard-Oeschger events by obliquity.
Quaternary Science Reviews, 192, 282–299.
Abstract: Astronomical forcing (obliquity and precession) has been thought to modulate Dansgaard-Oeschger (DO) events, yet the detailed quantification of such modulations has not been examined. In this study, we apply the novel Holo-Hilbert Spectral Analysis (HHSA) to five polar ice core records, quantifying astronomical forcing's time-varying amplitude modulation of DO events and identifying the preferred obliquity phases for large amplitude modulations. The unique advantages of HHSA over the widely used windowed Fourier spectral analysis for quantifying astronomical forcing's nonlinear modulations of DO events is first demonstrated with a synthetic data that closely resembles DO events recorded in Greenland ice cores (NGRIP, GRIP, and GISP2 cores on GICC05 modelext timescale). The analysis of paleoclimatic proxies show that statistically significantly more frequent DO events, with larger amplitude modulation in the Greenland region, tend to occur in the decreasing phase of obliquity, especially from its mean value to its minimum value. In the eastern Antarctic, although statistically significantly more DO events tend to occur in the decreasing obliquity phase in general, the preferred phase of obliquity for large amplitude modulation on DO events is a segment of the increasing phase near the maximum obliquity, implying that the physical mechanisms of DO events may be different for the two polar regions. Additionally, by using cross-spectrum and magnitude-squared analyses, Greenland DO mode at a timescale of about 1400 years leads the Antarctic DO mode at the same timescale by about 1000 years. (C) 2018 Elsevier Ltd. All rights reserved.
Deng, J., Wu, Z., Zhang, M., Huang, N. E., Wang, S., & Qiao, F. (2019). Data concerning statistical relation between obliquity and Dansgaard-Oeschger events.
Data Brief, 23.
Abstract: Data presented are related to the research article entitled “Using Holo-Hilbert spectral analysis to quantify the modulation of Dansgaard-Oeschger events by obliquity” (J. Deng et al., 2018). The datasets in Deng et al. (2018) are analyzed on the foundation of ensemble empirical mode decomposition (EEMD) (Z.H. Wu and N.E. Huang, 2009), and reveal more occurrences of Dansgaard-Oeschger (DO) events in the decreasing phase of obliquity. Here, we report the number of significant high Shannon entropy (SE) (C.E. Shannon and W. Weaver, 1949) of 95% significance level of DO events in the increasing and decreasing phases of obliquity, respectively. First, the proxy time series are filtered by EEMD to obtain DO events. Then, the time-varying SE of DO modes are calculated on the basis of principle of histogram. The 95% significance level is evaluated through surrogate data (T. Schreiber and A. Schmitz, 1996). Finally, a comparison between the numbers of SE values that are larger than 95% significance level in the increasing and decreasing phases of obliquity, respectively, is reported.
Fox-Kemper, B., Adcroft, A., Böning, C. W., Chassignet, E. P., Curchitser, E., Danabasoglu, G., et al. (2019). Challenges and Prospects in Ocean Circulation Models.
Front. Mar. Sci., 6.
Abstract: We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.
Griffies, S. M., Danabasoglu, G., Durack, P. J., Adcroft, A. J., Balaji, V., Böning, C. W., et al. (2016). OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project.
Geosci. Model Dev., 9(9), 3231–3296.
Wu, Z., Feng, J., Qiao, F., & Tan, Z. - M. (2016). Fast multidimensional ensemble empirical mode decomposition for the analysis of big spatio-temporal datasets.
Philos Trans A Math Phys Eng Sci, 374(2065), 20150197.
Abstract: In this big data era, it is more urgent than ever to solve two major issues: (i) fast data transmission methods that can facilitate access to data from non-local sources and (ii) fast and efficient data analysis methods that can reveal the key information from the available data for particular purposes. Although approaches in different fields to address these two questions may differ significantly, the common part must involve data compression techniques and a fast algorithm. This paper introduces the recently developed adaptive and spatio-temporally local analysis method, namely the fast multidimensional ensemble empirical mode decomposition (MEEMD), for the analysis of a large spatio-temporal dataset. The original MEEMD uses ensemble empirical mode decomposition to decompose time series at each spatial grid and then pieces together the temporal-spatial evolution of climate variability and change on naturally separated timescales, which is computationally expensive. By taking advantage of the high efficiency of the expression using principal component analysis/empirical orthogonal function analysis for spatio-temporally coherent data, we design a lossy compression method for climate data to facilitate its non-local transmission. We also explain the basic principles behind the fast MEEMD through decomposing principal components instead of original grid-wise time series to speed up computation of MEEMD. Using a typical climate dataset as an example, we demonstrate that our newly designed methods can (i) compress data with a compression rate of one to two orders; and (ii) speed-up the MEEMD algorithm by one to two orders.
Zhang, M., Wu, Z., & Qiao, F. (2018). Deep Atlantic Ocean Warming Facilitated by the Deep Western Boundary Current and Equatorial Kelvin Waves.
J. Climate, 31(20), 8541–8555.
Abstract: Increased heat storage in deep oceans has been proposed to account for the slowdown of global surface warming since the end of the twentieth century. How the imbalanced heat at the surface has been redistributed to deep oceans remains to be elucidated. Here, the evolution of deep Atlantic Ocean heat storage since 1950 on multidecadal or longer time scales is revealed. The anomalous heat in the deep Labrador Sea was transported southward by the shallower core of the deep western boundary current (DWBC). Upon reaching the equator around 1980, this heat transport route bifurcated into two, with one continuing southward along the DWBC and the other extending eastward along a narrow strip (about 4 degrees width) centered at the equator. In the 1990s and 2000s, meridional diffusion helped to spread warming in the tropics, making the eastward equatorial warming extension have a narrow head and wider tail. The deep Atlantic Ocean warming since 1950 had overlapping variability of approximately 60 years. The results suggest that the current basinwide Atlantic Ocean warming at depths of 1000-2000 m can be traced back to the subsurface warming in the Labrador Sea in the 1950s. An inference from these results is that the increased heat storage in the twenty-first century in the deep Atlantic Ocean is unlikely to partly account for the atmospheric radiative imbalance during the last two decades and to serve as an explanation for the current warming hiatus.
Zhang, M., Zhang, Y., Shu, Q., Zhao, C., Wang, G., Wu, Z., et al. (2018). Spatiotemporal evolution of the chlorophyll a trend in the North Atlantic Ocean.
Sci Total Environ, 612, 1141–1148.
Abstract: Analyses of the chlorophyll a concentration (chla) from satellite ocean color products have suggested the decadal-scale variability of chla linked to the climate change. The decadal-scale variability in chla is both spatially and temporally non-uniform. We need to understand the spatiotemporal evolution of chla in decadal or multi-decadal timescales to better evaluate its linkage to climate variability. Here, the spatiotemporal evolution of the chla trend in the North Atlantic Ocean for the period 1997-2016 is analyzed using the multidimensional ensemble empirical mode decomposition method. We find that this variable trend signal of chla shows a dipole pattern between the subpolar gyre and along the Gulf Stream path, and propagation along the opposite direction of the North Atlantic Current. This propagation signal has an overlapping variability of approximately twenty years. Our findings suggest that the spatiotemporal evolution of chla during the two most recent decades is part of the multidecadal variations and possibly regulated by the changes of Atlantic Meridional Overturning Circulation, whereas the mechanisms of such evolution patterns still need to be explored.