Search Results
You are looking at 1 - 10 of 14 items for
- Author or Editor: Kevin Raeder x
- Refine by Access: All Content x
Abstract
The authors propose a new procedure. designated the adjoint-based genesis diagnostic (AGD) procedure, for studying triggering mechanism and the subsequent genesis of the synoptic phenomena of interest. This procedure makes use of a numerical model sensitivity to initial conditions and the nonlinear evolution of the initial perturbations that are designed using this sensitivity. The model sensitivity is evaluated using the associated adjoint model. This study uses the dry version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System (MAMS) for the numerical experiments. The authors apply the AGD procedure to two cases of Alpine lee cyclogenesis that were observed during the Alpine Experiment special observation period. The results show that the sensitivity fields that are produced by the adjoint model and the associated initial perturbations are readily related to the probable triggering mechanisms for these cyclones. Additionally, the nonlinear evolution of these initial perturbations points toward the physical processes involved in the lee cyclone formation. The AGD experiments for a weak cyclone case indicate that the MAMS forecast model has an underrepresented topographic forcing due to the sigma vertical coordinate and that this model error can be compensated by adjustments in the initial conditions that are related to the triggering mechanism, which is not associated with the topographic blocking mechanism.
Abstract
The authors propose a new procedure. designated the adjoint-based genesis diagnostic (AGD) procedure, for studying triggering mechanism and the subsequent genesis of the synoptic phenomena of interest. This procedure makes use of a numerical model sensitivity to initial conditions and the nonlinear evolution of the initial perturbations that are designed using this sensitivity. The model sensitivity is evaluated using the associated adjoint model. This study uses the dry version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System (MAMS) for the numerical experiments. The authors apply the AGD procedure to two cases of Alpine lee cyclogenesis that were observed during the Alpine Experiment special observation period. The results show that the sensitivity fields that are produced by the adjoint model and the associated initial perturbations are readily related to the probable triggering mechanisms for these cyclones. Additionally, the nonlinear evolution of these initial perturbations points toward the physical processes involved in the lee cyclone formation. The AGD experiments for a weak cyclone case indicate that the MAMS forecast model has an underrepresented topographic forcing due to the sigma vertical coordinate and that this model error can be compensated by adjustments in the initial conditions that are related to the triggering mechanism, which is not associated with the topographic blocking mechanism.
Abstract
No abstract available
Abstract
No abstract available
Abstract
The importance of multivariate forecast error correlations between specific humidity, temperature, and surface pressure in perfect model assimilations of Global Positioning System radio occultation (RO) refractivity data is examined using the Ensemble Adjustment Filter (EAF) and the NCAR global Community Atmospheric Model, version 3. The goal is to explore whether inclusion of the multivariate forecast error correlations in the background term of 3D and 4D variational data assimilation systems (3DVAR and 4DVAR, respectively) is likely to improve RO data assimilation in the troposphere. It is not possible to explicitly neglect multivariate forecast error correlations with the EAF because they are not used directly in the algorithm. Instead, the filter only makes use of the forecast error correlations between observed quantities (RO here) and model state variables. However, because the forecast error correlations for RO observations are dominated by correlations with a subset of state variable types in certain regions, the importance of multivariate forecast error correlations between state variables can be indirectly assessed. This is done by setting the forecast error correlations of RO observations and some state variables (e.g., temperature) to zero in a set of assimilation experiments. Comparing these experiments to a control in which all state variables are impacted by RO observations allows an indirect assessment of the importance of multivariate correlations between state variables not impacted by the observations and those that are impacted. Results suggest that proper specification of the multivariate forecast error correlations in 3DVAR and 4DVAR systems should improve the analysis of specific humidity, surface pressure, and temperature in the troposphere when assimilating RO data.
Abstract
The importance of multivariate forecast error correlations between specific humidity, temperature, and surface pressure in perfect model assimilations of Global Positioning System radio occultation (RO) refractivity data is examined using the Ensemble Adjustment Filter (EAF) and the NCAR global Community Atmospheric Model, version 3. The goal is to explore whether inclusion of the multivariate forecast error correlations in the background term of 3D and 4D variational data assimilation systems (3DVAR and 4DVAR, respectively) is likely to improve RO data assimilation in the troposphere. It is not possible to explicitly neglect multivariate forecast error correlations with the EAF because they are not used directly in the algorithm. Instead, the filter only makes use of the forecast error correlations between observed quantities (RO here) and model state variables. However, because the forecast error correlations for RO observations are dominated by correlations with a subset of state variable types in certain regions, the importance of multivariate forecast error correlations between state variables can be indirectly assessed. This is done by setting the forecast error correlations of RO observations and some state variables (e.g., temperature) to zero in a set of assimilation experiments. Comparing these experiments to a control in which all state variables are impacted by RO observations allows an indirect assessment of the importance of multivariate correlations between state variables not impacted by the observations and those that are impacted. Results suggest that proper specification of the multivariate forecast error correlations in 3DVAR and 4DVAR systems should improve the analysis of specific humidity, surface pressure, and temperature in the troposphere when assimilating RO data.
Abstract
In this study, ensemble sensitivity analysis has been applied to examine the evolution of two extreme extratropical cyclones over the Pacific. Sensitivity using, as forecast metrics, forecast cyclone minimum pressure and location, as well as principal components (PCs) of the leading EOFs in forecast SLP variations near the cyclone center, has been computed for medium-range forecasts of up to 7.5 days. Results presented here show that coherent sensitivity patterns can be tracked from the forecast validation time back in time to at least day −6, with the sensitivity signal exhibiting downstream development characteristics in most cases. Comparing the different forecast metrics, sensitivity patterns derived from the PCs of the leading EOFs in forecast SLP variations are apparently more coherent than those derived from cyclone parameters.
To test whether the linear sensitivity analyses provide quantitatively accurate guidance under the highly nonlinear evolution of the atmospheric flow, perturbed initial condition experiments have been conducted using initial condition perturbations derived based on ensemble sensitivity. Results of this study suggest that in the medium range, perturbations derived from cyclone parameters are quite effective in modifying the evolution of the cyclones out to 5.5 days, but are largely ineffective for 7.5-day forecasts. On the other hand, perturbations derived based on the PCs of the leading EOFs are still quite effective in modifying forecast cyclone location out to 7.5 days. These results suggest that EOF-based sensitivities perform better than cyclone parameter-based sensitivities in the medium range.
Abstract
In this study, ensemble sensitivity analysis has been applied to examine the evolution of two extreme extratropical cyclones over the Pacific. Sensitivity using, as forecast metrics, forecast cyclone minimum pressure and location, as well as principal components (PCs) of the leading EOFs in forecast SLP variations near the cyclone center, has been computed for medium-range forecasts of up to 7.5 days. Results presented here show that coherent sensitivity patterns can be tracked from the forecast validation time back in time to at least day −6, with the sensitivity signal exhibiting downstream development characteristics in most cases. Comparing the different forecast metrics, sensitivity patterns derived from the PCs of the leading EOFs in forecast SLP variations are apparently more coherent than those derived from cyclone parameters.
To test whether the linear sensitivity analyses provide quantitatively accurate guidance under the highly nonlinear evolution of the atmospheric flow, perturbed initial condition experiments have been conducted using initial condition perturbations derived based on ensemble sensitivity. Results of this study suggest that in the medium range, perturbations derived from cyclone parameters are quite effective in modifying the evolution of the cyclones out to 5.5 days, but are largely ineffective for 7.5-day forecasts. On the other hand, perturbations derived based on the PCs of the leading EOFs are still quite effective in modifying forecast cyclone location out to 7.5 days. These results suggest that EOF-based sensitivities perform better than cyclone parameter-based sensitivities in the medium range.
Abstract
Finite-time growth of perturbations in the presence of moist physics (specifically, precipitation) is investigated using singular vectors (SVs) in the context of a primitive equation regional model. Two difficulties appear in the explicit consideration of the effect of moist physics when studying such optimal growth. First, the tangent-linear description of moist physics may not be as straightforward and accurate as for dry-adiabatic processes; second, because of the consideration of moisture, the design of an appropriate measure of growth (i.e., norm) is subject to even more ambiguity than in the dry situation.
In this study both of these problems are addressed in the context of the moist version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System, version 2, with emphasis on the second problem. Leading SVs are computed in an iterative fashion, using a Lanczos algorithm, for three norms over an optimization interval of 24 h; these norms are based on an expression related to (total) perturbation energy. The properties of these SVs are studied for a case of explosive cyclogenesis and a case of summer convection.
The consideration of moisture leads to faster growth of perturbations than in the dry situation, as well as to the appearance of new growing structures. Apparently, moist processes provide for new mechanisms of error growth and do not simply lead to a modulation of SVs obtained with the dry version of the model. Consequently, consideration of the linearized moist processes is essential for revealing all structures that might potentially grow in a moist primitive equation model. In the context of this investigation growth rates depend more on the choice of the basic state and linearized model (moist vs dry) than on the choice of the norm (moist vs dry total energy norm). A reference is cited that supports the validity of the moist tangent-linear SV perturbation growth studied here in the nonlinear regime.
Abstract
Finite-time growth of perturbations in the presence of moist physics (specifically, precipitation) is investigated using singular vectors (SVs) in the context of a primitive equation regional model. Two difficulties appear in the explicit consideration of the effect of moist physics when studying such optimal growth. First, the tangent-linear description of moist physics may not be as straightforward and accurate as for dry-adiabatic processes; second, because of the consideration of moisture, the design of an appropriate measure of growth (i.e., norm) is subject to even more ambiguity than in the dry situation.
In this study both of these problems are addressed in the context of the moist version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System, version 2, with emphasis on the second problem. Leading SVs are computed in an iterative fashion, using a Lanczos algorithm, for three norms over an optimization interval of 24 h; these norms are based on an expression related to (total) perturbation energy. The properties of these SVs are studied for a case of explosive cyclogenesis and a case of summer convection.
The consideration of moisture leads to faster growth of perturbations than in the dry situation, as well as to the appearance of new growing structures. Apparently, moist processes provide for new mechanisms of error growth and do not simply lead to a modulation of SVs obtained with the dry version of the model. Consequently, consideration of the linearized moist processes is essential for revealing all structures that might potentially grow in a moist primitive equation model. In the context of this investigation growth rates depend more on the choice of the basic state and linearized model (moist vs dry) than on the choice of the norm (moist vs dry total energy norm). A reference is cited that supports the validity of the moist tangent-linear SV perturbation growth studied here in the nonlinear regime.
Abstract
Sampling errors and model errors are major drawbacks from which ensemble Kalman filters suffer. Sampling errors arise because of the use of a limited ensemble size, while model errors are deficiencies in the dynamics and underlying parameterizations that may yield biases in the model’s prediction. In this study, we propose a new time-adaptive posterior inflation algorithm in which the analyzed ensemble anomalies are locally inflated. The proposed inflation strategy is computationally efficient and is aimed at restoring enough spread in the analysis ensemble after assimilating the observations. The performance of this scheme is tested against the relaxation to prior spread (RTPS) and adaptive prior inflation. For this purpose, two model are used: the three-variable Lorenz 63 system and the Community Atmosphere Model (CAM). In CAM, global refractivity, temperature, and wind observations from several sources are incorporated to perform a set of assimilation experiments using the Data Assimilation Research Testbed (DART). The proposed scheme is shown to yield better quality forecasts than the RTPS. Assimilation results further suggest that when model errors are small, both prior and posterior inflation are able to mitigate sampling errors with a slight advantage to posterior inflation. When large model errors, such as wind and temperature biases, are present, prior inflation is shown to be more accurate than posterior inflation. Densely observed regions as in the Northern Hemisphere present numerous challenges to the posterior inflation algorithm. A compelling enhancement to the performance of the filter is achieved by combining both adaptive inflation schemes.
Abstract
Sampling errors and model errors are major drawbacks from which ensemble Kalman filters suffer. Sampling errors arise because of the use of a limited ensemble size, while model errors are deficiencies in the dynamics and underlying parameterizations that may yield biases in the model’s prediction. In this study, we propose a new time-adaptive posterior inflation algorithm in which the analyzed ensemble anomalies are locally inflated. The proposed inflation strategy is computationally efficient and is aimed at restoring enough spread in the analysis ensemble after assimilating the observations. The performance of this scheme is tested against the relaxation to prior spread (RTPS) and adaptive prior inflation. For this purpose, two model are used: the three-variable Lorenz 63 system and the Community Atmosphere Model (CAM). In CAM, global refractivity, temperature, and wind observations from several sources are incorporated to perform a set of assimilation experiments using the Data Assimilation Research Testbed (DART). The proposed scheme is shown to yield better quality forecasts than the RTPS. Assimilation results further suggest that when model errors are small, both prior and posterior inflation are able to mitigate sampling errors with a slight advantage to posterior inflation. When large model errors, such as wind and temperature biases, are present, prior inflation is shown to be more accurate than posterior inflation. Densely observed regions as in the Northern Hemisphere present numerous challenges to the posterior inflation algorithm. A compelling enhancement to the performance of the filter is achieved by combining both adaptive inflation schemes.
Abstract
Global data assimilation systems for numerical weather prediction (NWP) are characterized by significant uncertainties in tropical analysis fields. Furthermore, the largest spread of global ensemble forecasts in the short range on all scales is in the tropics. The presented results suggest that these properties hold even in the perfect-model framework and the ensemble Kalman filter data assimilation with a globally homogeneous network of wind and temperature profiles. The reasons for this are discussed by using the modal analysis, which provides information about the scale dependency of analysis and forecast uncertainties and information about the efficiency of data assimilation to reduce the prior uncertainties in the balanced and inertio-gravity dynamics.
The scale-dependent representation of variance reduction of the prior ensemble by the data assimilation shows that the peak efficiency of data assimilation is on the synoptic scales in the midlatitudes that are associated with quasigeostrophic dynamics. In contrast, the variance associated with the inertia–gravity modes is less successfully reduced on all scales. A smaller information content of observations on planetary scales with respect to the synoptic scales is discussed in relation to the large-scale tropical uncertainties that current data assimilation methodologies do not address successfully. In addition, it is shown that a smaller reduction of the large-scale uncertainties in the prior state for NWP in the tropics than in the midlatitudes is influenced by the applied radius for the covariance localization.
Abstract
Global data assimilation systems for numerical weather prediction (NWP) are characterized by significant uncertainties in tropical analysis fields. Furthermore, the largest spread of global ensemble forecasts in the short range on all scales is in the tropics. The presented results suggest that these properties hold even in the perfect-model framework and the ensemble Kalman filter data assimilation with a globally homogeneous network of wind and temperature profiles. The reasons for this are discussed by using the modal analysis, which provides information about the scale dependency of analysis and forecast uncertainties and information about the efficiency of data assimilation to reduce the prior uncertainties in the balanced and inertio-gravity dynamics.
The scale-dependent representation of variance reduction of the prior ensemble by the data assimilation shows that the peak efficiency of data assimilation is on the synoptic scales in the midlatitudes that are associated with quasigeostrophic dynamics. In contrast, the variance associated with the inertia–gravity modes is less successfully reduced on all scales. A smaller information content of observations on planetary scales with respect to the synoptic scales is discussed in relation to the large-scale tropical uncertainties that current data assimilation methodologies do not address successfully. In addition, it is shown that a smaller reduction of the large-scale uncertainties in the prior state for NWP in the tropics than in the midlatitudes is influenced by the applied radius for the covariance localization.
The Data Assimilation Research Testbed (DART) is an open-source community facility for data assimilation education, research, and development. DART's ensemble data assimilation algorithms, careful software engineering, and diagnostic tools allow atmospheric scientists, oceanographers, hydrologists, chemists, and other geophysicists to build state-of-the-art data assimilation systems with unprecedented ease. For global numerical weather prediction, DART produces ensemble-mean analyses comparable to analyses from major centers while also providing initial conditions for ensemble predictions. In addition, DART supports more novel assimilation applications like parameter estimation, sensitivity analysis, observing system design, and smoothing. Implementing basic systems for large models requires only a few person-weeks; comprehensive systems have been built in a few months. Incorporating new observation types is also straightforward, requiring only a forward operator mapping between a model's state and an observation's expected value. Forward operators for standard, in situ observations and novel types, like GPS radio occultation soundings, are available. DART algorithms scale well on a variety of parallel architectures, allowing large data assimilation problems to be studied. DART also includes many low-order models and an ensemble assimilation tutorial appropriate for undergraduate and graduate instruction.
The Data Assimilation Research Testbed (DART) is an open-source community facility for data assimilation education, research, and development. DART's ensemble data assimilation algorithms, careful software engineering, and diagnostic tools allow atmospheric scientists, oceanographers, hydrologists, chemists, and other geophysicists to build state-of-the-art data assimilation systems with unprecedented ease. For global numerical weather prediction, DART produces ensemble-mean analyses comparable to analyses from major centers while also providing initial conditions for ensemble predictions. In addition, DART supports more novel assimilation applications like parameter estimation, sensitivity analysis, observing system design, and smoothing. Implementing basic systems for large models requires only a few person-weeks; comprehensive systems have been built in a few months. Incorporating new observation types is also straightforward, requiring only a forward operator mapping between a model's state and an observation's expected value. Forward operators for standard, in situ observations and novel types, like GPS radio occultation soundings, are available. DART algorithms scale well on a variety of parallel architectures, allowing large data assimilation problems to be studied. DART also includes many low-order models and an ensemble assimilation tutorial appropriate for undergraduate and graduate instruction.
Abstract
This paper explores the degree to which short-term forecasts with global models might be improved if clouds were fully included in a data assimilation system, so that observations of clouds affected all parts of the model state and cloud variables were adjusted during assimilation. The question is examined using a single ensemble data assimilation system coupled to two present-generation climate models with different treatments of clouds. “Perfect-model” experiments using synthetic observations, taken from a free run of the model used in subsequent assimilations, are used to circumvent complications associated with systematic model errors and observational challenges; these provide a rough upper bound on the utility of cloud observations with these models. A series of experiments is performed in which direct observations of the model’s cloud variables are added to the suite of observations being assimilated. In both models, observations of clouds reduce the 6-h forecast error, with much greater reductions in one model than in the other. Improvements are largest in regions where other observations are sparse. The two cloud schemes differ in their complexity and number of degrees of freedom; the model using the simpler scheme makes better use of the cloud observations because of the stronger correlations between cloud-related and dynamical variables (particularly temperature). This implies that the impact of real cloud observations will depend on both the strength of the instantaneous, linear relationships between clouds and other fields in the natural world, and how well each assimilating model’s cloud scheme represents those relationships.
Abstract
This paper explores the degree to which short-term forecasts with global models might be improved if clouds were fully included in a data assimilation system, so that observations of clouds affected all parts of the model state and cloud variables were adjusted during assimilation. The question is examined using a single ensemble data assimilation system coupled to two present-generation climate models with different treatments of clouds. “Perfect-model” experiments using synthetic observations, taken from a free run of the model used in subsequent assimilations, are used to circumvent complications associated with systematic model errors and observational challenges; these provide a rough upper bound on the utility of cloud observations with these models. A series of experiments is performed in which direct observations of the model’s cloud variables are added to the suite of observations being assimilated. In both models, observations of clouds reduce the 6-h forecast error, with much greater reductions in one model than in the other. Improvements are largest in regions where other observations are sparse. The two cloud schemes differ in their complexity and number of degrees of freedom; the model using the simpler scheme makes better use of the cloud observations because of the stronger correlations between cloud-related and dynamical variables (particularly temperature). This implies that the impact of real cloud observations will depend on both the strength of the instantaneous, linear relationships between clouds and other fields in the natural world, and how well each assimilating model’s cloud scheme represents those relationships.
Abstract
The authors report on the implementation and evaluation of a 48-member ensemble adjustment Kalman filter (EAKF) for the ocean component of the Community Climate System Model, version 4 (CCSM4). The ocean assimilation system described was developed to support the eventual generation of historical ocean-state estimates and ocean-initialized climate predictions with the CCSM4 and its next generation, the Community Earth System Model (CESM). In this initial configuration of the system, daily subsurface temperature and salinity data from the 2009 World Ocean Database are assimilated into the ocean model from 1 January 1998 to 31 December 2005. Each ensemble member of the ocean is forced by a member of an independently generated CCSM4 atmospheric EAKF analysis, making this a loosely coupled framework. Over most of the globe, the time-mean temperature and salinity fields are improved relative to an identically forced ocean model simulation without assimilation. This improvement is especially notable in strong frontal regions such as the western and eastern boundary currents. The assimilation system is most effective in the upper 1000 m of the ocean, where the vast majority of in situ observations are located. Because of the shortness of this experiment, ocean variability is not discussed. Challenges that arise from using an ocean model with strong regional biases, coarse resolution, and low internal variability to assimilate real observations are discussed, and areas of ongoing improvement for the assimilation system are outlined.
Abstract
The authors report on the implementation and evaluation of a 48-member ensemble adjustment Kalman filter (EAKF) for the ocean component of the Community Climate System Model, version 4 (CCSM4). The ocean assimilation system described was developed to support the eventual generation of historical ocean-state estimates and ocean-initialized climate predictions with the CCSM4 and its next generation, the Community Earth System Model (CESM). In this initial configuration of the system, daily subsurface temperature and salinity data from the 2009 World Ocean Database are assimilated into the ocean model from 1 January 1998 to 31 December 2005. Each ensemble member of the ocean is forced by a member of an independently generated CCSM4 atmospheric EAKF analysis, making this a loosely coupled framework. Over most of the globe, the time-mean temperature and salinity fields are improved relative to an identically forced ocean model simulation without assimilation. This improvement is especially notable in strong frontal regions such as the western and eastern boundary currents. The assimilation system is most effective in the upper 1000 m of the ocean, where the vast majority of in situ observations are located. Because of the shortness of this experiment, ocean variability is not discussed. Challenges that arise from using an ocean model with strong regional biases, coarse resolution, and low internal variability to assimilate real observations are discussed, and areas of ongoing improvement for the assimilation system are outlined.
