High-Frequency Orographically Forced Variability in a Single-Layer Model of the Martian Atmosphere

A shallow water model with realistic topography and idealized zonal wind forcing is used to investigate orographically forced modes in the Martian atmosphere. Locally, the model reproduces well the climatology at the sites of Viking Lander I and II (VLl and VL2) as inferred from the Viking Lander fall and spring observations. Its variability at those sites is dominated by a 3-sol (Martian solar day) oscillation in the region of VLl and by a 6-sol oscillation in that of VL2. These oscillations are forced by the zonal asymmetries of the Martian mountain field. It is suggested that they contribute to the observed variability by reinforcing the baroclinic oscillations with nearby periods identified in observational studies. The spatial variability associated with the orographically forced oscillations is studied by means of extended empirical orthogonal function analysis. The 3-sol VL1 oscillation corresponds to a tropical, eastward-traveling, zonal-wavenumber one pattern. The 6-sol VL2 oscillation is characterized by two midlatitude, eastward-traveling, mixed zonal-wavenumber one and two and zonal-wavenumber three and four patterns, with respective periods near 6.1 and 5.5 sols. The corresponding phase speeds arc in agreement with the conclusions drawn from the VL2 observations

An Ensemble Recentering Kalman Filter with an Application to Argo Temperature Data Assimilation into the NASA GEOS-5 Coupled Model

A two-step ensemble recentering Kalman filter (ERKF) analysis scheme is introduced. The algorithm consists of a recentering step followed by an ensemble Kalman filter (EnKF) analysis step. The recentering step is formulated such as to adjust the prior distribution of an ensemble of model states so that the deviations of individual samples from the sample mean are unchanged but the original sample mean is shifted to the prior position of the most likely particle, where the likelihood of each particle is measured in terms of closeness to a chosen subset of the observations. The computational cost of the ERKF is essentially the same as that of a same size EnKF. The ERKF is applied to the assimilation of Argo temperature profiles into the OGCM component of an ensemble of NASA GEOS-5 coupled models. Unassimilated Argo salt data are used for validation. A surprisingly small number (16) of model trajectories is sufficient to significantly improve model estimates of salinity over estimates from an ensemble run without assimilation. The two-step algorithm also performs better than the EnKF although its performance is degraded in poorly observed regions

Paralellized ensemble Kalman filter for hydraulic conductivity characterization

[EN] The ensemble Kalman filter (EnKF) is nowadays recognized as an excellent inverse method for hydraulic conductivity characterization using transient piezometric head data. Its implementation is well suited for a parallel computing environment. A parallel code has been designed that uses parallelization both in the forecast step and in the analysis step. In the forecast step, each member of the ensemble is sent to a different processor, while in the analysis step, the computations of the covariances are distributed between the different processors. An important aspect of the parallelization is to limit as much as possible the communication between the processors in order to maximize execution time reduction. Four tests are carried out to evaluate the performance of the parallelization with different ensemble and model sizes. The results show the savings provided by the parallel EnKF, especially for a large number of ensemble realizations. (c) 2012 Elsevier Ltd. All rights reserved.The first author acknowledges the financial support from China Scholarship Council (CSC). Financial support to carry out this work was also received from the Spanish Ministry of Science and Innovation through project CGL2011-23295, and from the Universitat Politecnica de Valencia through project PERFORA.Xu, T.; Gómez-Hernández, JJ.; Li ., L.; Zhou ., H. (2013). Paralellized ensemble Kalman filter for hydraulic conductivity characterization. Computers and Geosciences. 52:42-49. https://doi.org/10.1016/j.cageo.2012.10.007S42495

Completing the Feedback Loop: The Impact of Chlorophyll Data Assimilation on the Ocean State

In anticipation of the integration of a full biochemical model into the next generation GMAO coupled system, an intermediate solution has been implemented to estimate the penetration depth (1Kd_PAR) of ocean radiation based on the chlorophyll concentration. The chlorophyll is modeled as a tracer with sources-sinks coming from the assimilation of MODIS chlorophyll data. Two experiments were conducted with the coupled ocean-atmosphere model. In the first, climatological values of Kpar were used. In the second, retrieved daily chlorophyll concentrations were assimilated and Kd_PAR was derived according to Morel et al (2007). No other data was assimilated to isolate the effects of the time-evolving chlorophyll field. The daily MODIS Kd_PAR product was used to validate the skill of the penetration depth estimation and the MERRA-OCEAN re-analysis was used as a benchmark to study the sensitivity of the upper ocean heat content and vertical temperature distribution to the chlorophyll input. In the experiment with daily chlorophyll data assimilation, the penetration depth was estimated more accurately, especially in the tropics. As a result, the temperature bias of the model was reduced. A notably robust albeit small (2-5 percent) improvement was found across the equatorial Pacific ocean, which is a critical region for seasonal to inter-annual prediction

The GEOS-iODAS: Description and Evaluation

This report documents the GMAO's Goddard Earth Observing System sea ice and ocean data assimilation systems (GEOS iODAS) and their evolution from the first reanalysis test, through the implementation that was used to initialize the GMAO decadal forecasts, and to the current system that is used to initialize the GMAO seasonal forecasts. The iODAS assimilates a wide range of observations into the ocean and sea ice components: in-situ temperature and salinity profiles, sea level anomalies from satellite altimetry, analyzed SST, and sea-ice concentration. The climatological sea surface salinity is used to constrain the surface salinity prior to the Argo years. Climatological temperature and salinity gridded data sets from the 2009 version of the World Ocean Atlas (WOA09) are used to help constrain the analysis in data sparse areas. The latest analysis, GEOS ODAS5.2, is diagnosed through detailed studies of the statistics of the innovations and analysis departures, comparisons with independent data, and integrated values such as volume transport. Finally, the climatologies of temperature and salinity fields from the Argo era, 2002-2011, are presented and compared with the WOA09

Background Error Covariance Estimation Using Information from a Single Model Trajectory with Application to Ocean Data Assimilation

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory.SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications

Lessons Learned from Assimilating Altimeter Data into a Coupled General Circulation Model with the GMAO Augmented Ensemble Kalman Filter

Satellite altimetry measurements have provided global, evenly distributed observations of the ocean surface since 1993. However, the difficulties introduced by the presence of model biases and the requirement that data assimilation systems extrapolate the sea surface height (SSH) information to the subsurface in order to estimate the temperature, salinity and currents make it difficult to optimally exploit these measurements. This talk investigates the potential of the altimetry data assimilation once the biases are accounted for with an ad hoc bias estimation scheme. Either steady-state or state-dependent multivariate background-error covariances from an ensemble of model integrations are used to address the problem of extrapolating the information to the sub-surface. The GMAO ocean data assimilation system applied to an ensemble of coupled model instances using the GEOS-5 AGCM coupled to MOM4 is used in the investigation. To model the background error covariances, the system relies on a hybrid ensemble approach in which a small number of dynamically evolved model trajectories is augmented on the one hand with past instances of the state vector along each trajectory and, on the other, with a steady state ensemble of error estimates from a time series of short-term model forecasts. A state-dependent adaptive error-covariance localization and inflation algorithm controls how the SSH information is extrapolated to the sub-surface. A two-step predictor corrector approach is used to assimilate future information. Independent (not-assimilated) temperature and salinity observations from Argo floats are used to validate the assimilation. A two-step projection method in which the system first calculates a SSH increment and then projects this increment vertically onto the temperature, salt and current fields is found to be most effective in reconstructing the sub-surface information. The performance of the system in reconstructing the sub-surface fields is particularly impressive for temperature, but not as satisfactory for salt

Ensemble Kalman filter assimilation of temperature and altimeter data with bias correction and application to seasonal prediction

To compensate for a poorly known geoid, satellite altimeter data is usually analyzed in terms of anomalies from the time mean record. When such anomalies are assimilated into an ocean model, the bias between the climatologies of the model and data is problematic. An ensemble Kalman filter (EnKF) is modified to account for the presence of a forecast-model bias and applied to the assimilation of TOPEX/Poseidon (T/P) altimeter data. The online bias correction (OBC) algorithm uses the same ensemble of model state vectors to estimate biased-error and unbiased-error covariance matrices. Covariance localization is used but the bias covariances have different localization scales from the unbiased-error covariances, thereby accounting for the fact that the bias in a global ocean model could have much larger spatial scales than the random error.The method is applied to a 27-layer version of the Poseidon global ocean general circulation model with about 30-million state variables. Experiments in which T/P altimeter anomalies are assimilated show that the OBC reduces the RMS observation minus forecast difference for sea-surface height (SSH) over a similar EnKF run in which OBC is not used. Independent in situ temperature observations show that the temperature field is also improved. When the T/P data and in situ temperature data are assimilated in the same run and the configuration of the ensemble at the end of the run is used to initialize the ocean component of the GMAO coupled forecast model, seasonal SSH hindcasts made with the coupled model are generally better than those initialized with optimal interpolation of temperature observations without altimeter data. The analysis of the corresponding sea-surface temperature hindcasts is not as conclusive

Ensemble Data Assimilation Without Ensembles: Methodology and Application to Ocean Data Assimilation

Two methods to estimate background error covariances for data assimilation are introduced. While both share properties with the ensemble Kalman filter (EnKF), they differ from it in that they do not require the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The first method is referred-to as SAFE (Space Adaptive Forecast error Estimation) because it estimates error covariances from the spatial distribution of model variables within a single state vector. It can thus be thought of as sampling an ensemble in space. The second method, named FAST (Flow Adaptive error Statistics from a Time series), constructs an ensemble sampled from a moving window along a model trajectory. The underlying assumption in these methods is that forecast errors in data assimilation are primarily phase errors in space and/or time

(Co)constructing critical pedagogies: Expanding on our department’s approach to language teaching

In this report, we—the members of a curriculum working group (CWG) in Penn State’s German department—describe our efforts to decenter our German language sequence by integrating critical pedagogies into our department’s existing communicative language teaching (CLT) approach. We trace our process towards this goal, beginning with an exploration into and analysis of two critical pedagogies, namely Antiracist Pedagogy (ARP) and Social Justice Pedagogy (SJP). We ultimately adopt SJP because we find it to be a better fit for our purposes in German language instruction. We offer a framework to evaluate and didacticize existing as well as newly created course materials, guided by social justice (SJ) learning objectives. To illustrate our work, we describe the creation and implementation of an instructional unit in an intermediate German language course. Reflections from this course’s instructor and student reactions concerning this unit’s instruction—as well as SJP in the language classroom in general—make evident the importance of critical perspectives regarding curricular development in fostering equitable classrooms