Remote sensing-based field-scale evapotranspiration (ET) maps are useful for characterizing water use patterns and assessing crop performance. The relative impact of climate variability and water management decisions are better studied and quantified using historical data that are derived using a set of consistent datasets and methodology. Historical (1984–2014) Landsat-based ET maps were generated for major irrigation districts in California, i.e., Palo Verde and eight other sub-basins in parts of the middle and lower Central Valley. A total of 3396 Landsat images were processed using the Operational Simplified Surface Energy Balance (SSEBop) model that integrates weather and remotely sensed images to estimate monthly and annual ET within the study sites over the 31years. Model output evaluation and validation using gridded-flux data and water balance ET approaches indicated relatively good correspondence (R2 up to 0.88, root mean square error as low as 14mm/month) between SSEBop ET and validation datasets. In a pairwise comparison, annual variability of agro-hydrologic parameters of actual evapotranspiration (ETa), land surface temperature (Ts), and runoff (Q) were found to be more variable than their corresponding climatic counterparts of atmospheric water demand (ETo), air temperature (Ta), and precipitation (P), revealing process differences between regional climatic drivers and localized agro-hydrologic responses. However, only Ta showed a consistent increase (up to 1.2K) over study sites during the 31years, whereas other climate variables such as ETo and P showed a generally neutral trend. This study demonstrates a useful application of “Big Data” science where large volumes of historical Landsat and weather datasets were used to quantify and understand the relative importance of water management and climate variability in crop water use dynamics in regards to the linkages among water management decisions, hydrologic processes and economic transactions. Irrigation district-wide ETa estimates were used to compute historical crop water use volumes and monetary equivalents of water savings for the Palo Verde Irrigation District (PVID). During the peak crop fallowing year in PVID, the water saved reached a maximum of ~107,200acre-feet in 2011 with an estimated monetary payout value of $20.5 million. A significant decreasing trend in actual ET despite an increasing atmospheric demand in PVID highlights the role of management decisions in affecting local hydrologic processes. This study has importance for planning water resource allocation, managing water rights, sustaining agricultural production, and quantifying impacts of climate and land use/land cover changes on water resources. With increased computational efficiency, similar studies can be conducted in other parts of the world to help policy and decision makers understand and quantify various aspects of water resources management.