Research
Dr. Kavvas’ areas of specialization include assessment of climate change impact by mathematical modeling of the coupled hydrologic-atmospheric processes at watershed, regional and continental scales. The in-house-developed regional and watershed-scale coupled hydrologic-atmospheric numerical models are also used for the simulation and forecasting of hydrologic water balances and hydrologic extremes toward quantifying and estimating phenomena such as floods and droughts, and forecasting water supply at watershed and regional scales. The mathematical modeling of hydrologic processes at regional, watershed and hillslope scales are also used toward planning and management of water resources at these scales. Physics-based stochastic models are developed for investigating river flow, snow accumulation/melting, unsaturated soil water flow, unconfined and confined groundwater flow, and contaminant transport. Physical modeling of fish ecological hydraulics are performed at J. Amorocho Hydraulics Laboratory.
Major research areas of Dr. Kavvas are as follows:
- Climate change impact assessment at regional and watershed scales by means of coupled numerical hydrologic-atmospheric process modeling
Since mid-1990s Dr. Kavvas and his team of graduate students and postdocs have been working on the development of numerical computer models for the assessment of the impacts of climate change with respect to water balances and hydrologic extremes (floods and droughts) in California and over various regions of the world in collaboration with various government agencies.
- Within this framework, one of the first coupled numerical regional hydrologic-atmospheric models was developed by our group in collaboration with Japanese government scientists in mid-1990s in order to assess the impact of climate change on the climate and water resources of Japan’s main island, the Honshu Island Presentation of IRSHAM (Integrated Regional Scale Hydrologic-Atmospheric Model)-HRL1996; Application of a Coupled Regional-Scale Hydrological-Atmospheric Model to Japan for Climate Change Study. Journal of Hydrologic Engineering, 16(12), 2011. )
- In mid-1990s our team at UCD Hydrologic Research Laboratory (HRL) has developed the first numerical regional coupled hydrologic-atmospheric model of California to study the impact of climate change on the hydrologic and climatic conditions of California (CaliforniaClimateChangePresentation-2001.ppt ; Z.Q. Chen, M.L. Kavvas, L.Tan and S-T.Soong, “Development of a Regional Atmospheric-Hydrologic Model for the Study of Climate Change in California”, Proceedings of ASCE North American Water and Environment Congress, June, 1996)
- In 2000s our team collaborated with the Malaysian government scientists to assess the impact of climate change on the Peninsular Malaysia, and Sabah-Sarawak (Northern Borneo) water resources by means of coupled numerical hydrologic-atmospheric models we jointly developed for these regions (ClimChngProjection-Dynamical Approach-Malaysia -Presentation; Climate Change Projections for Sabah and Sarawak -presentation)
- In the last decade in California several climate change impact assessment studies were performed by our HRL team: i) in one study the water balances and drought conditions were estimated for Shasta Dam under the 21st century hydro-climate conditions (presentation: Impact of climate change on Shasta Dam Water Balances during 21st Century ) that were projected by our in-house-developed coupled hydrologic-atmospheric watershed hydro-climate model WEHY-HCM ( Presentation on WEHY-HCM Watershed Hydroclimate Model ); ii) in another study the 21st century hydro-climate conditions that were projected over the Upper Middle Fork subbasin and Sierra Valley aquifer of Feather river basin by means of our in-house WEHY-HCM model that was coupled to DWR’s groundwater model IWFM, were then used to quantify the water balances at the Sierra Valley Aquifer under the projected 21st century hydro-climate conditions and water use scenarios. This study is presented here: Climate Change Impact on the Sierra Valley Aquifer, California-2017 iii) In various other studies the estimation of future precipitation conditions over Northern California ( Trend analysis of watershed precipitation in Northern California under CMIP5 projections ), the estimation of future maximum precipitation at American River basin (Estimation of maximum precipitation at American River under climate change- Hydrol Processes- 2018.pdf ), the estimation of future snow conditions over Northern California (Impact of climate change on snow conditions at N California – STOTEN 2019.pdf ), and the estimation of future flood frequencies at Cache Creek basin under changing climate conditions, modeled by our in-house WEHY-HCM coupled atmospheric-hydrologic model (Estimation of Future Flood Frequencies under changing climate- JHE2016.pdf; presentation: Estimation of flood frequencies under changing climate), were performed. In these studies new methodologies for the assessment of climate change on water balances and hydrologic extremes were also developed by our HRL team.
- Our in-house developed coupled hydrologic-atmospheric model WEHY-HCM was also applied to climate change impact assessment of the hydrologic conditions in Vietnam ( Assessment of climate change effect on the hydrology of a Vietnamese basin.pdf) and Thailand ( Assessment of the Hydrologic Conditions under Climate Change at Ping River Basin, Thailand ).
- Physics-based quantification and assessment of extreme precipitation, floods, droughts and wildfires by means of numerical atmospheric-hydrologic modeling
In this focus area, starting in early 1990s and continuing to the present, we have been studying the physics of droughts, floods, maximum precipitation and wildfires over geographical regions and watersheds in order to understand their underlying causal mechanisms to be able to simulate such systems for their management toward the engineering planning and design of hydraulic structures.
- During the decade of 1990s our HRL (Hydrologic Research Laboratory) team has developed the first continental-scale coupled hydrologic-atmospheric model for the investigation and numerical simulation of droughts over Western USA (On the physics of droughts. I. A conceptual framework. Journal of Hydrology, Volume 129, 1991; On the physics of droughts. II. Analysis and simulation of the interaction of atmospheric and hydrologic processes during droughts. Journal of Hydrology, Volume 129, 1991). Then a Northern Hemisphere-scale energy-balance climate model ( Assessing Hydrologic Drought Risk Using Simplified Climate Model. J. Hydrol. Eng., 5(4). 2000) was developed to relate the effect of atmospheric high-pressure systems over the Pacific Ocean on the Western USA droughts. A summary presentation of this effort is here: On the Physics of Droughts Presentation
- In the last decade the drought conditions at Shasta Lake watershed (the main water resource of California) were investigated by numerical atmospheric-hydrologic modeling (Assessment of the drought conditions at Shasta Dam during 21st Century -presentation; Trinh, K.Ishida, M. L. Kavvas, A. Ercan, K. Carr, Science of the Total Environment, doi: 10.1016/ j.scitotenv.2017. 01.202; 586, 197-205, 2017).
- Starting in the beginning of 2000s and continuing to present, we have been developing new modeling techniques for the estimation of maximum precipitation by means of the regional atmospheric model WRF and applying them to:
- a) Atmospheric River-dominated regions over Western USA (Estimation of maximum precipitation over Feather-Yuba-American Watersheds-presentation-2015 ), and over Southern California (presentation is here: Maximum precipitation and maximum flood estimation over seven California Sierra Foothills Watersheds-Presentation-2019 )
- b) Tropical Cyclone/Hurricane-dominated regions over Eastern USA ( Presentation on the modeling of hurricanes/tropical cyclones in Eastern USA ; Mure-Ravaud, M. L. Kavvas, A. Dib, “Impact of increased atmospheric moisture on the precipitation depth caused by Hurricane Ivan (2004) over a target area”, Science of the Total Environment, 672; 916-926; 2019. https://doi.org/10.1016/j.scitotenv.2019.03.471; Mathieu Mure-Ravaud, M. Levent Kavvas, Alain Dib, “Investigation of intense precipitation from tropical cyclones during the 21st century by dynamical downscaling of CCSM4 RCP 4.5”, International Journal of Environmental Research and Public Health,16(5):687; 2019”.
- c) Maximization of precipitation fields from atmospheric rivers and extratropical cyclones over the 668,000 km2 Columbia River basin: The presentation for this 5-year research project, in collaboration with the US Army Corps of Engineers, is here : Maximum precipitation estimation for the Columbia River Basin Dams-2022
- During the last decade our HRL team has also been working on the development of new, physics-based methods that estimate maximum floods under current and changing climate conditions by means of numerical coupled atmospheric-hydrologic numerical modeling. Two representative presentations of this effort are : Estimation of Extreme Precipitation and Extreme Floods at Southern Sierra Foothills Watersheds of California ; Estimation of Extreme Floods under Changing Climate by Coupled Atmospheric-Hydrologic Modeling
- Recently, we have started to investigate the critical hydro-climatic conditions that result in wildfires over Western USA. The first product of this effort is reported in the paper: Hydroclimatic Conditions for Wildfires over Western USA.
- Numerical modeling of hydrologic and environmental processes at watershed scale by means of WEHY (Watershed Environmental Hydrology) Model for the Assessment of Flow and Environmental Conditions
Since the beginning of 2000s Dr. Kavvas and his graduate students and postdocs at Hydrologic Research Laboratory (HRL) and J. Amorocho Hydraulics Laboratory have been developing numerical modeling tools for simulating the hydrologic and environmental processes in time and space at watershed scale for the assessment of flow, sediment and contaminant (mercury) conditions at various watersheds within California.
- Besides the studies on water balance, floods and droughts, mentioned above, HRL has developed the physically-based, spatially-distributed Watershed Environmental Hydrology (WEHY) model in early 2000s ( Presentation of WEHY (Watershed Environmental Hydrology) Model);
- A few years later an environmental module to WEHY model was developed in order to be able to simulate and assess sediment and contaminant conditions at a watershed ( Environmental module of WEHY model );
- During the last decade the water and sediment inflows from Cache Creek Watershed to the Cache Creek Settling Basin under changing climate conditions were modelled and quantified by means of the hydrologic and environmental components of the WEHY model (a presentation of this effort is here: Estimation of water and sediment inflows from Cache Creek Watershed to Cache Creek Settling Basin under Changing Climate).
- Using the inflow hydrographs and sediment loads, corresponding to various return periods, the flood inundation and sediment bed variations were modelled by CCHE2D, a depth-integrated two-dimensional hydrodynamic and sediment transport model (Jia et al. 2013), over the City of Woodland and the Cache Creek Settling Basin (a presentation of this effort is here: Assessment of Flood Inundation and Sediment Bed Variation at City of Woodland and Cache Creek Settling Basin under Various Return Period Water and Sediment Inflows ).
- Modeling Hydrologic, Hydraulic and Environmental Processes Under Uncertainty
Hydrologic, hydraulic and environmental processes have uncertain parameters as function of spatial location within a specified modeling domain. For example, in soil water flow or in groundwater aquifer flow the hydraulic conductivity, the fundamental parameter of these processes, is uncertain as function of spatial location within a specified model domain. For surface flow processes, the surface roughness, which is a fundamental parameter of these processes, again has uncertain variation with spatial location within a model domain. The source/sink terms for these processes are again uncertain in time and space. For example, the recharge to a groundwater aquifer is an uncertain function of time and space. Similarly, the pumping rates from an aquifer are generally uncertain. Precipitation, which is the main source of overland flow, river flow, soil water flow and unconfined aquifer flow processes, is an uncertain function of time and space over a specified modeling domain. Furthermore, due to the scarcity or lack of data, the initial and boundary conditions of hydrologic-hydraulic-environmental processes are generally uncertain. In the mathematical modeling of these processes one solves their corresponding governing equations under specified initial and boundary conditions in time and space. Within the above framework, this so-called initial-boundary value problem solution must be performed under the appropriately quantified uncertainties with respect to their parameters, sources-sinks, initial values and boundary values. As such, these uncertain systems become stochastic initial-boundary value problems that need to be treated in general as stochastic partial differential equation systems, the governing equation of the process becoming a stochastic partial differential equation. Dr. Kavvas and his students and research associates at HRL and J. Amorocho Hydraulics Laboratory have been developing solutions to such uncertain systems since early 1990s. Below, these efforts will be described.
-
- The most comprehensive description of the time-space evolution of hydrologic-hydraulic-environmental processes is in terms of the time-space evolving probability distribution of the state variables of a particular process. The probability distribution of the state variables of a process describes the so-called ensemble (population) behavior of the process from which various statistics of the process (mean, variance, covariance, etc.) can be obtained. In order develop such probability distributions it is necessary to link the corresponding governing equation of a process to an equation that will describe its corresponding probability distribution which will be linked to the particular process. After about a decade of effort such a one-to-one correspondence between a specified governing equation of a process and the equation which describes the time-space evolution of the probability density of the state variable of the particular process was developed. The resulting equation that describes the evolution of the probability density of the state variables of a particular process is a generalization of the classical local Fokker-Planck-Kolmogorov equation to non-local Eulerian-Lagrangian framework. A presentation of this modeling framework is here: Physics-based modeling of hydrologic-hydraulic-environmental processes under uncertainty- presentation . The corresponding paper is here: Nonlinear stochastic hydrologic processes and their corresponding probability distributions-2003.
- The developed Eulerian-Lagrangian evolution equation for the probability density of any specified process (Kavvas, 2003) was then applied to solute transport in rivers under uncertainty ( Modeling solute transport by river flow under uncertainty-presentation ), modeling snow accumulation and melting under uncertainty (Modeling snow processes under uncertainty- presentation ), modeling of river flow under uncertainty ( Modeling kinematic wave open channel flow under uncertainty-presentation ; Modeling Saint-Venant Unsteady Open Channel Flow Under Uncertainty ), modeling of soil water flow under uncertainty ( Modeling soil water flow under uncertainty ), and modeling groundwater flow under uncertainty (Modeling two-dimensional groundwater flow under uncertainty-presentation ).
- Development of a New Scaling Theory for Modeling the Self-similarity of Hydrologic and Hydraulic Processes in Time and Space
Since the beginning of 1900’s engineers recognized that it is possible to scale the real-life (prototype) flow processes in reduced dimensions in a laboratory (model) in order to be able to investigate the behavior of such processes for the planning and design of hydraulic structures. Such scaling was possible due to Buckingham’s pi–theorem (Buckingham, 1914) which was based on dimensional analysis. By means of this theorem it is possible to scale any flow process, but not its initial and boundary conditions. Nevertheless, it was used successfully in hydraulics and fluid mechanics laboratories for more than a century. In the case of hydrologic processes, while the scaling of the process, analyzed in terms of its governing equation, is essential in order to investigate the process within a framework of different time-space dimensions, the initial and boundary conditions of the process are also very important in determining the behavior of the process. Accordingly, Dr. Kavvas and his team of graduate students and research associates have developed a new scaling theory that is based on mathematician Sophus Lie’s symmetry theory for transformations. This scaling theory does not only develop the scaling relationships for the process, in terms of its governing equation, but also the scaling relationships of its initial and boundary conditions. A simple presentation of this new scaling theory is here: A presentation of the scaling theory as applied to the scaling of hydrologic processes . The introductory paper on the scaling of various hydrologic processes is here: Scale invariance and self-similarity of hydrologic processes . After this initial development, the new scaling theory was applied to various hydrologic and hydraulic processes. For example, this new scaling theory was applied successfully to the scaled modeling of sediment transport while keeping the original material fixed ( Scaling and self-similarity in suspended sediment transport with fixed sediment material ). Development of the scaling relationships of various flow processes culminated in the development of the scaling relationships of the governing Navier-Stokes equations of fluid flow in multi-dimensional space and time ( Scaling relationships and self-similarity of the averaged Navier-Stokes equations ). This new scaling theory, coined as “Lie scaling theory” by our group, can be applied to any flow process in order to develop the scaling relationships for the process.
- Development of the fractional governing equations of hydrologic and hydraulic processes in fractional time and multi-fractional space in order to quantify the effect of long memory in time and space on these processes
The conventional governing equations of hydrologic and hydraulic processes are expressed as partial differential equations with integer power derivatives. As such, these conventional governing equations can only address finite memories in time and space within a process. As shown in this presentation on the fractional governing equations of unsteady open channel flow in fractional time-space: Fractional governing equations of unsteady open channel flow in fractional time-space – presentation , the conventional governing equations of hydrologic and hydraulic processes are essentially local, expressing the nearest-neighbour interactions among the grid cells of a modeling domain. However, as explained in the above presentation, in order to accommodate the long memories in time and/or in space in hydrologic-hydraulic processes it is necessary to have their governing equations as non-local. Such non-local structures for the governing equations of hydrologic-hydraulic processes can be achieved by developing these governing equations as partial differential equations with fractional powers (explained in the above presentation). When their fractional powers revert to integer powers, the developed fractional governing equations transform into the corresponding conventional governing equations of the processes. Besides the above-presented fractional governing equations of unsteady open channel flow, Dr. Kavvas and his team of graduate students and research associates have also developed fractional governing equations of confined groundwater flow ( Fractional governing equations of confined groundwater flow – presentation ; Fractional governing equations of confined groundwater flow in fractional time-space – paper ), unconfined groundwater flow ( Fractional governing equations of unconfined groundwater flow in fractional time-space – paper ) and soil water flow ( Fractional governing equations of soil water flow in fractional time-space – paper ) in fractional time-space. This effort has recently culminated in the development of fractional governing equations of fluid flow which are the generalizations of the Navier-Stokes-Euler governing equations of fluid flow ( Generalizations of the Navier-Stokes-Euler governing equations of fluid flow to fractional time-space- paper ).