A generic data assimilation framework for vehicle trajectory reconstruction on signalized urban arterials using particle filters

With trajectory data, a complete microscopic and macroscopic picture of traffic flow operations can be obtained. However, trajectory data are difficult to observe over large spatiotemporal regions—particularly in urban contexts—due to practical, technical and financial constraints. The next best thing is to estimate plausible trajectories from whatever data are available. This paper presents a generic data assimilation framework to reconstruct such plausible trajectories on signalized urban arterials using microscopic traffic flow models and data from loops (individual vehicle passages and thus vehicle counts); traffic control data; and (sparse) travel time measurements from whatever source available. The key problem we address is that loops suffer from miss- and over-counts, which result in unbounded errors in vehicle accumulations, rendering trajectory reconstruction highly problematic. Our framework solves this problem in two ways. First, we correct the systematic error in vehicle accumulation by fusing the counts with sparsely available travel times. Second, the proposed framework uses particle filtering and an innovative hierarchical resampling scheme, which effectively integrates over the remaining error distribution, resulting in plausible trajectories. The proposed data assimilation framework is tested and validated using simulated data. Experiments and an extensive sensitivity analysis show that the proposed method is robust to errors both in the model and in the measurements, and provides good estimations for vehicle accumulation and vehicle trajectories with moderate sensor quality. The framework does not impose restrictions on the type of microscopic models used and can be naturally extended to include and estimate additional trajectory attributes such as destination and path, given data are available for assimilation.     

Estimates of vertical mixing during a Lagrangian experiment off the Galician coast

An adjoint 1-D model was used to determine vertical diffusivity coefficients from temperature profiles collected within a filament escaping from the Galician coast following an upwelling event. The optimisation scheme ended with relatively high diffusivity values within the thermocline (9×10⁻⁵ m² s⁻¹). Such high values are relevant for biogeochemical exchanges between surface and deep waters in stratified areas.The optimised values were several orders of magnitude higher than the bulk of diffusivity measurements recorded with a free-falling device; however, the optimisation solution was consistent with the arithmetic mean of the measurements in the thermocline (7.7×10⁻⁵ m² s⁻¹), giving more weight to the few largest values. Below the thermocline, the data assimilation method failed because of the three-dimensional nature of the advective field of the upwelling system. Ignoring this advective forcing in the model led to estimates that were two orders of magnitude too high.The results suggest that turbulent mixing is a random process where a few intense events determine the average mixing that drives the long-term evolution of the water column structure. This statistical property is very important when one wants to use instantaneous diffusivity measurements for modelling purposes.     

Altimeter assimilation in the OCCAM global model: Part I: A twin experiment

An altimeter data assimilation scheme has been tested in the OCCAM (Ocean Circulation and Climate Advanced Modelling) global 1/4°, 36-level model using a twin experiment format. The Cooper and Haines displacement scheme is used. The method works well in most regions and depths. Currents and densities in the top 1000 m generally improve by over 50–70% after 5 months of sea level assimilation every 15 days. Below 1000 m, an error reduction of up to 50% is achieved. The errors remain low during a further 60-day run without assimilation. Diagnostics for the North Atlantic, the Tropical Pacific and the Antarctic Circumpolar Current are shown alongside the global averages.The main problems encountered were in weakly stratified regions of the Antarctic and Arctic seas. A scale selective filter is developed to avoid assimilating scales much larger than the local deformation radius, and this avoids the adverse assimilation effects in the southern oceans. A companion paper uses this scheme to assimilate TOPEX and ERS-1 altimeter maps.     

The influence of food quality on the nutritional acclimation of the copepod Acartia clausi

The influence of food quality on the nutritional metabolism of Acartia clausi was studied experimentally using four different diets: (1) diatom cells of Thalassiosira weissflogii, (2) detritus prepared from the same culture, (3) 50 : 50 mix on a protein basis of the two previous diets, and 4) dinoflagellate cells of Prorocentrum micans. For each trophic, ingestion, gut transit time, trypsin activity and K_m (half saturation constant) were measured at limiting and saturating concentration. Assimilation rates were also estimated for both pure diatoms and mixed live–detrital cell diets. Ingestion followed a Holling type 2 response for diets 1 and 4, a linear one for detritus and an intermediate response for diet 3. Gut transit time displayed different adaptive changes with food regime depending on protein concentration. Trypsin activity was lower for detrital food and trypsin K_m significantly decreased with increasing concentration of live diatoms. Assimilation rates were higher for live food than for mixed live–detrital food. Results illustrated that Acartia-type copepods optimize nitrogen or protein uptake. They suggest that besides chemoreception-mediated selectivity, internal controls by digestion and assimilation also regulate ingestion (feed-back).     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

美国海军业务化海洋预报系统综述

美国海军海洋预报系统最早始于苏联时期,自此之后,美国海军制定了一系列计划,构建海洋观测网络,开发海洋数值模式,发展军事海洋学,旨在建成一套完整、精确的全球海洋预报系统并实现业务化应用。美国海军海洋学一直走在世界前沿,本文将对美国海军海洋预报系统的发展历程、成果及其现状进行阐述,这对于我国军事海洋学以及海洋预报系统的发展具有重要借鉴作用。     

试析语流音变对词语发音和词形的影响

音位与音位结合的时候,一个音位由于受邻音的影响,或者由于说话的快慢、高低、强弱的不同,可能发生不同的变化．文章从音的同化、异化、强化（插入音或加音）、弱化（弱化和脱落）的角度深入分析了语流音变对词语发音和构形的影响．对于长期困扰语言学习者的一些常见语音变化规律和发音现象作出了较详尽的阐释,同时发音的改变最终必然反映到词形变化上．     

Gaussian anamorphosis extension of the DEnKF for combined state parameter estimation: Application to a 1D ocean ecosystem model

We consider the problem of combined state-parameter estimations in biased nonlinear models with non-Gaussian extensions of the Deterministic Ensemble Kalman Filter (DEnKF). We focus on the particular framework of ocean ecosystem models. Such models present important obstacles to the use of data assimilation methods based on Kalman filtering due to the non-linearity of the models, the constraints of positiveness that apply to the variables and parameters, and the non-Gaussian distribution of the variables in which they result.We present extensions of the DEnKF dealing with these difficulties by introducing a nonlinear change of variables (anamorphosis function) in order to execute the analysis step with Gaussian transformed variables and parameters. Several strategies to build the anamorphosis functions are investigated and compared within the framework of twin experiments realized in a simple 1D ocean ecosystem model. A solution to the problem of the specification of the observation error for transformed observations is suggested. The study highlights the inability of the plain DEnKF with a simple post-processing of the negative values to properly estimate parameters when constraints of positiveness apply to the variables. It goes on to show that the introduction of the Gaussian anamorphosis can remedy these assimilation biases.     

Quality control of ocean temperature and salinity profiles — Historical and real-time data

A new automated quality control system for oceanic temperature and temperature–salinity profiles is presented. Substantial development was needed for some of the quality control algorithms although the checks were based on documented procedures used elsewhere where appropriate. A new automated ship track check was developed: the results of an undetected position error can be very damaging to ocean analyses. Also important is a check against a gridded background, this can be a climatology but near the surface it is advantageous to use an estimate that is evolving over time. Bayesian probability theory is used in the background check and the associated check against nearby observations (buddy check). The system was used to process archive data for 1956–2004. As a by-product monthly model-free objective analyses for this period were produced. Versions of the system are used for near-real time ocean analysis and for initialising both short-range ocean forecasts and seasonal atmosphere–ocean forecasts. The main features of the oceanic observing systems are presented along with quality control statistics, examples of errors that can occur and some additional problematic cases.