Multi-objective optimization for free-phase LNAPL recovery using evolutionary computation algorithms

Abstract

A nonlinear, multi-objective optimization methodology is presented that seeks to maximize free product recovery of light non-aqueous phase liquids (LNAPLs) while minimizing operation cost, by introducing the novel concept of optimal alternating pumping and resting periods. This process allows more oil to flow towards the extraction wells, ensuring maximum free product removal at the end of the remediation period with minimum groundwater extraction. The methodology presented here combines FEHM (Finite Element Heat and Mass transfer code), a multiphase groundwater model that simulates LNAPL transport, with three evolutionary algorithms: the genetic algorithm (GA), the differential evolution (DE) algorithm and the particle swarm optimization (PSO) algorithm. The proposed optimal free-phase recovery strategy was tested using data from a field site, located near Athens, Greece. The PSO and DE solutions were very similar, while that provided by the GA was inferior, although the computation time was roughly the same for all algorithms. One of the most efficient algorithms (PSO) was chosen to approximate the optimal Pareto front, a method that provides multiple options to decision makers. When the optimal strategy is implemented, although a significant amount of LNAPL free product is captured, a spreading of the LNAPL plume occurs.

Editor Z.W. Kundzewicz; Associate editor L. See

Citation Dokou, Z. and Karatzas, G.P., 2013. Multi-objective optimization for free-phase LNAPL recovery using evolutionary computation algorithms. Hydrological Sciences Journal, 58 (3), 671–685.     

Hindcasting of tropical cyclone winds and waves

Ocean Dynamics - This paper describes a series of hindcast simulations of 17 tropical cyclones over the northwest shelf region of Australia. Tropical cyclone track and vortex details were obtained...     

Hindcasting multidecadal predevelopment groundwater levels in the Floridan aquifer

Establishing predevelopment benchmark groundwater conditions is challenging without long-term records to discern impacts of pumping and climate change on aquifer levels. Understanding periodic natural cycles and trends require 100 years or more data which rarely exist. Using limited records, we develop an approach to hindcast multidecadal levels and examine the temporal evolution of climatic and pumping impacts. The methodology includes a wavelet-aided statistical model, constrained by temporal scales of physical processes responsible for groundwater level variation, including rainfall, evapotranspiration and pumping stresses. The model and hindcasts are tested at three sites in Florida using traditional split calibration-verification methods for the period of record and with the documented historical drought and wet years for the period of no-record. The pumping impact is quantified over time and compared with regional groundwater models, revealing that withdrawals are responsible for 30 to 70% of the declines in levels since 1960s. Hindcasting yielding 110 years of monthly levels is used to assess the effect of climate change and pumping on the frequency of critical low levels. At all three sites, the frequencies of critical low levels increase significantly in the 1960 to 2015 period when compared to the 1904 to 1959 period. For example, at site 1, the return period of the critical low level is shortened by 3.9 years due to climate change and 2.2 years due to pumping.     

The importance of ocean temperature to global biogeochemistry

Variations in the mean temperature of the ocean, on time scales from millennial to millions of years, in the past and projected for the future, are large enough to impact the geochemistry of the carbon, oxygen, and methane geochemical systems. In each system, the time scale of the temperature perturbation is key. On time frames of 1-100 ky, atmospheric CO₂ is controlled by the ocean. CO₂ temperature-dependent solubility and greenhouse forcing combine to create an amplifying feedback with ocean temperature; the CaCO₃ cycle increases this effect somewhat on time scales longer than ∼5-10 ky. The CO₂/T feedback can be seen in the climate record from Vostok, and a model including the temperature feedback predicts that 10% of the fossil fuel CO₂ will reside in the atmosphere for longer than 100 ky. Timing is important for oxygen, as well; the atmosphere controls the ocean on short time scales, but ocean anoxia controls atmospheric pO₂ on million-year time scales and longer. Warming the ocean to Cretaceous temperatures might eventually increase pO₂ by approximately 25%, in the absence of other perturbations. The response of methane clathrate to climate change in the coming century will probably be small, but on longer time scales of 1-10 ky, there may be a positive feedback with ocean temperature, amplifying the long-term climate impact of anthropogenic CO₂ release.     

Spatiotemporal modes of global sea surface temperature variability

This study identifies the salient global modes of sea surface temperature variability based on 145 years of HadlSST data.Unlike the traditional mode identification by EOF analysis,a combination of wavelet and EOF analysis is used to extract the leading modes at distinct time scales.The spatial patterns of some well-known regional modes are recovered,with the global connection and frequency content of these modes being revealed.Our analysis indicates that,in terms of global influence,the Pacific Ocean is the major player,and the tropical Pacific is the center of action on various time scales.The Atlantic Ocean has its own outstanding modes,but their global impacts are not as strong as those from the Pacific.The Indian Ocean generally shows a passive response to the Pacific,with a basin-wide pattern in the tropics.Despite some preliminary theoretical attempts,how to elucidate the dynamics underlying the global modes of sea surface temperature variability is still an open question.     

An investigation on stream temperature analysis based on evolutionary computing

The data‐driven technique, evolutionary polynomial regression, has been tested and used for the study of water temperature behaviour in the River Barle (south‐west England). The study aimed to produce multiple models for forecasting water temperature, using air temperature as input. In addition, river discharge data were used to describe the hydrological regime of the study stream, even if they are not involved in the modelling phase. The availability of data sampled at hourly intervals allowed behaviour to be studied at several time scales, including short‐term lags between air temperature and water temperature. The approach to model building differs from previous studies in that the relationship between air temperature and water temperature is not evaluated on the basis of a multi‐parameter regression, nor does it identify particular structures; rather the evolutionary technique identifies the model by itself. In fact, the non‐linear relationship between air temperature and water temperature is investigated by an evolutionary search in the space of particular pseudo‐polynomials structures. Copyright © 2007 John Wiley & Sons, Ltd.     

Air temperature feedback and its contribution to global warming

Air temperature feedback results from the thermal-radiative coupling between the atmosphere and the surface and plays an important role in surface energy balance. This paper reveals the contribution of air temperature feedback to the global warming from 1980 to 2000. The air temperature feedback kernel, evaluated using the ERA-Interim reanalysis data, is used to discuss the physical mechanism for air temperature feedback, the dependency of the strength of air temperature feedback on the climatological spatial distributions of air temperature, water vapor and cloud content, and the contributions of air temperature feedback to rapid global warming. The coupling between temperature feedback and each of the external forcings and individual feedback processes will amplify the anomaly of direct energy flux convergence at the surface induced by the external forcings and individual processes. The air temperature feedback amplifies the initial surface warming due to the increase in CO₂ concentration, ice and snow melting, increase in water vapor, and change in ocean heat storage. It also amplifies the surface warming due to the longwave radiaitve forcing associated with the increase in cloud cover, which acts to suppress the cooling of the shortwave effect of cloud forcing. Overall, temperature feedback plays an important role in the global warming from 1980 to 2000, as the net positive contribution to the perturbation of global mean energy flux at the surface from the air temperature feedback is larger than the net negative contribution from external forcing and all non-temperature feedbacks.     

Traveltime computation by wavefront-orientated ray tracing

For multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts.     

Relationship between global mean sea-level and global mean temperature in a climate simulation of the past millennium

The possibility of using global mean near-surface temperature, its rate of change or the global mean ocean heat-flux as predictors to statistically estimate the change of global mean sea-level is explored in the context of a long climate simulation of the past millennium with the climate model ECHO-G. Such relationships have recently been proposed to by-pass the difficulty of estimating future sea-level changes based on simulations with coarse-resolution climate models. It is found that, in this simulation, a simple linear relationship between mean temperature and the rate of change of sea level does not exist. A regression parameter linking both variables, and estimated in sliding 120-year windows, varies widely along the simulation and, in some periods, even attains negative values. The ocean heat-flux and the rate-of-change of mean temperature seem to better capture the rate-of-change of sea level due to thermal expansion.     

The computation of the free barotropic oscillations of a global ocean model including friction and loading effects

A set of 61 normal modes with periods between 7.8 and 133.1 h has been calculated, using a 1° model of the global ocean, including the Arctic Ocean. The model explicitly considers frictional forces and ocean self-attraction and loading effects. The latter effects have generally been taken into account by parameterization, but for some modes the effects have also been considered fully. Due to friction, the computed eigenfrequencies are complex, exhibiting also the varying dissipative properties of the modes and their dependence on the distribution of potential and kinetic energies over the oceanic regions. In detail, gravity modes having periods less then 80 h and dominating the semi-diurnal and the diurnal tides, topographically controlled vorticity modes with periods longer than diurnal, and two planetary vorticity modes with periods of 96.8 and 119.4 h have been identified. These planetary vorticity modes have their energies distributed over Pacific, Atlantic, and Indian Oceans, while the other modes with periods longer than 80 h, as vorticity modes, have their energies concentrated on topographic structures of restricted extension. The modes are discussed with respect to their wave properties, e.g., concerning quasi-standing-wave resonances and to the appearance of Kelvin waves of different orders and trapped by different coastlines. In particular, the relevance of specific modes for the development of the fields of the most important semi-diurnal and diurnal tidal constituents is investigated.     

Multi-resolution analysis of global surface air temperature and solar activity relationship

In this work the surface temperature anomaly (dTG) and sunspot number (R_z) time series in the period 1880–2000 are studied with wavelet multi-resolution analysis. We found a very low correlation of 0.11 between dTG and R_z in the 11-yr-solar cycle band. A higher correlation of 0.66 is found in the ～22-yr-band with zero lag correlation coefficient between dTG and R_z. Furthermore, the long-term trend is markedly different between dTG and R_z. This might occurs because of the long-term warming on the last century, which is attributed mainly to anthropogenic effects.     

Changes in geomagnetic activity and global temperature during the past 40 years

The purpose of this study is to investigate the effect of geomagnetic activity (used as a measure of solar wind parameters) on the variability of large-scale climate patterns and on changes in the global temperature. We show that positive statistically significant correlations between global temperature and the distribution of surface temperature over Eurasia, the East and Equatorial Pacific and over the North Atlantic for the period 1966?C2009 correspond to large-scale climate patterns defined by climate indices. We found very similar positive correlations between geomagnetic activity and the distribution of surface temperature in the mentioned regions. As an effect of geomagnetic storms, energetic particles penetrate from the magnetosphere into the region of the stratospheric polar vortex. The increase of temperature and pressure can be observed over northern Canada. The vortex shifts towards Europe, rotates counter-clockwise and the wind blows from the polar region over Greenland southwards. It diverts the warm flow proceeding northward over the Atlantic, eastward along the deep Icelandic low extending as far as the Barents Sea and takes part in warming Eurasia. The strengthened zonal flow from Siberia cools the western Pacific with the impact on the warming of the equatorial and eastern Pacific when also a distinct 1976?C78 climate shift occurred. Processes in the Atlantic and Pacific play a significant role and a time delay (wind forcing over the previous 1?C4 yr) appears to be the most important for the relocation of the oceanic gyres. Results showing statistically significant relations between time series for geomagnetic activity, for the sum of climate indices and for the global temperature help to verify findings concerning the chain of processes from the magnetosphere to the troposphere.     

Role of the warming trend in global land surface air temperature variations

Global warming and its climatic and environmental effects have mainly been investigated in terms of the absolute warming rate. Little attention has been paid to the contribution of absolute warming rate to variability on various time scales of surface air temperature(SAT), which may be a more direct index for measuring the ecoclimatic effect of warming trend. The present study analyzed the role of secular warming trend in the variations of global land SAT for 1901–2016. Less than one-third of annual SAT variations were contributed by the warming trend over large parts of the globe generally. The ratios were up to two-thirds over eastern South America, parts of South Africa and the regions around the southwestern Mediterranean and Sunda islands where the absolute warming rate was moderate but the endemic species were undergoing exceptional loss of habitat. The ratios also exhibited smallest seasonal difference over these regions. Therefore, the ratio of the warming trend to the SAT variations may be a better measure compared to the absolute warming rate for the local ecoclimate. We should also pay more attention to the regions with high ratio, not only the regions with the high absolute warming rate.     

Limitations imposed by nonlinear phenomena on computation of acoustic paths

Approaches for estimating the influence of nonlinear effects on long-distance propagation of acoustic signals in an isothermal atmosphere are developed. It has been indicated that nonlinear effects can be ignored on paths up to 3000 km at beam angles smaller than 35° (nonlinear absorption). However, the spectrum will be pronouncedly distorted at these distances due to the generation of new harmonics (for initial Mach numbers about 10⁻⁴).     

CMIP6模式模拟的人类活动和自然强迫对全球地表气温多尺度变化的影响

量化人类活动在气候变化中的定量贡献是气候变化检测归因研究的核心科学问题,也是提高气候变化预测和预估水平的重要科学基础.本文基于最新的第6次国际耦合模式比较计划(CMIP6)多模式历史归因模拟试验(DAMIP),检测了人为因素(ANT)和自然因素(NAT)对近百年(1915—2014年)全球地表气温多尺度变化的影响,归因了温室气体(GHG)、气溶胶(AA)、土地利用(LU)等不同人为因素对全球地表气温长期变化的相对贡献及南北半球差异.结果显示,近百年来人为因素引起的全球陆地实际增温约为1.1℃(0.8℃~1.3℃),对南北半球的贡献则分别约为0.7℃和1.2℃;全球大多数区域人为排放GHG和AA的显著作用在1960—1980年期间就能够被检测到,其中北半球AA的冷却作用要超前于GHG的增温效应;气候系统内部自然变率是调制大多数区域气温年代际(10~30年)及多年代际变率(30~60年)的主导因子,而人为和自然外强迫在全球地表气温年代际变率中的方差贡献约为5%~20%,但二者在北半球尤其在东亚和欧洲中高纬度地区地表气温多年代际变率中的方差贡献可达50%.人为因素强迫可使近50年(1965—2014)极端高温事件的发生风险概率增加3倍,其中北半球发生风险要高于南半球.

     

Detecting global atmospheric oscillations by seismic instruments

Izvestiya, Physics of the Solid Earth - We briefly overview the periodic atmospheric motions in the period range from 1 h to a few dozen days. It is shown that many oscillations recorded by...     

Hydroclimatic teleconnection between global sea surface temperature and rainfall over India at subdivisional monthly scale

It is well established that sea surface temperature (SST) plays a significant role in the hydrologic cycle in which precipitation is the most important part. In this study, the influence of SST on Indian subdivisional monthly rainfall is investigated. Both spatial and temporal influences are investigated. The most influencing regions of sea surface are identified for different subdivisions and for different overlapping seasons in the year. The relative importance of SST, land surface temperature (LST) and ocean–land temperature contrast (OLTC) and their variation from subdivision to subdivision and from season to season are also studied. It is observed that LST does not show much similarity with rainfall series, but, in general, OLTC shows relatively higher influence in the pre‐monsoon and early monsoon periods, whereas SST plays a more important role in late‐ and post‐monsoon periods. The influence of OLTC is seen to be mostly confined to the Indian Ocean region, whereas the effect of SST indicates the climatic teleconnection between Indian regional rainfall and climate indices in Pacific and Atlantic Oceans. Copyright © 2006 John Wiley & Sons, Ltd.     

Long-term global sea surface temperature fluctuations and their possible influence on atmospheric CO2 concentrations

Changes of atmospheric CO₂ concentration since 1958 are shown to be related to sea surface temperature changes. The largest contribution to changes arises from the Pacific equatorial upwelling region, with the Indian Ocean and Atlantic contributing only small fractions to the variance. It is hypothesized that the observed relationship is related to the nutrients that are brought up by upwelling cold water, with photosynthesis contributing to a lowering of the partial pressure of CO₂ in the sea and thus to a greater tendency for a flux from the air to the sea.Possible longer term variations of sea temperature and CO₂ are discussed.