Search Results
You are looking at 11 - 20 of 24 items for
- Author or Editor: Ibrahim Hoteit x
- Refine by Access: All Content x
Abstract
The city of Jeddah, Saudi Arabia, is characterized by a hot and arid desert climate. On occasion, however, extreme precipitation events have led to flooding that caused extensive damage to human life and infrastructure. This study investigates the effect of incorporating an urban canopy model and urban land cover when simulating severe weather events over Jeddah using the Weather Research and Forecasting (WRF) Model at a convective-permitting scale (1.5-km resolution). Two experiments were conducted for 10 heavy rainfall events associated with the dominant large-scale patterns favoring convection over Jeddah: (i) an “urban” experiment that included the urban canopy model and modern-day land cover and (ii) a “desert” experiment that replaced the city area with its presettlement, natural land cover. The results suggest that urbanization plays an important role in modifying rainfall around city area. The urban experiment enhances the amount of rainfall by 26% on average over the Jeddah city area relative to the desert experiment in these extreme events. The changes in model-simulated precipitation are primarily tied to a nocturnal heat-island effect that modifies the planetary boundary layer and atmospheric instability of the convective events.
Abstract
The city of Jeddah, Saudi Arabia, is characterized by a hot and arid desert climate. On occasion, however, extreme precipitation events have led to flooding that caused extensive damage to human life and infrastructure. This study investigates the effect of incorporating an urban canopy model and urban land cover when simulating severe weather events over Jeddah using the Weather Research and Forecasting (WRF) Model at a convective-permitting scale (1.5-km resolution). Two experiments were conducted for 10 heavy rainfall events associated with the dominant large-scale patterns favoring convection over Jeddah: (i) an “urban” experiment that included the urban canopy model and modern-day land cover and (ii) a “desert” experiment that replaced the city area with its presettlement, natural land cover. The results suggest that urbanization plays an important role in modifying rainfall around city area. The urban experiment enhances the amount of rainfall by 26% on average over the Jeddah city area relative to the desert experiment in these extreme events. The changes in model-simulated precipitation are primarily tied to a nocturnal heat-island effect that modifies the planetary boundary layer and atmospheric instability of the convective events.
Abstract
This work combines two auxiliary techniques, namely the one-step-ahead (OSA) smoothing and the hybrid formulation, to boost the forecasting skills of a storm surge ensemble Kalman filter (EnKF) forecasting system. Bayesian filtering with OSA-smoothing enhances the robustness of the ensemble background statistics by exploiting the data twice: first to constrain the sampling of the forecast ensemble with the future observation, and then to update the resulting ensemble. This is expected to improve the behavior of EnKF-like schemes during the strongly nonlinear surges periods, but requires integrating the ensemble with the forecast model twice, which could be computationally demanding. The hybrid flow-dependent/static formulation of the EnKF background error covariance is then considered to enable the implementation of the filter with a small flow-dependent ensemble size, and thus less model runs. These two methods are combined within an ensemble transform Kalman filter (ETKF). The resulting hybrid ETKF with OSA smoothing is tested, based on twin experiments, using a realistic setting of the Advanced Circulation (ADCIRC) model configured for storm surge forecasting in the Gulf of Mexico and assimilating pseudo-observations of sea surface levels from a network of buoys. The results of our numerical experiments suggest that the proposed filtering system significantly enhances ADCIRC forecasting skills compared to the standard ETKF without increasing the computational cost.
Abstract
This work combines two auxiliary techniques, namely the one-step-ahead (OSA) smoothing and the hybrid formulation, to boost the forecasting skills of a storm surge ensemble Kalman filter (EnKF) forecasting system. Bayesian filtering with OSA-smoothing enhances the robustness of the ensemble background statistics by exploiting the data twice: first to constrain the sampling of the forecast ensemble with the future observation, and then to update the resulting ensemble. This is expected to improve the behavior of EnKF-like schemes during the strongly nonlinear surges periods, but requires integrating the ensemble with the forecast model twice, which could be computationally demanding. The hybrid flow-dependent/static formulation of the EnKF background error covariance is then considered to enable the implementation of the filter with a small flow-dependent ensemble size, and thus less model runs. These two methods are combined within an ensemble transform Kalman filter (ETKF). The resulting hybrid ETKF with OSA smoothing is tested, based on twin experiments, using a realistic setting of the Advanced Circulation (ADCIRC) model configured for storm surge forecasting in the Gulf of Mexico and assimilating pseudo-observations of sea surface levels from a network of buoys. The results of our numerical experiments suggest that the proposed filtering system significantly enhances ADCIRC forecasting skills compared to the standard ETKF without increasing the computational cost.
Abstract
The ensemble Kalman filter (EnKF) is widely used for sequential data assimilation. It operates as a succession of forecast and analysis steps. In realistic large-scale applications, EnKFs are implemented with small ensembles and poorly known model error statistics. This limits their representativeness of the background error covariances and, thus, their performance. This work explores the efficiency of the one-step-ahead (OSA) smoothing formulation of the Bayesian filtering problem to enhance the data assimilation performance of EnKFs. Filtering with OSA smoothing introduces an updated step with future observations, conditioning the ensemble sampling with more information. This should provide an improved background ensemble in the analysis step, which may help to mitigate the suboptimal character of EnKF-based methods. Here, the authors demonstrate the efficiency of a stochastic EnKF with OSA smoothing for state estimation. They then introduce a deterministic-like EnKF-OSA based on the singular evolutive interpolated ensemble Kalman (SEIK) filter. The authors show that the proposed SEIK-OSA outperforms both SEIK, as it efficiently exploits the data twice, and the stochastic EnKF-OSA, as it avoids observational error undersampling. They present extensive assimilation results from numerical experiments conducted with the Lorenz-96 model to demonstrate SEIK-OSA’s capabilities.
Abstract
The ensemble Kalman filter (EnKF) is widely used for sequential data assimilation. It operates as a succession of forecast and analysis steps. In realistic large-scale applications, EnKFs are implemented with small ensembles and poorly known model error statistics. This limits their representativeness of the background error covariances and, thus, their performance. This work explores the efficiency of the one-step-ahead (OSA) smoothing formulation of the Bayesian filtering problem to enhance the data assimilation performance of EnKFs. Filtering with OSA smoothing introduces an updated step with future observations, conditioning the ensemble sampling with more information. This should provide an improved background ensemble in the analysis step, which may help to mitigate the suboptimal character of EnKF-based methods. Here, the authors demonstrate the efficiency of a stochastic EnKF with OSA smoothing for state estimation. They then introduce a deterministic-like EnKF-OSA based on the singular evolutive interpolated ensemble Kalman (SEIK) filter. The authors show that the proposed SEIK-OSA outperforms both SEIK, as it efficiently exploits the data twice, and the stochastic EnKF-OSA, as it avoids observational error undersampling. They present extensive assimilation results from numerical experiments conducted with the Lorenz-96 model to demonstrate SEIK-OSA’s capabilities.
Abstract
A new hybrid ensemble Kalman filter/four-dimensional variational data assimilation (EnKF/4D-VAR) approach is introduced to mitigate background covariance limitations in the EnKF. The work is based on the adaptive EnKF (AEnKF) method, which bears a strong resemblance to the hybrid EnKF/three-dimensional variational data assimilation (3D-VAR) method. In the AEnKF, the representativeness of the EnKF ensemble is regularly enhanced with new members generated after back projection of the EnKF analysis residuals to state space using a 3D-VAR [or optimal interpolation (OI)] scheme with a preselected background covariance matrix. The idea here is to reformulate the transformation of the residuals as a 4D-VAR problem, constraining the new member with model dynamics and the previous observations. This should provide more information for the estimation of the new member and reduce dependence of the AEnKF on the assumed stationary background covariance matrix. This is done by integrating the analysis residuals backward in time with the adjoint model. Numerical experiments are performed with the Lorenz-96 model under different scenarios to test the new approach and to evaluate its performance with respect to the EnKF and the hybrid EnKF/3D-VAR. The new method leads to the least root-mean-square estimation errors as long as the linear assumption guaranteeing the stability of the adjoint model holds. It is also found to be less sensitive to choices of the assimilation system inputs and parameters.
Abstract
A new hybrid ensemble Kalman filter/four-dimensional variational data assimilation (EnKF/4D-VAR) approach is introduced to mitigate background covariance limitations in the EnKF. The work is based on the adaptive EnKF (AEnKF) method, which bears a strong resemblance to the hybrid EnKF/three-dimensional variational data assimilation (3D-VAR) method. In the AEnKF, the representativeness of the EnKF ensemble is regularly enhanced with new members generated after back projection of the EnKF analysis residuals to state space using a 3D-VAR [or optimal interpolation (OI)] scheme with a preselected background covariance matrix. The idea here is to reformulate the transformation of the residuals as a 4D-VAR problem, constraining the new member with model dynamics and the previous observations. This should provide more information for the estimation of the new member and reduce dependence of the AEnKF on the assumed stationary background covariance matrix. This is done by integrating the analysis residuals backward in time with the adjoint model. Numerical experiments are performed with the Lorenz-96 model under different scenarios to test the new approach and to evaluate its performance with respect to the EnKF and the hybrid EnKF/3D-VAR. The new method leads to the least root-mean-square estimation errors as long as the linear assumption guaranteeing the stability of the adjoint model holds. It is also found to be less sensitive to choices of the assimilation system inputs and parameters.
Abstract
The authors present a polynomial chaos (PC)–based Bayesian inference method for quantifying the uncertainties of the K-profile parameterization (KPP) within the MIT general circulation model (MITgcm) of the tropical Pacific. The inference of the uncertain parameters is based on a Markov chain Monte Carlo (MCMC) scheme that utilizes a newly formulated test statistic taking into account the different components representing the structures of turbulent mixing on both daily and seasonal time scales in addition to the data quality, and filters for the effects of parameter perturbations over those as a result of changes in the wind. To avoid the prohibitive computational cost of integrating the MITgcm model at each MCMC iteration, a surrogate model for the test statistic using the PC method is built. Because of the noise in the model predictions, a basis-pursuit-denoising (BPDN) compressed sensing approach is employed to determine the PC coefficients of a representative surrogate model. The PC surrogate is then used to evaluate the test statistic in the MCMC step for sampling the posterior of the uncertain parameters. Results of the posteriors indicate good agreement with the default values for two parameters of the KPP model, namely the critical bulk and gradient Richardson numbers; while the posteriors of the remaining parameters were barely informative.
Abstract
The authors present a polynomial chaos (PC)–based Bayesian inference method for quantifying the uncertainties of the K-profile parameterization (KPP) within the MIT general circulation model (MITgcm) of the tropical Pacific. The inference of the uncertain parameters is based on a Markov chain Monte Carlo (MCMC) scheme that utilizes a newly formulated test statistic taking into account the different components representing the structures of turbulent mixing on both daily and seasonal time scales in addition to the data quality, and filters for the effects of parameter perturbations over those as a result of changes in the wind. To avoid the prohibitive computational cost of integrating the MITgcm model at each MCMC iteration, a surrogate model for the test statistic using the PC method is built. Because of the noise in the model predictions, a basis-pursuit-denoising (BPDN) compressed sensing approach is employed to determine the PC coefficients of a representative surrogate model. The PC surrogate is then used to evaluate the test statistic in the MCMC step for sampling the posterior of the uncertain parameters. Results of the posteriors indicate good agreement with the default values for two parameters of the KPP model, namely the critical bulk and gradient Richardson numbers; while the posteriors of the remaining parameters were barely informative.
Abstract
Backward-in-time predictions are needed to better understand the underlying dynamics of physical fluid flows and improve future forecasts. However, integrating fluid flows backward in time is challenging because of numerical instabilities caused by the diffusive nature of the fluid systems and nonlinearities of the governing equations. Although this problem has been long addressed using a nonpositive diffusion coefficient when integrating backward, it is notoriously inaccurate. In this study, a physics-informed deep neural network (PI-DNN) is presented to predict past states of a dissipative dynamical system from snapshots of the subsequent evolution of the system state. The performance of the PI-DNN is investigated using several systematic numerical experiments and the accuracy of the backward-in-time predictions is evaluated in terms of different error metrics. The proposed PI-DNN can predict the previous state of the Rayleigh–Bénard convection with an 8-time-step average normalized
Abstract
Backward-in-time predictions are needed to better understand the underlying dynamics of physical fluid flows and improve future forecasts. However, integrating fluid flows backward in time is challenging because of numerical instabilities caused by the diffusive nature of the fluid systems and nonlinearities of the governing equations. Although this problem has been long addressed using a nonpositive diffusion coefficient when integrating backward, it is notoriously inaccurate. In this study, a physics-informed deep neural network (PI-DNN) is presented to predict past states of a dissipative dynamical system from snapshots of the subsequent evolution of the system state. The performance of the PI-DNN is investigated using several systematic numerical experiments and the accuracy of the backward-in-time predictions is evaluated in terms of different error metrics. The proposed PI-DNN can predict the previous state of the Rayleigh–Bénard convection with an 8-time-step average normalized
Abstract
This paper investigates the role of the linear analysis step of the ensemble Kalman filters (EnKF) in disrupting the balanced dynamics in a simple atmospheric model and compares it to a fully nonlinear particle-based filter (PF). The filters have a very similar forecast step but the analysis step of the PF solves the full Bayesian filtering problem while the EnKF analysis only applies to Gaussian distributions. The EnKF is compared to two flavors of the particle filter with different sampling strategies, the sequential importance resampling filter (SIRF) and the sequential kernel resampling filter (SKRF). The model admits a chaotic vortical mode coupled to a comparatively fast gravity wave mode. It can also be configured either to evolve on a so-called slow manifold, where the fast motion is suppressed, or such that the fast-varying variables are diagnosed from the slow-varying variables as slaved modes. Identical twin experiments show that EnKF and PF capture the variables on the slow manifold well as the dynamics is very stable. PFs, especially the SKRF, capture slaved modes better than the EnKF, implying that a full Bayesian analysis estimates the nonlinear model variables better. The PFs perform significantly better in the fully coupled nonlinear model where fast and slow variables modulate each other. This suggests that the analysis step in the PFs maintains the balance in both variables much better than the EnKF. It is also shown that increasing the ensemble size generally improves the performance of the PFs but has less impact on the EnKF after a sufficient number of members have been used.
Abstract
This paper investigates the role of the linear analysis step of the ensemble Kalman filters (EnKF) in disrupting the balanced dynamics in a simple atmospheric model and compares it to a fully nonlinear particle-based filter (PF). The filters have a very similar forecast step but the analysis step of the PF solves the full Bayesian filtering problem while the EnKF analysis only applies to Gaussian distributions. The EnKF is compared to two flavors of the particle filter with different sampling strategies, the sequential importance resampling filter (SIRF) and the sequential kernel resampling filter (SKRF). The model admits a chaotic vortical mode coupled to a comparatively fast gravity wave mode. It can also be configured either to evolve on a so-called slow manifold, where the fast motion is suppressed, or such that the fast-varying variables are diagnosed from the slow-varying variables as slaved modes. Identical twin experiments show that EnKF and PF capture the variables on the slow manifold well as the dynamics is very stable. PFs, especially the SKRF, capture slaved modes better than the EnKF, implying that a full Bayesian analysis estimates the nonlinear model variables better. The PFs perform significantly better in the fully coupled nonlinear model where fast and slow variables modulate each other. This suggests that the analysis step in the PFs maintains the balance in both variables much better than the EnKF. It is also shown that increasing the ensemble size generally improves the performance of the PFs but has less impact on the EnKF after a sufficient number of members have been used.
Abstract
The wind and wave climatology of the Red Sea is derived from a validated 30-yr high-resolution model simulation. After describing the relevant features of the basin, the main wind and wave systems are identified by using an innovative spectral partition technique to explain their genesis and characteristics. In the northern part of the sea, wind and waves of the same intensity are present throughout the year, while the central and southern zones are characterized by a marked seasonality. The partition technique allows the association of a general decrease in the energy of the different wave systems with a specific weather pattern. The most intense decrease is found in the northern storms, which are associated with meteorological pulses from the Mediterranean Sea.
Abstract
The wind and wave climatology of the Red Sea is derived from a validated 30-yr high-resolution model simulation. After describing the relevant features of the basin, the main wind and wave systems are identified by using an innovative spectral partition technique to explain their genesis and characteristics. In the northern part of the sea, wind and waves of the same intensity are present throughout the year, while the central and southern zones are characterized by a marked seasonality. The partition technique allows the association of a general decrease in the energy of the different wave systems with a specific weather pattern. The most intense decrease is found in the northern storms, which are associated with meteorological pulses from the Mediterranean Sea.
Abstract
The impacts of various climate modes on the Red Sea surface heat exchange are investigated using the MERRA reanalysis and the OAFlux satellite reanalysis datasets. Seasonality in the atmospheric forcing is also explored. Mode impacts peak during boreal winter [December–February (DJF)] with average anomalies of 12–18 W m−2 to be found in the northern Red Sea. The North Atlantic Oscillation (NAO), the east Atlantic–west Russia (EAWR) pattern, and the Indian monsoon index (IMI) exhibit the strongest influence on the air–sea heat exchange during the winter. In this season, the largest negative anomalies of about −30 W m−2 are associated with the EAWR pattern over the central part of the Red Sea. In other seasons, mode-related anomalies are considerably lower, especially during spring when the mode impacts are negligible. The mode impacts are strongest over the northern half of the Red Sea during winter and autumn. In summer, the southern half of the basin is strongly influenced by the multivariate ENSO index (MEI). The winter mode–related anomalies are determined mostly by the latent heat flux component, while in summer the shortwave flux is also important. The influence of the modes on the Red Sea is found to be generally weaker than on the neighboring Mediterranean basin.
Abstract
The impacts of various climate modes on the Red Sea surface heat exchange are investigated using the MERRA reanalysis and the OAFlux satellite reanalysis datasets. Seasonality in the atmospheric forcing is also explored. Mode impacts peak during boreal winter [December–February (DJF)] with average anomalies of 12–18 W m−2 to be found in the northern Red Sea. The North Atlantic Oscillation (NAO), the east Atlantic–west Russia (EAWR) pattern, and the Indian monsoon index (IMI) exhibit the strongest influence on the air–sea heat exchange during the winter. In this season, the largest negative anomalies of about −30 W m−2 are associated with the EAWR pattern over the central part of the Red Sea. In other seasons, mode-related anomalies are considerably lower, especially during spring when the mode impacts are negligible. The mode impacts are strongest over the northern half of the Red Sea during winter and autumn. In summer, the southern half of the basin is strongly influenced by the multivariate ENSO index (MEI). The winter mode–related anomalies are determined mostly by the latent heat flux component, while in summer the shortwave flux is also important. The influence of the modes on the Red Sea is found to be generally weaker than on the neighboring Mediterranean basin.
Abstract
The influence of the atmospheric circulation on the winter air–sea heat fluxes over the northern Red Sea is investigated during the period 1985–2011. The analysis based on daily heat flux values reveals that most of the net surface heat exchange variability depends on the behavior of the turbulent components of the surface flux (the sum of the latent and sensible heat). The large-scale composite sea level pressure (SLP) maps corresponding to turbulent flux minima and maxima show distinct atmospheric circulation patterns associated with each case. In general, extreme heat loss (with turbulent flux lower than −400 W m−2) over the northern Red Sea is observed when anticyclonic conditions prevail over an area extending from the Mediterranean Sea to eastern Asia along with a recession of the equatorial African lows system. Subcenters of high pressure associated with this pattern generate the required steep SLP gradient that enhances the wind magnitude and transfers cold and dry air masses from higher latitudes. Conversely, turbulent flux maxima (heat loss minimization with values from −100 to −50 W m−2) are associated with prevailing low pressures over the eastern Mediterranean and an extended equatorial African low that reaches the southern part of the Red Sea. In this case, a smooth SLP field over the northern Red Sea results in weak winds over the area that in turn reduce the surface heat loss. At the same time, southerlies blowing along the main axis of the Red Sea transfer warm and humid air northward, favoring heat flux maxima.
Abstract
The influence of the atmospheric circulation on the winter air–sea heat fluxes over the northern Red Sea is investigated during the period 1985–2011. The analysis based on daily heat flux values reveals that most of the net surface heat exchange variability depends on the behavior of the turbulent components of the surface flux (the sum of the latent and sensible heat). The large-scale composite sea level pressure (SLP) maps corresponding to turbulent flux minima and maxima show distinct atmospheric circulation patterns associated with each case. In general, extreme heat loss (with turbulent flux lower than −400 W m−2) over the northern Red Sea is observed when anticyclonic conditions prevail over an area extending from the Mediterranean Sea to eastern Asia along with a recession of the equatorial African lows system. Subcenters of high pressure associated with this pattern generate the required steep SLP gradient that enhances the wind magnitude and transfers cold and dry air masses from higher latitudes. Conversely, turbulent flux maxima (heat loss minimization with values from −100 to −50 W m−2) are associated with prevailing low pressures over the eastern Mediterranean and an extended equatorial African low that reaches the southern part of the Red Sea. In this case, a smooth SLP field over the northern Red Sea results in weak winds over the area that in turn reduce the surface heat loss. At the same time, southerlies blowing along the main axis of the Red Sea transfer warm and humid air northward, favoring heat flux maxima.
