SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
ENSO simulations are investigated in 30 year integrations of various Intermediate Coupled Models and compared with observed SST, wind and thermocline depth anomalies over the tropical Pacific.
Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyzed
The effort of this first year of research has been focused on the assimilation of TOPEX/Poseidon altimetric data into a primitive equation model of the Atlantic tropical/subtropical circulation. A reduced-rank, stationary Kalman filter has been constructed to assimilate the altimetric sea surface height anomaly (SHA) into the model. The goal is to assess how the inter-hemispheric transports between the Atlantic subtropics and tropics are affected by the assimilation and how the subsurface thermocline structure , and its variability ,is dynamically constrained by the SHA. The model is a reduced-gravity primitive equation GCM of the upper Atlantic Ocean between 30 S and 30 N. The assimilation scheme is an approximation to the extended Kalman filter in which the error covariances of the state estimates are calculated only in a reduced- dimension subspace. The subspace is defined by the leading empirical orthogonal functions calculated from an unconstrained model calculation. Both an identical twin experiment using simulated SHA observations and assimilation of the real TOPEX data were performed. Results from the twin experiments demonstrate the ability of the method to constrain the ocean circulation and the subsurface temperature structure. The impact on the subsurface temperature structure of TOPEX assimilation was assessed using data from expandable bathythermographs. This showed a substantial improvement in the estimated temperature variability only within 13 degrees in latitude around the equator. The impact of TOPEX SHA assimilation on zonally integrated meridional transport across different latitudes was also estimated. Again within 13 degrees from the equator both the mean amplitude and interannual variability of the surface and subsurface transports were significantly enhanced, while the transports were insensitive to the assimilation in the subtropics.
The paper, "Latitude-dependent vertical mixing and the tropical thermocline in a global OGCM", was revised and published in Geophysical Research Letters. It treats the new GISS mixing scheme which includes the latitudinal dependence of the interior ocean turbulence field reported by Gregg, Sanford & Winkel. When implemented in the 3x3 degree NCAR CSMl OGCM [NCOMl] the new mixing scheme produces an improved, sharper equatorial thermoclines in both the Atlantic and the Pacific while simultaneously maintaining the realistic meridional overturning and northward heat transports found already with the previous GISS scheme. Also the paper "Diagnostics of the oceanic thermohaline circulation in a coupled climate model" describing earlier work on the grany was published.
Phytoplankton primary productivity in the Subarctic Atlantic Ocean sustains valuable fisheries and exports anthropogenic atmospheric carbon dioxide, but multiple Earth system models show that this productivity declines substantially in 21st century global warming scenarios. Prior studies have investigated the spatial and temporal dynamics of the seasonal cycle of the phytoplankton, which exhibit a prominent but spatially variable spring bloom that is caused by the closely related seasonal cycles in solar radiation and vertical mixing that modulate both phytoplankton growth and loss rates. However, recent work indicates that multi-decade declines in this spring bloom with global warming in multiple Earth system models are driven by reductions in circulation of nutrients from the subtropical and tropical thermocline, rather than local changes in vertical mixing or shortwave radiation in the Subarctic. It is unclear if current models of how phytoplankton productivity varies seasonally or over a century in global warming can be interpolated to explain the interannual to decadal variability. In this poster, I will quantify the interannual to decadal variability of phytoplankton chlorophyll using historical satellite ocean color observations and model simulations with the Community Earth System Model. I will then use the observations and model simulations to elucidate the mechanisms of interannual to decadal variability and how that variability is restructured in 21st century global warming scenarios.
The breeding method has been implemented in the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Coupled General Circulation Model (CGCM) with the goal of improving operational seasonal to interannual climate predictions through ensemble forecasting and data assimilation. The coupled instability as cap'tured by the breeding method is the first attempt to isolate the evolving ENSO instability and its corresponding global atmospheric response in a fully coupled ocean-atmosphere GCM. Our results show that the growth rate of the coupled bred vectors (BV) peaks at about 3 months before a background ENSO event. The dominant growing BV modes are reminiscent of the background ENSO anomalies and show a strong tropical response with wind/SST/thermocline interrelated in a manner similar to the background ENSO mode. They exhibit larger amplitudes in the eastern tropical Pacific, reflecting the natural dynamical sensitivity associated with the presence of the shallow thermocline. Moreover, the extratropical perturbations associated with these coupled BV modes reveal the variations related to the atmospheric teleconnection patterns associated with background ENSO variability, e.g. over the North Pacific and North America. A similar experiment was carried out with the NCEP/CFS03 CGCM. Comparisons between bred vectors from the NSIPP CGCM and NCEP/CFS03 CGCM demonstrate the robustness of the results. Our results strongly suggest that the breeding method can serve as a natural filter to identify the slowly varying, coupled instabilities in a coupled GCM, which can be used to construct ensemble perturbations for ensemble forecasts and to estimate the coupled background error covariance for coupled data assimilation.
Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.
Much work has gone into revising and updating algorithms for converting satellite-measured radiances to useful ocean variables like sea surface salinity (e.g. SMOS - Boutin et al., 2017, SMAP - Fore et al., 2016 and Aquarius - Meissner et al., 2018). As part of our Ocean Salinity Science Team work, we utilize an intermediate-complexity air/sea coupled model as a transfer function to test if more mature satellite SSS model algorithms actually improve ENSO forecast skill. For initialization of the coupled forecast, we demonstrate that the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together provide the background state to amplify equatorial Kelvin waves and improve ENSO hindcasts (Hackert et al., 2019). Here we extend this work to compare the impact of various pairs of original and improved satellite SSS algorithms. For instance we compare SMAP V4.1 with the latest, SMAP V4.2, to see what impact algorithm improvements may have on ENSO forecasts. SSS observations are tested on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts by initializing twelve-month forecasts for each month of available data. All experiments assimilate satellite sea level (SL), sea surface temperature (SST), and in situ subsurface temperature and salinity observations (Tz, Sz). Additionally various satellite, blended, and in-situ SSS products are assimilated. We find that including satellite SSS significantly improves Niño3.4 sea surface temperature anomaly validation, more mature SSS model algorithms are generally improving ENSO forecasts over time, and more satellite SSS data coverage helps to extend useful forecasts.
The scheme for ocean vertical mixing presented in Canuto, Howard, Cheng & Dubovikov and Canuto, Howard, Hogan, Cheng, Dubovikov & Montenegro [GISS mixing scheme] was extended to include the latitudinal dependence reported by Gregg, Sanford & Winkel of the input to the interior ocean turbulence field of energy from internal waves. The resulting latitude dependence was implemented in our module of code for the GISS vertical mixing scheme and tested in the global NCAR-CSM ocean model with 3 degree X 3 degree resolution and 25 levels in the vertical. With the latitude dependence, the background diffusivity decreases from approx. 0.1 sq cm/sec at midlatitudes to approx. less than 0.01 sq cm/sec at the equator. A significant improvement was seen in the tropical ocean, with tropical thermoclines being sharpened in agreement with the requirements of observational studies and the needs of ENS0 modeling. At the same time global measures nf performance such as the meridional overturning and northward heat transport were not adversely affected. Results were presented at the 2004 Layered Ocean Modeling meeting and a paper has been prepared and submitted to Geophysical Research Letters.
Seasonal and interannual variations of the net surface heating F(sub NET) and sea surface temperature tendency (T(sub s)/dt) in the tropical eastern Indian and western Pacific Oceans are studied. The surface heat fluxes are derived from the Special Sensor Microwave/Imager and Japanese Geostationary Meteorological Satellite radiance measurements for the period October 1997-September 2000. It is found that the magnitude of solar heating is lager than that of evaporative cooling, but the spatial variation of the latter is significantly large than the former. As a result, the spatial variations of seasonal and interannual variability of F(sub NET), follow closely that of evaporative cooling. Seasonal variations of F(sub NET) and T(sub s)/dt are significantly correlated, except for the equatorial western Pacific. The high correlation is primarily attributable to high correlation between seasonal cycles of solar heating and T(sub s)/dt. The change of F(sub NET) between 1997-98 El Nino and 1998-99 La Nina is significantly larger in the tropical eastern Indian Ocean than tropical western Pacific. For the former region, the reduced evaporative cooling arising from weakened winds during the El Nino is generally associated with enhanced solar heating due to decreased cloudiness, and thus increases the interannual variability of F(sub NET). For the latter region, the reduced evaporative cooling due to weakened winds is generally associated with but exceeds the reduced solar heating arising from increased cloudiness, and vise versa. Thus the interannual variability of F(sub NET) is reduced due to this offsetting effect. Interannual variations of F(sub NET) and T(sub s)/dt have very low correlation. This is most likely related to interannual variability of ocean dynamics, which includes the variations of solar radiation penetrating through oceanic mixed layer, upwelling of cold thermocline water, Indonesian throughflow for transporting heat from the Pacific to Indian Ocean, and interhemispheric transport in the Indian Ocean.
We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO3.4 region. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.
We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO 3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10 degrees South latitude-5 degrees North latitude. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO 3.4 region. In addition, salting near the ITCZ (Intertropical Convergence Zone) leads to a deepening of the mixed layer and thermocline near 8 degrees North latitude. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.
Observational analyses and ocean general circulation model experiments were performed to understand the influence of oceanic processes on warm sea surface temperature anomalies (SSTAs) prior to convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden‐Julian Oscillation (MJO), in the tropical Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool (WP) region during the November–April season of the 2001–2012 period. Seventeen of them (44%) initiated over the Seychelles‐Chagos Thermocline Ridge (SCTR) before propagating eastward to the WP and Maritime Continent. Including only global‐scale MJOs, 71% (24%) of the remaining MJO events initiated from the SCTR (WP). Four (seven) SCTR (WP) ISO events were preceded by SSTAs that were strongly influenced by wind stress‐driven oceanic processes; all four (2/7) SCTR (WP) events were MJOs. Composite analyses of the four oceanic process influenced SCTR MJO events showed that wind stress‐driven oceanic processes contributed about 30% (between 0.1 and 0.2 °C) of the preconvection warming in the SCTR, which coincides with the Intertropical Convergence Zone with mean SST > 28 °C. Reduced upwelling and entrainment played a more important role than horizontal advection in elevating the +SSTAs. Case studies revealed that the two ocean process‐influenced primary MJO SCTR events were also associated with oceanic equatorial Rossby waves. Furthermore, results of a linear atmospheric mixed layer model indicated that boundary layer processes related to SSTAs contributed approximately half of the total convergence during MJO initiation, relative to the effect of the free troposphere, for the Rossby wave‐influenced events.