ON THE BREMMER SERIES DECOMPOSITION: EQUIVALENCE BETWEEN TWO DIFFERENT APPROACHES*

A Bremmer Series decomposition of the solution y(t) to the lossless wave equation in layered media is where the y_j(t) are physically meaningful constituents (i.e., y₁(t) are primaries, y₂(t) are secondaries, etc.). This paper reviews Mendel's state space models for generating the constituents; reviews Bremmer's integral equation models for generating the constituents; and demonstrates how Mendel's state space models can be obtained by a careful decomposition of Bremmer's integral equation models. It shows that Mendel's equations can be viewed as approximate numerical solutions of Bremmer's integral equations. In a lossless homogeneous medium, the approximations become exact.     

Global and temporal variations in hydrocarbons and nitriles in Titan's stratosphere for northern winter observed by Cassini/CIRS

Mid-infrared spectra measured by Cassini's Composite InfraRed Spectrometer (CIRS) between July 2004 and January 2007 (Ls=293°-328°) have been used to determine stratospheric temperature and abundances of C₂H₂, C₃H₄, C₄H₂, HCN, and HC₃N. Over 65,000 nadir spectra with spectral resolutions of 0.5 and 2.5 cm⁻¹ were used to probe spatial and temporal composition variations in Titan's stratosphere. Cassini's 180° orbital transfer in mid-2006 allowed low emission angle observations of the north polar region for the first time in the mission and allowed us to probe the full latitude range. We present the first measurements of composition variations within the polar vortex, which display increasing abundances right up to 90° N. The lack of a homogeneous abundance-latitude variation within the vortex indicates limited horizontal mixing and suggests that subsidence is greatest at the vortex core. Contrary to numerical model predictions and tropospheric cloud observations, we do not see any evidence for a secondary circulation cell near the south pole, which suggests a single Hadley-type circulation in the stratosphere at this epoch. This difference can be reconciled if the secondary cell is restricted to altitudes below 100 km, where there is no sensitivity in our data. Temporal variations in composition were observed in the south, with volatile species becoming less abundant as the season progressed. The observed variations are compared to numerical model predictions and observations from Voyager.     

Forecasting weather radar propagation conditions

Summary The increasing use of weather radar quantitative precipitation estimates, particularly in automatic applications such as operational hydrometeorological modelling or assimilation in numerical weather prediction (NWP) models, has promoted the development of quality control procedures on radar data. Anomalous propagation (AP) of the radar beam due to deviation from the standard refractivity vertical profile, is one of the factors that may affect seriously the quality of radar observations because of the increase in quantity and intensity of non-precipitating clutter echoes and consequent contamination of the estimated rainfall field. Another undesired effect of AP is the change in the expected radar echo height, which may be relevant when correcting for beam blockage in radar rainfall estimation in complex terrain. The aim of this paper is to study the use of NWP mesoscale forecasts to predict and monitor AP events. A nested 15-km grid resolution version of the MASS model has been used to retrieve refractivity profiles in the coastal area of Barcelona, near a weather radar and a radiosonde station. Using the refractivity profiles two different magnitudes were computed: the vertical refractivity profile of the lowest 1000 m layer and a ducting index which describes the existence and intensity of the most super-refractive layer contained in the lowest 3-km layer. A comparison between model forecasts and radiosonde diagnostics during a six-month period showed that the model tended to underestimate the degree of super-refraction, with a bias of 4 km⁻¹ and RMSE of 11 km⁻¹ in the 1-km vertical refractivity gradient. Further analysis of the data showed that a combination of previous observations and forecasts allowed to produce modified forecasts improving the original direct model output, decreasing substantially the bias, reducing the RMSE by 20% and improving the skill by 40%, beating also radiosonde observations persistence.     

An Overview of Dissolved and Suspended Matter Fluxes in the Loire River Basin: Natural and Anthropogenic Inputs

The spatial and temporal distributions of major elements were investigated in the surface waters and in associated suspended matter at two sites of the upper Loire basin (Orleans and Brehemont) between 1995 and 1998.According to geochemical and isotopic patterns, the dissolved load appears to result from a process of mixing rainwater inputs, weathering processes of carbonate and silicate bedrock, and agricultural and urban inputs. Natural inputs influence 60% of water chemical composition at both sites. Annual dissolved fluxes were estimated to be 1300 10³ t/y at Orleans and 1620 10³ t/y at Brehemont. Major elements are transported mainly in the dissolved fraction. After correcting for atmospheric and anthropogenic inputs, the silicate specific export rate was calculated to be 11 t/y/km² throughout the basin and the carbonate specific export rate to be from 47 t/y/km² at Orleans to 23 t/y/km² at Brehemont.The suspended load appears to result from at least two particle reservoirs: a silicate reservoir and a carbonate reservoir. The silicate reservoir has a detrital origin, mainly during periods of high flow, while the carbonate reservoir has a detrital origin during periods of high flow and an authigenic origin during periods of low flow. Of the total annual flow of suspended matter, this authigenic material represents 16% at Orleans, 25% at Brehemont and 37% in the fluvial part of the estuary. After correcting authigenic inputs, the specific export rate due to mechanical weathering was estimated to be 8 t/y/km² throughout the Loire basin.     

A distributed approach for calculating inundation height based on Dijkstra's algorithm

This research proposed a parallelized approach to scaling up the calculation of inundation height, the minimum sea‐level rise required to inundate a cell on a digital elevation model, which is based on Dijkstra's algorithm for shortest‐path calculations on a graph. Our approach is based on the concepts of spatial decomposition, calculate‐and‐correct, and a master/worker parallelization paradigm. The approach was tested using the U.S. Coastal Relief Model (CRM) dataset from the National Geophysical Data Center on a multicore desktop computer and various supercomputing resources through the U.S. Extreme Science and Engineering Discovery Environment (XSEDE) program. Our parallel implementation not only enables computations that were larger than previously possible, but also significantly outperforms serial implementations with respect to running time and memory footprint as the number of processing cores increases. The efficiency of the scalability seemed to be tied to tile size and flattened out at a certain number of workers.     

Trace elements in ocean ridge basalts

A correlary of sea floor spreading is that the production rate of ocean ridge basalts exceeds that of all other volcanic rocks on the earth combined. Basalts of the ocean ridges bring with them a continuous record in space and time of the chemical characteristics of the underlying mantle. The chemical record is once removed, due to chemical fractionation during partial melting. Chemical fractionations can be evaluated by assuming that peridotite melting has proceeded to an olivine-orthopyroxene stage, in which case the ratios of a number of magmaphile elements in the extracted melt closely match the ratios in the mantle. Comparison of ocean ridge basalts and chondritic meteorites reveals systematic patterns of element fractionation, and what is probably a double depletion in some elements. The first depletion is in volatile elements and is due to high accretion temperatures of a large percentage of the earth from the solar nebula. The second depletion is in the largest, most highly charged lithophile elements (“incompatible elements”), probably because the mantle source of the basalts was melted previously, and the melt, enriched in these elements, was removed. Migration of melt relative to solid under ocean ridges and oceanic plates, element fractionation at subduction zones, and fractional melting of amphibolite in the Precambrian are possible mechanisms for depleting the mantle in incompatible elements. Ratios of transition metals in the mantle source of ocean ridge basalts are close to chondritic, and contrast to the extreme depletion of refractory siderophile elements, the reason for which remains uncertain. Variation of ocean ridge basalt chemistry along the length of the ridge has been correlated with ridge elevation. Thus chemically anomalous ridge segments up to 1000 km long appear to broadly coincide with regions of high magma production (plumes, hot spots). Basalt heterogeneity at a single location indicates mantle heterogeneity on a smaller scale. Variation of ocean ridge basalt chemistry with time has not been established, in fact, criteria for recognizing old oceanic crust in ophiolite terrains are currently under debate. The similarity of rare earth element patterns in basalt from ocean ridges, back-arc basins, some young island arcs, and some continental flood basalts illustrates the dangers of tectonic labeling by rare earth element pattern.     

Analysis of reservoir water quality using fuzzy synthetic evaluation

A general methodology for fuzzy synthetic evaluation is developed and illustrated with a case study of trophic status assessment for Fei-Tsui Reservoir in Taiwan. The historical data base was collected from the management agency of Fei-Tsui Reservoir from 1987 to 1996. In fuzzy synthetic evaluation, the classification is determined by a matrix operation of the weighted vector with the fuzzy evaluation matrix. After all individual membership functions of evaluated factors have been determined, the fuzzy evaluation matrix can be established. The weighted vector is determined by the analytic hierarchy process method (AHP). The results of this investigation show that the long-term change of water quality and the overturn phenomena cannot be observed with the Carlson index from 1987 to 1992 but is expressed by fuzzy synthetic evaluation. Fuzzy synthetic evaluation is better suited than the Carlson index to rating the trophic status of self-sustaining lakes. Interpretation of the results can provide valuable information to decision makers and aid reservoir management.     

Road network detection from SPOT imagery for updating geographical information systems in the rural–urban fringe

Abstract

Visual interpretation of high-resolution satellite data has been useful for mapping linear features, such as roads and updating land-use changes. However, it would be beneficial to map new road networks digitally from satellite data to update digital databases using semi-automated techniques. In this paper, an algorithm called Gradient Direction Profile Analysis (GDPA) is used to extract road networks digitally from SPOT High Resolution Visible (HRV) panchromatic data. The roads generated are compared with a visual interpretation of the SPOT HRV multispectral and panchromatic data. The technique is most effective in areas where road development is relatively recent. This is due to the spectral consistency of new road networks. As new road networks are those of most interest to the land manager, this is a useful technique for updating digital road network files within a geographical information system of urban areas.     

SHORELINE CHANGE AT THOMPSON ISLAND,BOSTON HARBOR,MASSACHUSETTS, 1938-1977

A photogrammetric and sediment analysis is presented to illustrate the relationship between beach erosion and seacliff recession on Thompson Island, Boston Harbor, Massachusetts. Aerial photographs taken in 1938, 1952, 1963, and 1977 were measured to determine rates of shoreline change around the island. The 39-year average rate of beach erosion is 0.3 m/yr ± 2% with an average rate of cliff recession at 0.2 m/yr ± 2%. Rates of beach erosion between six orientations that reflect principal wave approach to the island were not found to be significantly different (0.05) for the 39-year period. The rank order correlation between beach erosion and cliff recession for these six orientations was moderately well correlated (0.63). To determine possible controls of cliff erosion, 58 sediment samples were collected from glacial cliffs along the shoreline. The textural composition was determined, and then tested with discriminant function analysis. Partial correlation analysis between beach erosion and cliff recession holding a surrogate for sediment size constant improved the rank order from 0.63 to 0.84. The results indicate that coarser-grained cliffs recede at faster rates, but with less erosion occurring on adjacent beaches. Alternatively, finer-grained cliffs recede at slower rates, but with greater erosion occurring on adjacent beaches.     

Predicting suitable habitat for the cold-water coral Lophelia pertusa (Scleractinia)

Ecological-niche factor analysis (ENFA) was applied to the reef framework-forming cold-water coral Lophelia pertusa. The environmental tolerances of this species were assessed using readily available oceanographic data, including physical, chemical, and biological variables. L. pertusa was found at mean depths of 468 and 480 m on the regional and global scales and occupied a niche that included higher than average current speed and productivity, supporting the theory that their limited food supply is locally enhanced by currents. Most records occurred in areas with a salinity of 35, mean temperatures of 6.2–6.7  °C and dissolved oxygen levels of 6.0–6.2 ml l⁻¹. The majority of records were found in areas that were saturated with aragonite but had low concentration of nutrients (silicate, phosphate, and nitrate). Suitable habitat for L. pertusa was predicted using ENFA on a global and a regional scale that incorporated the north-east Atlantic Ocean. Regional prediction was reliable due to numerous presence points throughout the area, whereas global prediction was less reliable due to the paucity of presence data outside of the north-east Atlantic. However, the species niche was supported at each spatial scale. Predicted maps at the global scale reinforced the general consensus that the North Atlantic Ocean is a key region in the worldwide distribution of L. pertusa. Predictive modelling is an approach that can be applied to cold-water coral species to locate areas of suitable habitat for further study. It may also prove a useful tool to assist spatial planning of offshore marine protected areas. However, issues with eco-geographical datasets, including their coarse resolution and limited geographical coverage, currently restrict the scope of this approach.