Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLIMATE MODELS”

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 1,441 records · Page 80

Projected U.S. drought extremes through the twenty-first century with vapor pressure deficit

Global warming is expected to enhance drought extremes in the United States throughout the twenty-first century. Projecting these changes can be complex in regions with large variability in atmospheric and soil moisture on small spatial scales. Vapor Pressure Deficit (VPD) is a valuable measure of evaporative demand as moisture moves from the surface into the atmosphere and a dynamic measure of drought. Here, VPD is used to identify short-term drought with the Standardized VPD Drought Index (SVDI); and used to characterize future extreme droughts using grid dependent stationary and non-stationary generalized extreme value (GEV) models, and a random sampling technique is developed to quantify multimodel uncertainties. The GEV analysis was performed with projections using the Weather Research and Forecasting model, downscaled from three Global Climate Models based on the Representative Concentration Pathway 8.5 for present, mid-century and late-century. Results show the VPD based index (SVDI) accurately identifies the timing and magnitude short-term droughts, and extreme VPD is increasing across the United States and by the end of the twenty-first century. The number of days VPD is above 9 kPa increases by 10 days along California’s coastline, 30–40 days in the northwest and Midwest, and 100 days in California’s Central Valley.

54 ENVIRONMENTAL SCIENCES↗

Exploring the role of clouds in offshore wind potential off the US West Coast in a changing climate

To meet US goals of deploying additional wind energy as part of the decarbonization strategy, wind plants are being planned for the deep water offshore the western US. The wind flow in that region is complex due to the proximity to the coast, cold water upwelling, and persistent stratiform clouds that interact with radiation in ways that have the potential to destabilize the atmosphere. That flow has the potential to change with a changing climate. To address these issues, we assess the flow and the clouds in that region using downscaled climate model data, under both historic climate (1975–2005) and projected future (2025–2055) conditions. We note that the climate simulations agree fairly well with the cloud patterns observed by satellite data in the nearshore and offshore regions. We then assess the projected changes in clouds, wind speed, and other important variables, noting that our simulations project that the predominant north/northwesterly low-level jet is expected to strengthen and clouds are likely to be commensurately enhanced, although projected changes are within about 10% of current conditions. Our examination of the dynamics associated with the changes in the climate simulations provides confidence in the dynamical consistency of these projected changes.

17 WIND ENERGY↗

TChem-atm v1.0

SAND2024-11300O TChem-atm is a software library that was developed to solve complex kinetic models for atmospheric chemistry applications. TChem-atm interface employs a hierarchical parallelism design to exploit the massive parallelism available from modern computing platforms. It also supports gas atmospheric chemistry applications, e.g., the energy exascale earth system model. TChem can be used as a box model or coupled with a climate model to compute the time evolution of gas tracer species. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Safta, Cosmin↗

Investigation of models for large-scale meteorological prediction experiments

The feasibility of long-range weather prediction through the use of global general circulation models (GCMs) was investigated. A climate model was developed to simulate the monthly mean state of the atmosphere from real global initial data at the beginning of the month. The model contains the same dynamic and physical ingredients as most numerical weather prediction models and GCMs. The model generates a one-day global simulation on the 8 x 10 grid in four minutes (on an IBM 360/95 computer), so that a 30 day forecast can be executed in two hours. The high speed of the model is achieved mainly at the price of its coarse resolution, which requires certain parameterizations of surface boundary conditions, as well as inherent filtering of smaller-scale features of the initial state.

Spar, J.↗

Climate change and the middle atmosphere. I - The doubled CO2 climate

The effect of doubling the atmospheric content of CO2 on the middle-atmosphere climate is investigated using the GISS global climate model. In the standard experiment, the CO2 concentration is doubled both in the stratosphere and troposphere, and the SSTs are increased to match those of the doubled CO2 run of the GISS model. Results show that the doubling of CO2 leads to higher temperatures in the troposphere, and lower temperatures in the stratosphere, with a net result being a decrease of static stability for the atmosphere as a whole. The middle atmosphere dynamical differences found were on the order of 10-20 percent of the model values for the current climate. These differences, along with the calculated temperature differences of up to about 10 C, may have a significant impact on the chemistry of the future atmosphere, including that of stratospheric ozone, the polar ozone 'hole', and basic atmospheric composition.

Rind, D.↗

A Stochastic Model of Space-Time Variability of Tropical Rainfall: I. Statistics of Spatial Averages

Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Kundu, Prasun K.↗

Top-of-Atmosphere Albedo Estimation from Angular Distribution Models using Scene Identification from Satellite Cloud Property Retrievals

The next generation of Earth radiation budget satellite instruments will routinely merge estimates of global top-of-atmosphere radiative fluxes with cloud properties. This information will offer many new opportunities for validating radiative transfer models and cloud parameterizations in climate models. In this study, five months of POLarization and Directionality of the Earth's Reflectances (POLDER) 670 nm radiance measurements are considered in order to examine how satellite cloud property retrievals can be used to define empirical Angular Distribution Models (ADMs) for estimating top-of-atmosphere (TOA) albedo. ADMs are defined for 19 scene types defined by satellite retrievals of cloud fraction and cloud optical depth. Two approaches are used to define the ADM scene types: The first assumes there are no biases in the retrieved cloud properties and defines ADMs for fixed discrete intervals of cloud fraction and cloud optical depth (fixed-tau approach). The second approach involves the same cloud fraction intervals, but uses percentile intervals of cloud optical depth instead (percentile-tau approach). Albedos generated using these methods are compared with albedos inferred directly from the mean observed reflectance field. Albedos based on ADMs that assume cloud properties are unbiased (fixed-tau approach) show a strong systematic dependence on viewing geometry. This dependence becomes more pronounced with increasing solar zenith angle, reaching approximately equals 12% (relative) between near-nadir and oblique viewing zenith angles for solar zenith angles between 60 deg and 70 deg. The cause for this bias is shown to be due to biases in the cloud optical depth retrievals. In contrast, albedos based on ADMs built using percentile intervals of cloud optical depth (percentile-tau approach) show very little viewing zenith angle dependence and are in good agreement with albedos obtained by direct integration of the mean observed reflectance field (less than 1% relative error). When the ADMs are applied separately to populations consisting of only liquid water and ice clouds, significant biases in albedo with viewing geometry are observed (particularly at low sun elevations), highlighting the need to account for cloud phase both in cloud optical depth retrievals and in defining ADM scene types. ADM-derived monthly mean albedos determined for all 5 deg x 5 deg latitude/longitude regions over ocean are in good agreement (regional RMS relative errors less than 2%) with those obtained by direct integration when ADM albedos inferred from specific angular bins are averaged together. Albedos inferred from near-nadir and oblique viewing zenith angles are the least accurate, with regional RMS errors reaching approximately 5-10% (relative). Compared to an earlier study involving ERBE ADMs, regional mean albedos based on the 19 scene types considered here show a factor of 4 reduction in bias error and a factor of 3 reduction in RMS error.

Loeb, N. G.↗

Receivers Gather Data for Climate, Weather Prediction

Signals from global positioning system (GPS) satellites are now being used for more than just location and navigation information. By looking at the radio waves from GPS satellites, a technology developed at NASA s Jet Propulsion Laboratory (JPL) not only precisely calculates its position, but can also use a technique known as radio occultation to help scientists study the Earth s atmosphere and gravity field to improve weather forecasts, monitor climate change, and enhance space weather research. The University Corporation for Atmospheric Research (UCAR), a nonprofit group of universities in Boulder, Colorado, compares radio occultation to the appearance of a pencil when viewed though a glass of water. The water molecules change the path of visible light waves so that the pencil appears bent, just like molecules in the air bend GPS radio signals as they pass through (or are occulted by) the atmosphere. Through measurements of the amount of bending in the signals, scientists can construct detailed images of the ionosphere (the energetic upper part of the atmosphere) and also gather information about atmospheric density, pressure, temperature, and moisture. Once collected, this data can be input into weather forecasting and climate models for weather prediction and climate studies. Traditionally, such information is obtained through the use of weather balloons. In 1998, JPL started developing a new class of GPS space science receivers, called Black Jack, that could take precise measurements of how GPS signals are distorted or delayed along their way to the receiver. By 2006, the first demonstration of a GPS radio occultation constellation was launched through a collaboration among Taiwan s National Science Council and National Space Organization, the U.S. National Science Foundation, NASA, the National Oceanic and Atmospheric Administration (NOAA), and other Federal entities. Called the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC), JPL was responsible for designing COSMIC s primary instrument - based on its revolutionary Black Jack receiver.

Source record↗

Climate Change Detection and Attribution of Infrared Spectrum Measurements

Climate change occurs when the Earth's energy budget changes due to natural or possibly anthropogenic forcings. These forcings cause the climate system to adjust resulting in a new climate state that is warmer or cooler than the original. The key question is how to detect and attribute climate change. The inference of infrared spectral signatures of climate change has been discussed in the literature for nearly 30 years. Pioneering work in the 1980s noted that distinct spectral signatures would be evident in changes in the infrared radiance emitted by the Earth and its atmosphere, and that these could be observed from orbiting satellites. Since then, a number of other studies have advanced the concepts of spectral signatures of climate change. Today the concept of using spectral signatures to identify and attribute atmospheric composition change is firmly accepted and is the foundation of the Climate Absolute Radiance and Refractivity Observatory (CLARREO) satellite mission being developed at NASA. In this work, we will present an overview of the current climate change detection concept using climate model calculations as surrogates for climate change. Any future research work improving the methodology to achieve this concept will be valuable to our society.

Phojanamongkolkij, Nipa↗

Stratospheric Ozone Response in Experiments G3 and G4 of the Geoengineering Model Intercomparison Project (GeoMIP)

Geoengineering with stratospheric sulfate aerosols has been proposed as a means of temporarily cooling the planet, alleviating some of the side effects of anthropogenic CO2 emissions. However, one of the known side effects of stratospheric injections of sulfate aerosols is a decrease in stratospheric ozone. Here we show results from two general circulation models and two coupled chemistry climate models that have simulated stratospheric sulfate aerosol geoengineering as part of the Geoengineering Model Intercomparison Project (GeoMIP). Changes in photolysis rates and upwelling of ozone-poor air in the tropics reduce stratospheric ozone, suppression of the NOx cycle increases stratospheric ozone, and an increase in available surfaces for heterogeneous chemistry modulates reductions in ozone. On average, the models show a factor 20-40 increase of the sulfate aerosol surface area density (SAD) at 50 hPa in the tropics with respect to unperturbed background conditions and a factor 3-10 increase at mid-high latitudes. The net effect for a tropical injection rate of 5 Tg SO2 per year is a decrease in globally averaged ozone by 1.1-2.1 DU in the years 2040-2050 for three models which include heterogeneous chemistry on the sulfate aerosol surfaces. GISS-E2-R, a fully coupled general circulation model, performed simulations with no heterogeneous chemistry and a smaller aerosol size; it showed a decrease in ozone by 9.7 DU. After the year 2050, suppression of the NOx cycle becomes more important than destruction of ozone by ClOx, causing an increase in total stratospheric ozone. Contribution of ozone changes in this experiment to radiative forcing is 0.23 W m-2 in GISS-E2-R and less than 0.1 W m-2 in the other three models. Polar ozone depletion, due to enhanced formation of both sulfate aerosol SAD and polar stratospheric clouds, results in an average 5 percent increase in calculated surface UV-B.

aerosol↗

Remote Sensing of Aerosol in the Terrestrial Atmosphere from Space: New Missions

The distribution and properties of atmospheric aerosols on a global scale are not well known in terms of determination of their effects on climate. This mostly is due to extreme variability of aerosol concentrations, properties, sources, and types. Aerosol climate impact is comparable to the effect of greenhouse gases, but its influence is more difficult to measure, especially with respect to aerosol microphysical properties and the evaluation of anthropogenic aerosol effect. There are many satellite missions studying aerosol distribution in the terrestrial atmosphere, such as MISR/Terra, OMI/Aura, AVHHR, MODIS/Terra and Aqua, CALIOP/CALIPSO. To improve the quality of data and climate models, and to reduce aerosol climate forcing uncertainties, several new missions are planned. The gap in orbital instruments for studying aerosol microphysics has arisen after the Glory mission failed during launch in 2011. In this review paper, we describe several planned aerosol space missions, including the Ukrainian project Aerosol-UA that obtains data using a multi-channel scanning polarimeter and wide-angle polarimetric camera. The project is designed for remote sensing of the aerosol microphysics and cloud properties on a global scale.

Aerosol↗

Computational Modeling of Atmospheric Processes at Texas Southern University

Texas Southern University (TSU) is strengthening its research program in atmospheric chemistry and physics with a climate science emphasis by leveraging partnerships with the U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) Facility, Brookhaven National Laboratory (BNL), and the Tracking Aerosol Convection Interactions ExpeRiment (TRACER). This RDPP-supported program focuses on secondary organic aerosols (SOAs) and reactive atmospheric species that influence cloud formation, precipitation processes, and radiative forcing. SOAs play a critical role in cloud microphysics and Earth’s energy balance, yet the chemical and physical mechanisms governing SOA–cloud interactions remain a significant source of uncertainty in predictive climate models. Through computational modeling, observational data analysis, and national laboratory collaboration, this program develops a skilled cohort of students trained in atmospheric science, environmental data analysis, and climate-relevant modeling. These research experiences build technical competencies that are transferable to careers in government laboratories, academia, and industry. By engaging students from historically underrepresented communities in high-impact climate research, TSU expands participation in the atmospheric sciences workforce while contributing meaningful scientific insights to DOE-supported ARM research activities. This partnership strengthens national capacity in climate science and supports the development of the next generation of atmospheric researchers.

54 ENVIRONMENTAL SCIENCES↗

Representation of Extratropical Cyclones, Blocking Anticyclones, and Alpine Circulation Types in Multiple Reanalyses and Model Simulations

Atmospheric circulation types, blockings, and cyclones are prime features of the extratropical flow and key to understanding the climate system. This study intercompares the representation of these features in 10 reanalyses and in an ensemble of 30 climate model simulations between 1980 and 2005. Both modern, full-input reanalyses and century-long, surface-input reanalyses are examined. Modern full-input reanalyses agree well on key statistics of blockings, cyclones, and circulation types. Yet, the intensity and depth of cyclones vary among them. Reanalyses with higher horizontal resolution show higher cyclone center densities and more intense cyclones. For blockings, no strict relationship is found between frequency or intensity and horizontal resolution. Full-input reanalyses contain more intense blocking, compared to surface-input reanalyses. Circulation-type classifications over central Europe show that both versions of the Twentieth Century Reanalysis dataset contain more easterlies and fewer westerlies than any other reanalysis, owing to their high pressure bias over northeast Europe. The temporal correlation of annual circulation types over central Europe and blocking frequencies over the North Atlantic–European domain between reanalyses is high (around 0.8). The ensemble simulations capture the main characteristics of midlatitudinal atmospheric circulation. Circulation types of westerlies to northerlies over central Europe are overrepresented. There are too few blockings in the higher latitudes and an excess of cyclones in the midlatitudes. Other characteristics, such as blocking amplitude and cyclone intensity, are realistically represented, making the ensemble simulations a rich dataset to assess changes in climate variability.

54 ENVIRONMENTAL SCIENCES↗

Shorebird Migration Patterns in Response to Climate Change: A Modeling Approach

The availability of satellite remote sensing observations at multiple spatial and temporal scales, coupled with advances in climate modeling and information technologies offer new opportunities for the application of mechanistic models to predict how continental scale bird migration patterns may change in response to environmental change. In earlier studies, we explored the phenotypic plasticity of a migratory population of Pectoral sandpipers by simulating the movement patterns of an ensemble of 10,000 individual birds in response to changes in stopover locations as an indicator of the impacts of wetland loss and inter-annual variability on the fitness of migratory shorebirds. We used an individual based, biophysical migration model, driven by remotely sensed land surface data, climate data, and biological field data. Mean stop-over durations and stop-over frequency with latitude predicted from our model for nominal cases were consistent with results reported in the literature and available field data. In this study, we take advantage of new computing capabilities enabled by recent GP-GPU computing paradigms and commodity hardware (general purchase computing on graphics processing units). Several aspects of our individual based (agent modeling) approach lend themselves well to GP-GPU computing. We have been able to allocate compute-intensive tasks to the graphics processing units, and now simulate ensembles of 400,000 birds at varying spatial resolutions along the central North American flyway. We are incorporating additional, species specific, mechanistic processes to better reflect the processes underlying bird phenotypic plasticity responses to different climate change scenarios in the central U.S.

Smith, James A.↗

Why Is Improvement of Earth System Models so Elusive? Challenges and Strategies from Dust Aerosol Modeling

Past decades have seen an accelerating increase in computing efficiency, while climate models are representing a rapidly widening set of physical processes. Yet simulations of some fundamental aspects of climate like precipitation or aerosol forcing remain highly uncertain and resistant to progress. Dust aerosol modeling of soil particles lofted by wind erosion has seen a similar conflict between increasing model sophistication and remaining uncertainty. Dust aerosols perturb the energy and water cycles by scattering radiation and acting as ice nuclei, while mediating atmospheric chemistry and marine photosynthesis (and thus the carbon cycle). These effects take place across scales from the dimensions of an ice crystal to the planetary-scale circulation that disperses dust far downwind of its parent soil. Representing this range leads to several modeling challenges. Should we limit complexity in our model, which consumes computer resources and inhibits interpretation? How do we decide if a process involving dust is worthy of inclusion within our model? Can we identify a minimal representation of a complex process that is efficient yet retains the physics relevant to climate? Answering these questions about the appropriate degree of representation is guided by model evaluation, which presents several more challenges. How do we proceed if the available observations do not directly constrain our process of interest? (This could result from competing processes that influence the observed variable and obscure the signature of our process of interest.) Examples will be presented from dust modeling, with lessons that might be more broadly applicable. The end result will either be clinical depression or there assuring promise of continued gainful employment as the community confronts these challenges.

Climate models↗

The Effects of Surface Longwave Spectral Emissivity on Atmospheric Circulation and Convection over the Sahara and Sahel

This study quantifies the impact of the inclusion of realistic surface spectral emissivity in the Sahara and Sahel on the simulated local climate and beyond. The surface emissivity in these regions can be as low as 0.6–0.7 over the infrared window band while close to unity in other spectral bands, but such spectral dependence has been ignored in current climate models. Realistic surface spectral emissivities over the Sahara and Sahel are incorporated into the Community Earth System Model (CESM) version 1.1.1, while treatments of surface emissivity for the rest of the globe remain unchanged. Both the modified and standard CESM are then forced with prescribed climatological SSTs and fixed present-day forcings for 35-yr simulations. The outputs from the last 30 years are analyzed. Compared to the standard CESM, the modified CESM has warmer surface air temperature, as well as a warmer and wetter planetary boundary layer over the Sahara and Sahel. The modified CESM thus favors more convection in these regions and has more convective rainfall, especially in the Sahara. The moisture convergence induced by such inclusion of surface spectral emissivity also contributes to the differences in simulated precipitation in the Sahel and the region south to it. Compared to observations, inclusion of surface spectral emissivity reduces surface temperature biases in the Sahara and precipitation biases in the Gulf of Guinea but exacerbates the wet biases in the Sahara. Such realistic representation of surface spectral emissivity can help unmask other factors contributing to regional biases in the original CESM.

54 ENVIRONMENTAL SCIENCES↗

Assessing the cloud radiative bias at Macquarie Island in the ACCESS-AM2 model

Abstract. As a long-standing problem in climate models, large positive shortwave radiation biases exist at the surface over the Southern Ocean, impacting the accurate simulation of sea surface temperature, atmospheric circulation, and precipitation. Underestimations of low-level cloud fraction and liquid water content are suggested to predominantly contribute to these radiation biases. Most model evaluations for radiation focus on summer and rely on satellite products, which have their own limitations. In this work, we use surface-based observations at Macquarie Island to provide the first long-term, seasonal evaluation of both downwelling surface shortwave and longwave radiation in the Australian Community Climate and Earth System Simulator Atmosphere-only Model version 2 (ACCESS-AM2) over the Southern Ocean. The capacity of the Clouds and the Earth’s Radiant Energy System (CERES) product to simulate radiation is also investigated. We utilize the novel lidar simulator, the Automatic Lidar and Ceilometer Framework (ALCF), and all-sky cloud camera observations of cloud fraction to investigate how radiation biases are influenced by cloud properties. Overall, we find an overestimation of +9.5±33.5 W m−2 for downwelling surface shortwave radiation fluxes and an underestimation of -2.3±13.5 W m−2 for downwelling surface longwave radiation in ACCESS-AM2 in all-sky conditions, with more pronounced shortwave biases of +25.0±48.0 W m−2 occurring in summer. CERES presents an overestimation of +8.0±18.0 W m−2 for the shortwave and an underestimation of -12.1±12.2 W m−2 for the longwave in all-sky conditions. For the cloud radiative effect (CRE) biases, there is an overestimation of +4.8±28.0 W m−2 in ACCESS-AM2 and an underestimation of -7.9±20.9 W m−2 in CERES. An overestimation of downwelling surface shortwave radiation is associated with an underestimated cloud fraction and low-level cloud occurrence. We suggest that modeled cloud phase is also having an impact on the radiation biases. Our results show that the ACCESS-AM2 model and CERES product require further development to reduce these radiation biases not just in shortwave and in all-sky conditions, but also in longwave and in clear-sky conditions.

54 ENVIRONMENTAL SCIENCES↗

Mineralogical Signatures of Cold and Icy Climates on Ancient and Modern Mars

Liquid water was abundant on early Mars, but whether the climate was warm and wet or cold and icy with punctuated periods of melting is still poorly understood. Modern climate models for Mars tend to predict a colder, icier early climate than previously imagined [e.g., 1]. However, any model for the early climate on Mars must be reconciled with the chemical record. We currently do not understand how alteration mineralogy formed in snow and ice dominated conditions compares to that of warmer climates, and it is unclear whether cold climate weathering could form all or any of the aqueous alteration phases expressed on early martian surfaces [2]. To help resolve this knowledge gap, we synthesize results from glacial Mars analog sites at the Three Sisters, Oregon and mafic regions of the Antarctic ice sheet, and compare them to the surface mineralogy of Mars. These sites provide the opportunity to investigate weathering in environments analogous to glacial environments on Mars throughout geologic time, including snowpacks or smaller wet-based or polythermal glaciers [3, 4] as well as the proposed extensive ice sheets of the late Noachian icy highlands model

Scudder, N. A.↗