Reduced-order modeling for thermal recovery processes

Thermal recovery can entail considerably higher costs than conventional oil recovery, so the use of computational optimization techniques in designing and operating these processes may be beneficial. Optimization, however, requires many simulations, which results in substantial computational cost. Here, we implement a model-order reduction technique that aims at large reductions in computational requirements. The technique considered, trajectory piecewise linearization (TPWL), entails the representation of new solutions in terms of linearizations around previously simulated (and saved) training solutions. The linearized representation is projected into a low-dimensional space, with the projection matrix constructed through proper orthogonal decomposition of solution “snapshots” generated in the training step. Two idealized problems are considered here: primary production of oil driven by downhole heaters and a simplified model for steam-assisted gravity drainage, where water and steam are treated as a single “effective” phase. The strong temperature dependence of oil viscosity is included in both cases. TPWL results for these systems demonstrate that the method can provide accurate predictions relative to full-order reference solutions. Observed runtime speedups are very substantial, over 2 orders of magnitude for the cases considered. The overhead associated with TPWL model construction is equivalent to the computation time for several full-order simulations (the precise overhead depends on the number of training runs), so the method is only applicable if many simulations are to be performed.     

Indoor Air Background Levels of Volatile Organic Compounds and Air-Phase Petroleum Hydrocarbons in Office Buildings and Schools

A background indoor air study has been completed which includes the collection of indoor air samples from office buildings and schools. The anonymous study was designed with input from the U.S. Environmental Protection Agency and the Massachusetts Department of Environmental Protection. The sampling was implemented in 2013, 2014, and 2015 and included the collection of 25 school building samples and 61 office building samples. The study generated 14,668 new indoor air background data points, with samples collected from buildings located in 26 cities in 18 states, including Arizona, California, Connecticut, Indiana, Kansas, Maine, Massachusetts, Minnesota, Montana, New Hampshire, New Jersey, New York, Nevada, North Carolina, Ohio, Texas, Utah, and Washington. Indoor air background concentrations of target compound volatile organic compounds (VOCs) ranged from less than the laboratory method reporting limit of 0.044 μg/m³ to concentrations up to 1190 μg/m³, with hydrocarbon ranges from less than the reporting method limit of 10 μg/m³ to concentrations up to 3000 μg/m³. Some VOCs were identified ubiquitously in indoor air background, and some were identified at concentrations which exceeded risk-based regulatory screening levels. These study results provide useful and updated information on indoor air background and air quality in offices and schools and can be used in future regulatory guidance update considerations, for further examination of relationships between these data and residential study data, in human health risk assessments and risk communication, and in planning future studies.     

Partition of arsenic in soils sediments and the origin of naturally elevated concentrations in groundwater of the southern pampa region (Argentina)

Excessive arsenic concentrations above the Argentinean and WHO guidelines for drinking water (10 μg L⁻¹) affects shallow aquifers of the southern Pampean Plain (Argentina) hosted in the Pampean and the Post Pampean formations (loess and reworked loess; Plio-Pleistocene–Holocene). Health problems related to high As concentrations in drinking waters are known as Endemic Regional Chronic Hydroarsenicism. Hydrochemistry of shallow groundwaters and soil geochemistry were investigated aiming to (1) understand the partition of As in the solid phase and its relationship with unacceptable As concentrations in waters, (2) identify the provision source of As to groundwaters. Only 5% of the samples had As concentrations <10 μg L⁻¹; in 27% As concentrations ranged from 10 to 50 μg L⁻¹ and in 58% it reached 60–500 μg L⁻¹. The coarse fraction (50–2,000 μm) hosts about 27% of the total As in the solid phase, being positively correlated to Ba (p < 0.01; r ² = 0.93). About 70% is included in the <2 μm fraction and had positive correlations of As–Fe (p < 0.05; r ² = 0.85) and As–Cr (p < 0.05; r ² = 0.68). Soils and sediment sand fractions of vadose zones are the primary sources of As in shallow groundwater while adsorption–desorption processes, codisolution–coprecipitation, and evaporation during the dry seasons raise As concentrations in waters exceeding the guideline value for drinking water.     

The Catalan seismic crisis (1427 and 1428; NE Iberian Peninsula): Geological sources and earthquake triggering

The Catalan seismic crisis of the years 1427 and 1428 is one of the most destructive seismic episodes that happened in the northeastern Iberian Peninsula in historical times. The main earthquakes of this crisis occurred on March 19th 1427 in the zone around Amer (I_EMS-98 = VIII), May 15th 1427 in the vicinity of Olot (I_EMS-98 = VIII) and on February 2nd 1428 in the area close to Camprodon (I_EMS-98 = IX). There is much evidence that the Amer fault produced the first two events of this crisis, but is still uncertain which fault generated the earthquake on February 2nd 1428. Using newly available macroseismic data, the earthquake area sources of the three main earthquakes of the crisis have been obtained and they corroborate that the Amer fault may be the origin of the first two events. However, the area source corresponding to the last earthquake of the crisis cannot be associated to a single fault and indicates three possible candidates: the Vallfogona and Ribes-Camprodon thrusts and the Amer normal fault. Modeling of the Coulomb failure stress transfer has been performed to help determine the best candidate responsible for the February event. The results of the modeling points to: (a) a triggering relationship between the three main events of the crisis and (b) the Amer fault, or a similar extensional fault close and parallel to it, as the most probable origin of the earthquake on February 2nd 1428.     

Coherence between atmospheric teleconnections,Great Lakes water levels,and regional climate

We investigated the frequency domain relationships between four atmospheric teleconnections (Trans-Niño Index TNI, Pacific Decadal Oscillation PDO, Northern Annular Mode/Arctic Oscillation Index NAM/AO, and Pacific/North American PNA pattern) and water levels in the Great Lakes from 1948 to 2002 by quantifying the coherence between these time series. The levels in all Great Lakes are significantly correlated with the TNI in the frequency range (3–7)⁻¹ cycles year⁻¹, and with the PDO in interdecadal frequencies. The levels in Lakes Superior, Michigan, and Erie are significantly correlated with the PNA pattern in interdecadal frequencies, and the levels in all Great Lakes are significantly correlated with the NAM/AO in interannual frequencies.     

Bioerosion by the sea urchin Diadema mexicanum along Eastern Tropical Pacific coral reefs

Bioerosion is a natural process in coral reefs. It is fundamental to the health of these ecosystems. In the Eastern Tropical Pacific (ETP) coral reefs, the most important bioeroders are sponges, bivalves, sea urchins and the fish Arothron meleagris. In the 1980s, El Niño caused high coral mortality and an increase in macroalgal growth. As a result, greater sea urchin bioerosion occurred. This weakened the reef framework. Considering the high vulnerability of the ETP coral reefs, the goal of this study was to determine the current bioerosion impact of the sea urchin Diadema mexicanum along the western coasts of Mexico, El Salvador, Costa Rica and Panamá. The balance between coral bioaccretion and sea urchin bioerosion was also calculated. Between 2009 and 2010, in 12 coral reefs localities, D. mexicanum density, bottom cover and rugosity were quantified along band transects. The daily bioerosion rate was obtained from the amount of carbonates evacuated by sea urchins per unit time. The rate of coral accretion was calculated by multiplying the coral growth rate of the dominant genus by the density of their skeleton and by their specific coral cover. The localities were dissimilar (R = 0.765, P < 0.001) in terms of live coral cover, crustose calcareous algae, turf cover, rugosity index, and density and size of D. mexicanum. At all sites, with the exception of Bahía Culebra (Costa Rica), coral bioerosion was less than coral bioaccretion. Diadema mexicanum plays a dominant role in the balance of carbonates in the ETP, but this depends on reef condition (protection, overfishing, eutrophication) and so the impacts can be either positive or negative.     

The formation of whitings on the Little Bahama Bank

Whitings (patches of suspended fine-grained calcium carbonate) on the Great Bahama Bank (GBB) have been studied since the 1940s. The source and cause of these whitings have been hotly debated for a number of years. Recent studies have shown that resuspension of underlying sediments act as seed crystals for the slow precipitation of CaCO₃. Similar studies have not been conducted on the nearby more northerly Little Bahama Bank (LBB) where satellite photographs indicate the presence of extensive whitings north of the Grand Bahama Island. Research cruises were made to the LBB in July 2003 and May 2005. On board measurements were used to examine the distribution of the components of the carbonic acid system in LBB waters. The distribution of suspended calcium carbonate was also determined and samples of the suspended material and underlying sediments were collected for analyses, including ¹⁴C age determinations.The maximum residence time for waters on the LBB is a little over half (～ 144 d) that on the GBB. A ～ 20% decrease was found in the salinity-normalized titration alkalinity for inner bank waters. This yields a carbonate precipitation rate similar to waters on the GBB. However, unlike the waters in the whitings on the GBB where no changes in carbonic acid system parameters have been detected, large changes in the carbonic acid system parameters were found in the whiting waters of the LBB. These results clearly indicate that active precipitation was occurring in these waters. These observations are interpreted as indicating that there is a much slower mixing of waters associated with whitings and surrounding bank waters on the LBB than on the GBB. The more closed nature of whitings on the LBB makes them better targets for study of the processes taking place in whitings.Suspended carbonate mineralogy and geochemistry closely matched those of the fine-grained portion of the underlying sediments. These results support a close relationship between carbonates found in whitings and those in underlying sediments. This is further evidenced by a similar age for whiting carbonates (～ 700 y) and easily resuspended sediment carbonates (～ 1000 y) and is what was previously observed on the GBB.These observations indicate that the processes leading to whiting formation and associated calcium carbonate precipitation along with changes in water chemistry are quite similar in these two shallow carbonate bank environments in the Bahamas.     

Efficient execution of the WRF model and other HPC applications in the cloud

There are many scientific applications that have high performance computing (HPC) demands. Such demands are traditionally supported by cluster- or Grid-based systems. Cloud computing, which has experienced a tremendous growth, emerged as an approach to provide on-demand access to computing resources. The cloud computing paradigm offers a number of advantages over other distributed platforms. For example, the access to resources is flexible and cost-effective since it is not necessary to invest a large amount of money on a computing infrastructure nor pay salaries for maintenance functions. Therefore, the possibility of using cloud computing for running high performance computing applications is attractive. However, it has been shown elsewhere that current cloud computing platforms are not suitable for running some of these kinds of applications since the performance offered is very poor. The reason is mainly the overhead from virtualisation which is extensively used by most cloud computing platforms as a means to optimise resource usage. Furthermore, running HPC applications in current cloud platforms is a complex task that in many cases requires configuring a cluster of virtual machines (VMs). In this paper, we present a lightweight virtualisation approach for efficiently running the Weather Research and Forecasting (WRF) model (a computing- and communication-intensive application) in a cloud computing environment. Our approach also provides a higher-level programming model that automates the process of configuring a cluster of VMs. We assume such a cloud environment can be shared with other types of HPC applications such as mpiBLAST (an embarrassingly parallel application), and MiniFE (a memory-intensive application). Our experimental results show that lightweight virtualisation imposes about 5 % overhead and it substantially outperforms traditional heavyweight virtualisation such as KVM.