Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing layer”

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.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Satellite Mapping of PM2.5 Episodes in the Wintertime San Joaquin Valley: A "Static" Model Using Column Water Vapor

The use of satellite Aerosol Optical Thickness (AOT) from imaging spectrometers has been successful in quantifying and mapping high PM2.5 (particulate matter mass <2.5µm diameter) episodes for pollution abatement and health studies. However, some regions have high PM2.5 but poor estimation success. The challenges in using Aerosol Optical Thickness (AOT) from imaging spectrometers to characterize PM2.5 worldwide was especially evident in the wintertime San Joaquin Valley (SJV). The SJV's attendant difficulties of high-albedo surfaces and very shallow, variable vertical mixing also occur in other significantly polluted regions around the world. We report on more accurate PM2.5 maps for the whole-winter period in the SJV, Nov 14, 2012?Dec 11, 2013. Intensive measurements by including NASA aircraft were made for several weeks in that winter, the DISCOVER-AQ California mission.We found success with a relatively simple method based on calibration and checking with surface monitors and a characterization of vertical mixing, and incorporating specific understandings of the region's climatology. We estimate PM2.5 to within ~7µgm?3 RMSE and with R values of ~0.9, based on remotely sensed MAIAC (Multi-Angle Implementation of Atmospheric Correction) observations, and that certain further work will improve that accuracy. Mapping is at 1km resolution. This allows a time sequence of mapped aerosols at 1km for cloud-free days. We describe our technique as a "static estimation". Estimation procedures like this one, not dependent on well-mapped source strengths or on transport error, should help full source-driven simulations by deconstructing processes. They also provide a rapid method to create a long-term climatology.Essential features of the technique are (a) daily calibration of the AOT to PM2.5 using available surface monitors, and (b) characterization of mixed-layer dilution using column water vapor (CWV, otherwise "precipitable water"). We noted that on multi-day timescales both water vapor and particles share near-surface sources and both fall to very low values with altitude; indeed, both are largely removed by precipitation. The existence of layers of H2O or aerosol not within the mixed layer adds complexity, but mixed-effects statistical regression captures essential proportionality of PM2.5 and the ratio variable (AOT/CWV). Accuracy is much higher than previous statistical models, and can be extended to the whole Aqua-satellite data record. The maps and time-series we show suggest a repeated pattern for large valleys like the SJV ? progressive stabilization of the mixing height after frontal passages: PM2.5 is somewhat more determined by day-by-day changes in mixing than it is by the progressive accumulation of pollutants (revealed as increasing AOT).

Chatfield, Robert B.↗
DOE OSTI2024

Convective Boundary-Layer Spatial Heterogeneity Experiment Field Campaign Report

Earth system models (ESMs) require accurate heat, mass, and momentum exchanges between components, which requires accurate observations and modeling of the atmospheric boundary layer. Over the land, the daytime convective atmospheric boundary layer (CBL, also called the convective mixing layer) develops and evolves daily, driven by solar surface heating. Doppler lidar (DL) measurements of vertical velocities provide an effective way to document the diurnal cycle of CBL (Chu et al. 2023). With the support of the U.S. Department of Energy (DOE) Atmospheric Research (ASR) program, we studied convective mixing-layer heights (MLH) across multiple DOE Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory sites based on DL measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Analysis for predicting adiabatic wall temperatures with single hole coolant injection into a low speed crossflow

Assuming the local adiabatic wall temperature equals the local total temperature in a low speed coolant mixing layer, integral conservation equations with and without the boundary layer effects are formulated for the mixing layer downstream of a single coolant injection hole oriented at a 30 degree angle to the crossflow. These equations are solved numerically to determine the center line local adiabatic wall temperature and the effective coolant coverage area. Comparison of the numerical results with an existing film cooling experiment indicates that the present analysis permits a simplified but reasonably accurate prediction of the centerline effectiveness and coolant coverage area downstream of a single hole crossflow streamwise injection at 30 degree inclination angle.

Wang, C. R.↗
NASA NTRS1981

Analysis for predicting adiabatic wall temperatures with single hole coolant injection into a low speed crossflow

Assuming the local adiabatic wall temperature equals the local total temperature in a low speed coolant mixing layer, integral conservation equations with and without the boundary layer effects are formulated for the mixing layer downstream of a single coolant injection hole oriented at a 30 degree angle to the crossflow. These equations are solved numerically to determine the center-line local adiabatic wall temperature and the effective coolant coverage area. Comparison of the numerical results with an existing film cooling experiment indicates that the present analysis permits a simplified but reasonably accurate prediction of the centerline effectiveness and coolant coverage area downstream of a single hole crossflow streamwise injection at 30-deg inclination angle.

Wang, C. R.↗
NASA NTRS1992

Magellan observations of extended impact crater related features on the surface of Venus

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber & Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Campbell, D. B.↗
NASA NTRS2012

Diagnostic modeling of dimethylsulfide production in coastal water west of the Antarctic Peninsula

The rate of gross biological dimethylsulfide (DMS) production at two coastal sites west of the Antarctic Peninsula, off Anvers Island, near Palmer Station, was estimated using a diagnostic approach that combined field measurements from 1 January 2006 through 1 March 2006 and a one-dimensional physical model of ocean mixing. The average DMS production rate in the upper water column (0-60 m) was estimated to be 3.1 +/- 0.6 nM/d at station B (closer to shore) and 2.7 +/- 0.6 nM/d1 at station E (further from shore). The estimated DMS replacement time was on the order of 1 d at both stations. DMS production was greater in the mixed layer than it was below the mixed layer. The average DMS production normalized to chlorophyll was 0.5 +/- nM/d)/(mg cubic m) at station B and 0.7 +/- 0.2 (nM/d)/(mg/cubic m3) at station E. When the diagnosed production rates were normalized to the observed concentrations of total dimethylsulfoniopropionate (DMSPt, the biogenic precursor of DMS), we found a remarkable similarity between our estimates at stations B and E (0.06 +/- 0.02 and 0.04 +/- 0.01 (nM DMS / d1)/(nM DMSP), respectively) and the results obtained in a previous study from a contrasting biogeochemical environment in the North Atlantic subtropical gyre (0.047 =/- 0.006 and 0.087 +/- 0.014 (nM DMS d1)/(nM DMSP) in a cyclonic and anticyclonic eddy, respectively).We propose that gross biological DMS production normalized to DMSPt might be relatively independent of the biogeochemical environment, and place our average estimate at 0.06 +/- 0.01 (nM DMS / d)/(nM DMSPt). The significance of this finding is that it can provide a means to use DMSPt measurements to extrapolate gross biological DMS production, which is extremely difficult to measure experimentally under realistic in situ conditions.

Dimethylsulfide↗
NASA NTRS1994

Highly ionized gas in the Galactic disk: The l = 295 deg, b,= 0 deg interarm direction

I present new high-resolution measurements of ultraviolet interstellar absorption lines toward three distant stars in the Galactic disk in the general direction l = 295 deg, b = 0 deg. The high ionization lines of Si IV and C IV toward HD 100276, HD 103779, and HD 104705 show substantial variations on angular scales of less than 1 deg, or about 50-100 pc at the distance of the stars. The Si IV and C IV midplane densities derived for HD 100276 and HD 103779 are N(sub o) (Si IV) approximately = 1.3 x 10(exp -9) cm(exp -3) and n(sub o)(C IV) approximately = 4.5 x 10(exp -9) cvm(exp -3), approximately a factor of 2-3 lower than the Galactic averages. Toward HD 104705, the Si IV midplane density is near the average value, but the C IV midplane density is a factor of 3 higher than average. The integrated C IV and Si IV ratios along the HD 100276 and HD 100379 sight lines are similar to the N(C IV)/N(Si IV) = 3.6 +/- 1.3 ratio for a sample of halo stars. For the HD 104705 sight line, N(C IV)/N(Si IV) = 11 and N(C IV)/N(N V) greater than 10. Turbulent mixing layers provide a possible source for the observed absorption along the three sight lines as does gas associated with cooling flows. If mixing layers are the primary source of the absorption along the disk sight lines, then the entrainment velocity of the hot gas is approximately 100 km/s and the post mixed gas temperature is approximately 2 x 10(exp 5) K; same differences may exist for these mixing layers compared to those inferred for gas toward the Galactic poles. Current models of cooling, hot gas in supernova remnants, and magnetized conduction fronts do not accurately describe the high ion properties toward the disk stars in this study.

Sembach, Kenneth R.↗
NASA NTRS1988

Boundary layer ozone - An airborne survey above the Amazon Basin

Ozone data obtained over the forest canopy of the Amazon Basin during July and August 1985 in the course of NASA's Amazon Boundary Layer Experiment 2A are discussed, and ozone profiles obtained during flights from Belem to Tabatinga, Brazil, are analyzed to determine any cross-basin effects. The analyses of ozone data indicate that the mixed layer of the Amazon Basin, for the conditions of undisturbed meteorology and in the absence of biomass burning, is a significant sink for tropospheric ozone. As the coast is approached, marine influences are noted at about 300 km inland, and a transition from a forest-controlled mixed layer to a marine-controlled mixed layer is noted.

Gregory, Gerald L.↗
NASA NTRS1984

Stratocumulus cloud deepening through entrainment

Under a substantial range of realistic conditions, stratocumulus cloud top entrainment is noted to either deepen an existing cloud layer or produce clouds in an unsaturated mixed layer, though the entrained air is warmer and drier than the mixed-layer air. These results, which apply irrespective of entrainment rate-determining mechanism, imply that the cloud top entrainment instability discussed by Randall (1980) and Deardorff (1980) does not necessarily destroy a layer cloud. Examples are given which include soundings, marine layer data, and simulation results produced by the UCLA general circulation model.

Randall, D. A.↗
DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Remote detection of ocean features in the Lesser Antilles using ERTS-1 data

The author has identified the following significant results. Photographic data received from the ERTS-1 satellite over the Lesser Antilles Islands shows distinct ocean features on the leeward side of each island. Attempts to relate these features to ocean eddy formations with the aid of ground truth data proved unsuccessful. However, analysis of surface and upper air wind data correlate extremely well with the size, shape, and downwind extent of the ocean features. Studies to date indicate strongly that these features result from horizontal differences in sea surface roughness due to the wind shadow effect of the islands. The results suggest that horizontal variations in the reflectance of the sea surface will make remote sensing of the ocean mixed layer more difficult than previously anticipated. The surface reflection seems to be large enough to mask the smaller variations in backscattered energy from the mixed layer. Efforts to limit the effect of surface reflectance by photographic differencing of two MSS bands were unsuccessful. A supplementary study to enhance the energy reflected from the mixed layer through numerical differencing of ERTS-1 digital data is planned.

Hanson, K. J.↗
NASA NTRS1993

Experiments on the enhancement of compressible mixing via streamwise vorticity. II - Vortex strength assessment and seed particle dynamics

The effect of streamwise vorticity on compressible axisymmetric mixing layers is examined using vortex strength assessment and seed particle dynamics analysis. Experimental results indicate that the particles faithfully represent the dynamics of the turbulent swirling flow. A comparison of the previously determined mixing layer growth rates with the present vortex strength data reveals that the increase of turbulent mixing up to 60 percent scales with the degree of swirl. The mixing enhancement appears to be independent of the compressibility level of the mixing layer.

Naughton, J. W.↗
NASA NTRS2016

Retrieval of Daily Maps of PM2.5 Aerosol in the Problematic California Valleys: Bright, Speckled Reflectances, Thin AOT, but High Pollution

The San Joaquin Valley suffers from severe episodes of respirable aerosol (PM2.5) in wintertime.We provide maps of aerosol episodes using daily snapshots of PM2.5 and its changing features despite numerous difficulties inherent to sampling the region. Linear relationships relating aerosol optical thickness (AOT) to PM2.5 give an explained variance of approximately 3.The GEO-CAPE mission has as a goal the provision of relevant measures of respirable aerosol to the community,but has not formulated a science goal beyond the limited goal of retrieval of AOT, bringing the usefulness of GEO-CAPE into doubt.Our special focus was on the DISCOVER-AQ period, Jan-Feb 2013, which had many supporting measurements.Both high pollution and retrieval difficulties tend to occur in many Mediterranean agricultural regions like the San Joauin. One difficulty is the relatively bright surfaces with considerable exposed soil. NASAs MAIAC and MODIS Deep Blue retrieval techniques are shown to have considerable skill even at low aerosol optical thickness (AOT) values, as evaluated by concurrent AERONET sunphotometer measurements.More significantly, these AOT values can correspond to high daytime PM2.5 since aerosol mixed layer depth is thin and variable, 200m 600 m. The thin layers derive from typical subsidence of dry air between more stormy periods. This situation provides an advantage: water vapor column is also almost completely limited to a similar mixed layer depth, and can thus serve as a measure of aerosol dilution.Using the MAIAC Water Vapor Column:In order to make the maps below, we used the MAIAC data but subtracted partial water-vapor columns estimated from MERRA Reanalysis Data availabe from the GSFC GMAO using kriging. We did not use the mixed-layer estimates from MERRA, since such analyses were found problematic during our forecasting exercises for DISCOVER-AQ. Observations from the aircraft soundings suggested that this overlying moisture was mostly due to larger scale flows, not ML venting.However, the specific humidity at the surface and a nearly well-mixed ML was analyzed by kriging from the surface network (MesoWest, University of Utah). These were thought to be truer, uninfluenced by physical process modeling that combines with data observations. (TBD: How different are they?) Procedure: Subtract overlying partial water columns from MAIAC column water and divide this by a surface value of water vapor. (MAIAC column is expressed in cm of water, i.e., water vapor at surface conditions.This method appears to bring out useful details in the distribution of submicron particles in the very problematic Wintertime San Joaquin Valley, and allow analysis of pollution episodes throughout the valley, rather than long-term averages.

san joaquin valley↗
DOE OSTI2026

Experiments and simulations on the Richtmyer-Meshkov instability with a thin intermediate layer

Experiments and simulations of the Richtmyer-Meshkov instability (RMI) in two- and three-layer configurations are presented. The two-layer case utilizes a light-over-heavy configuration and consists of air as the light gas and sulfur hexafluoride (SF 6 ) as the heavy gas. The three-layer case utilizes a light-intermediate-heavy configuration, with helium (He) as the light gas, air as the intermediate gas, and SF 6 as the heavy gas. Statistically significant differences in the mixing layer width of the lower interface are not observed between the two cases. This differs from the experiments of Schalles et al., where a small, though statistically significant, difference in mixing layer growth was observed between the two- and three-layer cases with a nominally two-dimensional, single mode perturbation. Notably, the perturbations on the lower interface in the present work do not grow large enough to significantly interact with the upper interface during the duration of the experiments. In conclusion, this suggests that the differences in mixing layer growth observed by Schalles et al. may be due to interactions of the perturbations on one interface with the other interface rather than being inherent to the three-layer problem.

42 ENGINEERING↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Demons in the North Atlantic: Variability of Deep Ocean Ventilation

Abstract Translation of atmospheric forcing variability into the ocean interior via ocean ventilation is an important aspect of transient climate change. On a seasonal timescale in the subtropics, this translation is mediated by a so‐called “Demon” that prevents access to all except late‐winter mixed‐layer water. Here, we use an eddy‐permitting numerical circulation model to investigate a similar process operating on longer (interannual) timescales in the subpolar North Atlantic. We find that variations in atmospheric forcing are mediated in their translation to the ocean interior, with year‐to‐year changes in the late‐winter mixed layer depth being the critical factor. The signature of persistent strong atmospheric forcing driving deep mixed layers is preferentially ventilated to the interior when the forcing is ceased. Susceptibility to this effect depends on the location and density of subduction—with the rate at which newly ventilated water escapes its region of subduction being the crucial factor.

MacGilchrist, G. A.↗
NASA NTRS2000

Estimating the Ocean Flow Field From Combined Sea Surface Temperature and Sea Surface Height Data

The primary focus of this project was on the estimation of the three-dimensional, absolute and time-evolving general circulation of the global ocean from a combined analysis of remotely sensed fields of sea surface temperature (SST) and sea surface height (SSH). The synthesis of those two fields was performed with other relevant physical data, and appropriate dynamical ocean models with emphasis on constraining ocean general circulation models by a combination of both SST and SSH data. This effort is directly related to an attempt to describe the mechanisms which give rise to observed SST and its variability on seasonal and inter-annual timescales, its relation to ocean-atmosphere interaction, and the dynamical coupling between the ocean mixed layer and the deep interior ocean. This is one of the fundamental climate related questions being pursued currently under the CLIVAR Program. Because of the strong turbulent mixing associated with atmospheric fluxes of momentum, heat and freshwater through the sea surface, the ocean forms a shallow surface boundary layer, the mixed layer which is largely homogeneous in its constituents. The relation between the temperature of the remotely sensed "skin" and the bulk of the mixed layer is largely understood (Reynolds and Smith 1994; Emery et al., 1995). However, because the surface mixed layer is effectively decoupled from the underlying ocean dynamics, an interpretation of satellite SST observations in isolation and in direct use for dynamical studies is very difficult. As a result, the impact of SST data on the understanding of ocean variability.

Stammer, Detlef↗
NASA NTRS1994

The roles of vertical mixing, solar radiation, and wind stress in a model simulation of the sea surface temperature seasonal cycle in the tropical Pacfic Ocean

The climatological seasonal cycle of sea surface temperature (SST) in the tropical Pacific is simulated using a newly developed upper ocean model. The roles of vertical mixing, solar radiation, and wind stress are investigated in a hierarchy of numerical experiments with various combinations of vertical mixing algorithms and surface-forcing products. It is found that the large SST annual cycle in the eastern equatorial Pacific is, to a large extent, controlled by the annually varying mixed layer depth which, in turn, is mainly determined by the competing effects of solar radiation and wind forcing. With the application of our hybrid vertical mixing scheme the model-simulated SST annual cycle is much improved in both amplitude and phase as compared to the case of a constant mixed layer depth. Beside the strong effects on vertical mixing, solar radiation is the primary heating term in the surface layer heat budget, and wind forcing influences SST by driving oceanic advective processes that redistribute heat in the upper ocean. For example, the SST seasonal cycle in the western Pacific basically follows the semiannual variation of solar heating, and the cycle in the central equatorial region is significantly affected by the zonal advective heat flux associated with the seasonally reversing South Equatorial Current. It has been shown in our experiments that the amount of heat flux modification needed to eliminate the annual mean SST errors in the model is, on average, no larger than the annual mean uncertainties among the various surface flux products used in this study. Whereas a bias correction is needed to account for remaining uncertainties in the annual mean heat flux, this study demonstrates that with proper treatment of mixed layer physics and realistic forcing functions the seasonal variability of SST is capable of being simulated successfully in response to external forcing without relying on a relaxation or damping formulation for the dominant surface heat flux contributions.

Chen, Dake↗