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.
In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific.
Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new generation GCMs predict too much liquid and too little ice in mixed-phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study (SOCRATES), this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOA) are the most important type of ice nucleating particles (INPs) over the SO with concentrations one order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice-ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than –20°C. SIP enhances the in-cloud ice number concentrations (Ni) by 1-3 orders of magnitude and predicts more mixed-phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. Finally, this study highlights the importance of accurately representing cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.
Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.
Abstract. The General Ocean Turbulence Model (GOTM) is a one-dimensional water column model, including a set of state-of-the-art turbulence closure models, and has widely been used in various applications in the ocean modeling community. Here, we extend GOTM to include a set of newly developed ocean surface vertical mixing parameterizations of Langmuir turbulence via coupling with the Community Vertical Mixing Project (CVMix). A Stokes drift module is also implemented in GOTM to provide the necessary ocean surface waves information to the Langmuir turbulence parameterizations, as well as to facilitate future development and evaluation of new Langmuir turbulence parameterizations. In addition, a streamlined workflow with Python and Jupyter notebooks is also described, enabled by the newly developed and more flexible configuration capability of GOTM. The newly implemented Langmuir turbulence parameterizations are evaluated against theoretical scalings and available observations in four test cases, including an idealized wind-driven entrainment case and three realistic cases at Ocean Station Papa, the northern North Sea, and the central Baltic Sea, and compared with the existing general length scale scheme in GOTM. The results are consistent with previous studies. This development extends the capability of GOTM towards including the effects of ocean surface waves and provides useful toolsets for the ocean modeling community to further study the effects of Langmuir turbulence in a broader scope.
Ocean angular momentum (OAM) calculations using forward model runs without any data constraints have, recently revealed the effects of OAM variability on the Earth's rotation. Here we use an ocean model and its adjoint to estimate OAM values by constraining the model to available oceanic data. The optimization procedure yields substantial changes in OAM, related to adjustments in both motion and mass fields, as well as in the wind stress torques acting on the ocean. Constrained and unconstrained OAM values are discussed in the context of closing the planet's angular momentum budget. The estimation procedure, yields noticeable improvements in the agreement with the observed Earth rotation parameters, particularly at the seasonal timescale. The comparison with Earth rotation measurements provides an independent consistency check on the estimated ocean state and underlines the importance of ocean state estimation for quantitative. studies of the variable large-scale oceanic mass and circulation fields, including studies of OAM.
A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.
Anomalies in ocean surface temperature have been identified as possible causes of variations in the climate of particular seasons or as a source of interannual climatic variability, and attempts have been made to forecast seasonal climate by using ocean temperatures as predictor variables. However, the seasonal atmospheric response to ocean temperature anomalies has not yet been systematically investigated with nonlinear models. The present investigation is concerned with ten-year integrations involving a model of intermediate complexity, the Held-Suarez climate model. The calculations have been performed to investigate the changes in seasonal climate which result from a fixed anomaly imposed on a seasonally varying, global ocean temperature field. Part I of the paper provides a report on the results of these decadal integrations. Attention is given to model properties, the experimental design, and the anomaly experiments.
Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.
Sensitivities of the oceanic biological pump within the GISS (Goddard Institute for Space Studies ) climate modeling system are explored here. Results are presented from twin control simulations of the air-sea CO2 gas exchange using two different ocean models coupled to the same atmosphere. The two ocean models (Russell ocean model and Hybrid Coordinate Ocean Model, HYCOM) use different vertical coordinate systems, and therefore different representations of column physics. Both variants of the GISS climate model are coupled to the same ocean biogeochemistry module (the NASA Ocean Biogeochemistry Model, NOBM), which computes prognostic distributions for biotic and abiotic fields that influence the air-sea flux of CO2 and the deep ocean carbon transport and storage. In particular, the model differences due to remineralization rate changes are compared to differences attributed to physical processes modeled differently in the two ocean models such as ventilation, mixing, eddy stirring and vertical advection. GISSEH(GISSER) is found to underestimate mixed layer depth compared to observations by about 55% (10 %) in the Southern Ocean and overestimate it by about 17% (underestimate by 2%) in the northern high latitudes. Everywhere else in the global ocean, the two models underestimate the surface mixing by about 12-34 %, which prevents deep nutrients from reaching the surface and promoting primary production there. Consequently, carbon export is reduced because of reduced production at the surface. Furthermore, carbon export is particularly sensitive to remineralization rate changes in the frontal regions of the subtropical gyres and at the Equator and this sensitivity in the model is much higher than the sensitivity to physical processes such as vertical mixing, vertical advection and mesoscale eddy transport. At depth, GISSER, which has a significant warm bias, remineralizes nutrients and carbon faster thereby producing more nutrients and carbon at depth, which eventually resurfaces with the global thermohaline circulation especially in the Southern Ocean. Because of the reduced primary production and carbon export in GISSEH compared to GISSER, the biological pump efficiency, i.e., the ratio of primary production and carbon export at 75 m, is half in the GISSEH of that in GISSER, The Southern Ocean emerges as a key region where the CO2 flux is as sensitive to biological parameterizations as it is to physical parameterizations. The fidelity of ocean mixing in the Southern Ocean compared to observations is shown to be a good indicator of the magnitude of the biological pump efficiency regardless of physical model choice.
An accurate long-term (84-month) climatology of net surface solar irradiance over the global oceans from Nimbus-7 earth radiation budget (ERB) wide-field-of-view planetary-albedo data is generated via an algorithm based on radiative transfer theory. Net surface solar irradiance is computed as the difference between the top-of-atmosphere incident solar irradiance (known) and the sum of the solar irradiance reflected back to space by the earth-atmosphere system (observed) and the solar irradiance absorbed by atmospheric constituents (modeled). It is shown that the effects of clouds and clear-atmosphere constituents can be decoupled on a monthly time scale, which makes it possible to directly apply the algorithm with monthly averages of ERB planetary-albedo data. Compared theoretically with the algorithm of Gautier et al. (1980), the present algorithm yields higher solar irradiance values in clear and thin cloud conditions and lower values in thick cloud conditions.
The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth's sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates.
Despite observational evidence that an El-Niño is likely to follow tropical volcanic eruptions (Predybaylo et al, 2017), modelling studies have so far reached no consensus on either the sign or physical mechanism of El Niño response to volcanism (Khodri et al, 2017) . A previous modelling study found that following the Pinatubo eruption, the tropical Pacific developed El Niño-like conditions that triggered a decrease in the upwelling of DIC and increased carbon uptake (Eddebbar et al., 2019). However, the carbon cycle response to volcanisms as it depends on ENSO phase remains unclear. Here we perform thirty ensemble control simulations using the GISS-E2-1-G model to sample the pre-industrial control ENSO cycle and thirty more to superimpose a Pinatubo-like volcanic eruption. We find that the model tends to favor the development of a weak El-Niño following the volcanic eruption only in the case of a volcano that erupts at the start of a neutral ENSO phase. In response to the weak El-Niño, some anomalies have the same sign as those that develop from a strong El-Niño, such as weaker winds, shallower mixed layer depths, less direct but more diffuse radiation at the surface of the ocean, less nitrate, and lower primary productivity. However, other anomalies show the opposite sign as those from a strong El-Niño, such as more precipitation, less surface DIC, and less surface ocean pCO 2 . In contrast to the neutral phase, the volcanic response during positive ENSO and negative ENSO phases are swamped by the model’s internal variability.
A coupled general ocean circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. The model is driven by climatological meteorological conditions, cloud cover, and sea surface temperature. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chorophytes, and picoplankton) and three nutrient groups (nitrate, ammonium, and silicate). Phytoplankton groups are initialized as homogeneous fields horizontally and vertically, and allowed to distribute themselves according to the prevailing conditions. Basin-scale model chlorophyll results are in very good agreement with CZCS pigments in virtually every global region. Seasonal variability observed in the CZCS is also well represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are also in good conformance, although occasional departures are apparent. Agreement of nitrate distributions with in situ data is even better, including seasonal dynamics, except for the equatorial Atlantic. The good agreement of the model with satellite and in situ data sources indicates that the model dynamics realistically simulate phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization, and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.
Seasonal variations in the speed of the Earth's rotation manifest themselves as fluctuations in the length of the day (LOD) with an amplitude of about 1000 microseconds. We know from previous work that at least 95% of these variations can be accounted for in terms of angular momentum exchanged between the atmosphere and the solid Earth. Here we examine the respective contributions of the Antarctic Circumpolar Current (ACC) and the global oceans to the Earth's seasonal angular momentum budget, using in situ data from the Drake Passage and results from both the oceanic regional model (Fine Resolution Antarctic Model -- FRAM) of Webb et al. (1991) and the global ocanic model of Maier-Reimer et al. (1993) as analyzed by Brosche et al. (1990). The estimated annual contribution of the ACC (2-4 microsec) is much smaller than the total variation in the oceanic models or the existing LOD-AAM residual (both approximately 15-20 microsec). The estimated semi-annual ACC contribution (3-8 microsec) is offset by counter-current further north in both oceanic models, which exhibit larger semi-annual variations in planetary angular momentum. Further refinements in the Earth's seasonal angular momentum budget, therefore, will require the full (planetary plus relative) contribution of the global oceans in addition to that of the ACC.
The Earth science community seeks to develop climate data records (CDRs) from satellite measurements of ocean color, a continuous data record that now exceeds 20 years. Space agencies will launch additional instruments in the coming decade that will continue this data record, including the NASA PACE spectrometer. Inherent optical properties (IOPs) quantitatively describe the absorbing and scattering constituents of seawater and can be estimated from satellite-observed spectroradiometric data using semi-analytical algorithms (SAAs). SAAs exploit the contrasting optical signatures of constituent matter at spectral bands observed by satellite sensors. SAA performance, therefore, depends on the spectral resolution of the satellite spectroradiometer. A CDR spanning SeaWiFS, MODIS, OLCI, and PACE, for example, would include IOPs derived using varied wavelength suites if all available wavelengths were considered. Here, we explored differences in derived IOPs that stem simply from the use of (eight) different wavelength suites of input radiometric measurements. Using synthesized data and SeaWiFS Level-3 mission-long composites, we demonstrated equivalent SAA performance for all wavelength suites, but that IOP retrievals vary by several percent across wavelength suites and as a function of water type. The differences equate to roughly ≤ 6, 12, and 7% for a(sub dg)(443), a(sub ph)(443), and b(sub bp)(443), respectively, for waters with C(sub a) ≤ 1 mg -cu m. These values shrink for sensors with similar wavelength suites (e.g., SeaWiFS, MODIS, and MERIS) and rise to substantially larger values for higher C(sub a) waters. Our results also indicate that including 400 nm (in the case of OLCI) influences the derived IOPs, using longer wavelengths (>600 nm) influences the derived IOPs when there is a red signal, and, including additional spectral information shows potential for improved IOP estimation, but not without revisiting SAA parameterizations and execution. While modest in scope, we believe this study contributes to the knowledge base for CDR development. The implication of ignoring such an analysis as CDRs continue to be developed is a prolonged inability to distinguish between algorithmic and environmental contributions to trends and anomalies in the IOP time-series.
The dependence of the depth and strength of the ocean's global meridional overturning cells (MOC) on the specification of mesoscale eddy diffusivity (K) is explored in two ocean models. The GISS and MIT ocean models are driven by the same prescribed forcing fields, configured in similar ways, spun up to equilibrium for a range of K 's and the resulting MOCs mapped and documented. Scaling laws implicit in modern theories of the MOC are used to rationalize the results. In all calculations the K used in the computation of eddy-induced circulation and that used in the representation of eddy stirring along neutral surfaces, is set to the same value but is changed across experiments. We are able to connect changes in the strength and depth of the Atlantic MOC, the southern ocean upwelling MOC, and the deep cell emanating from Antarctica, to changes in K.
Explore the source record for details and available documents.
The large scale seasonal variability is estimated jointly with the mean sea surface and the ocean tide signal not recovered by the Cartwright and Ray ocean tide model.