Estimation of water resources on continental surfaces by multi-sensor microwave remote sensing

L'estimació dels recursos hídrics de les superfícies continentals a escala regional i global és fonamental per a una bona gestió dels recursos hídrics. Aquesta estimació cobreix una àmplia gamma de temes i camps, incloent-hi la caracterització dels sòls i dels recursos hídrics a l’escala de la conca, la modelització hidrològica i la predicció i la cartografia d'inundacions. En aquest context, la caracterització dels estats de la superfície continental, per a obtenir millors paràmetres d’entrada als models hidrològics, és essencial per millorar la precisió en la simulació de cabals, sequeres i inundacions. L’estimació del contingut d’aigua en el sistema, incloses les diferents masses d’aigua i l’aigua lliure en el sòl, és especialment necessària per a una descripció precisa dels processos hidrològics i, en general, del cicle de l’aigua a les superfícies continentals. Per caracteritzar millor els processos hidrològics, les intervencions antropogèniques no es poden negligir. L'home influeix en el cicle de l'aigua, principalment mitjançant el reg i la construcció de preses, fet que s’ha de quantificar correctament. L’objectiu de la tesi és la millora de l’estimació remota dels recursos hídrics, incloent-hi la quantificació dels factors antròpics, mitjançant l’ús de diversos sensors llançats recentment, aprofitant recents desenvolupaments en la tecnologia de teledetecció. Amb l'arribada de les constel·lacions Sentinel (Sentinel-1, 2, 3), disposem de millors eines per estimar els recursos hídrics, incloent-hi els impactes humans, amb una major precisió i cobertura. Aquest treball de tesi consta principalment de dues línies de recerca on s’estimen les intervencions humanes en el cicle hidrològic: la cartografia del reg (com a aplicació en humitat del sòl), i el forçament d’embassaments en simulacions hidrològiques (com a aplicació de l’altimetria). En la primera linia s’estima la humitat del sòl a partir de l’anàlisi estadística de les dades SAR de Sentinel-1. Es desenvolupen dues metodologies per obtenir la humitat del sòl amb una resolució espacial de 100 m basant-se en la interpretació de les dades de Sentinel-1 obtingudes amb la polarització VV (vertical-vertical), que es combina amb dades òptiques Sentinel-2 per a l'anàlisi dels efectes de la vegetació. Com aplicació de la humitat del sòl, es cartografia el reg en diverses condicions meteorològiques, i amb una alta resolució espacial i temporal. Es proposa una metodologia per a la cartografia del reg mitjançant dades SAR obtingudes en polaritzacions VV (vertical-vertical) i VH (vertical-horitzontal). A partir de la sèrie temporal Sentinel-1, s’analitzen diferents estadístiques i mètriques, incloent-hi el valor mitjà, la variància del senyal, la longitud de la correlació i la dimensió fractal, a partir dels quals es classifiquen els arbres irrigats, els cultius irrigats i els cultius no irrigats. En la segona línia, s’estima el nivell dels embassaments a partir de les dades d’altimetria de Sentinel-3, amb l’altímetre SAR (SRAL), basant-se en diferents algorismes per millorar la precisió. Aquest estudi presenta tres algorismes especialitzats o retrackers destinats a obtenir el nivell de la superfície dels cossos d’aigua estudiats, minimitzant la contaminació de les formes d’ona degut al sòl que els envolta. Es compara el rendiment del mètode proposat de selecció de la porció d’ona amb tres retrackers, és a dir, un retracker de llindar, el retracker del centre de gravetat (OCOG) i un retracker de base física de dos passos. S’obtenen sèries temporals del nivell de la làmina d’aigua d’embassaments situats a la conca del riu Ebre (Espanya). Com aplicació, les sèries de nivell dels embassaments obtingudes s’utilitzen per a forçar els embassaments en simulacions hidrològiques.La estimación de los recursos hídricos de las superficies continentales a escala regional y global es fundamental para una buena gestión de los recursos hídricos. Esta estimación cubre una amplia gama de temas y campos, incluyendo la caracterización de los suelos y de los recursos hídricos a escala de cuenca, la modelización hidrológica y la predicción y la cartografía de inundaciones. En este contexto, la caracterización de los estados de la superficie continental, para obtener mejores parámetros de entrada para los modelos hidrológicos, es esencial para mejorar la precisión en la simulación de caudales, sequías e inundaciones. La estimación del contenido de agua en el sistema, incluidas las diferentes masas de agua y el agua libre en el suelo, es especialmente necesaria para una descripción precisa de los procesos hidrológicos y, en general, del ciclo del agua en las superficies continentales. Una caracterización precisa de los procesos hidrológicos requiere no descuidar las intervenciones humanas. El hombre influye en el ciclo del agua, principalmente mediante el riego y la construcción de embalses, lo que se debe cuantificar correctamente. El objetivo de la tesis es la mejora de la estimación remota de los recursos hídricos, incluyendo la cuantificación de los factores humanos, mediante el uso de varios sensores lanzados recientemente, aprovechando recientes desarrollos en la tecnología de teledetección. Con la llegada de las constelaciones Sentinel (Sentinel-1, 2, 3), disponemos de mejores herramientas para estimar los recursos hídricos, incluyendo los impactos humanos, con una mayor precisión y cobertura. Este trabajo de tesis consta principalmente en dos ejes de investigación donde se estiman las intervenciones humanas en el ciclo hidrológico: la cartografía del riego (como aplicación en humedad del suelo), y el forzamiento de embalses en simulaciones hidrológicas (como aplicación de la altimetría). En relación al primer eje, se estima la humedad del suelo a partir del análisis estadístico de los datos SAR de Sentinel-1. Se desarrollan dos metodologías para obtener la humedad del suelo con una resolución espacial de 100 m basándose en la interpretación de los datos de Sentinel-1 obtenidas con la polarización VV (vertical-vertical), que se combina con datos ópticas Sentinel-2 para el análisis de los efectos de la vegetación. Como aplicación de la humedad del suelo, se cartografía el riego en diversas condiciones meteorológicas, y con una alta resolución espacial y temporal. Se propone una metodología para la cartografía del riego mediante datos SAR obtenidos en polarizaciones VV (vertical-vertical) y VH (vertical-horizontal). A partir de la serie temporal Sentinel-1, se analizan diferentes estadísticas y métricas, incluyendo el valor medio, la varianza de la señal, la longitud de la correlación y la dimensión fractal, a partir de los cuales se clasifican los árboles irrigados, los cultivos irrigados y los cultivos no irrigados. En el segundo eje, se estima el nivel de los embalses a partir de los datos de altimetría de Sentinel-3, con el altímetro SAR (SRAL), basándose en diferentes algoritmos para mejorar la precisión. Este estudio presenta tres algoritmos especializados o retrackers destinados a obtener el nivel de la superficie de los cuerpos de agua estudiados, minimizando la contaminación de las formas de onda debido al suelo que los rodea. Se compara el rendimiento del método propuesto de selección de la porción de onda con tres retrackers, es decir, un retracker de umbral, el retracker del centro de gravedad (OCOG) y un retracker de base física de dos pasos. Se obtienen series temporales del nivel de la lámina de agua de embalses situados en la cuenca del río Ebro (España). Como aplicación, las series de nivel de los embalses obtenidas se utilizan para forzar los embalses en simulaciones hidrológicas.The estimation of the water resources of the continental surfaces at a regional and global scale is fundamental for good water resources management. This estimation covers a wide range of topics and fields, including the characterisation of soils and water resources at the basin scale, hydrological modelling and flood prediction and mapping. In this context, the characterisation of the states of the continental surface, to obtain better input parameters for hydrological models, is essential to improve the precision in the simulation of flows, droughts, and floods. The estimation of the water content in the system, including the different water bodies and the free water in the soil, is especially necessary for a precise description of the hydrological processes and, in general, of the water cycle on the continental surfaces. To better characterise hydrological processes, human interventions cannot be neglected. Humans influence the water cycle, mainly through irrigation and the construction of reservoirs, which must be correctly quantified. The objective of the thesis is the improvement of the remote estimation of water resources, including the quantification of human factors, using several sensors recently launched, taking advantage of recent developments in remote sensing technology. With the arrival of the Sentinel constellations (Sentinel-1, 2, 3), we have better tools to estimate water resources, including human impacts, with greater precision and coverage. This thesis consists mainly of two parts where human interventions in the water cycle are considered: irrigation cartography (as an application of soil moisture), and the forcing of reservoirs in hydrological simulations (as an application of altimetry). Firstly, soil moisture is estimated from the statistical analysis of Sentinel-1 SAR data. Two methodologies are developed to obtain soil moisture with a spatial resolution of 100 m based on the interpretation of Sentinel-1 data collected with the VV polarization (vertical-vertical), which is combined with optical data of Sentinel-2 for the analysis of the effects of vegetation. Secondly, irrigation is mapped under various meteorological conditions, including high spatial and temporal resolution. A methodology for irrigation mapping is proposed using SAR data obtained in VV (vertical-vertical) and VH (vertical-horizontal) polarizations. With Sentinel-1 time series, different statistics and metrics are analysed, including the mean value, the variance of the signal, the correlation length and the fractal dimension, based on which the classification of irrigated trees, irrigated crops, and non-irrigated crops are derived. Finally, the level of the reservoirs is estimated from the Sentinel-3 altimetry data, with the SAR altimeter (SRAL), based on different algorithms to improve the accuracy. This study presents three specialised algorithms or retrackers designed to obtain the level of the surface of the studied inland bodies of water, minimising the contamination of the waveforms due to the surrounding soil. The performance of the selection method of the proposed wave portion is compared with three retrackers, that is, the centre of gravity retracker (OCOG) and the two-step physical-based retracker. Temporal series of the water level of reservoirs located in the basin of the Ebro River (Spain) are obtained. As an application, the level series of the reservoirs obtained are used to force the reservoirs in hydrological simulations.L'estimation et le suivi des ressources en eau des surfaces continentales aux niveaux régional et global est essentielle pour la gestion du bilan hydrique, particulièrement dans le contexte des changements climatiques et anthropiques. Ils couvrent un large éventail de thèmes et de domaines, incluant la caractérisation des ressources en eau à l'échelle du bassin, la modélisation hydrologique ainsi que la prévision et la cartographie des inondations. Dans ce contexte, la caractérisation des états de surface, en tant que paramètres d’entrée dans les modèles hydrologiques, est essentielle pour obtenir une meilleure précision de la simulation, qui est liée à la précision prévisionnelle des débits des cours d’eau et le suivi des sécheresses et des inondations. L'estimation de la teneur en eau des surfaces continentales, incluant l’état hydrique du sol et les niveaux des surfaces couvertes d’eau, est particulièrement nécessaire pour une description précise des processus hydrologiques et plus généralement du cycle de l'eau sur les surfaces continentales. Afin de mieux comprendre les processus hydrologiques, l'influence humaine (l’effet anthropique) sur le cycle de l'eau nécessite une évaluation fine. Elle est particulièrement liée à la gestion de l’irrigation et la construction de barrages. L'objectif de la thèse était d'améliorer l'estimation des ressources en eau et une meilleure caractérisation des interventions anthropiques à travers l'utilisation de nouveaux capteurs satellitaires multi-configurations du programme européen Copernicus. Avec le développement de la technologie de télédétection spatiale, et plus particulièrement avec l’arrivée des constellations Sentinel (Sentinel-1, 2, 3) à haute résolution spatiale et temporelle, il existe un meilleur outil pour estimer les états des surfaces continentales. Ce travail de thèse comprend principalement deux priorités liées à des interventions humaines dans le cycle hydrologique:la cartographie de l'irrigation en tant que action humaine liée directement à l'humidité du sol et le forçage des barrages dans un modèle de simulation de rivière en tant qu'application liée à l’estimation du niveau de l'eau libre. Un premier axe de recherche a été basé sur une analyse statistique des données SAR Sentinel-1 pour caractériser l’état hydrique du sol. Deux méthodes ont été développées pour estimer ce paramètre avec une résolution spatiale de 100 m. Elles sont basées sur des approches de détection de changement à partir des données Sentinel-1 acquises en polarisation VV (verticale-verticale), combinées aux données optiques Sentinel-2 pour corriger les effets de la végétation. L’application consistait à cartographier l'irrigation, avec des résolutions spatiale et temporelle élevées. Une méthodologie de cartographie de l'irrigation utilisant des données SAR Sentinel-1 a été proposée. Elle estbasée sur les acquisitions en polarisations VV (vertical-vertical) et VH (vertical-horizontal). A partir de la série temporelle des mesures Sentinel-1, des paramètres statistiques tel que la valeur moyenne, la variance du signal, la longueur de corrélation temporelle et la dimension fractale, sont analysées, en fonction du type de culture; cultures annuelles irriguées, arbres irrigués et cultures pluviales. Des classifications supervisées utilisant les approches Random Forest et SVM sont testées. En deuxième axe, l'estimation de la hauteur de la surface de l'eau à partir des données altimétriques de Sentinel-3 avec l’altimètre SAR (SRAL) a été réalisée à l'aide de différents algorithmes afin d'améliorer la précision sur des petites surfaces. Cette étude présente trois algorithmes spécialisés (ou retrackers) dédiées à la minimisation de la contamination des sols par les formes d’ondes permettant de récupérer les niveaux d’eau à partir de données altimétriques SAR sur des masses d’eaux intérieures. Les performances de la méthode de sélection de portion de forme d'onde proposée avec trois retrackers, à savoir, le retracker à seuil, le retracker à centre de gravité décalé (OCOG) et le retracker à base physique à 2 étapes, sont comparées. Des séries chronologiques de niveaux d'eau sont extraites pour les masses d'eau du bassin de l'Èbre (Espagne). Une application des produits altimétriques est proposée. Le produit de niveau d’eau a été utilisé comme paramètre d’entrée pour analyser l’effet tampon des barrages dans les simulations de débits fluviaux

Desert roughness retrieval using CYGNSS GNSS-R data

The aim of this paper is to assess the potential use of data recorded by the Global Navigation Satellite System Reflectometry (GNSS-R) Cyclone Global Navigation Satellite System (CYGNSS) constellation to characterize desert surface roughness. The study is applied over the Sahara, the largest non-polar desert in the world. This is based on a spatio-temporal analysis of variations in Cyclone Global Navigation Satellite System (CYGNSS) data, expressed as changes in reflectivity (G). In general, the reflectivity of each type of land surface (reliefs, dunes, etc.) encountered at the studied site is found to have a high temporal stability. A grid of CYGNSS G measurements has been developed, at the relatively fine resolution of 0.03° x 0.03°, and the resulting map of average reflectivity, computed over a 2.5-year period, illustrates the potential of CYGNSS data for the characterization of the main types of desert land surface (dunes, reliefs, etc.). A discussion of the relationship between aerodynamic or geometric roughness and CYGNSS reflectivity is proposed. A high correlation is observed between these roughness parameters and reflectivity. The behaviors of the GNSS-R reflectivity and the Advanced Land Observing Satellite-2 (ALOS-2) Synthetic Aperture Radar (SAR) backscattering coeffcient are compared and found to be strongly correlated. An aerodynamic roughness (Z0) map of the Sahara is proposed, using four distinct classes of terrain roughness

Earth resources: A continuing bibliography with indexes (issue 61)

This bibliography lists 606 reports, articles, and other documents introduced into the NASA scientific and technical information system between January 1 and March 31, 1989. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, oceanography and marine resources, hydrology and water management, data processing and distribution systems, and instrumentation and sensors, and economic analysis

Summaries of the Fifth Annual JPL Airborne Earth Science Workshop. Volume 3: AIRSAR Workshop

This publication is the third containing summaries for the Fifth Annual JPL Airborne Earth Science Workshop, held in Pasadena, California, on January 23-26, 1995. The main workshop is divided into three smaller workshops as follows: (1) The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) workshop, on January 23-24. The summaries for this workshop appear in Volume 1; (2) The Airborne synthetic Aperture Radar (AIRSAR) workshop, on January 25-26. The summaries for this workshop appear in this volume; and (3) The Thermal Infrared Multispectral Scanner (TIMS) workshop, on January 26. The summaries for this workshop appear in Volume 2

Earth resources: A continuing bibliography with indexes (issue 62)

This bibliography lists 544 reports, articles, and other documents introduced into the NASA scientific and technical information system between April 1 and June 30, 1989. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

Earth resources: A continuing bibliography with indexes (issue 58)

This bibliography lists 500 reports, articles, and other documents introduced into the NASA scientific and technical information system between April 1 and June 30, 1988. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

Earth resources: A continuing bibliography with indexes (issue 51)

This bibliography lists 382 reports, articles and other documents introduced into the NASA scientific and technical information system between July 1 and September 30, 1986. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

GNSS reflectometry for land remote sensing applications

Soil moisture and vegetation biomass are two essential parameters from a scienti c and economical point of view. On one hand, they are key for the understanding of the hydrological and carbon cycle. On the other hand, soil moisture is essential for agricultural applications and water management, and vegetation biomass is crucial for regional development programs. Several remote sensing techniques have been used to measure these two parameters. However, retrieving soil moisture and vegetation biomass with the required accuracy, and the appropriate spatial and temporal resolutions still remains a major challenge. The use of Global Navigation Satellite Systems (GNSS) reflected signals as sources of opportunity for measuring soil moisture and vegetation biomass is assessed in this PhD Thesis. This technique, commonly known as GNSS-Reflectometry (GNSS-R), has gained increasing interest among the scienti c community during the last two decades due to its unique characteristics. Previous experimental works have already shown the capabilities of GNSS-R to sense small reflectivity changes on the surface. The use of the co- and cross-polarized reflected signals was also proposed to mitigate nuisance parameters, such as soil surface roughness, in the determination of soil moisture. However, experimental evidence of the suitability of that technique could not be demonstrated. This work analyses from a theoretical and an experimental point of view the capabilities of polarimetric observations of GNSS reflected signals for monitoring soil moisture and vegetation biomass. The Thesis is structured in four main parts. The fi rst part examines the fundamental aspects of the technique and provides a detailed review of the GNSS-R state of the art for soil moisture and vegetation monitoring. The second part deals with the scattering models from land surfaces. A comprehensive description of the formation of scattered signals from rough surfaces is provided. Simulations with current state of the art models for bare and vegetated soils were performed in order to analyze the scattering components of GNSS reflected signals. A simpli ed scattering model was also developed in order to relate in a straightforward way experimental measurements to soil bio-geophysical parameters. The third part reviews the experimental work performed within this research. The development of a GNSS-R instrument for land applications is described, together with the three experimental campaigns carried out in the frame of this PhD Thesis. The analysis of the GNSS-R and ground truth data is also discussed within this part. As predicted by models, it was observed that GNSS scattered signals from natural surfaces are a combination of a coherent and an incoherent scattering components. A data analysis technique was proposed to separate both scattering contributions. The use of polarimetric observations for the determination of soil moisture was demonstrated to be useful under most soil conditions. It was also observed that forests with high levels of biomass could be observed with GNSS reflected signals. The fourth and last part of the Thesis provides an analysis of the technology perspectives. A GNSS-R End-to-End simulator was used to determine the capabilities of the technique to observe di erent soil reflectivity conditions from a low Earth orbiting satellite. It was determined that high accuracy in the estimation of reflectivity could be achieved within reasonable on-ground resolution, as the coherent scattering component is expected to be the predominant one in a spaceborne scenario. The results obtained in this PhD Thesis show the promising potential of GNSS-R measurements for land remote sensing applications, which could represent an excellent complementary observation for a wide range of Earth Observation missions such as SMOS, SMAP, and the recently approved ESA Earth Explorer Mission Biomass.La humedad del suelo y la biomasa de la vegetaci on son dos parametros clave desde un punto de vista tanto cient co como econ omico. Por una parte son esenciales para el estudio del ciclo del agua y del carbono. Por otra parte, la humedad del suelo es esencial para la gesti on de las cosechas y los recursos h dricos, mientras que la biomasa es un par ametro fundamental para ciertos programas de desarrollo. Varias formas de teledetección se han utilizado para la observaci on remota de estos par ametros, sin embargo, su monitorizaci on con la precisi on y resoluci on necesarias es todav a un importante reto tecnol ogico. Esta Tesis evalua la capacidad de medir humedad del suelo y biomasa de la vegetaci on con señales de Sistemas Satelitales de Posicionamiento Global (GNSS, en sus siglas en ingl es) reflejadas sobre la Tierra. La t ecnica se conoce como Reflectometr í a GNSS (GNSS-R), la cual ha ganado un creciente inter es dentro de la comunidad científ ca durante las dos ultimas d ecadas. Experimentos previos a este trabajo ya demostraron la capacidad de observar cambios en la reflectividad del terreno con GNSS-R. El uso de la componente copolar y contrapolar de la señal reflejada fue propuesto para independizar la medida de humedad del suelo de otros par ametros como la rugosidad del terreno. Sin embargo, no se pudo demostrar una evidencia experimental de la viabilidad de la t ecnica. En este trabajo se analiza desde un punto de vista te orico y experimental el uso de la informaci on polarim etrica de la señales GNSS reflejadas sobre el suelo para la determinaci on de humedad y biomasa de la vegetaci on. La Tesis se estructura en cuatro partes principales. En la primera parte se eval uan los aspectos fundamentales de la t ecnica y se da una revisi on detallada del estado del arte para la observaci on de humedad y vegetaci on. En la segunda parte se discuten los modelos de dispersi on electromagn etica sobre el suelo. Simulaciones con estos modelos fueron realizadas para analizar las componentes coherente e incoherente de la dispersi on de la señal reflejada sobre distintos tipos de terreno. Durante este trabajo se desarroll o un modelo de reflexi on simpli cado para poder relacionar de forma directa las observaciones con los par ametros geof sicos del suelo. La tercera parte describe las campañas experimentales realizadas durante este trabajo y discute el an alisis y la comparaci on de los datos GNSS-R con las mediciones in-situ. Como se predice por los modelos, se comprob o experimentalmente que la señal reflejada est a formada por una componente coherente y otra incoherente. Una t ecnica de an alisis de datos se propuso para la separacióon de estas dos contribuciones. Con los datos de las campañas experimentales se demonstr o el bene cio del uso de la informaci on polarim etrica en las señales GNSS reflejadas para la medici on de humedad del suelo, para la mayor a de las condiciones de rugosidad observadas. Tambi en se demostr o la capacidad de este tipo de observaciones para medir zonas boscosas densamente pobladas. La cuarta parte de la tesis analiza la capacidad de la t ecnica para observar cambios en la reflectividad del suelo desde un sat elite en orbita baja. Los resultados obtenidos muestran que la reflectividad del terreno podr a medirse con gran precisi on ya que la componente coherente del scattering ser a la predominante en ese tipo de escenarios. En este trabajo de doctorado se muestran la potencialidades de la t ecnica GNSS-R para observar remotamente par ametros del suelo tan importantes como la humedad del suelo y la biomasa de la vegetaci on. Este tipo de medidas pueden complementar un amplio rango de misiones de observaci on de la Tierra como SMOS, SMAP, y Biomass, esta ultima recientemente aprobada para la siguiente misi on Earth Explorer de la ESA

Earth Resources: A continuing bibliography with indexes, issue 1

This bibliography lists 616 reports, articles, and other documents introduced into the NASA scientific and technical information system between January 1974 and March 1974. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory, natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, oceanography and marine resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis