Search NASA⌕ Search

SEARCH · Search NASA

Results for “wind advection”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Sensitivity of calculated odd nitrogen distributions to the diabatic wind fields

This paper compares atmospheric total odd nitrogen and ozone computed with two different advective wind fields, one using climatological averages of ozone and temperature to obtain monthly averaged horizontal and vertical winds and the other using measurements from the LIMS instrument. Calculations using the former data show stronger poleward and downward motion in the winter season compared to those using the LIMS data. This leads to NO(y) mixing ratios in the lower stratosphere that are about 20 percent larger in the polar regions of both hemispheres and about 40 percent higher in the equatorial region for climatological transport fields compared to those derived from LIMS data. Consequently, the NO(y) distributions calculated with the LIMS advective field show worse agreement with the NO(y) values inferred from the LIMS measurements than similar results obtained with the climatological wind field.

Boughner, R. E.↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2025

Short-Term Energy and Meteorological Impacts on Thanksgiving CO2 in Salt Lake City

Abstract Long-term, high-frequency atmospheric CO2 measurements at multiple sites in the Salt Lake City (SLC), Utah, reveal that annual and monthly CO2 variability aligns with a priori estimates of emissions from anthropogenic and biological sources. In this study, we investigate whether short-term fluctuations in anthropogenic emissions, as captured in the Vulcan3 dataset for the United States, can be detected in atmospheric CO2 observations. Specifically, we focus on Thanksgiving holidays, when traffic and energy usage patterns differ from the rest of November. Onroad CO2 emissions exhibit a double peak during weekday morning and evening rush hours but remain relatively low on weekends and Thanksgiving. Interestingly, CO2 mole fractions during Thanksgiving were higher than the rest of November at all SLC monitoring sites, particularly from 2008 to 2013. This increase is partially attributed to elevated energy-related emissions — especially residential sources — and meteorological factors such as weak wind speeds, cold temperature, and a low planetary boundary layer height (PBLH).

 While CO₂ emissions and mole fraction patterns align over time, notable spatial differences exist. For instance, the near-highway site in Murray shows the highest CO₂ mole fractions despite low local emissions, suggesting pollution transport via highways and wind advection. Random Forest model-based SHapley Additive exPlanations (SHAP) analysis reveals that onroad emissions dominate CO2 contributions on weekdays and weekends, while energy-related emissions play a larger role during Thanksgiving, alongside meteorological drivers such as wind speed and PBLH. Across six urban cities, CO2 emissions display a consistent pattern: residential and commercial (onroad) emissions peak during Thanksgiving (weekday) with substantial (minimal) year-to-year variability. These findings highlight that urban CO₂ variability is driven by the combined influence of emissions and meteorology, underscoring the need for integrated mitigation strategies. Additionally, multi-site measurements are essential for accurate source attribution and the development of effective policy interventions. 

Ryoo, Ju-Mee (ORCID:0000000234256296)↗
NASA NTRS2012

3D Drop Size Distribution Extrapolation Algorithm Using a Single Disdrometer

Determining the Z-R relationship (where Z is the radar reflectivity factor and R is rainfall rate) from disdrometer data has been and is a common goal of cloud physicists and radar meteorology researchers. The usefulness of this quantity has traditionally been limited since radar represents a volume measurement, while a disdrometer corresponds to a point measurement. To solve that problem, a 3D-DSD (drop-size distribution) method of determining an equivalent 3D Z-R was developed at the University of Central Florida and tested at the Kennedy Space Center, FL. Unfortunately, that method required a minimum of three disdrometers clustered together within a microscale network (.1-km separation). Since most commercial disdrometers used by the radar meteorology/cloud physics community are high-cost instruments, three disdrometers located within a microscale area is generally not a practical strategy due to the limitations of these kinds of research budgets. A relatively simple modification to the 3D-DSD algorithm provides an estimate of the 3D-DSD and therefore, a 3D Z-R measurement using a single disdrometer. The basis of the horizontal extrapolation is mass conservation of a drop size increment, employing the mass conservation equation. For vertical extrapolation, convolution of a drop size increment using raindrop terminal velocity is used. Together, these two independent extrapolation techniques provide a complete 3DDSD estimate in a volume around and above a single disdrometer. The estimation error is lowest along a vertical plane intersecting the disdrometer position in the direction of wind advection. This work demonstrates that multiple sensors are not required for successful implementation of the 3D interpolation/extrapolation algorithm. This is a great benefit since it is seldom that multiple sensors in the required spatial arrangement are available for this type of analysis. The original software (developed at the University of Central Florida, 1998.- 2000) has also been modified to read standardized disdrometer data format (Joss-Waldvogel format). Other modifications to the software involve accounting for vertical ambient wind motion, as well as evaporation of the raindrop during its flight time.

Lane, John↗
DOE OSTI2021

Sunlight-absorbing aerosol amplifies the seasonal cycle in low-cloud fraction over the southeast Atlantic

Abstract. The mean altitude of the smoke loading over the southeast Atlantic moves from the boundary layer in July to the free troposphere by October. This study details the month-by-month changes in cloud properties and the large-scale environment as a function of the biomass burning aerosol loading at Ascension Island (8∘ S, 14.5∘ W) from July to October, based on island measurements, satellite retrievals, and reanalysis. In July and August, the smoke loading predominantly varies within the boundary layer. During both months, the low-cloud fraction is less and is increasingly cumuliform when more smoke is present, with the exception of a late morning boundary layer deepening that encourages a short-lived cloud development. The meteorology varies little, suggesting aerosol–cloud interactions explain the cloudiness changes. September marks a transition month during which midlatitude disturbances can intrude into the Atlantic subtropics, constraining the free tropospheric aerosol closer to the African coast. Stronger boundary layer winds on cleaner days help deepen, dry, and cool much of the marine boundary layer compared to that on days with high smoke loadings, with stratocumulus reducing everywhere but at the northern deck edge. The September free troposphere is better mixed on smoky days compared to October. Longwave cooling rates, generated by a sharp water vapor gradient at the aerosol layer top, encourage a small-scale vertical mixing that could help maintain the well-mixed smoky September free troposphere. The October meteorology primarily varies as a function of the strength of the free tropospheric winds advecting aerosol offshore. The free tropospheric aerosol loading is less than in September, and the moisture variability is greater. Low-level clouds increase and are more stratiform in October when the smoke loadings are higher. The increased free tropospheric moisture can help sustain the clouds through a reduction in evaporative drying during cloud-top entrainment. Enhanced subsidence above the coastal upwelling region, increasing cloud droplet number concentrations, may further prolong cloud lifetime through microphysical interactions. Reduced subsidence underneath stronger free tropospheric winds at Ascension Island supports slightly higher cloud tops during smokier conditions. Overall, the monthly changes in the large-scale aerosol and moisture vertical structure act to amplify the seasonal cycle in low-cloud amount and morphology. This is climatically important, as cloudiness changes dominate changes in the top-of-atmosphere radiation budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Atmospheric Variability Drives Anomalies in the Bering Sea Air–Sea Heat Exchange

Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance Statement A long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change.

Hayden, Emily E. [a College of Earth, Ocean, and A↗
DOE OSTI2023

Insensitivity to propagation timing in a preview-enabled wind turbine control experiment

Lidar scanners are capable of taking measurements of a wind field upstream of a wind turbine. The wind turbine controller can use these measurements as a “preview” of future disturbances impacting the turbine. Such preview-enabled (or feedforward) controllers show superior performance to standard wind turbine control configurations based purely on a feedback architecture. To capitalize on the performance improvements that preview wind measurements can provide, feedforward control actions should be timed to coincide with the arrival of the wind field at the wind turbine location. However, the time of propagation of the wind field between the lidar measurement location and the wind turbine is not perfectly known. Moreover, the best time to take feedforward control action may not perfectly coincide with the true arrival time of the wind disturbance. This contribution presents results from an experiment where preview-enabled model predictive control was deployed on a fully-actuated, scaled model wind turbine operating in a wind tunnel testbed. In the study, we investigate the sensitivity of the controller performance to the assumed propagation delay using a range of wind input sequences. We find that the preview-enabled controller outperforms the feedback only case across a wide range of assumed propagation delays, demonstrating a level of robustness to the time alignment of the incoming disturbances.

17 WIND ENERGY↗
NASA NTRS1977

Saturn - Long-term variation of H2 and CH4 absorptions

Results are presented for observations of various H2 quadrupole lines and CH4 absorption bands in spectra of Saturn obtained with a coude scanner at each apparition since 1969. The data are found to indicate a long-term variation in Saturn's H2 equivalent widths, a seasonal dependence of H2 line strength, and a fairly steady increase in CH4 absorption during the period from 1973 to 1976. A possible correlation between planetary shading and the strength of H2 absorption features is discussed along with seasonal variations in Saturn's NH3 cloud deck and atmospheric models that correspond to the observed H2 equivalent widths. It is suggested that: (1) significant shadowing by the rings may cause dynamical atmospheric activity due to the occurrence of a temperature gradient in the penumbral boundaries of the ring shadow in the planet's upper atmosphere; (2) cooling in the shadow may be a source of a high-altitude CH4 mist; and (3) spreading of this mist over the globe by advective winds would reduce the H2 and CH4 absorptions, as observed.

Trafton, L.↗
NASA NTRS1985

The evolution of Ertel's potential vorticity during stratospheric sudden warmings

In the winter stratosphere of the Northern Hemisphere, the disruption of the westerly vortex and associated warming of polar latitudes is a well known phenomenon. It has become apparent that some important dynamical processes in the stratosphere are highly nonlinear and are best thought of locally rather than in terms of the interaction between the zonal-mean flow and harmonic waves around latitude circles. The importance of nonlinear processes was suggested by McIntyre and Palmer (1983, 1984) who used isentropic maps of Ertel's Potential Vorticity to show that during disturbed episodes material lines may become strongly and irreversibly deformed in certain places. They adopted the term planetary wave breaking to describe this process. Isentropic maps of Q are used to follow the evolution of a Canadian warming in November - December 1981 and a particularly strong warming in January 1982. The advection of Q over large distances on isentropic surfaces was a striking feature of the flow during each event. This could be identified because of our ability to follow the movement of material lines due to the approximate conservation of Q over several days. The advection of Q was a nonlinear process because its changing distribution affected the advecting wind field. The Canadian warming did not lead to a permanent change in the structure of the westerly vortex, as defined by the coarse-grain field of Q, whereas the January event was accompanied by a substantial loss of resolved Q which was never fully recovered.

Fairlie, T. D. A.↗
NASA NTRS1993

Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint

An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.

Li, Y.↗
NASA NTRS1994

Ozone formation during an episode over Europe: A 3-D chemical/transport model simulation

A 3-D regional photochemical tracer/transport model for Europe and the Eastern Atlantic has been developed based on the NASA/GISS CTM. The model resolution is 4x5 degrees latitude and longitude with 9 layers in the vertical (7 in the troposphere). Advective winds, convection statistics and other meteorological data from the NASA/GISS GCM are used. An extensive gas-phase chemical scheme based on the scheme used in our global 2D model has been incorporated in the 3D model. In this work ozone formation in the troposphere is studied with the 3D model during a 5 day period starting June 30. Extensive local ozone production is found and the relationship between the source regions and the downwind areas are discussed. Variations in local ozone formation as a function of total emission rate, as well as the composition of the emissions (HC/NO(x)) ratio and isoprene emissions) are elucidated. An important vertical transport process in the troposphere is by convective clouds. The 3D model includes an explicit parameterization of this process. It is shown that this process has significant influence on the calculated surface ozone concentrations.

Berntsen, Terje↗
NASA NTRS1994

Quasi-biennial oscillation and tropical waves in total ozone

Westward and eastward propagating tropical waves in total ozone are investigated in 13 years (1979-1991) of version 6 total column ozone data from the Nimbus 7 total ozone mapping spectrometer (TOMS) satellite instrument. A clear synchronization between the stratospheric quasi-biennial osciallation (QBO) zonal winds and the fast (periods less than 15 days) propagating waves in tropical TOMS data is detailed. Largest total ozone wave amplitudes (about 3-6 Dobson units) occur when their phase propagation direction is primarily opposite the Singapore QBO lower-stratospheric winds. This effect is most apparent in meridionally symmetric components. Examination of specific episodes, including cross-spectral calculations with Singapore rawinsonde wind data (10-70 hPa), reveals signatures of tropically confined eastward propagating Kelvin waves of zonal wavenumbers 1-2 during the descending eastward QBO phase, consistent with acceleration of that QBO phase by Kelvin waves. The TOMS results are also consistent with possible forcing of the westward QBO wind phase by episodes of both meridionally symmetric and asymmetric westward waves. However, in contrast to the case of eastward (Kelvin) waves the strongest westward events appear to be filtered by, rather than forcing, the westward phase of the stratospheric QBO wind. These dominant westward episodes are interpreted as meridionally symmetric westward global normal modes and tropically confined equatorial-Rossby waves 2-6. The events exhibit phase and group speeds characteristic of wave dynamics rather than simple wind advection. These results underscore the utility of the long time series and excellent horizontal coverage of TOMS data for dynamical investigations in the relatively observation-poor tropical stratosphere.

Ziemke, J. R.↗
NASA NTRS1994

Incorporation of an evaporative cooling scheme into a dynamic model of orographic precipitation

A simple evaporative cooling scheme was incorporated into a dynamic model to estimate orographic precipitation in mountainous regions. The orographic precipitation model is based on the transport of atmospheric moisture and the quantification of preciptable water across a 3D representation of the terrain from the surface up to 250 hPa. Advective wind fields are computed independently and boundary conditions are extracted from radiosonde data. Precipitation rates are obtained through calibration of a spatially distributed precipitation efficiency parameter. The model was applied to the central Sierra Nevada. Results show a gain of the order of 20% in threat-score coefficients designed to measure the forecast ability of the model. Accuracy gains are largest at high elevations and during intense storms associated with warm air masses.

Barros, Ana Paula↗
DOE OSTI2021

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part I: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LESs. We derive these quantities from measurements and reanalysis obtained for two cold-air outbreak (CAO) events during Phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) in February–March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from ERA5 data agree reasonably well with those derived from ACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time-evolving heat fluxes, wind, and advective tendencies profiles from ERA5 data to drive the LES. We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LESs. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path simulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol–cloud–meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Large-Eddy Simulations of Marine Boundary-Layer Clouds Associated With Cold Air Outbreaks During the Activate Campaign – Part 1: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary-layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LES simulations. We derive these quantities from measurements and reanalysis obtained for two cold air outbreak (CAO) events during Phase I of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) in February-March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from the reanalysis data ERA5 agree reasonablywell with those derived fromACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time evolving heat fluxes, wind and advective tendencies profiles from ERA5 reanalysis data to drive the LES.We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 reanalysis data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LES simulations. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path imulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol-cloud-meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

Large-eddy simulation↗
NASA NTRS2001

The DROPPS/MIDAS Campaign Neutral Atmosphere Measurements and the Occurrence of PMSE and NLC

Measurements of the neutral atmosphere and their relationship to electrodynamic conditions in the mesosphere have been of interest for many years. Inflatable falling sphere measurements along with electrodynamic measurements were obtained in conjunction with the occurrence of PMSE and NLC during the DROPPS/MIDAS Campaign conducted in July 1999 from Andenes Rocket Range, Norway. The inflatable failing sphere measurements in conjunction with a PMSE event on 5-6 July and with a NLC event on 14 July are used to infer thermal advection and its influence on the clouds' maintenance. Hodograph analysis, an early tropospheric tool used by analyst and forecasters, will be used to determine the magnitude and direction of thermal advection from measured wind data. Analysis of the wind structure through the use of hodographs and some assumptions can determine thermal advection, wind shear, and possible vertical motion. Changes in the temperature structure between allied observations were subtle which may be explained by advection. Because of meteorological instabilities in the mesosphere it is possible that hodograph analysis may not fully work. It is our intention to show that such analysis has value and has a place in the mesosphere.

Schmidlin, F. J.↗
NASA NTRS2014

Sea Ice Concentration and Extent

Among the most seasonal and most dynamic parameters on the surface of the Earth is sea ice which at any one time covers about 3-6% of the planet. In the Northern Hemisphere, sea ice grows in extent from about 6 x 10(exp 6) sq km to 16 x 10(exp 6) sq km, while in the Southern Hemisphere, it grows from about 3 x 10(exp 6) sq km to about 19 x 10(exp 6) sq km (Comiso, 2010; Gloersen et al., 1992). Sea ice is up to about 2-3 m thick in the Northern Hemisphere and about 1 m thick in the Southern Hemisphere (Wadhams, 2002), and compared to the average ocean depth of about 3 km, it is a relatively thin, fragile sheet that can break due to waves and winds or melt due to upwelling of warm water. Being constantly advected by winds, waves, and currents, sea ice is very dynamic and usually follows the directions of the many gyres in the polar regions. Despite its vast expanse, the sea ice cover was previously left largely unstudied and it was only in recent years that we have understood its true impact and significance as related to the Earths climate, the oceans, and marine life.

sea ice extent↗
NASA NTRS2018

The NASA Eulerian Snow on Sea Ice Model (NESOSIM) v1.0: initial model development and analysis

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a new, open-source snow budget model that is currentlyconfigured to produce daily estimates of the depth and density of snow on sea ice across the Arctic Ocean throughthe accumulation season. NESOSIM has been developed in a three-dimensional Eulerian framework and includes two(vertical) snow layers and several simple parameterizations (accumulation, wind packing, advection_divergence, blowingsnow lost to leads) to represent key sources and sinks of snow on sea ice. The model is forced with daily inputs of snowfall and near-surface winds (from reanalyses), sea ice concentration (from satellite passive microwave data) and sea ice drift (from satellite feature tracking) during the accumulation season (August through April). In this study, we present the NESOSIM formulation, calibration efforts, sensitivity studies and validation efforts across an Arctic Ocean domain (100 km horizontal resolution). The simulated snow depth and density are calibrated with in situ data collected on drifting ice stations during the 1980s. NESOSIM shows strong agreement with the in situ seasonal cycles of snow depth and density, and shows good (moderate) agreement with the regional snow depth (density) distributions. NESOSIM is run for a contemporary period (2000 to 2015), with the results showing strong sensitivity to the reanalysis derived snowfall forcing data, with the Modern-Era Retrospective analysis for Research and Applications (MERRA) and the Japanese Meteorological Agency 55-year reanalysis (JRA-55) forced snow depths generally higher than ERA-Interim,and the Arctic System Reanalysis (ASR) generally lower. We also generate and force NESOSIM with a consensus "median" daily snowfall dataset from these reanalyses.The results are compared against snow depth estimates derived from NASA's Operation IceBridge (OIB) snow radar data from 2009 to 2015, showing moderate_strong correlations and root mean squared errors of 10 cm depending on the OIB snow depth product analyzed, similar to the comparisons between OIB snow depths and the commonly used modified Warren snow depth climatology. Potential improvements to this initial NESOSIM formulation are discussed in the hopes of improving the accuracy and reliability of these simulated snow depths and densities.

Petty, Alek A.↗