A Model for Solving the Optimal Water Allocation Problem in River Basins with Network Flow Programming When Introducing Non-Linearities

[EN] The allocation of water resources between different users is a traditional problem in many river basins. The objective is to obtain the optimal resource distribution and the associated circulating flows through the system. Network flow programming is a common technique for solving this problem. This optimisation procedure has been used many times for developing applications for concrete water systems, as well as for developing complete decision support systems. As long as many aspects of a river basin are not purely linear, the study of non-linearities will also be of great importance in water resources systems optimisation. This paper presents a generalised model for solving the optimal allocation of water resources in schemes where the objectives are minimising the demand deficits, complying with the required flows in the river and storing water in reservoirs. Evaporation from reservoirs and returns from demands are considered, and an iterative methodology is followed to solve these two non-network constraints. The model was applied to the Duero River basin (Spain). Three different network flow algorithms (Out-of-Kilter, RELAX-IVand NETFLO) were used to solve the allocation problem. Certain convergence issues were detected during the iterative process. There is a need to relate the data from the studied systems with the convergence criterion to be able to find the convergence criterion which yields the best results possible without requiring a long calculation time.We thank the Spanish Ministry of Economy and Competitivity (Comision Interministerial de Ciencia y Tecnologia, CICYT) for funding the projects INTEGRAME (contract CGL2009-11798) and SCARCE (program Consolider-Ingenio 2010, project CSD2009-00065). We also thank the European Commission (Directorate-General for Research & Innovation) for funding the project DROUGHT-R&SPI (program FP7-ENV-2011, project 282769). And last, but not least, to the Fundacion Instituto Euromediterraneo del Agua with the project "Estudio de Adaptaciones varias del modelo de optimizacion de gestiones de recursos hidricos Optiges".Haro Monteagudo, D.; Paredes Arquiola, J.; Solera Solera, A.; Andreu Álvarez, J. (2012). A Model for Solving the Optimal Water Allocation Problem in River Basins with Network Flow Programming When Introducing Non-Linearities. Water Resources Management. 26(14):4059-4071. https://doi.org/10.1007/s11269-012-0129-7S405940712614Ahuja R, Magnanti T, Orlin J (1993) Network flows: theory, algorithms and applications. Prentice Hall, New YorkAndreu J, Capilla J, Sanchís E (1996) AQUATOOL, a generalized decision-support system for water resources planning and operational management. J Hydrol 177:269–291Bersetkas D (1985) A unified framework for primal-dual methods in minimum cost network flows problems. Math Program 32:125–145Bersetkas D, Tseng P (1988) The relax codes for linear minimum cost network flow problems. Ann Oper Res 13:125–190Bersetkas D, Tseng P (1994) RELAX-IV: A faster version of the RELAX code for solving minimum cost flow problems. Completion Report under NSFGrant CCR-9103804. Dept. of Electrical Engineering and Computer Science, MIT, BostonChou F, Wu C, Lin C (2006) Simulating multi-reservoir operation rules by network flow model. ASCE Conf Proc 212:33Chung F, Archer M, DeVries J (1989) Network flow algorithm applied to California aqueduct simulation. J Water Resour Plan Manag 115:131–147Ford L, Fulkerson D (1962) Flows in networks. Princeton University Press, PrincetonFredericks J, Labadie J, Altenhofen J (1998) Decision support system for conjunctive stream-aquifer management. J Water Resour Plan Manag 124:69–78Harou JJ, Medellín-Azuara J, Zhu T et al (2010) Economic consequences of optimized water management for a prolonged, severe drought in California. Water Resour Res 46:W05522Hsu N, Cheng K (2002) Network Flow Optimization Model for Basin-Scale Water Supply Planning. J Water Resour Plan Manag 128:102–112Ilich N (1993) Improvement of the return flow allocation in the Water Resources Management Model of Alberta Environment. Can J Civ Eng 20:613–621Ilich N (2009) Limitations of network flow algorithms in river basin modeling. J Water Resour Plan Manag 135:48–55Kennington JL, Helgason RV (1980) Algorithms for network programming. John Wiley and Sons, New YorkKhaliquzzaman, Chander S (1997) Network flow programming model for multireservoir sizing. J Water Resour Plan Manag 123:15–21Kuczera G (1989) Fast Multireservoir Mulltiperiod Linear Programming Models. Water Resour Res 25:169–176Kuczera G (1993) Network linear programming codes for water-supply headworks modeling. J Water Resour Plan Manag 119:412–417Labadie J (2004) Optimal operation of multireservoir systems: state-of-the-art review. J Water Resour Plan Manag 130:93–111Labadie J (2006) MODSIM: river basin management decision support system. In: Singh W, Frevert D (eds) Watershed models. CRC, Boca Raton, pp 569–592Labadie J, Baldo M, Larson R (2000) MODSIM: decision support system for river basin management. Documentation and user manual. Dept. Of Civil Engineering, CSU, Fort CollinsManca A, Sechi G, Zuddas P (2010) Water supply network optimisation using equal flow algorithms. Water Resour Manag 24:3665–3678MMA (2000) Libro blanco del agua en España. Ministerio de Medio Ambiente, Secretaría general Técnica, Centro de PublicacionesMMA (2008) Confederación Hidrográfica del Duero. Memoria 2008. http://www.chduero.es/Inicio/Publicaciones/tabid/159/Default.aspx . Last accessed 25 June 2012Perera B, James B, Kularathna M (2005) computer software tool REALM for sustainable water allocation and management. J Environ Manag 77:291–300Rani D, Moreira M (2010) Simulation-optimization modeling: a survey and potential application in reservoir systems operation. Water Resour Manag 24:1107–1138Reca J, Roldán J, Alcaide M, López R, Camacho E (2001a) Optimisation model for water allocation in deficit irrigation systems I. Description of the model. Agric Water Manag 48:103–116Reca J, Roldán J, Alcaide M, López R, Camacho E (2001b) Optimisation model for water allocation in deficit irrigation systems II. Application to the Bembézar irrigation system. Agric Water Manag 48:117–132Sechi G, Zuddas P (2008) Multiperiod hypergraph models for water systems optimization. Water Resour Manag 22:307–320Sun H, Yeh W, Hsu N, Louie P (1995) Generalized network algorithm for water-supply-system optimization. J Water Resour Plan Manag 121:392–398Wurbs R (1993) Reservoir-system simulation and optimization models. J Water Resour Plan Manag 119:455–472Wurbs R (2005) Modeling river/reservoir system management, water allocation, and supply reliability. J Hydrol 300:100–113Yamout G, El-Fadel M (2005) An optimization approach for multi-sectoral water supply management in the greater Beirut area. Water Resour Manag 19:791–812Yates D, Sieber J, Purkey D, Hubert-Lee A (2005) WEAP21 – a demand-, priority-, and preference-driven water planning model. Part 1: model characteristics. Water Int 30:487–500Zoltay V, Vogel R, Kirshen P, Westphal K (2010) Integrated watershed management modeling: generic optimization model applied to the Ipswich river basin. J Water Resour Plan Manag 136:566–57

KELT-20b: A Giant Planet With A Period Of P ~ 3.5 Days Transiting The V ~ 7.6 Early A Star HD 185603

We report the discovery of KELT-20b, a hot Jupiter transiting a early A star, HD 185603, with an orbital period of days. Archival and follow-up photometry, Gaia parallax, radial velocities, Doppler tomography, and AO imaging were used to confirm the planetary nature of KELT-20b and characterize the system. From global modeling we infer that KELT-20 is a rapidly rotating ( ) A2V star with an effective temperature of K, mass of , radius of , surface gravity of , and age of . The planetary companion has a radius of , a semimajor axis of au, and a linear ephemeris of . We place a upper limit of on the mass of the planet. Doppler tomographic measurements indicate that the planetary orbit normal is well aligned with the projected spin axis of the star ( ). The inclination of the star is constrained to , implying a three-dimensional spin–orbit alignment of . KELT-20b receives an insolation flux of , implying an equilibrium temperature of of ∼2250 K, assuming zero albedo and complete heat redistribution. Due to the high stellar , KELT-20b also receives an ultraviolet (wavelength nm) insolation flux of , possibly indicating significant atmospheric ablation. Together with WASP-33, Kepler-13 A, HAT-P-57, KELT-17, and KELT-9, KELT-20 is the sixth A star host of a transiting giant planet, and the third-brightest host (in V ) of a transiting planet

Ultraviolet Spectropolarimetry with Polstar: on the origin of rapidly rotating B stars

The proposed MIDEX mission, Polstar, will provide high resolution UV spectroscopy and spectropolarimetry and offers a unique opportunity to study massive stars in this wavelength range with unprecedented detail. We demonstrate that these observations will provide critical new knowledge of several types of massive stars (specifically the B-emission stars and the Bn stars). We will determine accurate stellar parameters including their rotation rates and variation of surface temperature associated with stellar oblateness. Our work will allow us to detect binary companions and determine binary orbital properties. Binary population synthesis predictions will allow us to determine the fraction of these stars that are spun up due to binary interaction compared to single star evolution. These rapidly rotating stars have the potential to probe mass loss and the mixing of chemical elements which affects their evolution and ultimately the evolution of their surroundings.Comment: white paper, 24 page

A Giant Planet Undergoing Extreme-Ultraviolet Irradiation By Its Hot Massive-Star Host

The amount of ultraviolet irradiation and ablation experienced by a planet depends strongly on the temperature of its host star. Of the thousands of extrasolar planets now known, only six have been found that transit hot, A-type stars (with temperatures of 7,300–10,000 kelvin), and no planets are known to transit the even hotter B-type stars. For example, WASP-33 is an A-type star with a temperature of about 7,430 kelvin, which hosts the hottest known transiting planet, WASP-33b (ref. 1); the planet is itself as hot as a red dwarf star of type M (ref. 2). WASP-33b displays a large heat differential between its dayside and nightside2, and is highly inflated–traits that have been linked to high insolation3,4. However, even at the temperature of its dayside, its atmosphere probably resembles the molecule-dominated atmospheres of other planets and, given the level of ultraviolet irradiation it experiences, its atmosphere is unlikely to be substantially ablated over the lifetime of its star. Here we report observations of the bright star HD 195689 (also known as KELT-9), which reveal a close-in (orbital period of about 1.48 days) transiting giant planet, KELT-9b. At approximately 10,170 kelvin, the host star is at the dividing line between stars of type A and B, and we measure the dayside temperature of KELT-9b to be about 4,600 kelvin. This is as hot as stars of stellar type K4 (ref. 5). The molecules in K stars are entirely dissociated, and so the primary sources of opacity in the dayside atmosphere of KELT-9b are probably atomic metals. Furthermore, KELT-9b receives 700 times more extreme-ultraviolet radiation (that is, with wavelengths shorter than 91.2 nanometres) than WASP-33b, leading to a predicted range of mass-loss rates that could leave the planet largely stripped of its envelope during the main-sequence lifetime of the host star (ref. 6)

KELT-11b: A Highly Inflated Sub-Saturn Exoplanet Transiting the V=8 Subgiant HD 93396

We report the discovery of a transiting exoplanet, KELT-11b, orbiting the bright ($V=8.0$) subgiant HD 93396. A global analysis of the system shows that the host star is an evolved subgiant star with $T_{\rm eff} = 5370\pm51$ K, $M_{*} = 1.438_{-0.052}^{+0.061} M_{\odot}$, $R_{*} = 2.72_{-0.17}^{+0.21} R_{\odot}$, log $g_*= 3.727_{-0.046}^{+0.040}$, and [Fe/H]$= 0.180\pm0.075$. The planet is a low-mass gas giant in a $P = 4.736529\pm0.00006$ day orbit, with $M_{P} = 0.195\pm0.018 M_J$, $R_{P}= 1.37_{-0.12}^{+0.15} R_J$, $\rho_{P} = 0.093_{-0.024}^{+0.028}$ g cm$^{-3}$, surface gravity log ${g_{P}} = 2.407_{-0.086}^{+0.080}$, and equilibrium temperature $T_{eq} = 1712_{-46}^{+51}$ K. KELT-11 is the brightest known transiting exoplanet host in the southern hemisphere by more than a magnitude, and is the 6th brightest transit host to date. The planet is one of the most inflated planets known, with an exceptionally large atmospheric scale height (2763 km), and an associated size of the expected atmospheric transmission signal of 5.6%. These attributes make the KELT-11 system a valuable target for follow-up and atmospheric characterization, and it promises to become one of the benchmark systems for the study of inflated exoplanets.Comment: 15 pages, Submitted to AAS Journal

KELT-18b: Puffy Planet, Hot Host, Probably Perturbed

We report the discovery of KELT-18b, a transiting hot Jupiter in a 2.87-day orbit around the bright ( V = 10.1), hot, F4V star BD+60 1538 (TYC 3865-1173-1). We present follow-up photometry, spectroscopy, and adaptive optics imaging that allow a detailed characterization of the system. Our preferred model fits yield a host stellar temperature of K and a mass of , situating it as one of only a handful of known transiting planets with hosts that are as hot, massive, and bright. The planet has a mass of , a radius of , and a density of , making it one of the most inflated planets known around a hot star. We argue that KELT-18b’s high temperature and low surface gravity, which yield an estimated ∼600 km atmospheric scale height, combined with its hot, bright host, make it an excellent candidate for observations aimed at atmospheric characterization. We also present evidence for a bound stellar companion at a projected separation of ∼1100 au, and speculate that it may have contributed to the strong misalignment we suspect between KELT-18\u27s spin axis and its planet’s orbital axis. The inferior conjunction time is 2457542.524998 ± 0.000416 (BJD TDB ) and the orbital period is 2.8717510 ± 0.0000029 days. We encourage Rossiter–McLaughlin measurements in the near future to confirm the suspected spin–orbit misalignment of this system

Dynamical masses of the primary Be star and the secondary sdB star in the single-lined binary kappa Dra (B6 IIIe)

Because many classical Be stars may owe their nature to mass and angular-momentum transfer in a close binary, the present masses, temperatures, and radii of their components are of high interest for comparison to stellar evolution models. Kappa Dra is a 61.5-day single-lined binary with a B6 IIIe primary. With the CHARA Array instruments MIRC/MIRC-X and MYSTIC, we detected the secondary at (approximately photospheric) flux ratios of 1.49 +- 0.10% and 1.63 +- 0.09% in the H and K band, respectively. From a large and diverse optical spectroscopic database only the radial velocity curve of the Be star could be extracted. However, employing the parallaxes from Hipparcos and Gaia, which agree within their nominal 1-sigma errors, we could derive the total mass and found component masses of 3.65 +- 0.48 Msun and 0.426 +- 0.043 Msun for the Be star and the companion, respectively. Previous cross-correlation of the observed far-UV spectrum with sdO spectral model templates had not detected a companion belonging to the hot O-type subdwarf (sdO) population known from ~20 earlier-type Be stars. Guided by our full 3D orbital solution, we found a strong cross-correlation signal for a stripped subdwarf B-type companion (far-UV flux ratio of 2.3 +- 0.5%), enabling the first firm characterization of such a star, and making kappa Dra the first mid- to late-type Be star with a directly-observed subdwarf companion.Comment: Accepted to ApJ. Figure sets available on reques

The KELT Follow-Up Network And Transit False-Positive Catalog: Pre-Vetted False Positives For TESS

The Kilodegree Extremely Little Telescope (KELT) project has been conducting a photometric survey of transiting planets orbiting bright stars for over 10 years. The KELT images have a pixel scale of ~23\u27\u27 pixel⁻¹—very similar to that of NASA\u27s Transiting Exoplanet Survey Satellite (TESS)—as well as a large point-spread function, and the KELT reduction pipeline uses a weighted photometric aperture with radius 3\u27. At this angular scale, multiple stars are typically blended in the photometric apertures. In order to identify false positives and confirm transiting exoplanets, we have assembled a follow-up network (KELT-FUN) to conduct imaging with spatial resolution, cadence, and photometric precision higher than the KELT telescopes, as well as spectroscopic observations of the candidate host stars. The KELT-FUN team has followed-up over 1600 planet candidates since 2011, resulting in more than 20 planet discoveries. Excluding ~450 false alarms of non-astrophysical origin (i.e., instrumental noise or systematics), we present an all-sky catalog of the 1128 bright stars (6 \u3c V \u3c 13) that show transit-like features in the KELT light curves, but which were subsequently determined to be astrophysical false positives (FPs) after photometric and/or spectroscopic follow-up observations. The KELT-FUN team continues to pursue KELT and other planet candidates and will eventually follow up certain classes of TESS candidates. The KELT FP catalog will help minimize the duplication of follow-up observations by current and future transit surveys such as TESS