A multiobjective steepest descent method with applications to optimal well control

Multiobjective optimization deals with mathematical optimization problems where two or more objective functions (cost functions) are to be optimized (maximized or minimized) simultaneously. In most cases of interest, the objective functions are in conflict, i.e., there does not exist a decision (design) vector (vector of optimization variables) at which every objective function takes on its optimal value. The solution of a multiobjective problem is commonly defined as a Pareto front, and any decision vector which maps to a point on the Pareto front is said to be Pareto optimal. We present an original derivation of an analytical expression for the steepest descent direction for multiobjective optimization for the case of two objectives. This leads to an algorithm which can be applied to obtain Pareto optimal points or, equivalently, points on the Pareto front when the problem is the minimization of two conflicting objectives. The method is in effect a generalization of the steepest descent algorithm for minimizing a single objective function. The steepest-descent multiobjective optimization algorithm is applied to obtain optimal well controls for two example problems where the two conflicting objectives are the maximization of the life-cycle (long-term) net-present-value (NPV) and the maximization of the short-term NPV. The results strongly suggest the multiobjective steepest-descent (MOSD) algorithm is more efficient than competing multiobjective optimization algorithms.     

Assessing Winter Wheat Responses to Climate Change Scenarios: A Simulation Study in the U.S. Great Plains

Hard red winter wheat (Triticum aestivum L.) is a major crop in the Great Plains region of the U.S. The goal of this assessment effort was to investigate the influence of two contrasting global climate change projections (U.K. Hadley Center for Climate Prediction and Research and Canadian Centre for Climate Modelling and Analysis) on the yield and percent kernel nitrogen content of winter wheat at three locations in Nebraska. These three locations represent sub-humid and semi arid areas and the transition between these areas and are also representative of major portions of the winter wheat growing areas of the central Great Plains. Climate scenarios based on each of the projections for each location were developed using the LARS-WG weather generator along with data from automated weather stations. CERES-Wheat was used to simulate the responses for two contrasting cultivars of wheat using two sowing dates. The first sowing date represented current sowing dates appropriate for each location. The second sowing date was later and represents the approximate date when the mean air temperature from the climate scenarios is the same as the mean air temperature from the actual climate data at the current sowing dates. The yield and percent kernel nitrogen content using the two climate scenarios generally decrease going from the sub-humid eastern to the semi arid western parts of Nebraska. Results from these simulations indicate that yield and percent kernel nitrogen content using the two climate scenarios could not both be maintained at levels currently simulated. Protein content (directly related to kernel nitrogen content) and end-use quality are the primary determinants for the use of hard red winter wheat in baked goods. Nitrogen management and new cultivars, which can enhance the uptake and translocation of nitrogen, will be proactive steps to meet the challenges of global climate change as represented by these climate scenarios.     

Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes,Boundary Conditions and Spin-Up

We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about \(2\) K) and substantially underestimates it at night (about \(4\) K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( \(\approx \) 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.     

Uranium speciation and stability after reductive immobilization in aquifer sediments

It has generally been assumed that the bioreduction of hexavalent uranium in groundwater systems will result in the precipitation of immobile uraninite (UO₂). In order to explore the form and stability of uranium immobilized under these conditions, we introduced lactate (15 mM for 3 months) into flow-through columns containing sediments derived from a former uranium-processing site at Old Rifle, CO. This resulted in metal-reducing conditions as evidenced by concurrent uranium uptake and iron release. Despite initial augmentation with Shewanella oneidensis, bacteria belonging to the phylum Firmicutes dominated the biostimulated columns. The immobilization of uranium (∼1 mmol U per kg sediment) enabled analysis by X-ray absorption spectroscopy (XAS). Tetravalent uranium associated with these sediments did not have spectroscopic signatures representative of U-U shells or crystalline UO₂. Analysis by microfocused XAS revealed concentrated micrometer regions of solid U(IV) that had spectroscopic signatures consistent with bulk analyses and a poor proximal correlation (μm scale resolution) between U and Fe. A plausible explanation, supported by biogeochemical conditions and spectral interpretations, is uranium association with phosphoryl moieties found in biomass; hence implicating direct enzymatic uranium reduction. After the immobilization phase, two months of in situ exposure to oxic influent did not result in substantial uranium remobilization. Ex situ flow-through experiments demonstrated more rapid uranium mobilization than observed in column oxidation studies and indicated that sediment-associated U(IV) is more mobile than biogenic UO₂. This work suggests that in situ uranium bioimmobilization studies and subsurface modeling parameters should be expanded to account for non-uraninite U(IV) species associated with biomass.     

Statewide groundwater‐vulnerability assessment in nebraska using the drastic/GIS model

The paper summarizes a technique for implementing the “DRASTIC” groundwater‐vulnerability model within the context of an automated raster‐based geographic information system. Discussion focuses on a methodological development and a statewide project completed recently in Nebraska. The final products, a comprehensive flow chart illustrating procedures and a map of calculated potential pollution hazard, are presented. The methodology can be executed successfully with minimal training and experience. Areas of Nebraska considered vulnerable to groundwater pollution are identified.     

Solving Kepler's equation with high efficiency and accuracy

We present a method for solving Kepler's equation for elliptical orbits that represents a gain in efficiency and accuracy compared with those currently in use. The gain is obtained through a starter algorithm which uses Mikkola's ideas in a critical range, and less costly methods elsewhere. A higher-order Newton method is used thereafter. Our method requires two trigonometric evaluations.     

The Taconite Inlet Lakes Project: a systems approach to paleoclimatic reconstruction

A comprehensive study of meteorological, hydrological, limnological and sedimentological conditions in the watersheds of density-stratified (meromictic) lakes around Taconite Inlet, Northern Ellesmere Island, N.W.T., Canada was carried out from 1990–1992. Lakes C1 and C2 contain seawater trapped by isostatic uplift as the former embayments became isolated from the sea. These lakes, and Lake C3, contain varved sediments which provide an annually resolvable paleoclimatic record. By studing the major systems influencing sedimentation in one of these lakes (Lake C2) a better understanding of the climatic controls on varve formation, and hence on the paleoclimatic signal in the varved sediment record, was obtained. The varves of Lake C2 provide a proxy record of summer temperature for the region.This is the first in a series of papers published in this issue on the Taconite Inlet Lakes Project. These papers were collected by Dr R. S. Bradley.     

Response of low‐angle dunes to variable flow

Current understanding of bedform dynamics is largely based on field and laboratory observations of bedforms in steady flow environments. There are relatively few investigations of bedforms in flows dominated by unsteadiness associated with rapidly changing flows or tides. As a consequence, the ability to predict bedform response to variable flow is rudimentary. Using high‐resolution multibeam bathymetric data, this study explores the dynamics of a dune field developed by tidally modulated, fluvially dominated flow in the Fraser River Estuary, British Columbia, Canada. The dunes were dominantly low lee angle features characteristic of large, deep river channels. Data were collected over a field ca 1·0 km long and 0·5 km wide through a complete diurnal tidal cycle during the rising limb of the hydrograph immediately prior to peak freshet, yielding the most comprehensive characterization of low‐angle dunes ever reported. The data show that bedform height and lee angle slope respond to variable flow by declining as the tide ebbs, then increasing as the tide rises and the flow velocities decrease. Bedform lengths do not appear to respond to the changes in velocity caused by the tides. Changes in the bedform height and lee angle have a counterclockwise hysteresis with mean flow velocity, indicating that changes in the bedform geometry lag changes in the flow. The data reveal that lee angle slope responds directly to suspended sediment concentration, supporting previous speculation that low‐angle dune morphology is maintained by erosion of the dune stoss and crest at high flow, and deposition of that material in the dune trough.     

Tendances et ruptures au niveau des pluies journalières extrêmes en Côte d'Ivoire (Afrique de l'Ouest)

Résumé

Le changement climatique est une réalité qui affecte plusieurs variables climatiques dont les précipitations. Néanmoins, son impact sur les évènements extrêmes et en particulier sur les pluies journalières extrêmes n'est pas encore certain car peu de travaux y ont été consacrés en Afrique de l'Ouest. Dans ce contexte, il a été proposé de détecter d'éventuels tendances et ruptures dans les propriétés statistiques (moyenne, variance) des pluies journalières extrêmes à l'aide de tests statistiques locaux et régionaux. Pour détecter ces changements, les indices caractérisant la pluie maximale journalière annuelle (PJmaxan), le nombre annuel de jours de pluie dépassant 50 mm (NJsup50) et la contribution des pluies dépassant 50 mm dans les cumuls annuels (R(PJsup50/Pan)) ont été définis. L'analyse de 44 postes pluviométriques en Côte d'Ivoire sur la période 1942–2002 ne montre pas de changement généralisé ni en moyenne, ni en variance. Toutefois, en subdivisant la Côte d'Ivoire en régions climatiques homogènes, des tendances à la baisse ont été observées dans les régions IV (au Nord) et II (au Sud-Est).

Editeur Z.W. Kundzewicz; Assistant editeur G. Mahé

Citation Goula, A.B.T., Gneneyougo Soro, E., Kouassi, W. et Srohourou, B., 2012. Tendances et ruptures au niveau des pluies journalières extrêmes en Côte d'Ivoire (Afrique de l'Ouest). Hydrological Sciences Journal, 57 (6), 1067–1080.     

Sedimentation in Lake Biel,an eutrophic,hard-water lake in northwestern Switzerland

A 2-year study of sedimentation in Lake Biel reveals that three major factors can account for the spatial and temporal patterns of tripton sedimentation. Allochthonous inputs of particulate matter mostly derived from the Aare river during spring snowmelt and other high-water periods, provide about 50% of the carbonate and 50% of the siliceous matter collected in sediment traps. Phytoplankton activity in this eutrophic, hard-water lake accounts for the rest of the carbonate and siliceous matter, the latter being mostly diatom frustules. Sedimentation rates are thus highest during periods of maximum river flooding, maximum phytoplankton activity and following the breakdown of the thermocline, at which time trapped particles settle out. Resuspension of sediment is important in at least one basin. Resolution of carbonate appears to be minor.