Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “land surface model”

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,153 records · Page 64

Soil Moisture‐Cloud‐Precipitation Feedback in the Lower Atmosphere From Functional Decomposition of Satellite Observations

Abstract The feedback of topsoil moisture (SM) content on convective clouds and precipitation is not well understood and represented in the current generation of weather and climate models. Here, we use functional decomposition of satellite‐derived SM and cloud vertical profiles (CVP) to quantify the relationship between SM and the vertical distribution of cloud water in the central US. High‐dimensional model representation is used to disentangle the contributions of SM and other land‐surface and atmospheric variables to the CVP. Results show that the sign and strength of the SM‐cloud‐precipitation feedback varies with cloud height and time lag and displays a large spatial variability. Positive anomalies in antecedent 7‐hr SM and land‐surface temperature enhance cloud reflectivity up to 4 dBZ in the lower atmosphere about 1–3 km above the surface. Our approach presents new insights into the SM‐cloud‐precipitation feedback and aids in the diagnosis of land‐atmosphere interactions simulated by weather and climate models.

54 ENVIRONMENTAL SCIENCES↗

Estimation of Mesoscale Atmospheric Latent Heating Profiles from TRMM Rain Statistics Utilizing a Simple One-Dimensional Model

In this study, a model is developed to estimate mesoscale-resolution atmospheric latent heating (ALH) profiles. It utilizes rain statistics deduced from Tropical Rainfall Measuring Mission (TRMM) data, and cloud vertical velocity profiles and regional surface thermodynamic climatologies derived from other available data sources. From several rain events observed over tropical ocean and land, ALH profiles retrieved by this model in convective rain regions reveal strong warming throughout most of the troposphere, while in stratiform rain regions they usually show slight cooling below the freezing level and significant warming above. The mesoscale-average, or total, ALH profiles reveal a dominant stratiform character, because stratiform rain areas are usually much larger than convective rain areas. Sensitivity tests of the model show that total ALH at a given tropospheric level varies by less than +/- 10 % when convective and stratiform rain rates and mesoscale fractional rain areas are perturbed individually by +/- 15 %. This is also found when the non-uniform convective vertical velocity profiles are replaced by one that is uniform. Larger variability of the total ALH profiles arises when climatological ocean- and land-surface temperatures (water vapor mixing ratios) are independently perturbed by +/- 1.0 K (+/- 5%) and +/- 5.0 K (+/- 15%), respectively. At a given tropospheric level, such perturbations can cause a +/- 25% variation of total ALH over ocean, and a factor-of-two sensitivity over land. This sensitivity is reduced substantially if perturbations of surface thermodynamic variables do not change surface relative humidity, or are not extended throughout the entire model evaporation layer. The ALH profiles retrieved in this study agree qualitatively with tropical total diabatic heating profiles deduced in earlier studies. Also, from January and July 1999 ALH-profile climatologies generated separately with TRMM Microwave Imager and Precipitation Radar rain statistics, it is shown that ALH profiles can be retrieved utilizing diverse satellite-derived rain products that offer convective and stratiform discrimination. Therefore, the ALH retrieval model developed in this study can be used to make regional estimates of total diabatic heating profiles in the future Global Precipitation Measurement mission, and to assimilate these profiles into numerical weather forecast and climate models.

Iacovazzi, Robert A., Jr.↗

Estimation of Mesoscale Atmospheric Latent Heating Profiles from TRMM Rain Statistics Utilizing a Simple One-Dimensional Model

In this study, a model is developed to estimate mesoscale-resolution atmospheric latent heating (ALH) profiles. It utilizes rain statistics deduced from Tropical Rainfall Measuring Mission (TRMM) data, and cloud vertical velocity profiles and regional surface thermodynamic climatologies derived from other available data sources. From several rain events observed over tropical ocean and land, ALH profiles retrieved by this model in convective rain regions reveal strong warming throughout most of the troposphere, while in stratiform rain regions they usually show slight cooling below the freezing level and significant warming above. The mesoscale-average, or total, ALH profiles reveal a dominant stratiform character, because stratiform rain areas are usually much larger than convective rain areas. Sensitivity tests of the model show that total ALH at a given tropospheric level varies by less than +/- 10 % when convective and stratiform rain rates and mesoscale fractional rain areas are perturbed individually by 1 15 %. This is also found when the non-uniform convective vertical velocity profiles are replaced by one that is uniform. Larger variability of the total ALH profiles arises when climatological ocean- and land-surface temperatures (water vapor mixing ratios) are independently perturbed by +/- 1.0 K (+/- 5 %) and +/- 5.0 K (+/- 15 %), respectively. At a given tropospheric level, such perturbations can cause a +/- 25 % variation of total ALH over ocean, and a factor-of-two sensitivity over land. This sensitivity is reduced substantially if perturbations of surface thermodynamic variables do not change surface relative humidity, or are not extended throughout the entire model evaporation layer. The ALH profiles retrieved in this study agree qualitatively with tropical total diabatic heating profiles deduced in earlier studies. Also, from January and July 1999 ALH-profile climatologies generated separately with TRMM Microwave Imager and Precipitation Radar rain statistics, it is shown that ALH profiles can be retrieved utilizing diverse satellite-derived rain products that offer convective and stratiform discrimination. Therefore, the ALH retrieval model developed in this study can be used to make regional estimates of total diabatic heating profiles in the future Global Precipitation Measurement mission, and to assimilate these profiles into numerical weather forecast and climate models.

Iacovazzi, Robert A., Jr.↗

A Numerical Testbed for Remote Sensing of Aerosols, and its Demonstration for Evaluating Retrieval Synergy from a Geostationary Satellite Constellation of GEO-CAPE and GOES-R

We present a numerical testbed for remote sensing of aerosols, together with a demonstration for evaluating retrieval synergy from a geostationary satellite constellation. The testbed combines inverse (optimal-estimation) software with a forward model containing linearized code for computing particle scattering (for both spherical and non-spherical particles), a kernel-based (land and ocean) surface bi-directional reflectance facility, and a linearized radiative transfer model for polarized radiance. Calculation of gas absorption spectra uses the HITRAN (HIgh-resolution TRANsmission molecular absorption) database of spectroscopic line parameters and other trace species cross-sections. The outputs of the testbed include not only the Stokes 4-vector elements and their sensitivities (Jacobians) with respect to the aerosol single scattering and physical parameters (such as size and shape parameters, refractive index, and plume height), but also DFS (Degree of Freedom for Signal) values for retrieval of these parameters. This testbed can be used as a tool to provide an objective assessment of aerosol information content that can be retrieved for any constellation of (planned or real) satellite sensors and for any combination of algorithm design factors (in terms of wavelengths, viewing angles, radiance and/or polarization to be measured or used). We summarize the components of the testbed, including the derivation and validation of analytical formulae for Jacobian calculations. Benchmark calculations from the forward model are documented. In the context of NASA's Decadal Survey Mission GEOCAPE (GEOstationary Coastal and Air Pollution Events), we demonstrate the use of the testbed to conduct a feasibility study of using polarization measurements in and around the O2 A band for the retrieval of aerosol height information from space, as well as an to assess potential improvement in the retrieval of aerosol fine and coarse mode aerosol optical depth (AOD) through the synergic use of two future geostationary satellites, GOES-R (Geostationary Operational Environmental Satellite R-series) and TEMPO (Tropospheric Emissions: Monitoring of Pollution). Strong synergy between GEOS-R and TEMPO are found especially in their characterization of surface bi-directional reflectance, and thereby, can potentially improve the AOD retrieval to the accuracy required by GEO-CAPE.

Aerosols↗

Upstream surface roughness and terrain are strong drivers of contrast in tornado potential between North and South America

Central North America is the global hotspot for tornadoes, fueled by elevated terrain of the Rockies to the west and a source of warm, moist air from equatorward oceans. This conventional wisdom argues that central South America, with the Andes to the west and Amazon basin to the north, should have a “tornado alley” at least as active as central North America. Central South America has frequent severe thunderstorms yet relatively few tornadoes. Here, we show that conventional wisdom is missing an important ingredient specific to tornadoes: a smooth, flat ocean-like upstream surface. Using global climate model experiments, we show that central South American tornado potential substantially increases if its equatorward land surface is smoothed and flattened to be ocean-like. Similarly, we show that central North American tornado potential substantially decreases if its equatorward ocean surface is roughened to values comparable to forested land. A rough upstream surface suppresses the formation of tornadic environments principally by weakening the poleward low-level winds, characterized by a weakened low-level jet east of the mountain range. Results are shown to be robust for any midlatitude landmass using idealized experiments with a simplified continent and mountain range. Our findings indicate that large-scale upstream surface roughness is likely a first-order driver of the strong contrast in tornado potential between North and South America.

54 ENVIRONMENTAL SCIENCES↗

Asian Monsoons: Variability, Predictability, and Sensitivity to External Forcing

In this study, we have addressed the interannual variations of Asian monsoons including both broad-scale and regional monsoon components. Particular attention is devoted to the identities of the South China Sea monsoon and Indian monsoon. We use CPC Merged Analysis of Precipitation and NCEP reanalyses to define regional monsoon indices and to depict the various monsoons. Parallel modeling studies have also been carried out to assess the role of boundary forcing and the potential predictability of the monsoons. Each monsoon is characterized by its unique features. While the South Asian monsoon represents a classical monsoon in which anomalous circulation is governed by Rossby-wave dynamics, the Southeast Asian monsoon symbolizes a "hybrid" monsoon that features multi-cellular meridional circulation over eastern Asia. The broad-scale Asian monsoon links to the basin-wide atmospheric circulation over the Indian-Pacific oceans. Both SST and land surface processes are important for determining the variations of all monsoons. For the broad-scale monsoon, SST anomalies are more important than land surface processes. For regional monsoons, however, land surface processes may become equally important. Both observation and model shows that the broad-scale monsoon is potentially more predictable than regional monsoons, and that the Southeast Asian monsoon may possess higher predictability than the South Asian monsoon.

Yang, Song↗

Surface Flow Visualization of the High-Lift Common Research Model

A 10% scale version of the High-Lift Common Research Model (CRM-HL) was tested in the NASA Langley 14- by 22-Foot Subsonic Tunnel (14x22) in support of the NASA Advanced Air Transport Technology (AATT) Project. The CRM-HL experiment included various configurations such as conventional and simple-hinged flaps, with and without engine nacelle/pylon, with and without nacelle chine, different Active Flow Control (AFC) methods (sweeping jets, alternating pulsed jets, and preconditioned boundary layer blowing), and their various parameters. This particular study is focused on the surface flow visualization of the conventional CRM-HL model at landing configuration. The conventional CRM-HL model with the single-slotted Fowler flap system serves as a baseline for the AFC-enabled simplified high-lift configuration as well as a high-lift technology development platform due to its publicly open geometry. Surface flow visualizations were performed using fluorescent minitufts, which were found to be nonintrusive to the aerodynamic performance. Tuft flow visualizations are supplemented with the relevant pressure and force measurements in order to understand the flow characteristics developed on the conventional CRM- HL model. In addition, three dimensional, unsteady, compressible Computational Fluid Dynamic (CFD) simulations were performed for selective cases. The surface streamlines and transverse velocity fluctuations obtained by the CFD simulations are qualitatively compared to the tuft direction and tuft unsteadiness, respectively. Force measurements of the CRM-HL model show performance degradation at higher angles of attack. Surface flow visualizations revealed the performance loss due to the nacelle/pylon wake that grows with angle of attack and eventually promotes flow separation over the inboard wing. This performance loss was successfully recovered by placing a chine on the engine nacelle.

Koklu, Mehti↗

Asian Monsoons: Variability, Predictability, and Sensitivity to External Forcing

In this study, we have addressed the interannual variations of Asian monsoons including both broad-scale and regional monsoon components. Particular attention is devoted to the identities of the South China Sea monsoon and Indian monsoon. We use CPC Merged Analysis of Precipitation and NCEP reanalyses to define regional monsoon indices and to depict the various monsoons. Parallel modeling studies have also been carried out to assess the potential predictability of the broad-scale and regional monsoons. Each monsoon is characterized by its unique features. While the South Asian monsoon represents a classical monsoon in which anomalous circulation is governed by Rossby-wave dynamics, the Southeast Asian monsoon symbolizes a "hybrid" monsoon that features multi-cellular meridional circulation over eastern Asia. The broad-scale Asian monsoon links to the basin-wide atmospheric circulation over the Indian-Pacific oceans. Both Sea Surface Temperatures (SST) and land surface processes are important for determining the variations of all monsoons. For the broad-scale monsoon, SST anomalies are more important than land surface processes. However, for regional monsoons, land surface processes may become equally important. Both observation and model shows that the broad-scale monsoon is potentially more predictable than regional monsoons, and that the Southeast Asian monsoon may possess higher predictability than the South Asian monsoon.

Yang, Song↗

Improving the Estimation of the Atmospheric Water Vapor Pressure Using Interpretable Long Short-Term Memory Networks: Dataset, Python code, and trained models

Atmospheric water vapor pressure is an essential meteorological control on land surface and hydrologic processes. It is not as frequently observed as other meteorologic conditions, but often inferred through the August–Roche–Magnus formula by simply assuming dew point and daily minimum temperatures are equivalent or by empirically correlating the two temperatures using an aridity correction. The performance of both methods varies considerably across different regions and during different time periods; obtaining consistently accurate estimates across space and time remains a great challenge. We applied an interpretable Long Short-Term Memory (iLSTM) network conditioned on static, location specific attributes to estimate daily vapor pressure for 83 FLUXNET sites in the United States and Canada. This data package includes all raw data of the 83 FLUXNET sites, input data for model training/validation/test, trained models and results, and python codes for the manuscript "Improving the Estimation of the Atmospheric Water Vapor Pressure Using an Interpretable Long Short-term Memory Network". Specifically, it consists of five parts. - First, "1_Daymet_data_83sites.zip" includes raw data downloaded from Daymet for the 83 sites used in the paper according to their longitude and latitude, in which vapor pressure is used. It also includes a pre-processed CSV data file combining all data from the 83 sites which is specifically used for the paper. - Second, "2_Fluxnet2015_data_83sites.zip" includes raw half hourly data of the 83 sites downloaded from FLUXNET2015 data portal, pre-processed daily data of the 83 sites, a CSV file including combined pre-processed daily data of the 83 sites, and a CSV file including the information (site ID, site name, latitude, longitude, data available period) of the 83 sites. - Third, "3_MODIS_LAI_data_83sites_raw.zip" includes raw leaf area index (LAI) data downloaded from the AppEEARs data portal. - Fourth, "4_Scripts.zip" includes all scripts related to model training and post-processing of a trained model, and a jupyter notebook showing an example for model post-processing. Two typo errors in files titled "run2get_args.py" and "postprocess.py" were corrected on March 27, 2024 to avoid confusions. - Finally, "Trained_models_and_results.zip" includes three folders and three files with suffix ".npy", and each folder corresponds to one file with suffix ".npy" with the same title. Each of the three folders include all trained models associated with one iLSTM model configuration (35 models for each configuration, details are described in the paper). Each file with suffix ".npy" includes the post-processed results of the corresponding 35 models under one iLSTM model configuration.

54 ENVIRONMENTAL SCIENCES↗

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗

Constraining effects of aerosol-cloud interaction by accounting for coupling between cloud and land surface

Aerosol-cloud interactions (ACIs) are vital for regulating Earth’s climate by influencing energy and water cycles. Yet, effects of ACI bear large uncertainties, evidenced by systematic discrepancies between observed and modeled estimates. This study quantifies a major bias in ACI determinations, stemming from conventional surface or space measurements that fail to capture aerosol at the cloud level unless the cloud is coupled with land surface. We introduce an advanced approach to determine radiative forcing of ACI by accounting for cloud-surface coupling. By integrating field observations, satellite data, and model simulations, this approach reveals a drastic alteration in aerosol vertical transport and ACI effects caused by cloud coupling. In coupled regimes, aerosols enhance cloud droplet number concentration across the boundary layer more homogeneously than in decoupled conditions, under which aerosols from the free atmosphere predominantly affect cloud properties, leading to marked cooling effects. Our findings spotlight cloud-surface coupling as a key factor for ACI quantification, hinting at potential underassessments in traditional estimates.

54 ENVIRONMENTAL SCIENCES↗

Hyper‐Local Temperature Prediction Using Detailed Urban Climate Informatics

The accurate modeling of urban microclimate is a challenging task given the high surface heterogeneity of urban land cover and the vertical structure of street morphology. Recent years have witnessed significant efforts in numerical modeling and data collection of the urban environment. Nonetheless, it is difficult for the physical‐based models to fully utilize the high‐resolution data under the constraints of computing resources. The advancement in machine learning (ML) techniques offers the computational strength to handle the massive volume of data. In this study, we proposed a modeling framework that uses ML approach to estimate point‐scale street‐level air temperature from the urban‐resolving meso‐scale climate model and a suite of hyper‐resolution urban geospatial data sets, including three‐dimensional urban morphology, parcel‐level land use inventory, and weather observations from a sensor network. We implemented this approach in the City of Chicago as a case study to demonstrate the capability of the framework. The proposed approach vastly improves the resolution of temperature predictions in cities, which will help the city with walkability, drivability, and heat‐related behavioral studies. Moreover, we tested the model's reliability on out‐of‐sample locations to investigate the modeling uncertainties and the application potentials to the other areas. This study aims to gain insights into next‐gen urban climate modeling and guide the observation efforts in cities to build the strength for the holistic understanding of urban microclimate dynamics.

54 ENVIRONMENTAL SCIENCES↗

Modelling planetary bidirectional reflectance over land

Simulations of the bidirectional reflectance distribution functions at the top of the atmosphere (TOA) over land surfaces were performed and compared with observations derived from Nimbus-7. The two methods produce similar patterns even though the model computations were performed under the assumptions of a Lambertian surface. These results raise a question regarding the need to account for surface anisotropy in simulations intended for the interpretation of low-resolution satellite observations.

Pinker, R. T.↗

Estimating Evaporative Fraction From Readily Obtainable Variables in Mangrove Forests of the Everglades, U.S.A.

A remote-sensing-based model to estimate evaporative fraction (EF) the ratio of latent heat (LE; energy equivalent of evapotranspiration -ET-) to total available energy from easily obtainable remotely-sensed and meteorological parameters is presented. This research specifically addresses the shortcomings of existing ET retrieval methods such as calibration requirements of extensive accurate in situ micro-meteorological and flux tower observations, or of a large set of coarse-resolution or model-derived input datasets. The trapezoid model is capable of generating spatially varying EF maps from standard products such as land surface temperature [T(sub s)] normalized difference vegetation index (NDVI)and daily maximum air temperature [T(sub a)]. The 2009 model results were validated at an eddy-covariance tower (Fluxnet ID: US-Skr) in the Everglades using T(sub s) and NDVI products from Landsat as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors. Results indicate that the model accuracy is within the range of instrument uncertainty, and is dependent on the spatial resolution and selection of end-members (i.e. wet/dry edge). The most accurate results were achieved with the T(sub s) from Landsat relative to the T(sub s) from the MODIS flown on the Terra and Aqua platforms due to the fine spatial resolution of Landsat (30 m). The bias, mean absolute percentage error and root mean square percentage error were as low as 2.9% (3.0%), 9.8% (13.3%), and 12.1% (16.1%) for Landsat-based (MODIS-based) EF estimates, respectively. Overall, this methodology shows promise for bridging the gap between temporally limited ET estimates at Landsat scales and more complex and difficult to constrain global ET remote-sensing models.

trapezoid↗

A method for coupling a parameterization of the planetary boundary layer with a hydrologic model

Deardorff's parameterization of the planetary boundary layer is adapted to drive a hydrologic model. The method converts the atmospheric conditions measured at the anemometer height at one site to the mean values in the planetary boundary layer; it then uses the planetary boundary layer parameterization and the hydrologic variables to calculate the fluxes of momentum, heat and moisture at the atmosphere-land interface for a different site. A simplified hydrologic model is used for a simulation study of soil moisture and ground temperature on three different land surface covers. The results indicate that this method can be used to drive a spatially distributed hydrologic model by using observed data available at a meteorological station located on or nearby the site.

Lin, J. D.↗

Sarasota Climate: Monitoring Heat and Assessing Heat Vulnerability to Identify Locations for Heat Mitigation Efforts in Sarasota, Florida

The coastal county of Sarasota, Florida, is located within the humid subtropical climate region and experiences an average of 250 days of sunshine every year. The county has large, urbanized communities which are vulnerable to urban heat island (UHI) effects. These rapidly growing communities contribute to the increasing surface temperatures and placing more residents at risk of heat-related illness. In partnership with Sarasota County Sustainability, the Sarasota Climate team utilized Earth observation data from NASA Landsat 8 and Landsat 9’s Thermal Infrared Sensors (TIRS), and the International Space Station’s Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) to model UHI effects within the county during the Summer for the last five years, 2019 to 2023. Data analysis with the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) and Urban Heat Exposure Assessment Tempe 1.0 (UHEAT 1.0) models and within ModelBuilder produced maps that identified the land surface temperature (LST) variance within the county, the regions that are most susceptible to extreme heat, and areas least capable of mitigating the effects of UHI. The results revealed that heat intensity varies significantly across Sarasota County with the highest temperatures in the more developed western part of the county. Additionally, the team identified that there are at least three vulnerable communities that exist in high-heat regions, including North Sarasota, Venice, and North Port. These regions have an overlap between socioeconomic sensitivity and environmental hazard that indicate a high priority in future heat mitigation efforts.

Remote sensing↗

Desert Peak Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to the Desert Peak Geothermal Field. It includes all input and output files used in the project. The files include data categories of raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs including Radar, SWIR, Thermal, Geophysics, Geology, and Wells. The files for the Desert Peak Geothermal Site are used with the Geothermal Exploration Artificial Intelligence to identify indicators of blind geothermal systems. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Desert Peak Geothermal Field.

15 GEOTHERMAL ENERGY↗

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗