Slivinski
Laura
Slivinski
Laura
No Thumbnail Available
Search Results
Now showing
1 - 4 of 4
-
ArticleUncertainties in ocean latent heat flux variations over recent decades in satellite-based estimates and reduced observation reanalyses(American Meteorological Society, 2020-08-31) Robertson, Franklin R. ; Roberts, Jason B. ; Bosilovich, Michael G. ; Bentamy, Abderrahim ; Clayson, Carol A. ; Fennig, Karsten ; Schröder, Marc ; Tomita, Hiroyuki ; Compo, Gilbert P. ; Gutenstein, Marloes ; Hersbach, Hans ; Kobayashi, Chiaki ; Ricciardulli, Lucrezia ; Sardeshmukh, Prashant ; Slivinski, LauraFour state-of-the-art satellite-based estimates of ocean surface latent heat fluxes (LHFs) extending over three decades are analyzed, focusing on the interannual variability and trends of near-global averages and regional patterns. Detailed intercomparisons are made with other datasets including 1) reduced observation reanalyses (RedObs) whose exclusion of satellite data renders them an important independent diagnostic tool; 2) a moisture budget residual LHF estimate using reanalysis moisture transport, atmospheric storage, and satellite precipitation; 3) the ECMWF Reanalysis 5 (ERA5); 4) Remote Sensing Systems (RSS) single-sensor passive microwave and scatterometer wind speed retrievals; and 5) several sea surface temperature (SST) datasets. Large disparities remain in near-global satellite LHF trends and their regional expression over the 1990–2010 period, during which time the interdecadal Pacific oscillation changed sign. The budget residual diagnostics support the smaller RedObs LHF trends. The satellites, ERA5, and RedObs are reasonably consistent in identifying contributions by the 10-m wind speed variations to the LHF trend patterns. However, contributions by the near-surface vertical humidity gradient from satellites and ERA5 trend upward in time with respect to the RedObs ensemble and show less agreement in trend patterns. Problems with wind speed retrievals from Special Sensor Microwave Imager/Sounder satellite sensors, excessive upward trends in trends in Optimal Interpolation Sea Surface Temperature (OISST AVHRR-Only) data used in most satellite LHF estimates, and uncertainties associated with poor satellite coverage before the mid-1990s are noted. Possibly erroneous trends are also identified in ERA5 LHF associated with the onset of scatterometer wind data assimilation in the early 1990s.
-
ArticleA hybrid particle–ensemble Kalman filter for Lagrangian data assimilation(American Meteorological Society, 2015-01) Slivinski, Laura ; Spiller, Elaine ; Apte, Amit ; Sandstede, BjornLagrangian measurements from passive ocean instruments provide a useful source of data for estimating and forecasting the ocean’s state (velocity field, salinity field, etc.). However, trajectories from these instruments are often highly nonlinear, leading to difficulties with widely used data assimilation algorithms such as the ensemble Kalman filter (EnKF). Additionally, the velocity field is often modeled as a high-dimensional variable, which precludes the use of more accurate methods such as the particle filter (PF). Here, a hybrid particle–ensemble Kalman filter is developed that applies the EnKF update to the potentially high-dimensional velocity variables, and the PF update to the relatively low-dimensional, highly nonlinear drifter position variable. This algorithm is tested with twin experiments on the linear shallow water equations. In experiments with infrequent observations, the hybrid filter consistently outperformed the EnKF, both by better capturing the Bayesian posterior and by better tracking the truth.
-
ArticleExploring practical estimates of the ensemble size necessary for particle filters(American Meteorological Society, 2015-11-11) Slivinski, Laura ; Snyder, ChrisParticle filtering methods for data assimilation may suffer from the “curse of dimensionality,” where the required ensemble size grows rapidly as the dimension increases. It would, therefore, be useful to know a priori whether a particle filter is feasible to implement in a given system. Previous work provides an asymptotic relation between the necessary ensemble size and an exponential function of , a statistic that depends on observation-space quantities and that is related to the system dimension when the number of observations is large; for linear, Gaussian systems, the statistic can be computed from eigenvalues of an appropriately normalized covariance matrix. Tests with a low-dimensional system show that these asymptotic results remain useful when the system is nonlinear, with either the standard or optimal proposal implementation of the particle filter. This study explores approximations to the covariance matrices that facilitate computation in high-dimensional systems, as well as different methods to estimate the accumulated system noise covariance for the optimal proposal. Since may be approximated using an ensemble from a simpler data assimilation scheme, such as the ensemble Kalman filter, the asymptotic relations thus allow an estimate of the ensemble size required for a particle filter before its implementation. Finally, the improved performance of particle filters with the optimal proposal, relative to those using the standard proposal, in the same low-dimensional system is demonstrated.
-
PreprintAssimilating Lagrangian data for parameter estimation in a multiple-inlet system( 2017-04) Slivinski, Laura ; Pratt, Lawrence J. ; Rypina, Irina I. ; Orescanin, Mara M. ; Raubenheimer, Britt ; MacMahan, Jamie ; Elgar, SteveNumerical models of ocean circulation often depend on parameters that must be tuned to match either results from laboratory experiments or field observations. This study demonstrates that an initial, suboptimal estimate of a parameter in a model of a small bay can be improved by assimilating observations of trajectories of passive drifters. The parameter of interest is the Manning's n coefficient of friction in a small inlet of the bay, which had been tuned to match velocity observations from 2011. In 2013, the geometry of the inlet had changed, and the friction parameter was no longer optimal. Results from synthetic experiments demonstrate that assimilation of drifter trajectories improves the estimate of n, both when the drifters are located in the same region as the parameter of interest and when the drifters are located in a different region of the bay. Real drifter trajectories from field experiments in 2013 also are assimilated, and results are compared with velocity observations. When the real drifters are located away from the region of interest, the results depend on the time interval (with respect to the full available trajectories) over which assimilation is performed. When the drifters are in the same region as the parameter of interest, the value of n estimated with assimilation yields improved estimates of velocity throughout the bay.