Search NASA⌕ Search

SEARCH · Search NASA

Results for “Upper atmosphere”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stellar wind impact on early atmospheres around unmagnetized Earth-like planets

ABSTRACT Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability.

Astronomy & Astrophysics↗
DOE OSTI2020

ACCESS: Confirmation of No Potassium in the Atmosphere of WASP-31b

We present a new optical (400–950 nm) transmission spectrum of the hot Jupiter WASP-31b (M = 0.48 M {sub J}; R = 1.54 R {sub J}; P = 3.41 days), obtained by combining four transit observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet, as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na i and K i, with particular focus on K i, for which there have been two previously published disagreeing results. Observations with Hubble Space Telescope (HST)/STIS detected K i, whereas ground-based low- and high-resolution observations did not. We use equilibrium and nonequilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (α=5.3{sub −3.1}{sup +2.9}), high-altitude clouds at a log cloud top pressure of −3.6{sub −2.1}{sup +2.7} bars, and possible muted H{sub 2}O features. We find that our observations support other ground-based claims of no K i. Clouds are likely why signals like H{sub 2}O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (α = 8.3 ± 1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Effects of Global and Regional Dust Storms on the Martian Hot O Corona and Photochemical Loss

Here, we examine for the first time the effects of both global and regional dust storms on the formation of the Martian hot O corona and associated photochemical loss of O. Our study is conducted by utilizing our integrated model framework, which couples our Martian hot O corona model with a multifluid magnetohydrodynamic model for Mars for the dusty and clear atmospheric condition cases. We present our results with the most up-to-date cross sections for the O( 3 P )-CO 2 collisions. The main effect of dust storms on the ionosphere is the upward shift of the ionosphere on the dayside, which results in an increase in production of hot O at all altitudes above the ionospheric peak. However, the dust-induced inflation of the neutral upper atmosphere results in an enhancement in collisional loss of hot O and thus effectively suppresses the hot O density, reducing the global photochemical loss rate by ~28% for the global dust storm scenario. The relative density structure of the hot O corona does not show any significant changes, while its magnitude decreases at all altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

EMISSA (Exploring Millimeter Indicators of Solar-Stellar Activity): I. The initial millimeter–centimeter main-sequence star sample

Due to their wide wavelength coverage across the millimeter to centimeter (mm–cm) range and their increased sensitivity, modern interferometric arrays facilitate observations of the thermal and non-thermal radiation that is emitted from different layers in the outer atmospheres of stars. We study the spectral energy distribution (S obs (ν)) of main-sequence stars based on archival observations in the mm–cm range with the aim to study their atmospheric stratification as a function of stellar type. The main-sequence stars with significant detection in mm bands were identified in the ALMA Science Archive. These data were then complemented with spectral flux data in the extreme ultraviolet to cm range as compiled from various catalogues and observatory archives. We compared the resultant S obs (ν) of each star with a photospheric emission model (S mod (ν)) calculated with the PHOENIX code. The departures of S obs (ν) from S mod (ν) were quantified in terms of a spectral flux excess parameter (ΔS/S mod ) and studied as a function of stellar type. The initial sample consists of 12 main-sequence stars across a broad range of spectral types from A1 to M3.5 and the Sun-as-a-star as reference. The stars with T eff = 3000–7000 K (F–M type) showed a systematically higher S obs (ν) than S mod (ν) in the mm–cm range. Their ΔS/S mod exhibits a monotonic rise with decreasing frequency. The steepness of this rise is higher for cooler stars in the T eff = 3000–7000 K range, although the single fully convective star (T eff ~ 3000 K) in the sample deviates from this trend. Meanwhile, Sobs(ν) of the A-type stars agrees with S mod (ν) within errors. The systematically high ΔS/S mod in F–M stars indicates hotter upper atmospheric layers, that is, a chromosphere and corona in these stars, like for the Sun. The mm–cm ΔS/S mod spectrum offers a way to estimate the efficiency of the heating mechanisms across various outer atmospheric layers in main-sequence stars, and thereby to understand their structure and activity. We emphasise the need for dedicated surveys of main-sequence stars in the mm–cm range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2022

Meso-Scale Electrodynamic Coupling of the Earth Magnetosphere-Ionosphere System

Within the fully integrated magnetosphere-ionosphere system, many electrodynamic processes interact with each other. We review recent advances in understanding three major meso-scale coupling processes within the system: the transient field-aligned currents (FACs), mid-latitude plasma convection, and auroral particle precipitation. (1) Transient FACs arise due to disturbances from either dayside or nightside magnetosphere. As the interplanetary shocks suddenly compress the dayside magnetosphere, short-lived FACs are induced at high latitudes with their polarity successively changing. Magnetotail dynamics, such as substorm injections, can also disturb the current structures, leading to the formation of substorm current wedges and ring current disruption. (2) The mid-latitude plasma convection is closely associated with electric fields in the system. Recent studies have unraveled some important features and mechanisms of subauroral fast flows. (3) Charged particles, while drifting around the Earth, often experience precipitating loss down to the upper atmosphere, enhancing the auroral conductivity. Recent studies have been devoted to developing more self-consistent geospace circulation models by including a better representation of the auroral conductance. It is expected that including these new advances in geospace circulation models could promisingly strengthen their forecasting capability in space weather applications. The remaining challenges especially in the global modeling of the circulation system are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Effects of Field Line Curvature (FLC) Scattering on Ring Current Dynamics and Isotropic Boundary

Abstract In the ring current dynamics, various loss mechanisms contribute to the ring current decay, including losses to the upper atmosphere through particle precipitation. This study implements the field line curvature (FLC) scattering mechanism in a kinetic ring current model and investigates its role in precipitating ions into the ionosphere during the 17 March 2013 storm. Simulation results indicate that (1) the FLC scattering process exerts on energetic ions on the nightside where the magnetospheric configuration is more stretching. It is more effective on heavy ions (e.g., O + ). These ion losses thereafter lead to a faster recovery of the ring current. (2) The FLC‐associated ion precipitation mainly occurs in the outer region (L > 5 for protons and L > 4.5 for oxygen ions) on the nightside. The O + precipitation takes places in a wider region than protons although its intensity is much lower. Comparisons with POES observations suggest that more proton precipitation is needed in the inner region. This is probably caused by the less stretched configuration in the simulation that prevents more precipitation. It may also imply that other loss process is required in the model such as wave‐particle interactions. (3) The storm time precipitating proton flux of tens of keV due to the FLC scattering sometimes becomes comparable to that of electrons on the nightside, although electrons usually dominate the ionospheric energy deposition from the midnight eastward toward the dayside. (4) The FLC scattering process seems to be capable of explaining the formation of isotropic boundary in the ionosphere during the investigated event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

First Analysis of Earth's Stratospheric γ ‐Ray Emission in Geographical Coordinates With Fermi LAT

Abstract From low‐Earth orbit, the Earth is an extremely bright γ ‐ray source compared to distant astrophysical objects. The high‐energy γ ‐ray emission, from the Earth's upper atmosphere, is produced by the interactions between cosmic rays (CRs) and air molecules. As CRs are predominantly charged particles, the geomagnetic field acts as a spectrometer, forbidding CRs below a certain energy to arrive and interact with the atmosphere at a given geomagnetic location from a given direction. We use photon data between 0.2 and 20 GeV as observed by the Fermi Large Area Telescope (LAT) to create, for the first time, stratospheric γ ‐ray emission maps of the Earth in the geographical coordinate system, for different viewing directions from LAT altitude. These maps confirm a strong association between measured γ ‐ray fluxes and geomagnetic cutoffs for CRs. We also report the spectrum of the Earth's stratospheric γ ‐ray emission in cardinal directions and use the geomagnetic spectrometer effect to infer the atmosphere's yield function of γ rays per incoming CR proton as a function of CR rigidity at the altitude of the Fermi spacecraft. This work is relevant to the study of near‐Earth radiation exposure for astronauts and spacecraft, modeling the geomagnetic field, and CR‐air interactions.

Madlee, S.↗
DOE OSTI2022

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Informed Multi‐Scale Approach Applied to the British Columbia Fires of Late Summer 2017

Abstract Pyrocumulonimbus (PyroCb) clouds have a complex origin dependent on fire dynamics and meteorological conditions. When a pyrocumulonimbus cloud develops and is maintained over a period of time, it can inject significant aerosol into the troposphere and lower stratosphere, resulting in a longer‐term (months to years) occurrence of aerosol in the stratosphere. In this work, we investigate the British Columbia wildfires on 12–13 August 2017 using a multi‐scale simulation framework. We use the output of a physics‐based wildfire model (FIRETEC) with parameterized energy, particle, and gas emissions to drive the upper atmospheric aerosol mass injection within a regional cloud resolving model (HIGRAD). We demonstrate that vertical motions produced by latent heat release of the condensation of ice and cloud particles within the PyroCbs induce another 5 km of lifting of the simulated aerosol plume. Primary black carbon and organic aerosols (OAs) alone may not be enough to explain the observed aerosol burden, thus we show that secondary OA produced via condensation of gases by the fires, ash, and possibly dust can enhance lofted aerosol mass. A simulation with all emission mechanisms active, driven by the observed fuel load and environmental conditions, reasonably reproduces an aerosol profile inferred from observational data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Substorm Driven Chorus Waves: Decay Timescales and Implications for Pulsating Aurora

Abstract Energetic electron precipitation (EEP) associated with pulsating aurora can transfer greater than 30 keV electrons from the outer radiation belt region into the upper atmosphere and can deplete atmospheric ozone via collisions that produce NOx and HOx molecules. Our knowledge of exactly how EEP occurs is incomplete. Previous studies have shown that pitch angle scattering between electrons and lower‐band chorus waves can cause pulsating aurora associated with EEP and that substorms play an important role. In this work, we quantify the timescale of chorus wave decay following substorms and compare that to previously determined timescales. We find that the chorus decay e‐folding time varies based on magnetic local time (MLT), magnetic latitude, and wave frequency. The shortest timescales occur for lower‐band chorus in the 21 to 9 MLT region and compares, within uncertainty, to the energetic pulsating aurora timescale of Troyer et al. (2022, https://doi.org/10.3389/fspas.2022.1032552 ) for energetic pulsating aurora. We are able to further support this connection by modeling our findings in a quasi‐linear diffusion simulation. These results provide observations of how chorus waves behave after substorms and add additional statistical evidence linking energetic pulsating aurora to substorm driven lower‐band chorus waves.

Troyer, Riley N.↗
DOE OSTI2024

Exploring the Relative Importance of the MJO and ENSO to North Pacific Subseasonal Predictability

Abstract Here we explore the relative contribution of the Madden‐Julian Oscillation (MJO) and El Niño Southern Oscillation (ENSO) to midlatitude subseasonal predictive skill of upper atmospheric circulation over the North Pacific, using an inherently interpretable neural network applied to pre‐industrial control runs of the Community Earth System Model version 2. We find that this interpretable network generally favors the state of ENSO, rather than the MJO, to make correct predictions on a range of subseasonal lead times and predictand averaging windows. Moreover, the predictability of positive circulation anomalies over the North Pacific is comparatively lower than that of their negative counterparts, especially evident when the ENSO state is important. However, when ENSO is in a neutral state, our findings indicate that the MJO provides some predictive information, particularly for positive anomalies. We identify three distinct evolutions of these MJO states, offering fresh insights into opportune forecasting windows for MJO teleconnections.

58 GEOSCIENCES↗
DOE OSTI2021

Characteristics of thermospheric infrasound predicted using ray tracing and weakly non-linear waveform analyses

Infrasonic signals refracted by thermal gradients in the rarefied upper atmosphere are modeled using a combination of ray tracing and weak shock theory to develop an understanding of thermospheric infrasound signals produced by energetic, transient sources. Canonical arrival structures in the form of u-wave signatures are identified for returns refracted at lower altitudes within the thermosphere, and possible multi-pathing produced by effective sound speed inflections are investigated to elucidate more complex arrival structures, which are found to be spatially localized. Furthermore, variability in the source characteristics is investigated and it is found that whereas some waveform phase information is lost due to finite amplitude effects, arrival characteristics are strongly dependent on the peak overpressure near the source. Variability in the propagation path is considered using archived atmospheric specifications and implies that despite uncertainties related to the dynamic and sparsely sampled nature of the atmosphere, thermospheric signatures might be useful in estimating the yield for explosive sources. Last, thermospheric arrivals from a failed rocket launch, as well as those from several large chemical explosions, are analyzed and it is found that qualitative trends match those predicted, and analyses here provide additional insight into such signatures.

42 ENGINEERING↗
DOE OSTI2025

Measuring Local Turbulence Along the Optical Path: Multi-Beam Optical Seeing Sensor

Deflection of light along the optical path is a major source of image degradation for ground-based telescopes. Methods have been developed to measure upper atmospheric seeing based on models of the turbulence in the atmosphere, but due to boundary conditions, transmission within telescope enclosures is more complex. The Multi-beam Optical Seeing Sensor (MOSS) directly measures the component of the image quality degradation from inhomogeneity of the index of refraction within the telescope dome. MOSS outputs four near-parallel beams of light that travel along the optical path and are imaged by the telescope’s detector, landing like starlight on the telescope’s focal plane. By using a strobed light source, we can ‘freeze’ the instantaneous index variations transverse to the optical path. This system captures both ‘dome’ and ‘mirror’ seeing. Through plotting the standard deviation of differential motion between pairs of beams, MOSS enables characterization of the length scale of turbulence within the dome. The temporal coherence of temperature gradients can be probed with different pulse lengths, and the spatial coherence by comparing pairs at different separations across the aperture of the telescope. Optical path turbulence measurements, alongside other telemetry metrics, will guide thermal and airflow management to optimize image quality. A MOSS prototype was installed in the 1.2[Formula: see text]m Auxiliary Telescope (AuxTel) at the Vera C. Rubin Observatory in Chile, and preliminary data constrain the optical path turbulence with a lower bound of 1.4 arcsec. The optical path turbulence varied throughout the night of observing.

Astronomical seeing↗
DOE OSTI2024

Enhanced material identification via momentum-integrated muon scattering tomography

Cosmic ray muons, originating from interactions in the upper atmosphere, possess high energy and unique penetrative capabilities suitable for non-traditional radiographic inspection. This study explores their application in various fields such as nuclear fuel cask monitoring, nuclear reactor imaging, and archaeology, leveraging the principle of multiple Coulomb scattering for imaging dense materials. While muon scattering tomography has shown promise, accurately measuring muon momentum remains challenging. This research introduces the Momentum Integrated Point-of-Closest Approach (mPoCA) algorithm, integrating muon momentum data into the traditional Point-of-Closest Approach (PoCA) framework. Utilizing the Cherenkov muon spectrometer, renowned for precise muon momentum estimation, the mPoCA algorithm offers a novel imaging approach. Simulations conducted with GEANT4 evaluate the mPoCA algorithm’s performance against the standard PoCA method, demonstrating superior image resolution and enhanced material identification capabilities, particularly in distinguishing materials like uranium and lead. These findings underscore the potential of the mPoCA algorithm for advancing muon scattering tomography applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimizing photonic ring-resonator filters for OH-suppressed near-infrared astronomy

Near-infrared wavelength observations are crucial for understanding numerous fields of astrophysics, such as supernova cosmology and positronium annihilation detection. However, current ground-based observations suffer from an enormous background due to OH emission in the upper atmosphere. One promising way to solve this problem is to use ring-resonator filters to suppress OH emission lines. Here, we discuss our optimization of ring-resonator filter performance from five perspectives: resonance wavelength matching, polarization-independent operation, low insertion loss, low-loss coupling to astronomical instruments, and broadband operation. In the end, we discuss next steps needed for reliable supernova and positronium observations, thus providing a roadmap for future advances in near-infrared astronomy.

42 ENGINEERING↗
DOE OSTI2020

Geophysical and Planetary Acoustics on Balloons

Natural events and human activity often generate acoustic waves capable of traveling tens to tens of thousands of kilometers across the globe. Ground-based acoustic sensors are limited to dry land and often suffer from wind noise. In contrast, balloon borne acoustic sensors can cross oceans, polar ice caps, and other inhospitable areas, greatly expanding sensor coverage. Since they move with the mean wind speed, their background noise levels are exceptionally low. In the last six years, such sensors have recorded sounds from colliding ocean waves, surface and buried chemical explosions, thunder, wind/mountain interactions, wind turbines, aircraft, and possibly meteors and the aurora. These results have led to new insights on acoustic heating of the upper atmosphere, the detectability of underground explosions, and directional sound fields generated by ocean waves.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Puerto Rico Historical Climate Analysis: A closer look at complex tropical terrain

Small tropical islands like Puerto Rico are especially vulnerable to climate change impacts, yet are often underrepresented in most datasets because the spatial resolution is too coarse to sufficiently cover their complex terrain. This study was done to better understand how the climate has changed in the various regions of Puerto Rico, which will support a transition to renewable energy and aid in projecting future climate impacts on the residents and land they live on. This study used multiple datasets including Daymet, a daily surface weather and climatological summary; European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5); and station observations. The climate variables examined here include minimum and maximum near-surface air temperature, daily precipitation, incident shortwave radiation, and large-scale upper atmosphere conditions. The time periods cover the years 1950–2019. Results in both Daymet and station observations show long-term warming trends throughout every season. In particular, the minimum winter temperatures increased the most compared with maximum and mean air temperature. Northeast Puerto Rico, where the majority of the island’s population is located, experienced the largest warming. Summertime precipitation has decreased over time, and spring and fall precipitation has increased. The largest increase in fall is because there has been more available precipitable water, lower sea-level pressure, and favorable cyclonic circulations from the trade winds over the past seven decades. Overall, wet days are getting wetter, while dry days are getting drier. While there are good solar energy potential, we see a long-term decreasing trend in shortwave radiation over time in all seasons across entire island, possibly due to the increase of cloud cover.

54 ENVIRONMENTAL SCIENCES↗