Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied meteorology”

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 523 records · Page 29

Multiple dipole eddies in the Alaska Coastal Current detected with Landsat thematic mapper data

Seventeen dipole eddies, including five large, well-formed ones, three second-generation eddies, and two double eddies, were observed in the Alaska Coastal Current near Kayak Island in one single scene of the Landsat thematic mapper (TM) on April 22, l985. The digital Landsat TM satellite data were computer analyzed to extract details in the near coastal circulation in the northern Gulf of Alaska. Enhancement techniques were applied to the visible and thermal IR bands. The features are evident only in the visible bands because of the ability of these bands to detect the distribution of sediments in the near surface. These eddies did not have a significant thermal signature.The sources of these sediments are the glacial streams found throughout the Gulf of Alaska coast. Eddies of this configuration and frequency have never been observed here previously. However, the oceanographic and meteorological conditions are typical for this time of year. These eddies should be important to the cross-shelf mixing processes in the Alaska Coastal Current and are an indicator that the flow here can be unstable at certain times of the year.

Ahlnas, Kristina↗

Ka-Band ARM Zenith Radar Corrections (KAZRCOR, KAZRCFRCOR) Value-Added Products

The Ka-Band Atmospheric Radiation Measurement (ARM) Zenith Radar Corrections (KAZRCOR) and KA-Band ARM Zenith Radar CF-Radial, Corrected (KAZRCFRCOR) value-added products (VAPs) perform several corrections to the ingested KAZR moments and also create a significant detection mask for each radar mode. The VAPs compute gaseous attenuation as a function of time and radial distance from the radar antenna, based on ambient meteorological observations, and correct observed reflectivities for that effect. Mean Doppler velocities are dealiased to correct velocities whose magnitudes exceed the radar’s Nyquist velocity. Input KAZR data fields are passed through to the KAZRCOR or KAZRCFRCOR output files, in their native time and range coordinates. Complementary corrected reflectivity and velocity fields are provided, along with a mask of significant detections and a number of data quality flags. This report covers the KAZRCOR VAP as applied to the original KAZR radars and the upgraded KAZR2 radars. Originally, prior to late 2019, there were two separate code bases for the different radar versions. Following the harmonization of KAZR and KAZR2 data formats in 2019, only a single code base is required. The new combined KAZR and KAZR2 code base is called the KAZRCFRCOR VAP. The ‘cfr portion of the VAP name refers to the Radial Climate and Forecasting data format standards that are used in these data sets. Throughout this report, references to ‘KAZRCOR’ should be taken to apply to the ‘KAZRCFRCOR’ VAP as well unless there is an explicit statement to the contrary.

54 ENVIRONMENTAL SCIENCES↗

Retrieval of sizes and total masses of particles in volcanic clouds using AVHRR bands 4 and 5

The advanced very high resolution radiometer (AVHRR) sensor on polar orbiting NOAA satellites can discriminate between volcanic clouds and meteorological clouds using two-band data in the thermal infrared. This paper is aimed at developing a retrieval of the particle sizes, optical depth, and total masses of particles from AVHRR two-band data of volcanic clouds. Radiative transfer calculations are used with a semi-transparent cloud model that is based on assumptions of spherical particle shape, a homogeneous underlying surface, and a simple thin cloud parallel to the surface. The model is applied to observed AVHRR data from a 13-hour old drifting cloud from the August 19, 1992, eruption of Crater Peak/Spurr Volcano, Alaska. The AVHRR data fit in the range of results calculated by the model, which supports its credibility. According to the model results, the average of effective particle radius in the test frame of this cloud is in the range of 2 to 2.5 micrometer, the optical depth at 12 micrometer is about 0.60 - 0.65. The total estimated mass of ash in the air amounts to 0.24 - 0.31 x 10(exp 6) tons, which is about 0.7-0.9% of the mass measured in the ashfall blanket. Sensitivity tests show that the mass estimate is more sensitive to the assumed ash size distribution than it is to the ash composition.

Wen, Shiming↗

A New Operational Snow Retrieval Algorithm Applied to Historical AMSR-E Brightness Temperatures

Snow is a key element of the water and energy cycles and the knowledge of spatio-temporal distribution of snow depth and snow water equivalent (SWE) is fundamental for hydrological and climatological applications. SWE and snow depth estimates can be obtained from spaceborne microwave brightness temperatures at global scale and high temporal resolution (daily). In this regard, the data recorded by the Advanced Microwave Scanning Radiometer-Earth Orbiting System (EOS) (AMSR-E) onboard the National Aeronautics and Space Administration's (NASA) AQUA spacecraft have been used to generate operational estimates of SWE and snow depth, complementing estimates generated with other microwave sensors flying on other platforms. In this study, we report the results concerning the development and assessment of a new operational algorithm applied to historical AMSR-E data. The new algorithm here proposed makes use of climatological data, electromagnetic modeling and artificial neural networks for estimating snow depth as well as a spatio-temporal dynamic density scheme to convert snow depth to SWE. The outputs of the new algorithm are compared with those of the current AMSR-E operational algorithm as well as in-situ measurements and other operational snow products, specifically the Canadian Meteorological Center (CMC) and GlobSnow datasets. Our results show that the AMSR-E algorithm here proposed generally performs better than the operational one and addresses some major issues identified in the spatial distribution of snow depth fields associated with the evolution of effective grain size.

passive microwave↗

Physically Interpretable Neural Networks for the Geosciences: Applications to Earth System Variability

Abstract Neural networks have become increasingly prevalent within the geosciences, although a common limitation of their usage has been a lack of methods to interpret what the networks learn and how they make decisions. As such, neural networks have often been used within the geosciences to most accurately identify a desired output given a set of inputs, with the interpretation of what the network learns used as a secondary metric to ensure the network is making the right decision for the right reason. Neural network interpretation techniques have become more advanced in recent years, however, and we therefore propose that the ultimate objective of using a neural network can also be the interpretation of what the network has learned rather than the output itself. We show that the interpretation of neural networks can enable the discovery of scientifically meaningful connections within geoscientific data. In particular, we use two methods for neural network interpretation called backward optimization and layerwise relevance propagation, both of which project the decision pathways of a network back onto the original input dimensions. To the best of our knowledge, LRP has not yet been applied to geoscientific research, and we believe it has great potential in this area. We show how these interpretation techniques can be used to reliably infer scientifically meaningful information from neural networks by applying them to common climate patterns. These results suggest that combining interpretable neural networks with novel scientific hypotheses will open the door to many new avenues in neural network‐related geoscience research.

58 GEOSCIENCES↗

Understanding the Transport of Patagonian Dust and Its Influence on Marine Biological Activity in the South Atlantic Ocean

Modeling and remote sensing techniques were applied to examine the horizontal and vertical transport pathways of Patagonian dust and quantify the effect of soluble-iron- laden mineral dust deposition on marine primary productivity in the South Atlantic Ocean (SAO) surface waters. The global chemistry transport model GEOS-Chem, implemented with an iron dissolution scheme, was applied to evaluate the atmospheric transport and deposition of mineral dust and bioavailable iron during two dust outbreaks originating in the source regions of Patagonia. In addition to this "rapidly released" iron, offline calculations were also carried out to estimate the amount of bioavailable iron leached during the residence time of dust in the ocean mixed layer. Model simulations showed that the horizontal and vertical transport pathways of Patagonian dust plumes were largely influenced by the synoptic meteorological patterns of high and low pressure systems. Model-predicted horizontal and vertical transport pathways of Patagonian dust over the SAO were in reasonable agreement with remotely-sensed data. Comparison between remotely-sensed and offline calculated ocean surface chlorophyll-a concentrations indicated that, for the two dust outbreaks examined in this study, the deposition of bioavailable iron in the SAO through atmospheric pathways was insignificant. As the two dust transport episodes examined here represent typical outflows of mineral dust from South American sources, our study suggests that the atmospheric deposition of mineral dust is unlikely to induce large scale marine primary productivity and carbon sequestration in the South Atlantic sector of the Southern Ocean.

Johnson, Matthew↗

The Kinematic and Microphysical Control of Storm Integrated Lightning Flash Extent

Objective: To investigate the kinematic and microphysical control of lightning properties, particularly those that may govern the production of nitrogen oxides (NOx) in thunderstorms, such as flash rate, type (intracloud [IC] vs. cloud-to-ground [CG] ) and extent. Data and Methodology: a) NASA MSFC Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection Network(TradeMark) (NLDN) observations following ordinary convective cells through their lifecycle. b) LNOM provides estimates of flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles (Koshak et al. 2012). c) LNOM lightning characteristics are compared to the evolution of updraft and precipitation properties inferred from dual-Doppler (DD) and polarimetric radar analyses of UAHuntsville Advanced Radar for Meteorological and Operational Research (ARMOR, Cband, polarimetric) and KHTX (S-band, Doppler).

Carey, Lawrence D.↗

Timing the Flames: Geostationary Satellite Detection of Diurnally Shifting Stubble Burning in Northwestern India

Post-monsoon open-field stubble burning in northwestern (NW) India—a key agricultural region known as the “breadbasket”—is a longstanding practice used to clear fields. Satellite observations spanning over two decades have revealed significant upward trends in crop production, vegetative greenness, and the frequency of post-harvest fires, with this last contributing to hazardous air quality during the peak burning season (mid-October to mid-November). Since 2022, thermal anomaly data from Aqua-MODIS and SNPP-VIIRS sensors have shown a sharp decline in reported fire events—an observation that contrasts starkly with the concurrent rise in regional aerosol loading detected from space. This apparent discrepancy became particularly pronounced in 2024–2025, prompting a closer examination using high-temporal-resolution imagery from the Advanced Meteorological Imager (AMI) on the geostationary satellite GEO-KOMPSAT-2A. These observations revealed a clear spike in fire-related signals occurring around and after 4:00 p.m. local time, i.e., outside the typical noon to 2:00 p.m. detection window of the MODIS and VIIRS. A fire detection algorithm exploiting the fire-sensitive shortwave-infrared 3.8 μm signal and its contrast to 11.2 μm infrared observations is designed to adopt AMI observations and applied to its multi-year observations (2019–2025). The resulting fire dataset unambiguously shows a gradual shift in stubble burning activity toward the late afternoon hours beginning in 2022 which is underreported by polar-orbiting satellites. The orbital drift of NASA’s MODIS sensor on the Aqua platform allows detection of some of the gradually shifting fires during afternoon hours, but the MODIS still misses a large number of fires occurring around and after 4 p.m. The AMI’s relatively coarse spatial resolution (~4 km), a consequence of its slant viewing geometry over NW India, imposes inherent limitations on quantifying the full extent of fire occurrences. The operational air quality forecasting models currently assimilate satellite fire detections predominantly captured during early afternoon overpasses of the MODIS and VIIRS. The temporal shift in fire activity complicates such forecast, leading to a substantial underestimation of emissions. Intense stubble burning and the resulting air pollution highlight the need for effective crop residue management practices for mitigating the frequency of open biomass burning and thereby reducing episodic degradation of air quality and its associated public health and economic impacts.

post-monsoon stubble burning; northwestern India;↗

Ice Nucleating Particles, Aerosols and Clouds over the Higher Latitude Southern Ocean

Improved understanding of aerosol-cloud-precipitation interactions is a key goal of the DOE-ASR program. Ice Nucleating Particles (INPs) are rare aerosol particles that trigger ice formation in clouds. By doing so they may initiate precipitation in clouds, and this impacts cloud lifetime. Nowhere may the role of INPs be manifested so dramatically as in the Southern Ocean (SO). This vast region was the focus of the DOE-ARM MARCUS (Measurement of Aerosols, Radiation and CloUds over the Southern Oceans) and MICRE (Macquarie Island Cloud and Radiation Experiment) campaigns. In the SO, a deficit of INPs has been a leading hypothesis for why liquid clouds persist, which underlays the bias of global climate models in predicting excess shortwave radiation reaching the ocean surface between 55°S and Antarctica. In this study, we used the 6-month acquisition of INP data on four MARCUS cruises from Hobart, Tasmania, to Antarctica, and the single-location full annual INP cycle observed in MICRE, to gain a holistic, quantitative picture of the INP number concentrations as a function of temperature and their variability over SO latitudes from 43°S to Antarctica. Augmenting initial processing of filter collections of aerosol particles during MARCUS and MICRE, additional freezing studies following thermal treatment (to removes INPs associated with microbes) and peroxide treatment (to remove all organics), as well as ionic chemistry and total organic carbon analysis, were applied to samples to improve time and space resolution, and differentiate marine from terrestrial contributions. We also performed Next Generation Sequencing to determine that bacterial composition can tag periods of enhanced and degraded marine organic INPs and the relative lack of occurrence of land-sourced bio-particles. Through these and other analyses, INP data were categorized into representative source types (e.g., marine versus terrestrially-influenced). The products of our analyses are being used to parameterize ice formation for use in numerical modeling studies to determine if, and when, ocean-derived INPs control the microphysical composition and radiative balance of SO clouds in a manner not presently captured by global models. The specific objectives and approaches proposed were largely accomplished, including: 1) Additional and value-added chemical and biological analyses of archived samples of aerosols were completed to categorize and quantify INP types and concentrations over the SO. 2) Analyses were completed to place INP data in meteorological and aerosol context. 3) Sea spray aerosol INPs were demonstrated to dominate the SO region in all but episodic events. 4) New parameterizations for INP sources relevant to the SO were constructed. 5) Collaborative modeling studies demonstrating the crucial role of INPs in determining cloud radiative and precipitation properties were conducted with MARCUS and MICRE partners. Collaborative publication submissions in this regard ensued within the two year study. This work will advance the science of interactions of aerosols, clouds and precipitation, to improve their representation in regional and global climate models. Results have informed representation of primary ice crystal nucleation, give inference to the role of secondary processes, and improve investigations of aerosol influences on clouds via ice nucleation over SO high latitudes and similar vast ocean regions. Application will improve representation of such interactions for clouds in both regional and global climate models. The data base and analyses methods applied will continue to serve research studies in the future.

58 GEOSCIENCES↗

DownScaleBench for developing and applying a deep learning based urban climate downscaling- first results for high-resolution urban precipitation climatology over Austin, Texas

Abstract Cities need climate information to develop resilient infrastructure and for adaptation decisions. The information desired is at the order of magnitudes finer scales relative to what is typically available from climate analysis and future projections. Urban downscaling refers to developing such climate information at the city (order of 1 – 10 km) and neighborhood (order of 0.1 – 1 km) resolutions from coarser climate products. Developing these higher resolution (finer grid spacing) data needed for assessments typically covering multiyear climatology of past data and future projections is complex and computationally expensive for traditional physics-based dynamical models. In this study, we develop and adopt a novel approach for urban downscaling by generating a general-purpose operator using deep learning. This ‘DownScaleBench’ tool can aid the process of downscaling to any location. The DownScaleBench has been generalized for both in situ (ground- based) and satellite or reanalysis gridded data. The algorithm employs an iterative super-resolution convolutional neural network (Iterative SRCNN) over the city. We apply this for the development of a high-resolution gridded precipitation product (300 m) from a relatively coarse (10 km) satellite-based product (JAXA GsMAP). The high-resolution gridded precipitation datasets is compared against insitu observations for past heavy rain events over Austin, Texas, and shows marked improvement relative to the coarser datasets relative to cubic interpolation as a baseline. The creation of this Downscaling Bench has implications for generating high-resolution gridded urban meteorological datasets and aiding the planning process for climate-ready cities.

Singh, Manmeet (ORCID:0000000233747149)↗

Yonkers Urban Development II: Leveraging NASA Earth Observations to Support Modeling Urban Cooling Interventions and Urban Heat Vulnerability in Yonkers, New York

The City of Yonkers, New York, located in Westchester County, is experiencing rising temperatures which are a growing threat to the health and safety of its residents. Furthermore, the risk of heat-related illnesses and mortality disproportionately affects neighborhoods in Yonkers historically subjected to race-based housing segregation. To better understand these inequities, Groundwork Hudson Valley and NASA DEVELOP collaborated for a second term to evaluate community-level heat vulnerability, landcover distribution, street-level thermal comfort, and modeled urban cooling interventions. This team applied 2019 5-year American Community Survey (ACS) data and social and biophysical heat vulnerability variables established by the New York State Department of Health (NYSDOH), along with land surface temperature (LST) data collected from Landsat 8 Thermal Infrared Sensor (TIRS), and ISS ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) to identify communities in Yonkers in need of prioritized cooling intervention at the census tract level. Data from the Real-Time Mesoscale Analysis (RTMA) provided relevant meteorological data for the ENVI-met model to conduct street-level thermal observations and model tree canopy cooling interventions in the Yonkers neighborhoods of Kimball and Old 7th Ward. The project results will support the prioritization and equitable distribution of cooling infrastructure in identified neighborhoods. Additionally, Groundwork Hudson Valley will use the analyses as a heat literacy tool to improve advocacy efforts and inform both residents and officials about how investment in deliberate modification to tree canopy cover improves the city’s thermal environment and helps mitigate extreme heat.

Tamara Barbakova↗

Datasets for Widespread Residential Space Heating Electrification in Texas

In this experiment, we explore long term patterns in electricity demand driven by the dual effects of full electrification of space heating in Texas (by adoption of electric heat pumps), and climate change. We use a predictive model of electricity demand, climate projections, and an open source nodal power system (DC Optimal Power Flow) model of the Electric Reliability Council of Texas (ERCOT) system. Heat pumps are a more energy efficient way of providing space heating and cooling in homes. We attempt to exhaustively investigate the impacts of full residential space heating electrification by adoption of heat pumps for the segment of Texas households that currently rely on fossil fuels (about 40%), while simultaneously incorporating climate change meteorological variables. We explore a range of scenarios of heat pump efficiency and climate uncertainty over a long period of future years (2020-2099). In total, the simulation experiment generates 1,280 simulation years of hourly data. We report and analyze results in form of impacts on residential load, total load, peak load, seasonality of peaking, and reliability measured by occurrence and frequency of loss of load events. While the experiment is for ERCOT, the insights and approach can be applied to other regions. The results from the analysis can inform system planners on a range of potential capacity requirements/ reliability implications and/or risks of full space heating electrification via the adoption of electric heat pumps, given the uncertainty in the scenarios/ climate futures. The dataset includes model output for residential, non residential and total load, and the results from the GO ERCOT model runs for 4 RCP Scenarios (RCP 4.5 Cooler, RCP 4.5 Hotter, RCP 8.5 Cooler, RCP 8.5 Hotter), 4 Heating electrification Scenarios (Base , Standard Efficiency HP, High Efficiency HP, Ultra-High Efficiency HP) over 80 years (2020-2099). The metrological variables at BA scale were weighted weighted using population projections consistent with the SSP3 scenario.

Climate Change↗

Revision of the Dobson total ozone series at Hohenpeissenberg

Total ozone measurements with the Dobson No. 104 (D 104) have been performed at the Meteorological Observatory Hohenpeissenberg since 1967. A critical review of this time series and the comparison with other instruments like TOMS or Brewer spectrophotometer revealed some intervals with uncertainties. Especially in the early eighties a monthly mean bias of about minus 3 percent to TOMS-data with annual variations depending on the mean sun-height exists. An extreme amplitude of 5.6 percent occurs in 1980 with minus 0.76 percent (February) and minus 6.36 percent (July). Two different methods were applied to reprocess the Dobson data set. A comparison of the differently recalculated data showed, that the application of N-corrections by means of the standard-lamp tests starting from the reference values of the Arosa Intercomparison 1986 yields better results than the N-corrections based on a Langley-pilot of the Arosa Intercomparison 1978. The extreme amplitude of the year 1980 is now reduced to 3.02 percent. There is still a slight drift in the monthly and yearly mean differences between TOMS and revised Dobson data. It cannot be excluded, that the satellite data may be responsible for the trend.

Koehler, U.↗

Stratospheric and mesospheric pressure-temperature profiles from rotational analysis of CO2 lines in atmospheric trace molecule spectroscopy/ATLAS 1 infrared solar occultation spectra

A simple, classical, and expedient method for the retrieval of atmospheric pressure-temperature profiles has been applied to the high-resolution infrared solar absorption spectra obtained with the atmospheric trace molecule spectroscopy (ATMOS) instrument. The basis for this method is a rotational analysis of retrieved apparent abundances from CO2 rovibrational absorption lines, employing existing constituent concentration retrieval software used in the analysis of data returned by ATMOS. Pressure-temperature profiles derived from spectra acquired during the ATLAS 1 space shuttle mission of March-April 1992 are quantitatively evaluated and compared with climatological and meteorological data as a means of assessing the validity of this approach.

Stiller, G. P.↗

A Multisensor Investigation of Convection During HyMeX SOP1 IOP13

A multisensor analysis of the convective precipitation event occurred over Rome during the IOP13 (October 15th, 2012) of the HyMeX (Hydrological cycle in the Mediterranean eXperiment) Special Observation Period (SOP) 1 is presented. Thanks to the cooperation among Italian meteorological services and scientific community and a specific agreement with NASA-GSFC, different types of devices for meteorological measurements were made available during the HyMeX SOP.1. For investigating this event, used are the 3-D lightning data provided by the LINET, the CNR ISAC dual-pol C-band radar (Polar 55C), located in Rome, the Drop Size Distributions (DSD) collected by the 2D Video Disdrometer (2DVD) and the collocated Micro Rain Radar (MRR) installed at the Radio Meteorology Lab. of "Sapienza" University of Rome, located 14 km from the Polar 55C radar. The relation between microphysical structure and electrical activity during the convective phase of the event was investigated using LINET lightning data and Polar 55C (working both in PPI and RHI scanning mode) observations. Location of regions of high horizontal reflectivity (Zh) values ( > 50 dBz), indicating convective precipitation, were found to be associated to a high number of LINET strokes. In addition, an hydrometeor classification scheme applied to the Polar 55C scans was used to detect graupel and to identify a relation between number of LINET strokes and integrated IWC of graupel along the event. Properties of DSDs measured by the 2DVD and vertical DSD profiles estimated by MRR and their relation with the lighting activity registered by LINET were investigated with specific focus on the transition from convective to stratiform regimes. A good agreement was found between convection detected by these instruments and the number of strokes detected by LINET.

Roberto, N.↗

Estimation of Snow Parameters Based on Passive Microwave Remote Sensing and Meteorological Information

A method to incorporate passive microwave remote sensing measurements within a spatially distributed snow hydrology model to provide estimates of the spatial distribution of Snow Water Equivalent (SWE) as a function of time is implemented. The passive microwave remote sensing measurements are at 25 km resolution. However, in mountain regions the spatial variability of SWE over a 25 km footprint is large due to topographic influences. On the other hand, the snow hydrology model has built-in topographic information and the capability to estimate SWE at a 1 km resolution. In our work, the snow hydrology SWE estimates are updated and corrected using SSM/I passive microwave remote sensing measurements. The method is applied to the Upper Rio Grande River Basin in the mountains of Colorado. The change in prediction of SWE from hydrology modeling with and without updating is compared with measurements from two SNOTEL sites in and near the basin. The results indicate that the method incorporating the remote sensing measurements into the hydrology model is able to more closely estimate the temporal evolution of the measured values of SWE as a function of time.

Tsang, Leung↗

Hacking Kilometer-Scale Models: A Participative Model for Climate Information

In May 2025, nearly 700 participants from all around the world coalesced at 10 regional nodes and a few satellite nodes to take part in a global hackathon of kilometer-scale (horizontal grid spacing < 10 km) regional and global Earth system models. Exciting science is emerging from these efforts, ranging across novel model analysis, new ways of integrating with satellite data, and emulation with machine learning. New technologies were trialed that enable the community to work in new and complementary ways to democratize access to global information at a local scale from a set of the world’s highest-resolution climate models. The hackathon demonstrated how exascale data can be organized to be accessible to anyone. Fundamentally, the community could apply these techniques and technologies to move toward more participative models for coproduction and delivery of diverse sources of climate information for climate scientists and citizens alike.

Climate models↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗