HOUSING INTEGRATION AND RENT SUPPLEMENTS TO EXISTING HOUSING*

This work examines the use of an annually updated metropolitan data system as a tool for evaluating benign racial location policies during the first year of the Section 8 Existing Rent Supplement program. Case-study evidence shows that the rent supplement program (with benign locational policies) is ineffective in reducing the isolation faced by minority groups, and that the information used to monitor compliance of the program with benign locational policies lags badly behind the true state of affairs. As a result, new areas of racial concentration have been encouraged and augmented.     

Improved magnetic anomalies of the Antarctic lithosphere from satellite and near-surface data

The Antarctic magnetic anomaly map compiled marine and airborne surveys collected south of 60°S through 1999 and used Magsat data to help fill in the regional gaps between the surveys. Ørsted and CHAMP satellite magnetic observations with greatly improved measurement accuracies and temporal and spatial coverage of the Antarctic, have now supplanted the Magsat data. We combined the new satellite observations with the near-surface survey data for an improved magnetic anomaly map of the Antarctic lithosphere. Specifically, we separated the crustal from the core and external field components in the satellite data using crustal thickness variations estimated from the terrain and the satellite-derived free-air gravity observations. Regional gaps in the near-surface surveys were then filled with predictions from crustal magnetization models that jointly satisfied the near-surface and satellite crustal anomalies. Comparisons in some of the regional gaps that also considered newly acquired aeromagnetic data demonstrated the enhanced anomaly estimation capabilities of the predictions over those from conventional minimum curvature and spherical harmonic geomagnetic field models. We also noted that the growing number of regional and world magnetic survey compilations involve coverage gaps where these procedures can contribute effective near-surface crustal anomaly estimates.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II - Performance and application to the California Current System

The Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation systems have been systematically applied to the mesoscale circulation environment of the California Current to demonstrate the performance and practical utility of the various components of ROMS 4D-Var. In particular, we present a comparison of three approaches to 4D-Var, namely: the primal formulation of the incremental strong constraint approach; the dual formulation “physical-space statistical analysis system”; and the dual formulation indirect representer approach. In agreement with theoretical considerations all three approaches converge to the same ocean circulation estimate when using the same observations and prior information. However, the rate of convergence of the dual formulation was found to be inferior to that of the primal formulation. Other aspects of the 4D-Var performance that relate to the use of multiple outer-loops, preconditioning, and the weak constraint are also explored. A systematic evaluation of the impact of the various components of the 4D-Var control vector (i.e. the initial conditions, surface forcing and open boundary conditions) is also presented. It is shown that correcting for uncertainties in the model initial conditions exerts the largest influence on the ability of the model to fit the available observations. Various important diagnostics of 4D-Var are also examined, including estimates of the posterior error, the information content of the observation array, and innovation-based consistency checks on the prior error assumptions. Using these diagnostic tools, we find that more than 90% of the observations assimilated into the model provide redundant information. This is a symptom of the large percentage of satellite data that are used and to some extent the nature of the data processing employed. This is the second in a series of three papers describing the ROMS 4D-Var systems.     

Spatiotemporal aspects of silica buffering in restored tidal marshes

Losses of pelagic diatom production resulting from silica limitation have not only been blamed for toxic algal blooms, but for the reduction in ability of coastal food webs to support higher trophic levels. Recent research has shown the importance of advective seepage water fluxes of dissolved silica (DSi) from freshwater marshes to pelagic waters during moments of riverine Si-limitation. In this study, we investigated the potential impact of recently installed new tidal areas along the Schelde estuary, located in former polder areas and characterized by so-called controlled reduced tidal regimes (CRT). Nine mass-balance studies were conducted in a newly constructed CRT in the freshwater Schelde estuary. During complete tidal cycles both DSi and amorphous silica (ASi) concentrations were monitored at the entrance culverts and in different habitats in the marsh. A swift DSi-delivery capacity was observed despite the shifted spatiotemporal frame of exchange processes compared to reference marshes. As silica-accumulating vegetation is not yet present, and difference with reference marshes’ deliveries is surprisingly small, we indicate diatomaceous debris and phytoliths to be the main silica source. Although further research is necessary on the driving forces of the different processes involved, restoration of former agricultural areas under CRT-regime provide the potential to buffer silica in the estuary.     

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.     

Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.     

Shallow-water longshore drift-fed submarine fan deposition (Moisie River Delta, Eastern Canada)

Submarine canyons and associated submarine fans are in some cases located at the end of a littoral cell where they act as conduits for the transfer of eroded terrigenous sediments to the marine environment. Such fans are generally found in deep-water settings at >500 m water depth. Offshore the Moisie River Delta (NW Gulf of St. Lawrence, Eastern Canada), high-resolution multibeam bathymetry and seismic data led to the discovery of an unusually shallow submarine fan (≤60 m) located at the end of a littoral cell. Sediment is transported westward on the shallow coastal shelf, as demonstrated by the downcurrent displacement of oblique nearshore sandbars where the shelf narrows to less than 1 km. The steep slope near the end of the littoral cell is incised by a channel that feeds a submarine fan composed of smaller channels and depositional lobes. According to existing Holocene evolution models for the region, the fan formed within the last 5,000 years. Its evolution is largely due to the transport of sediment by longshore drift. Multibeam echosounder and seismic data also reveal that the gravity-driven accretion of the submarine fan is characterized mainly by two processes, i.e., frequent small-scale, downslope migration of sandwaves on the slope, and more episodic slumping/turbidity-current activity in the deeper part of the fan. This study documents that, besides their common deep-water location, smaller-scale submarine fans can occur also in very shallow water, implying that they could be more frequent than previously thought both in modern environments and in the rock record.     

Modeled Sverdrup flow in the North Atlantic from 11 different wind stress climatologies

In studies of large-scale ocean dynamics, often quoted values of Sverdrup transport are computed using the Hellerman–Rosenstein wind stress climatology. The Sverdrup solution varies, however, depending on the wind set used. We examine the differences in the large-scale upper ocean response to different surface momentum forcing fields for the North Atlantic Ocean by comparing the different Sverdrup interior/Munk western boundary layer solutions produced by a 1/16° linear numerical ocean model forced by 11 different wind stress climatologies. Significant differences in the results underscore the importance of careful selection of a wind set for Sverdrup transport calculation and for driving nonlinear models. This high-resolution modeling approach to solving the linear wind-driven ocean circulation problem is a convenient way to discern details of the Sverdrup flow and Munk western boundary layers in areas of complicated geometry such as the Caribbean and Bahamas. In addition, the linear solutions from a large number of wind sets provide a well-understood baseline oceanic response to wind stress forcing and thus, (1) insight into the dynamics of observed circulation features, by themselves and in conjunction with nonlinear models, and (2) insight into nonlinear model sensitivity to the choice of wind-forcing product.The wind stress products are evaluated and insight into the linear dynamics of specific ocean features is obtained by examining wind stress curl patterns in relation to the corresponding high-resolution linear solutions in conjunction with observational knowledge of the ocean circulation. In the Sverdrup/Munk solutions, the Gulf Stream pathway consists of two branches. One separates from the coast at the observed separation point, but penetrates due east in an unrealistic manner. The other, which overshoots the separation point at Cape Hatteras and continues to flow northward along the continental boundary, is required to balance the Sverdrup interior transport. A similar depiction of the Gulf Stream is commonly seen in the mean flow of nonlinear, eddy-resolving basin-scale models of the North Atlantic Ocean. An O(1) change from linear dynamics is required for realistic simulation of the Gulf Stream pathway. Nine of the eleven Sverdrup solutions have a C-shaped subtropical gyre, similar to what is seen in dynamic height contours derived from observations. Three mechanisms are identified that can contribute to this pattern in the Sverdrup transport contours. Along 27°N, several wind sets drive realistic total western boundary current transport (within 10% of observed) when a 14 Sv global thermohaline contribution is added (COADS, ECMWF 10 m re-analysis and operational, Hellerman–Rosenstein and National Centers for Environmental Prediction (NCEP) surface stress re-analysis), a few drive transport that is substantially too high (ECMWF 1000 mb re-analysis and operational and Isemer–Hasse) and Fleet Numerical Meteorology and Oceanography Center (FNMOC) surface stresses give linear transport that is slightly weaker than observed. However, higher order dynamics are required to explain the partitioning of this transport between the Florida Straits and just east of the Bahamas (minimal in the linear solutions vs. 5 Sv observed east of the Bahamas). Part of the Azores Current transport is explained by Sverdrup dynamics. So are the basic path of the North Atlantic Current (NAC) and the circulation features within the Intra-Americas Sea (IAS), when a linear rendition of the northward upper ocean return flow of the global thermohaline circulation is added in the form of a Munk western boundary layer.