Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wind forecasting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

The core mass growth and stellar lifetime of thermally pulsing asymptotic giant branch stars

We establish new constraints on the intermediate-mass range of the initial-final mass relation, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). These constraints derive from newly discovered (bright) white dwarfs in the nearby Hyades and Praesepe star clusters, including a total of 18 high signal-to-noise ratio measurements with progenitor masses of M {sub initial} = 2.8-3.8 M {sub ☉}. We also include a new analysis of existing white dwarfs in the older NGC 6819 and NGC 7789 star clusters, M {sub initial} = 1.6 and 2.0 M {sub ☉}. Over this range of initial masses, stellar evolutionary models for metallicity Z {sub initial} = 0.02 predict the maximum growth of the core of TP-AGB stars. By comparing the newly measured remnant masses to the robust prediction of the core mass at the first thermal pulse on the AGB (i.e., from stellar interior models), we establish several findings. First, we show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with M {sub initial} = 1.6 to 2.0 M {sub ☉}. At larger masses, the core-mass growth decreases steadily to ∼10% at M {sub initial} = 3.4 M {sub ☉}, after which there is a small hint of a upturn out to M {sub initial} = 3.8 M {sub ☉}. These observations are in excellent agreement with predictions from the latest TP-AGB evolutionary models in Marigo et al. We also compare to models with varying efficiencies of the third dredge-up and mass loss, and demonstrate that the process governing the growth of the core is largely the stellar wind, while the third dredge-up plays a secondary, but non-negligible role. Based on the new white dwarf measurements, we perform an exploratory calibration of the most popular mass-loss prescriptions in the literature, as well as of the third dredge-up efficiency as a function of the stellar mass. Finally, we estimate the lifetime and the integrated luminosity of stars on the TP-AGB to peak at t ∼ 3 Myr and E = 1.2 × 10{sup 10} L {sub ☉} yr for M {sub initial} ∼ 2 M {sub ☉} (t ∼ 2 Myr for luminosities brighter than the red giant branch tip at log (L/L {sub ☉}) > 3.4), decreasing to t = 0.4 Myr and E = 6.1 × 10{sup 9} L {sub ☉} yr for stars with M {sub initial} ∼ 3.5 M {sub ☉}. The implications of these results are discussed, especially with respect to general studies aimed at characterizing the integrated light output of TP-AGB stars in population synthesis models.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of the Rapid Refresh Numerical Weather Prediction Model over Arctic Alaska

Abstract Despite a need for accurate weather forecasts for societal and economic interests in the U.S. Arctic, thorough evaluations of operational numerical weather prediction in the region have been limited. In particular, the Rapid Refresh Model (RAP), which plays a key role in short-term forecasting and decision-making, has seen very limited assessment in northern Alaska, with most evaluation efforts focused on lower latitudes. In the present study, we verify forecasts from version 4 of the RAP against radiosonde, surface meteorological, and radiative flux observations from two Arctic sites on the northern Alaskan coastline, with a focus on boundary layer thermodynamic and dynamic biases, model representation of surface inversions, and cloud characteristics. We find persistent seasonal thermodynamic biases near the surface that vary with wind direction, and may be related to the RAP’s handling of sea ice and ocean interactions. These biases seem to have diminished in the latest version of the RAP (version 5), which includes refined handling of sea ice, among other improvements. In addition, we find that despite capturing boundary layer temperature profiles well overall, the RAP struggles to consistently represent strong, shallow surface inversions. Further, while the RAP seems to forecast the presence of clouds accurately in most cases, there are errors in the simulated characteristics of these clouds, which we hypothesize may be related to the RAP’s treatment of mixed-phase clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unified Modeling Architecture for Load Management in Extreme Heat: The New York City Case

Integration of renewable resources to meet growing energy demand is becoming a global priority under decarbonization mandates. This study contributes to ongoing efforts on this key subject by assessing the feasibility of using coastal-urban renewable energy resources, namely, offshore wind and rooftop photovoltaic systems, to meet electricity demand of New York City during the intense recent heat wave period of June 2025. A unified modeling framework, based on the urbanized weather research and forecasting model, is used to simulate climate, renewable resources, and energy demand variables. Findings show significant energy load mismatch of approximately 1150 GWh over the month, between the demand and the combined renewable generation outcome. Three storage integration scenarios are analyzed to mitigate the deficits, reducing said deficits by a minimum of approximately 9% over the duration of the month. This study provides a transferable modeling framework tool for evaluating renewable integration in dense urban environments that can be used by grid operators to support grid resilience during extreme heat events.

54 ENVIRONMENTAL SCIENCES↗

The Pale Green Dot: A Method to Characterize Proxima Centauri b Using Exo-Aurorae

We examine the feasibility of detecting auroral emission from the potentially habitable exoplanet Proxima Centauri b. Detection of aurorae would yield an independent confirmation of the planet’s existence, constrain the presence and composition of its atmosphere, and determine the planet’s eccentricity and inclination, thereby breaking the mass-inclination degeneracy. If Proxima Centauri b is a terrestrial world with an Earth-like atmosphere and magnetic field, we estimate that the power at the 5577 Å O i auroral line is on the order of 0.1 TW under steady-state stellar wind, or ∼100× stronger than that on Earth. This corresponds to a planet–star contrast ratio of 10{sup −6}--10{sup −7} in a narrow band about the 5577 Å line, though higher contrast (10{sup −4}--10{sup −5}) may be possible during periods of strong magnetospheric disturbance (auroral power 1–10 TW). We searched the Proxima Centauri b HARPS data for the 5577 Å line and for other prominent oxygen and nitrogen lines, but find no signal, indicating that the O i auroral line contrast must be lower than 2×10{sup −2} (with power ≲3000 TW), consistent with our predictions. We find that observations of 0.1 TW auroral emission lines are likely infeasible with current and planned telescopes. However, future observations with a space-based coronagraphic telescope or a ground-based extremely large telescope (ELT) with a coronagraph could push sensitivity down to terawatt oxygen aurorae (contrast 7×10{sup −6}) with exposure times of ∼1 day. If a coronagraph design contrast of 10{sup −7} can be achieved with negligible instrumental noise, a future concept ELT could observe steady-state auroral emission in a few nights.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of an atmospheric transport model for annual inverse estimates of California greenhouse gas emissions

Abstract Atmospheric inverse estimates of gas emissions depend on transport model predictions, hence driving a need to assess uncertainties in the transport model. In this study we assess the uncertainty in WRF‐STILT (Weather Research and Forecasting and Stochastic Time‐Inverted Lagrangian Transport) model predictions using a combination of meteorological and carbon monoxide (CO) measurements. WRF configurations were selected to minimize meteorological biases using meteorological measurements of winds and boundary layer depths from surface stations and radar wind profiler sites across California. We compare model predictions with CO measurements from four tower sites in California from June 2013 through May 2014 to assess the seasonal biases and random errors in predicted CO mixing ratios. In general, the seasonal mean biases in boundary layer wind speed (< ~ 0.5 m/s), direction (< ~ 15°), and boundary layer height (< ~ 200 m) were small. However, random errors were large (~1.5–3.0 m/s for wind speed, ~ 40–60° for wind direction, and ~ 300–500 m for boundary layer height). Regression analysis of predicted and measured CO yielded near‐unity slopes (i.e., within 1.0 ± 0.20) for the majority of sites and seasons, though a subset of sites and seasons exhibit larger (~30%) uncertainty, particularly when weak winds combined with complex terrain in the South Central Valley of California. Looking across sites and seasons, these results suggest that WRF‐STILT simulations are sufficient to estimate emissions of CO to up to 15% on annual time scales across California.

Bagley, Justin E.↗

Quantifying the Thermodynamic Impacts on the Atmospheric Boundary Layer due to the Sea Breeze in the Coastal Houston Region

The atmospheric boundary layer (ABL) is unique in coastal regions because of kinematic and thermodynamic influences from continental and marine environments. Sea-breeze (SB) circulations act to equilibrate the land–sea temperature gradient through advecting marine air onshore. The strength of the SB varies in terms of stability, temperature, and moisture advection and influences air quality and weather forecasts. The Tracking Aerosol Convection Interactions Experiment (TRACER) collected a wealth of data on coastal boundary layer evolution, including observations from uncrewed aerial systems (UASs). Vertical profiles of temperature, humidity, and winds were collected by the OU CopterSonde UAS from June to September in the coastal region of Houston. These profiles offer 5-m vertical resolution, on average, every 30 min through diurnal transitions, SB events, and nearby deep convection. During the campaign, CopterSonde observations were gathered through 17 SB events, six of which led to convection initiation. The UAS data can resolve the thermodynamic evolution and interactions between the SB and the preexisting convective boundary layer. Results show large variability across observed SBs and their impacts on temperature and moisture. The intensity of thermodynamic changes depends on the time of sea-breeze passage and influence from the Galveston Bay Breeze, a secondary marine circulation commonly observed in this region. In quantifying the spectrum of SB impacts, equivalent potential temperature θ e is used to contextualize its role in convection initiation and evolution. In conclusion, while all SBs tend to increase θ e from moisture advection, the rate and timing of the θ e rise can distinguish convective from nonconvective cases.

54 ENVIRONMENTAL SCIENCES↗

Mariner Venus-Mercury 1973 mission solar proton environment: fluence and dose

A derivation is presented for proton fluence over the duration of the Mariner Venus-Mercury 1973 (MVM73) mission in terms of a relatively constant, low energy component, the solar wind, and a probabilistic high energy component from discrete solar events. An updated correlation of yearly energetic proton fluence with yearly average sunspot number is presented. This correlation and sunspot cycle forecasts for the period of the MVM73 mission (late 1973 through early 1975) form the basis for the high energy proton fluence estimates with various confidence levels. Uncertainties in the probability estimates and in calculation of the scaling with distance from the Sun are discussed. Selection of a particular 95-percentile model as the design constraint is recommended, and reasons are presented for not using the worst-case model. Interior fluences were calculated and expressed in terms useful for spacecraft design.

Thomas, J. R.↗

NASA's program in lidar remote sensing

Several major NASA research efforts in lidar remote sensing are reviewed, with attention given to hardware and key sensor issues along with results and expectations. The discussion covers temperature and pressure measurements, measurements methods and instrumentation, pressure data, wind field measurements, atmospheric backscatter measurement, aerosol and cloud measurements, and water vapor measurement. Consideration is also given to the applicability of lidar measurements to problems of operational weather analysis and forecasting, climate studies, mesoscale and severe storm analysis and forecastig, and studies of atmosphere/surface interface.

Theon, John S.↗

Impact of the ongoing Amazonian deforestation on local precipitation: A GCM simulation study

Numerical simulation experiments were conducted to delineate the influence of in situ deforestation data on episodic rainfall by comparing two ensembles of five 5-day integrations performed with a recent version of the Goddard Laboratory for Atmospheres General Circulation Model (GCM) that has a simple biosphere model (SiB). The first set, called control cases, used the standard SiB vegetation cover (comprising 12 biomes) and assumed a fully forested Amazonia, while the second set, called deforestation cases, distinguished the partially deforested regions of Amazonia as savanna. Except for this difference, all other initial and prescribed boundary conditions were kept identical in both sets of integrations. The differential analyses of these five cases show the following local effects of deforestation. (1) A discernible decrease in evapotranspiration of about 0.80 mm/d (roughly 18%) that is quite robust in the averages for 1-, 2-, and 5-day forecasts. (2) A decrease in precipitation of about 1.18 mm/d (roughly 8%) that begins to emerge even in 1-2 day averages and exhibits complex evolution that extends downstream with the winds. (3) A significant decrease in the surface drag force (as a consequence of reduced surface roughness of deforested regions) that, in turn, affects the dynamical structure of moisture convergence and circulation. The surface winds increase significantly during the first day, and thereafter the increase is well maintained even in the 2- and 5-day averages.

Walker, G. K.↗

Convective Systems Over the Japan Sea: Cloud-Resolving Model Simulations

Wintertime observations of MCSs (Mesoscale Convective Systems) over the Sea of Japan - 2001 (WMO-01) were collected from January 12 to February 1, 2001. One of the major objectives is to better understand and forecast snow systems and accompanying disturbances and the associated key physical processes involved in the formation and development of these disturbances. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, wind profilers, radiometers, etc.) during WMO-01 provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with winter storms over the Sea of Japan region. WMO-01 also provided estimates of the apparent heat source (Q1) and apparent moisture sink (Q2). The vertical integrals of Q1 and Q2 are equal to the surface precipitation rates. The horizontal and vertical adjective components of Q1 and Q2 can be used as large-scale forcing for the Cloud Resolving Models (CRMs). The Goddard Cumulus Ensemble (GCE) model is a CRM (typically run with a 1-km grid size). The GCE model has sophisticated microphysics and allows explicit interactions between clouds, radiation, and surface processes. It will be used to understand and quantify precipitation processes associated with wintertime convective systems over the Sea of Japan (using data collected during the WMO-01). This is the first cloud-resolving model used to simulate precipitation processes in this particular region. The GCE model-simulated WMO-01 results will also be compared to other GCE model-simulated weather systems that developed during other field campaigns (i.e., South China Sea, west Pacific warm pool region, eastern Atlantic region and central USA).

Tao, Wei-Kuo↗

Gravity Waves and Mesospheric Clouds in the Summer Middle Atmosphere: A Comparison of Lidar Measurements and Ray Modeling of Gravity Waves Over Sondrestrom, Greenland

We conducted gravity wave ray-tracing experiments within an atmospheric region centered near the ARCLITE lidar system at Sondrestrom, Greenland (67N, 310 deg E), in efforts to understand lidar observations of both upper stratospheric gravity wave activity and mesospheric clouds during August 1996 and the summer of 2001. The ray model was used to trace gravity waves through realistic three-dimensional daily-varying background atmospheres in the region, based on forecasts and analyses in the troposphere and stratosphere and climatologies higher up. Reverse ray tracing based on upper stratospheric lidar observations at Sondrestrom was also used to try to objectively identify wave source regions in the troposphere. A source spectrum specified by reverse ray tracing experiments in early August 1996 (when atmospheric flow patterns produced enhanced transmission of waves into the upper stratosphere) yielded model results throughout the remainder of August 1996 that agreed best with the lidar observations. The model also simulated increased vertical group propagation of waves between 40 km and 80 km due to intensifying mean easterlies, which allowed many of the gravity waves observed at 40 km over Sondrestrom to propagate quasi-vertically from 40-80 km and then interact with any mesospheric clouds at 80 km near Sondrestrom, supporting earlier experimentally-inferred correlations between upper stratospheric gravity wave activity and mesospheric cloud backscatter from Sondrestrom lidar observations. A pilot experiment of real-time runs with the model in 2001 using weather forecast data as a low-level background produced less agreement with lidar observations. We believe this is due to limitations in our specified tropospheric source spectrum, the use of climatological winds and temperatures in the upper stratosphere and mesosphere, and missing lidar data from important time periods.

LIDAR (LIGHT DETECTION AND RANGING)↗

AE Geomagnetic Index Predictability for High Speed Solar Wind Streams: A Wavelet Decomposition Technique

High speed solar wind streams cause geomagnetic activity at Earth. In this study we have applied a wavelet interactive filtering and reconstruction technique on the solar wind magnetic field components and AE index series to allowed us to investigate the relationship between the two. The IMF Bz component was found as the most significant solar wind parameter responsible by the control of the AE activity. Assuming magnetic reconnection associated to southward directed Bz is the main mechanism transferring energy into the magnetosphere, we adjust parameters to forecast the AE index. The adjusted routine is able to forecast AE, based only on the Bz measured at the L1 Lagrangian point. This gives a prediction approximately 30-70 minutes in advance of the actual geomagnetic activity. The correlation coefficient between the observed AE data and the forecasted series reached values higher than 0.90. In some cases the forecast reproduced particularities observed in the signal very well.The high correlation values observed and the high efficacy of the forecasting can be taken as a confirmation that reconnection is the main physical mechanism responsible for the energy transfer during HILDCAAs. The study also shows that the IMF Bz component low frequencies are most important for AE prediction.

Space weather↗

Earth Global Reference Atmospheric Model (GRAM) Overview and Updates: DOLWG Meeting

What is Earth-GRAM (Global Reference Atmospheric Model): Provides monthly mean and standard deviation for any point in atmosphere - Monthly, Geographic, and Altitude Variation; Earth-GRAM is a C++ software package - Currently distributed as Earth-GRAM 2016; Atmospheric variables included: pressure, density, temperature, horizontal and vertical winds, speed of sound, and atmospheric constituents; Used by engineering community because of ability to create dispersions in atmosphere at a rapid runtime - Often embedded in trajectory simulation software; Not a forecast model; Does not readily capture localized atmospheric effects.

White, Patrick↗

Building a standardized Observing System Simulation Experiment (OSSE) framework for Mars

We advocate that the Decadal Survey recommends the NASA Science Mission Directorate to develop a rigorous Observing System Simulation Experiment (OSSE) framework for Mars, to optimize future atmospheric observations. Atmospheric conditions on Mars are a potential hazard source for landing missions. Errors in the estimates of atmospheric density profiles, inadequate knowledge of wind vertical structure and dust concentration as a function of height are likely causes of uncertainty at the landing site on the order of kilometers. An operational real-time weather forecasting capability for Mars would reduce such uncertainties, carrying enormous benefits to future robotic missions, and would be an invaluable prerequisite for human missions.A real-time forecasting capability relies upon three fundamental components: a critical mass of observing systems, a data assimilation system (DAS), and a global forecast model. The DAS allows the model to ingest the data effectively, optimizing the observational information content,and transforming them into a gridded representation of the atmosphere at a given time, called an ‘analysis’. The analysis is the best estimate of the atmospheric state for that time, and also represents a set of ‘initial conditions’ from which a global model can be initialized, to predict a future state of the atmosphere. The connection between analysis and forecast represents the foundation of modern weather forecasting. However, from the point of view of a forecast system,not all observations are equally impactful, partially because of the problem of “observational error correlation”, one important research topic in data assimilation development. For the Earth, partly due to the spontaneous and deregulated development of observations and forecast capabilities worldwide for more than half a century,the use of observations in contemporary operational forecast systems is suboptimal, with many potentially useful data being underutilized. On the contrary, Mars atmospheric scientists are in the unique situation of designing the next-generation observing systems by learning from the experience gathered on the Earth, so as to assure that the future instruments are specifically optimized to give the maximum benefit to a future weather forecast capability.An immensely powerful tool that has been firmly established by atmospheric scientists on the Earth is represented by a properly designed OSSE framework. A realistic OSSE framework cannot only quantify the benefit of future data types, be them surface based or space borne, but can also help design and optimize an entire observational network. Furthermore, OSSEs can provide deep insights into an atmosphere’s behavior, by addressing conceptual problems of its intrinsic predictability and delineating the regions or features of the atmosphere which are more sensitive to additional data and would benefit from a denser sampling. The difficulties posed by OSSEs are fundamentally different for Earth and Mars. For Earth, the enormous data volume imposes a tremendous constraint on any innovation in the observing systems: it is very hard for a single sensor to impact the skill. For Mars, the problem is the opposite: almost any additional instrument will exert some impact. However, OSSEs can help to evaluate the cost/benefit for every sensor and suggest optimal data configuration and density.The purpose of this white paper is to provide an introduction to a rigorously designed OSSE framework, explain the underlying problems and challenges, and engage the Mars community to collaborate with Earth Atmospheric scientists in order to develop a joint-OSSE framework for Mars with the largest consensual basis possible. An OSSE infrastructure would increase the understanding of the Martian atmosphere, would help NASA to optimize instrument specifications and orbit choice, providing the maximium benefit for a given expenditure of resources, and could even help establishing a roadmap for a future real-time weather forecasting capability.

Oreste Reale↗

From the Solar Limb and Out: Results from the Wide-Field EUV Image Campaigns with GOES/SUVI

Traditional approaches to tracking solar outflows for space weather forecasting rely primarily on coronagraph images, which generally observe the solar corona above a minimum height of about 2.5 solar radii. EUV images have been widely used to characterize features on the solar disk, but the limited fields of view of most current EUV imagers have prevented their use for tracking outflows through the inner and middle coronae. A series of off-point campaigns with the GOES 16-18 Solar Ultraviolet Imager (SUVI) between 2018 and 2022 from three Flight Models have provided an opportunity to assess the value of extended EUV images for space weather forecasting applications. These new results demonstrate that wide field-of-view EUV images are useful for characterizing the early onset of eruptive events and tracking smaller outflow into the solar wind. They also reveal the origins of shocks that are known to accelerate particles and drive solar energetic particle (SEP) events. Because CMEs generally experience the bulk of their acceleration below the height of white light coronagraphic observations, these images provide information about the origins of these events that has not been available traditionally. Together with coronagraphic measurements, EUV images provide the continuous views needed to connect CMEs back to their source regions. Here, we present these new SUVI observations and discuss their potential use in space weather operations.

SUVI↗

Automatic Detection of Large-scale Flux Ropes and Their Geoeffectiveness with a Machine-learning Approach

Detecting large-scale flux ropes (FRs) embedded in interplanetary coronal mass ejections (ICMEs) and assessing their geoeffectiveness are essential, since they can drive severe space weather. At 1 au, these FRs have an average duration of 1 day. Their most common magnetic features are large, smoothly rotating magnetic fields. Their manual detection has become a relatively common practice over decades, although visual detection can be time-consuming and subject to observer bias. Our study proposes a pipeline that utilizes two supervised binary classification machine-learning models trained with solar wind magnetic properties to automatically detect large-scale FRs and additionally determine their geoeffectiveness. The first model is used to generate a list of autodetected FRs. Using the properties of the southward magnetic field, the second model determines the geoeffectiveness of FRs. Our method identifies 88.6% and 80% of large-scale ICMEs (duration 1day) observed at 1au by the Wind and the Solar TErrestrial RElations Observatory missions, respectively. While testing with continuous solar wind data obtained from Wind, our pipeline detected 56 of the 64 large-scale ICMEs during the 2008–2014 period (recall = 0.875), but also many false positives (precision = 0.56), as we do not take into account any additional solar wind properties other than the magnetic properties. We find an accuracy of 0.88 when estimating the geoeffectiveness of the autodetected FRs using our method. Thus, in space-weather nowcasting and forecasting at L1 or any planetary missions, our pipeline can be utilized to offer a first-order detection of large-scale FRs and their geoeffectiveness.

Sanchita Pal↗

Environmental controls on simulated deep moist convection initiation occurring during RELAMPAGO-CACTI

This study synthesizes the results of 13 high-resolution simulations of deep convective updrafts forming over idealized terrain using environments observed during the RELAMPAGO and CACTI field projects. Using composite soundings from multiple observed cases, and variations upon them, we explore the sensitivity of updraft properties (e.g., size, buoyancy, and vertical pressure gradient forces) to influences of environmental relative humidity, wind shear, and mesoscale orographic forcing that support or suppress deep convection initiation (CI). Emphasis is placed on differentiating physical processes affecting the development of updrafts (e.g., entrainment-driven dilution of updrafts) in environments typifying observed successful and null (i.e., no CI despite affirmative operational forecasts) CI events. Thermally-induced mesoscale orographic lift favors the production of deep updrafts originating from ~1–2-km-wide boundary layer thermals. Simulations without terrain forcing required much larger ( ~5-km-wide) thermals to yield precipitating convection. CI outcome was quite sensitive to environmental relative humidity; updrafts with increased buoyancy, depth, and intensity thrived in otherwise inhospitable environments by simply increasing the free tropospheric relative humidity. This implicates the entrainment of free-tropospheric air into updrafts as a prominent governor of CI, consistent with previous studies. Sensitivity of CI to the environmental wind is manifested by: 1) low-level flow affecting the strength and depth of mesoscale convergence along the terrain, and 2) clouds encountering updraft-suppressing pressure gradient forces while interacting with vertical wind shear in the free-troposphere. Among the ensemble of thermals occurring in each simulation, the widest deep updrafts in each simulation were the most sensitive to environmental influences.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning for Daily Forecasts of Arctic Sea Ice Motion: An Attribution Assessment of Model Predictive Skill

Physics-based simulations of Arctic sea ice are highly complex, involving transport between different phases, length scales, and time scales. Resultantly, numerical simulations of sea ice dynamics have a high computational cost and model uncertainty. We employ data-driven machine learning (ML) to make predictions of sea ice motion. The ML models are built to predict present-day sea ice velocity given present-day wind velocity and previous-day sea ice concentration and velocity. Models are trained using reanalysis winds and satellite-derived sea ice properties. We compare the predictions of three different models: persistence (PS), linear regression (LR), and a convolutional neural network (CNN). We quantify the spatiotemporal variability of the correlation between observations and the statistical model predictions. Additionally, we analyze model performance in comparison to variability in properties related to ice motion (wind velocity, ice velocity, ice concentration, distance from coast, bathymetric depth) to understand the processes related to decreases in model performance. Results indicate that a CNN makes skillful predictions of daily sea ice velocity with a correlation up to 0.81 between predicted and observed sea ice velocity, while the LR and PS implementations exhibit correlations of 0.78 and 0.69, respectively. The correlation varies spatially and seasonally: lower values occur in shallow coastal regions and during times of minimum sea ice extent. LR parameter analysis indicates that wind velocity plays the largest role in predicting sea ice velocity on 1-day time scales, particularly in the central Arctic. Regions where wind velocity has the largest LR parameter are regions where the CNN has higher predictive skill than the LR.

54 ENVIRONMENTAL SCIENCES↗