Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wind forecasting”

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 811 records · Page 45

GC13I-0857: Designing a Frost Forecasting Service for Small Scale Tea Farmers in East Africa

Kenya is the third largest tea exporter in the world, producing 10% of the world's black tea. Sixty percent of this production occurs largely by small scale tea holders, with an average farm size of 1.04 acres, and an annual net income of $1,075. According to a recent evaluation, a typical frost event in the tea growing region causes about $200 dollars in losses which can be catastrophic for a small holder farm. A 72-hour frost forecast would provide these small-scale tea farmers with enough notice to reduce losses by approximately 80 USD annually. With this knowledge, SERVIR, a joint NASA-USAID initiative that brings Earth observations for improved decision making in developing countries, sought to design a frost monitoring and forecasting service that would provide farmers with enough lead time to react to and protect against a forecasted frost occurrence on their farm. SERVIR Eastern and Southern Africa, through its implementing partner, the Regional Centre for Mapping of Resources for Development (RCMRD), designed a service that included multiple stakeholder engagement events whereby stakeholders from the tea industry value chain were invited to share their experiences so that the exact needs and flow of information could be identified. This unique event allowed enabled the design of a service that fit the specifications of the stakeholders. The monitoring service component uses the MODIS Land Surface Temperature product to identify frost occurrences in near-real time. The prediction component, currently under testing, uses the 2-m air temperature, relative humidity, and 10-m wind speed from a series of high-resolution Weather Research and Forecasting (WRF) numerical weather prediction model runs over eastern Kenya as inputs into a frost prediction algorithm. Accuracy and sensitivity of the algorithm is being assessed with observations collected from the farmers using a smart phone app developed specifically to report frost occurrences, and from data shared through our partner network developed at the stakeholder engagement meeting. This presentation will illustrate the efficacy of our frost forecasting algorithm, and a way forward for incorporating these forecasts in a meaningful way to the key decision makers - the small-scale farmers of East Africa.

frost↗

Structure of the Highly Sheared Tropical Storm Chantal During CAMEX -4

On 20 August 2001 during the Convection and Moisture Experiment 4 (CAMEX-4) and NOAA Hurricane Field Program (HFP2001), the NASA high-altitude ER-2 and medium-altitude DC-8, and lower-altitude NOAA P3 aircraft conducted a coordinated Quantitative Precipitation Estimation (QPE) mission focused on convection in Tropical Storm Chantal. This storm first became a depression on 14 August, a tropical storm on 17 August, and it maintained maximum winds of about 65-70 mph during 19-20 August with minimum pressures ranging from 1008 mb on 19 August to 1001 mb late on 20 August. The storm was westward moving and was forecasted to intensify and landfall near the Yucatan-Belize border late on 20 August. Chanter failed to intensify and instead exhibited a highly sheared structure with an open low-level circulation and intense convection well to the northeast of this circulation center. The NASA ER-2 and DC-8 aircraft were closely coordinated with the NOAA P3 (NOAA-42). The NASA aircraft collected remote sensing and in situ data sets, while the P3 collected lower level in situ and radar data; both the DC-8 and P3 released 7 and 24 dropsondes, respectively. These aircraft measurements provided a unique opportunity to examine the structure of a sheared system and why it did not develop as forecasted a few days earlier. This paper will describe a preliminary study of the precipitation and wind structure provided by the NASA aircraft within the context of the NOAA P3 measurements.

Heymsfield, Gerald M.↗

Validation of turbulence intensity as simulated by the Weather Research and Forecasting model off the US northeast coast

We report turbulence intensity (TI) is often used to quantify the strength of turbulence in wind energy applications and serves as the basis of standards in wind turbine design. Thus, accurately characterizing the spatiotemporal variability in TI should lead to improved predictions of power production. Nevertheless, turbulence measurements over the ocean are far less prevalent than over land due to challenges in instrumental deployment, maintenance, and operation. Atmospheric models such as mesoscale (weather prediction) and large-eddy simulation (LES) models are commonly used in the wind energy industry to assess the spatial variability of a given site. However, the TI derivation from atmospheric models has not been well examined. An algorithm is proposed in this study to realize online calculation of TI in the Weather Research and Forecasting (WRF) model. Simulated TI is divided into two components depending on scale, including sub-grid (parameterized based on turbulence kinetic energy (TKE)) and grid resolved. The sensitivity of sea surface temperature (SST) on simulated TI is also tested. An assessment is performed by using observations collected during a field campaign conducted from February to June 2020 near the Woods Hole Oceanographic Institution Martha's Vineyard Coastal Observatory. Results show that while simulated TKE is generally smaller than the lidar-observed value, wind speed bias is usually small. Overall, this leads to a slight underestimation in sub-grid-scale estimated TI. Improved SST representation subsequently reduces model biases in atmospheric stability as well as wind speed and sub-grid TI near the hub height. Large TI events in conjunction with mesoscale weather systems observed during the studied period pose a challenge to accurately estimating TI from models. Due to notable uncertainty in accurately simulating those events, this suggests summing up sub-grid and resolved TI may not be an ideal solution. Efforts in further improving skills in simulating mesoscale flow and cloud systems are necessary as the next steps.

17 WIND ENERGY↗

Three Western Pacific Typhoons Strengthened Fire Weather in the Recent Northwest U.S. Conflagration

Abstract Strong winds that accentuated a fire outbreak in the western United States in early September of 2020 resulted from an atmospheric wave train that spanned the Pacific Ocean. Days before the atmospheric waves developed in the United States, three western Pacific tropical cyclones (typhoons) underwent an extratropical transition over Korea within an unprecedentedly short span of 12 days. Using a climate diagnostic approach and historical forecast data, it was found that the amplitude of the atmospheric waves accompanying the western U.S. fire weather would not have been so profound if not for the influence of these typhoons. Together, the recurving typhoons provided a significant source of wave activity flux directed toward North America – amplifying the ridge over the U.S. west coast while deepening the trough in central Canada. This anomalous circulation produced the severe frontal system that caused extreme winds in western Oregon, Washington and California – rapidly spreading fire.

58 GEOSCIENCES↗

A Public-Private-Acadmic Partnership to Advance Solar Power Forecasting

The National Center for Atmospheric Research (NCAR) is pleased to have led a partnership to advance the state-of-the-science of solar power forecasting by designing, developing, building, deploying, testing, and assessing the SunCast™ Solar Power Forecasting System. The project has included cutting edge research, testing in several geographically- and climatologically-diverse high penetration solar utilities and Independent System Operators (ISOs), and wide dissemination of the research results to raise the bar on solar power forecasting technology. The partners include three other national laboratories, six universities, and industry partners. This public-private-academic team has worked in concert to perform use-inspired research to advance solar power forecasting through cutting-edge research to advance both the necessary forecasting technologies and the metrics for evaluating them. The project has culminated in a year-long, full-scale demonstration of provide irradiance and power forecasts to utilities and ISOs to use in their operations. The project focused on providing elements of a value chain, beginning with the weather that causes a deviation from clear sky irradiance and progresses through monitoring and observations, modeling, forecasting, dissemination and communication of the forecasts, interpretation of the forecast, and through decision-making, which produces outcomes that have an economic value. The system has been evaluated using metrics developed specifically for this project, which has provided rich information on model and system performance. Research was accomplished on the very short range (0-6 hours) Nowcasting system as well as on the longer term (6-72 hour) forecasting system. The shortest range forecasts are based on observations in the field. The shortest range system, built by Brookhaven National Laboratory (BNL) and based on Total Sky Imagers (TSIs) is TSICast, which operates on the shortest time scale with a latency of only a few minutes and forecasts that currently go out to about 15 min. This project has facilitated research in improving the hardware and software so that the new high definition cameras deployed at multiple nearby locations allow discernment of the clouds at varying levels and advection according to the winds observed at those levels. Improvements over “smart persistence” are about 29% for even these very short forecasts. StatCast is based on pyranometer data measured at the site as well as concurrent meteorological observations and forecasts. StatCast is based on regime-dependent artificial intelligence forecasting techniques and has been shown to improve on “smart persistence” forecasts by 15-50%. A second category of short-range forecasting systems employ satellite imagery and use that information to discern clouds and their motion, allowing them to project the clouds, and the resulting blockage of irradiance, in time. CIRACast (the system produced by the Cooperative Institute for Atmospheric Research [CIRA] at Colorado State University) was already one of the more advanced cloud motion systems, which is the reason that team was brought to this project. During the project timeframe, the CIRA team was able to advance cloud shadowing, parallax removal, and implementation of better advecting winds at different altitudes. CIRACast shows generally a 25-40% improvement over Smart Persistence between sunrise and approximately 1600 UTC (Coordinated Universal Time) . A second satellite-based system, MADCast (Multi-sensor Advective Diffusive foreCast system), assimilates data from multiple satellite imagers and profilers to assimilate a fully three-dimensional picture of the cloud into the dynamic core of WRF. During 2015, MADCast (provided at least 70% improvement over Smart Persistence, with most of that skill being derived during partly cloudy conditions. That allows advection of the clouds via the Weather Research and Forecasting (WRF) model dynamics directly. After WRF-Solar™ showed initial success, it was also deployed in nowcasting mode with coarser runs out to 6 hours made hourly. It provided improvements on the order of 50-60% over Smart Persistence for forecasts up to 1600 UTC. The advantages of WRF-Solar-Nowcasting and MADCast were then blended to develop the new MAD-WRF model that incorporates the most important features of each of those models, both assimilating satellite cloud fields and using WRF-So far physics to develop and dissipate clouds. MAE improvements for MAD-WRF for forecasts from 3-6 hours are improved over WRF-Solar-Now by 20%. While all the Nowcasting system components by themselves provide improvement over Smart Persistence, the largest benefit is derived when they are smartly blended together by the Nowcasting Integrator to produce an integrated forecast. The development of WRF-Solar™ under this project has provided the first numerical weather prediction (NWP) model specifically designed to meet the needs of irradiance forecasting. The first augmentation improved the solar tracking algorithm to account for deviations associated with the eccentricity of the Earth’s orbit and the obliquity of the Earth. Second, WRF-Solar™ added the direct normal irradiance (DNI) and diffuse (DIF) components from the radiation parameterization to the model output. Third, efficient parameterizations were implemented to either interpolate the irradiance in between calls to the expensive radiative transfer parameterization, or to use a fast radiative transfer code that avoids computing three-dimensional heating rates but provides the surface irradiance. Fourth, a new parameterization was developed to improve the representation of absorption and scattering of radiation by aerosols (aerosol direct effect). A fifth advance is that the aerosols now interact with the cloud microphysics, altering the cloud evolution and radiative properties, an effect that has been traditionally only implemented in atmospheric computationally costly chemistry models. A sixth development accounts for the feedbacks that sub-grid scale clouds produce in shortwave irradiance as implemented in a shallow cumulus parameterization Finally, WRF-Solar™ also allows assimilation of infrared irradiances from satellites to determine the three dimensional cloud field, allowing for an improved initialization of the cloud field that increases the performance of short-range forecasts. We find that WRF-Solar™ can improve clear sky irradiance prediction by 15-80% over a standard version of WRF, depending on location and cloud conditions. In a formal comparison to the NAM baseline, WRF-Solar™ showed improvements in the Day-Ahead forecast of 22-42%. The SunCast™ system requires substantial software engineering to blend all of the new model components as well as existing publically available NWP model runs. To do this we use an expert system for the Nowcasting blender and the Dynamic Integrated foreCast (DICast®) system for the NWP models. These two systems are then blended, we use an empirical power conversion method to convert the irradiance predictions to power, then apply an analog ensemble (AnEn) approach to further tune the forecast as well as to estimate its uncertainty. The AnEn module decreased RMSE (root mean squared error) by 17% over the blended SunCast™ power forecasts and provided skill in the probabilistic forecast with a Brier Skill Score of 0.55. In addition, we have also developed a Gridded Atmospheric Forecast System (GRAFS) in parallel, leveraging cost share funds. An economic evaluation based on Production Cost Modeling in the Public Service Company of Colorado showed that the observed 50% improvement in forecast accuracy will save their customers $819,200 with the projected MW deployment for 2024. Using econometrics, NCAR has scaled this savings to a national level and shown that an annual expected savings for this 50% forecast error reduction ranges from $11M in 2015 to $43M expected in 2040 with increased solar deployment. This amounts to a $455M discounted savings over the 26 year period of analysis.

14 SOLAR ENERGY↗

Chapter 15: Reliability of Wind Turbines

The global wind industry has witnessed exciting developments in recent years. The future will be even brighter with further reductions in capital and operation and maintenance costs, which can be accomplished with improved turbine reliability, especially when turbines are installed offshore. One opportunity for the industry to improve wind turbine reliability is through the exploration of reliability engineering life data analysis based on readily available data or maintenance records collected at typical wind plants. If adopted and conducted appropriately, these analyses can quickly save operation and maintenance costs in a potentially impactful manner. This chapter discusses wind turbine reliability by highlighting the methodology of reliability engineering life data analysis. It first briefly discusses fundamentals for wind turbine reliability and the current industry status. Then, the reliability engineering method for life analysis, including data collection, model development, and forecasting, is presented in detail and illustrated through two case studies. The chapter concludes with some remarks on potential opportunities to improve wind turbine reliability. An owner and operator's perspective is taken and mechanical components are used to exemplify the potential benefits of reliability engineering analysis to improve wind turbine reliability and availability.

17 WIND ENERGY↗

Chapter 14 - Reliability of Wind Turbines

The global wind energy industry has grown at a fast pace during the past half-decade. Advancements from design and manufacturing to operation and maintenance have led to reduced capital and maintenance costs, which make wind power an indispensable source for a comprehensive solution to global electricity needs. Once wind turbines are installed, the opportunity to lower wind power costs is mainly through improved operation and maintenance practices. Modern wind turbines are equipped with tens or hundreds of measurement channels and are generating an abundance of data, with lot of efforts being put into data analysis by both the research community and the industry. One type of analysis is through the exploration of reliability engineering methods based on readily available data or maintenance records collected at typical wind power plants. If adopted and conducted appropriately, these analyses can quickly save operation and maintenance costs in a potentially impactful manner. The wind industry has adopted this discipline more broadly in recent years. This chapter discusses wind turbine reliability by highlighting the methodology of reliability engineering life data analysis. It first briefly discusses the fundamentals of wind turbine reliability and the current industry status. Then, the reliability engineering method for life analysis, including data collection, model development, and forecasting, is presented in detail and illustrated through two case studies. The chapter concludes with some remarks on potential opportunities to improve wind turbine reliability. An owner and operator's perspective is taken and mechanical components are used to exemplify the potential benefits of reliability engineering analysis to improve wind turbine reliability and availability.

database↗

POD-based constrained sensor placement and field reconstruction from noisy wind measurements: A perturbation study

Sensor placement at the extrema of Proper Orthogonal Decomposition (POD) is efficient and leads to accurate reconstruction of the wind field from a limited number of measure- ments. In this paper we extend this approach of sensor placement and take into account measurement errors and detect possible malfunctioning sensors. We use the 48 hourly spa- tial wind field simulation data sets simulated using the Weather Research an Forecasting (WRF) model applied to the Maine Bay to evaluate the performances of our methods. Specifically, we use an exclusion disk strategy to distribute sensors when the extrema of POD modes are close. It turns out that this strategy can also reduce the error of recon- struction from noise measurements. Also, by a cross-validation technique, we successfully locate the malfunctioning sensors.

54 ENVIRONMENTAL SCIENCES↗

Aircraft measurements and analysis of severe storms: 1976 field experiment

Severe storm aircraft measurements are documented, as well as the instrumentation and operational features of aircraft mobility capabilities. The measurements and data analyses indicate that the concept of a highly mobile research aircraft capability for obtaining detailed measurements of wind, temperature, moisture, spherics, etc., near and within severe storm systems, forecast 48 hours in advance in a 1000 nm operating radius, is feasible, and was successfully demonstrated. The measurements and analyses reveal several severe storm features and insights with respect to storm air flow circulations and inflow-outflow orientation. Precipitation downdraft air is recirculated back into the updraft core below the scud cloud in both back and front feeder type storms. In a back feeder type storm, the downdraft outflow air ahead of the storm is also recirculated back into the updraft region near cloud base.

Sinclair, P. C.↗

The effects of tropical wind data on the prediction of ultralong waves

The influence of tropical wind data on the prediction of planetary waves were studied. Two assimilation experiments were performed, one with and one without FGGE tropical winds. The planetary wave error was then analyzed in 72 h forecasts from the initial conditions provided by the two assimilations.

Baker, W. E.↗

The impact of Seasat-A scatterometer data on the numerical prediction of the Queen Elizabeth II storm

The development of the Queen Elizabeth II storm from 1200 UT on September 9, 1978, to 0000 UT on September 11, 1978, is forecast with a limited-area, fine resolution (100 km) numerical weather prediction model. Forecasts from initial conditions that include and exclude Seasat scatterometer surface winds are performed. When the Seasat winds are allowed to influence the upper levels of the atmosphere, a large positive impact is found.

Duffy, D. G.↗

Simulation studies of the impact of future observing systems on weather prediction

The features and preliminary results from a simulation system being implemented to develop realistic estimates of the impacts future data acquisition systems will have on large-scale numerical weather simulation are described. The new instruments may include advanced passive IR and microwave satellite sensors, as well as active scatterometer and lidar sounders. A main goal of the impact study is to identify those sensor systems which will provide the most benefit. The realism of the study is being enhanced by assimilating as much real-world data as possible and generating global weather maps for comparison with maps generated on the bases on the projected new, higher resolution data. Early results have indicated a preference for higher resolution wind data than for temperature data for making 1-5 day forecasts. The prime instrument candidate for collecting the data is lidar, provided the sensor resolution design goals are met.

Atlas, R.↗

The atmospheric excitation of earth orientation changes during MERIT

Geodetic estimates of earth orientation variations detected during the MERIT campaign (September 1983-November 1984) are compared with corresponding meteorological data. The geodetic data were obtained using VLBI, lunar laser ranging, and satellite laser ranging and the meteorological data were from the NMC in the U.S. and the European Centre for Medium Range Weather Forecasting (ECMRWF) in the U.K. The effects of changes in pressure, wind, and the inverted barometer ocean response on the excitation of the polar motion and the length of day are examined. The comparison between the meteorological and geodetic data reveals that the equatorial vector component of the atmospheric angular momentum (AAM) has significant annual and semiannual variations. Good correlation is detected between the geodetic polar motion data and the NMC pressure inverted barometer data combined with the ECMRWF wind estimates and semiannual agreement with the NMC pressure data is observed. There is also good correlation between the ECMRWF and NMC polar vector component of the AAM data and the geodetic length of day estimates, and good semiannual agreement with NMC pressure data is noted.

Eubanks, T. M.↗

An atlas of monthly mean distributions of GEOSAT sea surface height, SSMI surface wind speed, AVHRR/2 sea surface temperature, and ECMWF surface wind components during 1988

Monthly mean global distributions for 1988 are presented with a common color scale and geographical map. Distributions are included for sea surface height variation estimated from GEOSAT; surface wind speed estimated from the Special Sensor Microwave Imager on the Defense Meteorological Satellite Program spacecraft; sea surface temperature estimated from the Advanced Very High Resolution Radiometer on NOAA spacecrafts; and the Cartesian components of the 10m height wind vector computed by the European Center for Medium Range Weather Forecasting. Charts of monthly mean value, sampling distribution, and standard deviation value are displayed. Annual mean distributions are displayed.

Halpern, D.↗

Application of a linear spectral model to the study of Amazonian squall lines during GTE/ABLE 2B

A linear nonhydrostatic spectral model is run with the basic state, or large scale, vertical profiles of temperature and wind observed prior to convective development along the northern coast of South America during the GTE/ABLE 2B. The model produces unstable modes with mesoscale wavelength and propagation speed comparable to observed Amazonian squall lines. Several tests with different vertical profiles of low-level winds lead to the conclusion that a shallow and/or weak low-level jet either does not produce a scale selection or, if it does, the selected mode is stationary, indicating the absence of a propagating disturbance. A 700-mbar jet of 13 m/s, with a 600-mbar wind speed greater or equal to 10 m/s, is enough to produce unstable modes with propagating features resembling those of observed Amazonian squall lines. However, a deep layer of moderate winds (about 10 m/s) may produce similar results even in the absence of a low-level wind maximum. The implications in terms of short-term weather forecasting are discussed.

Silva Dias, Maria A. F.↗

Principal modes of atmospheric circulation anomalies associated with global angular momentum fluctuations

This paper provides a description of the variability of global atmospheric angular momentum (GAM) and its relationship with principal modes of three-dimensional atmospheric circulation anomalies. The data used are 5-day mean global wind fields from the European Centre for Medium-Range Weather Forecasts initialized dataset for 1980-1989. Significant seasonal variation of GAM is observed with maxima in April and November and a minimum during late July. The amplitude of the annual cycle is largest in the upper troposphere and decreases toward the surface. Although the lower tropospheric contribution to the total angular momentum is relatively small, its annual cycle is out of phase with those of the upper atmosphere and GAM. Also identified is a distinct semiannual component, with double peaks appearing in April and November. This signal is most noticeable in the upper troposphere above the 300-mb level. The principal modes of zonal-mean angular momentum and meridional circulation anomalies and their coupled modes are obtained by using empirical orthogonal function analysis and singular value decomposition. It is shown that the leading modes of the angular momentum and meridional circulation are coupled with each other and are responsible for much of the variability in GAM. The coupled modes represent fluctuations of upper-level subtropical zonal flow, which are linked to the modulation of Hadley circulation intensity in both hemispheres. It is found that GAM is highly correlated with the first eigenvector of upper-level streamfunction anomalies, which consists of a superrotational flow in the tropics and subtropics, except over the central Pacific where a 'blocked' flow with two subtropical anticyclonic circulation cells straddling the equator is found. Much of the blocked flow is due to the establishment of dipole anomalies in the velocity potential with centers over the central Pacific and the Maritime Continent on the interannual time scale. On the intraseasonal time scale, GAM fluctuation is dominated by superrotational flow in the tropics, with the blocked flow present to a much lesser extent. The associated velocity potential anomaly has a weak dipole structure with centers over the Indian Ocean and the eastern Pacific. The implications of the above results on the total angular momentum balance of the earth-atmosphere system are also discussed.

Kang, In-Sik↗

TOGA COARE Satellite data summaries available on the World Wide Web

Satellite data summary images and analysis plots from the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE), which were initially prepared in the field at the Honiara Operations Center, are now available on the Internet via World Wide Web browsers such as Mosaic. These satellite data summaries consist of products derived from the Japanese Geosynchronous Meteorological Satellite IR data: a time-size series of the distribution of contiguous cold cloudiness areas, weekly percent high cloudiness (PHC) maps, and a five-month time-longitudinal diagram illustrating the zonal motion of large areas of cold cloudiness. The weekly PHC maps are overlaid with weekly mean 850-hPa wind calculated from the European Centre for Medium-Range Weather Forecasts (ECMWF) global analysis field and can be viewed as an animation loop. These satellite summaries provide an overview of spatial and temporal variabilities of the cloud population and a large-scale context for studies concerning specific processes of various components of TOGA COARE.

Chen, S. S.↗

Sea-air of CO2 in the North Pacific using shipboard and satellite data

A method has been developed to produce high-resolution maps of pCO2 in surface water for the North Pacific using satellite sea surface temperature (SST) data and statistical relationships between measured pCO2 and temperature. In the subtropical North Pacific the pCO in seawater is controlled primarily by temperature. Accordingly, pCO2 values that are calculated from the satellite SST data have good agreement with the measured values (rms deviation of +/- microatm). In the northwestern subpolar region the pCO2 is controlled not only by temperature, but also by significant seasonal changes in the total CO2 concentration, which are caused by seasonal changes in primary production, mixing with subsurface waters and sea-air exchange. Consequently, the parameterization of oceanic p CO2 based on SST data alone is not totally successful in the northwestern region (rms deviation of +/- 40 microatm). The use of additional satellite products, such as wind and ocean color data, as planned for a future study, is considered necessary to account for the pCO2 variability caused by seasonal changes in the total CO2 concentration. The net CO2 flux for the area of the North Pacific included in this study (north of 10 deg N) has been calculated using the monthly pCO2 distributions computed, and monthly wind speeds from the European Centre for Medium-Range Weather Forecasts. The region is found to be a net source to the atmosphere of 1.9 x 10(exp 12) to 5.8 x 10(exp 12) moles of CO2 per year (or 0.02-0.07 Gt C/yr), most of the outflux occurring in the subtropics.

Stephens, Mark P.↗