Hybrid-optimization algorithm for the management of a conjunctive-use project and well field design

Hi-Desert Water District (HDWD), the primary water-management agency in the Warren Groundwater Basin, California, plans to construct a waste water treatment plant to reduce future septic-tank effluent from reaching the groundwater system. The treated waste water will be reclaimed by recharging the groundwater basin via recharge ponds as part of a larger conjunctive-use strategy. HDWD wishes to identify the least-cost conjunctive-use strategies for managing imported surface water, reclaimed water, and local groundwater. As formulated, the mixed-integer nonlinear programming (MINLP) groundwater-management problem seeks to minimize water-delivery costs subject to constraints including potential locations of the new pumping wells, California State regulations, groundwater-level constraints, water-supply demand, available imported water, and pump/recharge capacities. In this study, a hybrid-optimization algorithm, which couples a genetic algorithm and successive-linear programming, is developed to solve the MINLP problem. The algorithm was tested by comparing results to the enumerative solution for a simplified version of the HDWD groundwater-management problem. The results indicate that the hybrid-optimization algorithm can identify the global optimum. The hybrid-optimization algorithm is then applied to solve a complex groundwater-management problem. Sensitivity analyses were also performed to assess the impact of varying the new recharge pond orientation, varying the mixing ratio of reclaimed water and pumped water, and varying the amount of imported water available. The developed conjunctive management model can provide HDWD water managers with information that will improve their ability to manage their surface water, reclaimed water, and groundwater resources.     

Arsenic Management Through Well Modification and Simulation

Arsenic concentrations can be managed with a relatively simple strategy of grouting instead of completely destroying a selected interval of well. The strategy of selective grouting was investigated in Antelope Valley, California, where groundwater supplies most of the water demand. Naturally occurring arsenic typically exceeds concentrations of 10 µg/L in the water produced from these long-screened wells. The vertical distributions of arsenic concentrations in intervals of the aquifer contributing water to selected supply wells were characterized with depth-dependent water-quality sampling and flow logs. Arsenic primarily entered the lower half of the wells where lacustrine clay deposits and a deeper aquifer occurred. Five wells were modified by grouting from below the top of the lacustrine clay deposits to the bottom of the well, which reduced produced arsenic concentrations to less than 2 µg/L in four of the five wells. Long-term viability of well modification and reduction of specific capacity was assessed for well 4-54 with AnalyzeHOLE, which creates and uses axisymmetric, radial MODFLOW models. Two radial models were calibrated to observed borehole flows, drawdowns, and transmissivity by estimating hydraulic-conductivity values in the aquifer system and gravel packs of the original and modified wells. Lithology also constrained hydraulic-conductivity estimates as regularization observations. Well encrustations caused as much as 2 µg/L increase in simulated arsenic concentration by reducing the contribution of flow from the aquifer system above the lacustrine clay deposits. Simulated arsenic concentrations in the modified well remained less than 3 µg/L over a 20-year period.     

Analysis of dust events in 2008 and 2009 using the lidar network, surface observations and the CFORS model

The Asian dust events in 2008 (May 24–June 4 in 2008) and in 2009 (March 12–25, October 13–26, and December 15–28 in 2009) were analyzed with the lidar network observations, surface observations in China, Korea, Japan, and Mongolia, and with the chemical transport model CFORS. Transport of Asian dust and mixing of dust with air pollution aerosols were studied. The event of May 24 to June 4 in 2008 was a significant event unusually late in the spring dust season. The dust event of March 12–25, 2009 was an interesting example of elevated dust layer, and transport of dust from the elevated dust layer to the ground by the boundary layer activity was observed with the lidars and surface observations in Japan. The concentration of air pollution aerosols was relatively high during the dust event, and the results suggest that vertical structure as well as transport path is important for the mixing of dust and air pollution aerosols. The dust events in October and December 2009 were examples of dust events in autumn and winter. The online mode CFORS reproduced the observation data generally well, except for the event of May 24 to June 4 in 2008. The results of the fourdimensional variational assimilation of the lidar network data reproduced the dust concentration in Korea and Japan reasonably in that event.     

A model of ‘disparitions brusques’ (sudden disappearance of eruptive prominences) as an instability driven by mhd waves

A model of ‘disparitions brusques’ (sudden disappearence of eruptive prominences) is discussed based on the Kippenhahn ans Schlüter configuration. It is shown that Kippenhahn and Schlüter's current sheet is very weakly unstable against magnetic reconnecting modes during the lifetime of quiescent prominences. Disturbances in the form of fast magnetosonic waves originating from nearby active regions or the changes of whole magnetic configuration due to newly emerged magnetic flux may trigger a rapid growing instability associated with magnetic field reconnection. This instability gives rise to disruptions of quiescent prominences and also generates high energy particles.     

Origin of the weakening of EUV emission lines formed in the chromosphere-corona transition zone

EUV spectroheliograms of the quiet Sun obtained with the Harvard experiment on Skylab are analyzed to identify the structure causing the weakening of the EUV line emission due to Lyman continuum absorption. The weakening at the network boundaries can be explained by overlapping of several spicules each of which being wrapped in an EUV emitting sheath. Part of the cell interiors show moderate weakening, this has the shape of a belt surrounding the network boundaries. There are a number of patches showing intense weakening near network boundaries and in cell interiors; the weakening at the points cannot be explained by overlapping of chromospheric structures with an EUV emitting sheath. A possible explanation is that the intense weakening is caused by cool chromospheric clouds or moving blobs over the EUV emitting sources in the cell interiors. Some of the points showing intense weakening are associated with an enhancement of the EUV emission. These points have lifetime shorter than the time interval of 5.5 min between successive observations, which stresses again that the chromosphere-corona transition zone is in a dynamic state.     

Excitation of whistler waves driven by an electron temperature anisotropy

A 3-D particle simulation of excitation of whistler waves driven by an electron temperature anisotropy (T > T ) is presented. Results show that whistler waves can have appreciable growth driven by the anisotropy. The maximum intensity of the excited whistler waves increases as a quadratic function of the anisotropy. Due to the presence of a threshold, one needs a relatively large electron temperature anisotropy above threshold to generate large-amplitude whistler waves. The average amplitude of turbulence in the context of whistler waves is up to as large as about 1% of the ambient magnetic field when T /T . The total energy density of the whistler turbulence is adequate for production of relativistic electrons in solar flares through stochastic acceleration.     

Biomass, Abundance, and Vertical Distribution of Micronekton and Large Gelatinous Zooplankton in the Subarctic Pacific and the Bering Sea during the Summer of 1997

The biomass, abundance, and vertical distribution of micronekton, including enidarians, mysids, euphausiids, decapods, thaliaceans, and fishes, were studied on the basis of samples collected with an 8-m² opening-closing rectangular midwater trawl (RMT-8, mesh size: 4.5 mm) at three stations in the subarctic Pacific (the western subarctic gyre, the central Subarctic, and the Gulf of Alaska) and one station in the oceanic Bering Sea. The total biomass in the 0–1000 m water column ranged from 2.9 to 5.1 gDW m^–2. Except for primary consumers that showed highly variable biomass (thaliaceans and euphausiids), biomass was highest in the oceanic Bering Sea followed by the central (boundary between eastern and western gyres), western gyre, and eastern Gulf of Alaska. The biomass compositions by higher taxa were basically similar between regions: fishes were most dominant, followed by enidarians at all stations, except for the marked predominance of thaliaceans in the Gulf of Alaska. High biomasses of gelatinous animals (31% of overall dry weight), occasionally comparable to those of fishes and crustaceans, suggest their potential importance in the subarctic Pacific. Characteristics in vertical patterns of micronekton biomass common in all stations were: (1) a mesopelagic peak around 500–600 m both day and night, (2) a layer of low biomass in the cold intermediate water and/or in the upper mesopelagic zone, (3) a nighttime shift of biomass to upper layers, and (4) an highly variable biomass of epipelagic/interzonal migrants (euphausiids and thaliaceans).     

Nutrient and Phytoplankton Dynamics in Harima-Nada, Eastern Seto Inland Sea, Japan During a 35-Year Period from 1973 to 2007

Long-term monitoring of water quality and phytoplankton was conducted at 19 sampling stations in Harima-Nada, eastern Seto Inland Sea, Japan for 35 years from 1973 to 2007. There were two significant long-term changes, an increase in winter water temperatures of 0.042°C year^?1, and a decrease in dissolved inorganic nitrogen (DIN) from about 10 μM in the 1970s to ~5 μM in the late 1990s due to the reduction in nutrient inputs. DIN concentrations and total phytoplankton cell density were both higher during the 1970s to the early 1980s and then exhibited a significant decrease in the mid 1980s and remained relatively constant thereafter. Diatoms were the dominant phytoplankton group (>90%) over the 35-year period, and there was a dramatic shift from Skeletonema dominance (~70%) to Chaetoceros in the mid 1980s. This shift in diatom species may be attributed to differences in the life cycle of Skeletonema and Chaetoceros and the response to the decrease in DIN concentration.     

3-D Rpic Simulations of Relativistic Jets: Particle Acceleration, Magnetic Field Generation, and Emission

We have applied numerical simulations and modeling to the particle acceleration, magnetic field generation, and emission from relativistic shocks. We investigate the nonlinear stage of theWeibel instability and compare our simulations with the observed gamma-ray burst emission. In collisionless shocks, plasma waves and their associated instabilities (e.g., the Weibel, Buneman and other two-stream instabilities) are responsible for particle (electron, positron, and ion) acceleration and magnetic field generation. 3-D relativistic electromagnetic particle (REMP) simulations with three different electron-positron jet velocity distributions and also with an electron-ion plasma have been performed and show shock processes including spatial and temporal evolution of shocks in unmagnetized ambient plasmas. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-positron shocks are larger for smaller jet Lorentz factor. This comes from the fact that the growth time of the Weibel instability is proportional to the square of the jet Lorentz factor. We have performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which is assumed to be created by photon annihilation. Simulation results with this broad distribution show that the Weibel instability is excited continuously by the wide-range of jet Lorentz factor from lower to higher values. In all simulations the Weibel instability is responsible for generating and amplifying magnetic fields perpendicular to the jet propagation direction, and contributes to the electron’s (positron’s) transverse deflection behind the jet head. This small scale magnetic field structure contributes to the generation of “jitter” radiation from deflected electrons (positrons), which is different from synchrotron radiation in uniform magnetic fields. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks. The detailed studies of shock microscopic process evolution may provide some insights into early and later GRB afterglows.