Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean model comparisons”

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 199 records · Page 11

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions

Simulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere-ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25deg42.5degN, 125deg-105degW). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes in exogenous forcings. Although the ECHO-G AOGCM is capable of simulating megadroughts through persistent La Nina-like conditions in the tropical Pacific, other mechanisms can produce similarly extreme NASW moisture anomalies in the model. In particular, the lack of low-frequency coherence between NASW soil moisture and simulated modes of climate variability like the El Nino-Southern Oscillation, Pacific decadal oscillation, and Atlantic multidecadal oscillation during identified drought periods suggests that stochastic atmospheric variability can contribute significantly to the occurrence of simulated megadroughts in the NASW. These findings indicate that either an expanded paradigm is needed to understand multidecadal hydroclimate variability in the NASW or AOGCMs may incorrectly simulate the strength and/or dynamics of the connection between NASW hydroclimate variability and the tropical Pacific.

variability↗

Magnetization modeling in the north and equatorial Atlantic Ocean using MAGSAT data

Magsat 2 x 2-deg scalar anomalous-magnetic-field data (Langel et al., 1982) for the northern and equatorial Atlantic are inverted and combined with physiographic data and laboratory results on the magnetization of oceanic rocks and the oceanic crust to construct models explaining the shorter-wavelength component of the anomalies. An annihilator is applied to the inversion results to eliminate reverse-magnetized sources and facilitate comparisons of areas inverted separately, and a latitude effect on source spacing is tentatively attributed to greater noise contamination at lower latitudes. It is found that remanent magnetization combined with considerable crustal thickening can best explain the high intensity levels observed, although viscous magnetization or contamination of the data by noncrustal sources must also be considered.

Hayling, K. L.↗

Unforced decadal fluctuations in a coupled model of the atmosphere and ocean mixed layer

Global average temperature in a 100-year control run of a model used for greenhouse gas response simulations showed low-frequency natural variability comparable in magnitude to that observed over the last 100 years. The model variability was found to be barotropic in the atmosphere, and located in the tropical strip with largest values near the equator in the Pacific. The model variations were traced to complex, low-frequency interactions between the meridional sea surface temperature gradients in the eastern equatorial Pacific, clouds at both high and low levels, and features of the tropical atmospheric circulation. The variations in these and other model parameters appear to oscillate between two limiting climate states. The physical scenario accounting for the oscillations on decadal time scales is almost certainly not found in the real world on shorter time scales due to limited resolution and the omission of key physics (e.g., equatorial ocean dynamics) in the model. The real message is that models with dynamical limitations can still produce significant long-term variability. Only a thorough physical diagnosis of such simulations and comparisons with decadal-length data sets will allow one to decide if faith in the model results is, or is not, warranted.

Barnett, T. P.↗

Use of global positioning system measurements to determine geocentric coordinates and variations in Earth orientation

Geocentric tracking station coordinates and short-period Earth-orientation variations can be measured with Global Positioning System (GPS) measurements. Unless calibrated, geocentric coordinate errors and changes in Earth orientation can lead to significant deep-space tracking errors. Ground-based GPS estimates of daily and subdaily changes in Earth orientation presently show centimeter-level precision. Comparison between GPS-estimated Earth-rotation variations, which are the differences between Universal Time 1 and Universal Coordinated Time (UT1-UTC), and those calculated from ocean tide models suggests that observed subdaily variations in Earth rotation are dominated by oceanic tidal effects. Preliminary GPS estimates for the geocenter location (from a 3-week experiment) agree with independent satellite laser-ranging estimates to better than 10 cm. Covariance analysis predicts that temporal resolution of GPS estimates for Earth orientation and geocenter improves significantly when data collected from low Earth-orbiting satellites as well as from ground sites are combined. The low Earth GPS tracking data enhance the accuracy and resolution for measuring high-frequency global geodynamical signals over time scales of less than 1 day.

R P Malla↗

Efficient Flowline Simulations of Ice Shelf-Ocean Interactions: Sensitivity Studies with a Fully Coupled Model

Thermodynamic flowline and plume models for the ice shelf-ocean system simplify the ice and ocean dynamics sufficiently to allow extensive exploration of parameters affecting ice-sheet stability while including key physical processes. Comparison between geophysically and laboratory-based treatments of ice-ocean interface thermodynamics shows reasonable agreement between calculated melt rates, except where steep basal slopes and relatively high ocean temperatures are present. Results are especially sensitive to the poorly known drag coefficient, highlighting the need for additional field experiments to constrain its value. These experiments also suggest that if the ice-ocean interface near the grounding line is steeper than some threshold, further steepening of the slope may drive higher entrainment that limits buoyancy, slowing the plume and reducing melting; if confirmed, this will provide a stabilizing feedback on ice sheets under some circumstances.

polar oceanography↗

Preliminary Results From GMAO’s JEDI Based Ocean/Sea-ice Data Assimilation System

Preparations for GMAO’s transition to next generation coupled DA framework includes MOM6 for the ocean and CICE6 for sea ice, alongside assembling the data assimilation systems. The coupled DA system will primarily rely on the Joint Effort for Data Assimilation Integration (JEDI) and involves the use of the Sea-ice Ocean and Coupled Analysis (SOCA) package. The early stages of the work involved developing a coupled, cycling (ocean-only) DA workflow named “Swell”. Some other highlights from this transdisciplinary work: Collaborated with NOAA/EMC to ensure direct comparisons using the same observations and time (July -September 2021) during model development. Conducted in-kind work with the JCSDA to improve various JEDI components. Swell coupled DA framework will be used to assimilate observations from NASA’s PACE and SWOT missions

Dorukhan Ardag↗

A high-resolution OGCM simulation of the Tropical Pacific Ocean during the 1985-1994 TOGA period. Part I: Long equatorial waves

A high resolution oceanic general circulation model of the three topical oceans is used to investigate long equatorial wave activity in the Pacific Ocean during the 1985-1994 TOGA period. Zonal wind stress forcing and simulated dynamic height are interpreted using techniques previously applied to data. Kelvin and first Rossby waves are observed propagating during all the period. A seasonal cycle and interannual anomalies are computed for each long equatorial wave. The east Pacific basin is mainly dominated by seasonal cycle variations while strong interannual anomalies are observed west of the dateline. Long wave interannual anomalies are then compared to wave coefficients simulated by a simple wind-forced model. Our results outline the major role played by wind forcing on interannual time scales in generating long equatorial waves. However, near both eastern and western boundaries, some differences can be attributed to long wave reflections. A comparison to wave coefficients calculated from GEOSAT sea-level data gives some insight of the model behavior.

Long Equatorial Waves↗

Modeling the Optical Properties of Biomass Burning Aerosols: Young Smoke Aerosols From Savanna Fires and Comparisons to Observations from SAFARI 2000

Annually, farmers in southern Africa manage their land resources and prepare their fields for cultivation by burning crop residual debris, with a peak in the burning season occurring during August and September. The emissions from these fires in southern Africa are among the greatest from fires worldwide, and the gases and aerosol particles produced adversely affect air quality large distances from their source regions, and can even be tracked in satellite imagery as they cross the Atlantic and Pacific Ocean basins. During August and September 2000 an international group of researchers participating in the Southern African Regional Science Initiate field experiment (SAFARI 2000) made extensive ground-based, airborne, and satellite measurements of these gases and aerosols in order to quantify their amounts and effects on Earth's atmosphere. In this study we interpreted the measurements of smoke aerosol particles made during SAFARI 2000 in order to better represent these particles in a numerical model simulating their transport and fate. Typically, smoke aerosols emitted from fires are concentrated by mass in particles about 0.3 micrometers in diameter (1,000,000 micrometers = 1 meter, about 3 feet); for comparison, the thickness of a human hair is about 50 micrometers, almost 200 times as great. Because of the size of these particles, at the surface they can be easily inhaled into the lungs, and in high concentrations have deleterious health effects on humans. Additionally, these particles reflect and absorb sunlight, impacting both visibility and the balance of sunlight reaching -Earth's surface, and ultimately play a role in modulating Earth's climate. Because of these important effects, it is important that numerical models used to estimate Earth's climate response to changes in atmospheric composition accurately represent the quantity and evolution of smoke particles. In our model, called the Community Aerosol and Radiation Model for Atmospheres (CARMA) we used estimates of smoke emissions based on field studies and observations made with the NASA Terra and TRMM satellites. The meteorology used to calculate the transport was based on an assimilation of observed meteorological conditions provided by the National Center for Atmospheric Research.

Matichuk, R. I.↗

Unique Observational Constraints on the Seasonal and Longitudinal Variability of the Earth’s Planetary Albedo and Cloud Distribution Inferred From EPIC Measurements

Thorough comparison to observations is key to developing a credible climate model forecasting capability. Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation provide a unique observational perspective that permits a more reliable model/data comparison than is possible with the otherwise available satellite data. The uniqueness is in the DSCOVR satellite’s viewing geometry, which enables continuous viewing of the Earth’s sunlit hemisphere from its Lissajous orbit around the Lagrangian L1 point. The key instrument is the Earth Polychromatic Imaging Camera (EPIC), which views the Earth’s sunlit hemisphere with 1024-by-1024-pixel imagery in 10 narrow spectral bands from 317 to 780 nm, acquiring up to 22 high spatial resolution images per day. The additional feature is that the frequency of EPIC image acquisition is nearly identical to that of the climate GCM data generation scheme where climate data for the entire globe are ‘instantaneously’ calculated at 1-h radiation time-step intervals. Implementation of the SHS (Sunlit Hemisphere Sampling) EPIC-view geometry for the in-line GCM output data sampling establishes a precise self-consistency in the space-time data sampling between EPIC observational and GCM output data generation and sampling. The remaining problem is that the GCM generated data are radiative fluxes, while the EPIC measurements are backscatter-dependent radiances. Radiance to flux conversion is a complex problem with no simple way to convert GCM radiative fluxes into spectral radiances. The more expedient approach is to convert the EPIC spectral radiances into broadband radiances by MODIS/CERES-based regression relationships and then into solar radiative fluxes using the CERES angular distribution models. Averaging over the sunlit hemisphere suppresses the meteorological weather noise, but preserves the intra-seasonal larger scale variability. Longitudinal slicing by the Earth’s rotation permits a self-consistent model/data comparison of the longitudinal model/data differences in the variability of the reflected solar radiation. Ancillary EPIC Composite data provide additional cloud property information for climate model diagnostics. Comparison of EPIC-derived seasonal and longitudinal variability of the Earth’s planetary albedo with the GISS ModelE2 results shows systematic overestimate of cloud reflectivity over the Pacific Ocean with corresponding underestimates over continental land areas.

DSCOVR↗

On-Orbit Radiometric Characterization of OLI (Landsat 8) for Applications in Aquatic Remote Sensing

Landsat-8 carries two separate sensors, namely the Operational Land Imager (OLI) and the Thermal Infrared Radiometer Suite (TIRS), that image the earth surface throughout the visible and thermal portions of the spectrum. Compared to Landsat heritage sensors, the OLI has enhanced features, which include its 12-bit radiometric resolution and the addition of a band centered at 443 nm. The dramatically improved data quality/quantity expands existing applications of Landsat imagery in aquatic sciences from the retrieval of bio-geochemical properties, such as near-surface concentrations of chlorophyll-a (CHL) and total suspended solids (TSS), to benthic mapping. This study offers analysis of OLI's absolute radiometric performance over bodies ofwater using benchmark observations, namely the top-of-atmosphere (TOA) ocean color observations and marine in situ radiometric measurements. Sensor-to-sensor comparisons are performed to derive gain factors (g1) from near-concurrent observations in TOA radiance and reflectance domains. The gains in the radiance domain were further validated/adjusted by determining a second set of gains (g2) via analysis of OLI-derived water-leaving radiance, i.e., L(sub w) (gamma), against in situ measurements made either at the Ocean Color AErosol RObotic NETwork (AERONETOC) sites or during field campaigns. The analyses yield the OLI calibration uncertainties that need to be accounted forwhen studying aquatic environments. Itwas found that, for the visible and near-infrared channels, the OLI radiometric responses, on average, are well in agreement (b 2 % discrepancies) with the TOA radiances estimated by ocean color satellites or those predicted by models based onmeasurements of aquatic and atmospheric properties. However, the TOA radiance at the new 443-nm band is found to be, on average, 3.4 % larger than the reference observations. The inter-sensor comparisons in the reflectance domain, however, indicated slightly different results with the OLI responses being low in the blue bands. To enhance the retrieval accuracy of aquatic-science products from OLI datasets, sets of temporally averaged gains (radiance and reflectance) are derived and recommended for use prior to the retrieval of in-water products.

Remote Sensing↗

Simulating PACE Global Ocean Radiances

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P < 0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and a CDOC (r =0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250-800 nm. These unassimilated radiances were within 0.074 mW/sq cm/micron/sr of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of 10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PACE, ocean color, water-leaving radiances, biogeochemical model, radiative transfer model

PIC↗

Compression of oceanic lithosphere - An analysis of intraplate deformation in the Central Indian Basin

The development of intraplate structure in the Central Indian Basin is examined using models in which deformation is due to flexural buckling and the hydrodynamic growth of instabilities. Comparison of the models reveal that in a strong viscous lithosphere, deformation of the layer occurs by flexural folding at a wavelength which corresponds to the flexural buckling theory; in a lithosphere of intermediate strength, the layer deforms by folding characterized by thickening which localizes beneath topographic heights; and in a relatively weak lithosphere, the layer incurs an even greater amount of localized thickening and deforms in the symmetric or pinch-and-swell mode by inverse boudinage. It is noted that the models in which the layer folds either flexurally or with periodic thickening correspond with observed depth distribution of seismicity in the Central Indian Basin and with experimental rock rheological data.

Zuber, Maria T.↗

Radiative Flux Changes by Aerosols from North America, Europe, and Africa over the Atlantic Ocean: Measurements and Calculations from TARFOX and ACE-2

Aerosol effects on atmospheric radiative fluxes provide a forcing function that is a major source of uncertainty in understanding the past climate and predicting climate change. To help reduce this uncertainty, the 1996 Tropospheric Aerosol Radiative Forcing Experiment (TARFOX) and the 1997 second Aerosol Characterization Experiment (ACE-2) measured the properties and radiative effects of American, European, and African aerosols over the Atlantic. In TARFOX, radiative fluxes and microphysics of the American aerosol were measured from the UK C-130 while optical depth spectra, aerosol composition, and other properties were measured by the University of Washington C-131A and the CIRPAS Pelican. Closure studies show that the measured flux changes agree with those derived from the aerosol measurements using several modelling approaches. The best-fit midvisible single-scatter albedos (approx. 0.89 to 0.93) obtained from the TARFOX flux comparisons are in accord with values derived by independent techniques. In ACE-2 we measured optical depth and extinction spectra for both European urban-marine aerosols and free-tropospheric African dust aerosols, using sunphotometers on the R/V Vodyanitskiy and the Pelican. Preliminary values for the radiative flux sensitivities (Delta Flux / Delta Optical depth) computed for ACE-2 aerosols (boundary layer and African dust) over ocean are similar to those found in TARFOX. Combining a satellite-derived optical depth climatology with the aerosol optical model validated for flux sensitivities in TARFOX provides first-cut estimates of aerosol-induced flux changes over the Atlantic Ocean.

Russell, P. B.↗

A comparison between general circulation model simulations using two sea surface temperature datasets for January 1979

Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large-scale circulation anomalies with a pattern that resembles that obtained using COADS SSTs as well as EOF 1 of the control simulation; the warm pool in the Atlantic does not. These results suggest that the differences obtained with COADS SSTs and HIRS SSTs are mostly due to the differences in the datasets over the northern Pacific. There was a blocking episode near Greenland in late January 1979. Both simulations with warm SST anomalies over the northwest and northeast Pacific show a tendency toward increased incidence of North Atlantic blocking; the simulation with warm SST anomalies over the northwest Atlantic shows a tendency toward decreased incidence. These results suggest that features in both SST datasets that do not have a counterpart in the other dataset contribute signficantly to the differences between the simulated and observed fields. The results of this study imply that uncertainties in current SST distributions for the world oceans can be as important as the SST anomalies themselves in terms of their impact on the atmospheric circulation. Caution should be exercised, therefore, when linking anomalous circulation and SST patterns, especially in long-range prediction.

Ose, Tomoaki↗

A comparative study on synthetic aperture radar and optical imagery of ocean waves

The ocean imaging capability of synthetic aperture radar (SAR) is examined in terms of an intercomparison between simultaneous, spatially collocated images of a coastal ocean scene, as recorded by aerial photography and aircraft-borne SAR. Based on a detailed study of these images, this paper examines various local features of oceanographic interest as they are imaged by the two methods. This provides information not available in the usual comparisons of normalized SAR and single point buoy spectra. In addition, aspects of the intensity wavenumber spectra of the two images are examined and discussed. An interesting and surprising conclusion arising from the comparison study is the ability of the SAR to image accurately non-uniform (possibly nonlinear) wave features which are imaged by the aerial photography. This result should have important implications on the development of SAR imaging models and the interpretation of the effect of ocean surface motion on the SAR images.

Banner, M. L.↗

Evaluating the Dominant Components of Warming in Pliocene Climate Simulations

The Pliocene Model Intercomparison Project (PlioMIP) is the first coordinated climate model comparison for a warmer palaeoclimate with atmospheric CO2 significantly higher than pre-industrial concentrations. The simulations of the mid-Pliocene warm period show global warming of between 1.8 and 3.6 C above pre-industrial surface air temperatures, with significant polar amplification. Here we perform energy balance calculations on all eight of the coupled ocean-atmosphere simulations within PlioMIP Experiment 2 to evaluate the causes of the increased temperatures and differences between the models. In the tropics simulated warming is dominated by greenhouse gas increases, with the cloud component of planetary albedo enhancing the warming in most of the models, but by widely varying amounts. The responses to mid-Pliocene climate forcing in the Northern Hemisphere midlatitudes are substantially different between the climate models, with the only consistent response being a warming due to increased greenhouse gases. In the high latitudes all the energy balance components become important, but the dominant warming influence comes from the clear sky albedo, only partially offset by the increases in the cooling impact of cloud albedo. This demonstrates the importance of specified ice sheet and high latitude vegetation boundary conditions and simulated sea ice and snow albedo feedbacks. The largest components in the overall uncertainty are associated with clouds in the tropics and polar clear sky albedo, particularly in sea ice regions. These simulations show that albedo feedbacks, particularly those of sea ice and ice sheets, provide the most significant enhancements to high latitude warming in the Pliocene.

climate models↗

Hydrodynamics of internal solitons and a comparison of SIR-A and SIR-B data with ocean measurements

Large internal solitary waves have been observed by Shuttle SIR-A and SIR-B at locations in the Andaman Sea and the New York Bight. Satellite imagery and oceanographic measurements are used in conjunction with hydrodynamic interaction and electromagnetic scattering models to estimate the expected SAR image intensity modulations associated with the internal waves. There is reasonable agreement between the predicted and observed internal wave signatures.

Apel, J. R.↗

Air-Sea Interactions in a High-Resolution Ocean-Atmosphere Simulation

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. We present here results from a coupled GEOS-MIT "nature run" simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere simulation to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the coupled versus uncoupled differences in interactions between Sea Surface Temperature (SST) and ocean surface wind. We discuss, in particular, a several-day mode of temporal variability in the SST-wind cycle and how it is represented in the different model simulations and in observationally-based products. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗