Seamless unconfined groundwater heads for coastal California groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (i.e. land surface less than approximately 10 m above mean sea level) areas. In areas where coastal elevations increase rapidly (e.g., bluff stretches), the model boundary was set approximately 1 kilometer inland of the present-day shoreline. Steady-state MODFLOW groundwater flow models were used to obtain detailed (10-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0 to 2 meters (m) in 0.25 m increments, 2.5 m, 3 m, and 5 m) using a range of horizontal hydraulic conductivity (Kh) scenarios (0.1, 1, and 10 m/day). For each SLR/Kh combination, results are provided for two marine boundary conditions, local mean sea level (LMSL) and mean higher-high water (MHHW), and two model versions. In the first model version, groundwater reaching the land surface is removed from the model, simulating loss via natural drainage. In the second model version, groundwater reaching the land surface is retained, simulating the worst-case "linear" response of groundwater head to sea-level rise. Additional details about the groundwater model and data sources are outlined in Befus and others (2020) and in Groundwater_model_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5b8ef008e4b0702d0e7ec72b?name=Groundwater_model_methods.pdf). Methods specific to groundwater head and water table depth products are outlined in Groundwater_head_and_water_table_depth_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bda1563e4b0b3fc5cec39b4?name=Groundwater_head _and_water_table_depth_methods.pdf). Please read the model details, data sources and methods summaries and inspect model output carefully. Data are complete for the information presented. Users should note that while the metadata Spatial Reference Information/UTM Zone Number in this document is 10, some files in southern California are in UTM Zone 11, as noted in the Format Specification for individual downloadable files. As a result users may need to modify the metadata for automated import and display of Zone 11 datafiles.

Seamless unconfined groundwater heads for coastal California groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (i.e. land surface less than approximately 10 m above mean sea level) areas. In areas where coastal elevations increase rapidly (e.g., bluff stretches), the model boundary was set approximately 1 kilometer inland of the present-day shoreline. Steady-state MODFLOW groundwater flow models were used to obtain detailed (10-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0 to 2 meters (m) in 0.25 m increments, 2.5 m, 3 m, and 5 m) using a range of horizontal hydraulic conductivity (Kh) scenarios (0.1, 1, and 10 m/day). For each SLR/Kh combination, results are provided for two marine boundary conditions, local mean sea level (LMSL) and mean higher-high water (MHHW), and two model versions. In the first model version, groundwater reaching the land surface is removed from the model, simulating loss via natural drainage. In the second model version, groundwater reaching the land surface is retained, simulating the worst-case "linear" response of groundwater head to sea-level rise. Modeled groundwater heads were then subtracted from high-resolution topographic digital elevation model (DEM) data to obtain the water table depths, which are represented as polygons for specific depth ranges in this dataset. Additional details about the groundwater model and data sources are outlined in Befus and others (2020) and in Groundwater_model_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5b8ef008e4b0702d0e7ec72b?name=Groundwater_model_methods.pdf). Methods specific to groundwater head and water table depth products are outlined in Groundwater_head_and_water_table_depth_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bda1563e4b0b3fc5cec39b4?name=Groundwater_head _and_water_table_depth_methods.pdf). Methods specific to groundwater emergence and shoaling products are outlined in Groundwater_emergence_and_shoaling_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bd9f318e4b0b3fc5cec20ed?name=Groundwater_emergence_and_shoaling_methods.pdf). Please read the model details, data sources and methods summaries and inspect model output carefully. Data are complete for the information presented. Users should note that while the metadata Spatial Reference Information/UTM Zone Number in this document is 10, some files in southern California are in UTM Zone 11, as noted in the Format Specification for individual downloadable files. As a result users may need to modify the metadata for automated import and display of Zone 11 datafiles.

Seamless unconfined groundwater heads for coastal California groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (i.e. land surface less than approximately 10 m above mean sea level) areas. In areas where coastal elevations increase rapidly (e.g., bluff stretches), the model boundary was set approximately 1 kilometer inland of the present-day shoreline. Steady-state MODFLOW groundwater flow models were used to obtain detailed (10-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0 to 2 meters (m) in 0.25 m increments, 2.5 m, 3 m, and 5 m) using a range of horizontal hydraulic conductivity (Kh) scenarios (0.1, 1, and 10 m/day). For each SLR/Kh combination, results are provided for two marine boundary conditions, local mean sea level (LMSL) and mean higher-high water (MHHW), and two model versions. In the first model version, groundwater reaching the land surface is removed from the model, simulating loss via natural drainage. In the second model version, groundwater reaching the land surface is retained, simulating the worst-case "linear" response of groundwater head to sea-level rise. Modeled groundwater heads were then subtracted from high-resolution topographic digital elevation model (DEM) data to obtain the water table depths, which are represented as polygons for specific depth ranges in this dataset. Additional details about the groundwater model and data sources are outlined in Befus and others (2020) and in Groundwater_model_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5b8ef008e4b0702d0e7ec72b?name=Groundwater_model_methods.pdf). Methods specific to groundwater head and water table depth products are outlined in Groundwater_head_and_water_table_depth_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bda1563e4b0b3fc5cec39b4?name=Groundwater_head _and_water_table_depth_methods.pdf). Methods specific to groundwater emergence and shoaling products are outlined in Groundwater_emergence_and_shoaling_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bd9f318e4b0b3fc5cec20ed?name=Groundwater_emergence_and_shoaling_methods.pdf). Please read the model details, data sources and methods summaries and inspect model output carefully. Data are complete for the information presented. Users should note that while the metadata Spatial Reference Information/UTM Zone Number in this document is 10, some files in southern California are in UTM Zone 11, as noted in the Format Specification for individual downloadable files. As a result users may need to modify the metadata for automated import and display of Zone 11 datafiles.

Seamless unconfined groundwater heads for coastal California groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (i.e. land surface less than approximately 10 m above mean sea level) areas. In areas where coastal elevations increase rapidly (e.g., bluff stretches), the model boundary was set approximately 1 kilometer inland of the present-day shoreline. Steady-state MODFLOW groundwater flow models were used to obtain detailed (10-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0 to 2 meters (m) in 0.25 m increments, 2.5 m, 3 m, and 5 m) using a range of horizontal hydraulic conductivity (Kh) scenarios (0.1, 1, and 10 m/day). For each SLR/Kh combination, results are provided for two marine boundary conditions, local mean sea level (LMSL) and mean higher-high water (MHHW), and two model versions. In the first model version, groundwater reaching the land surface is removed from the model, simulating loss via natural drainage. In the second model version, groundwater reaching the land surface is retained, simulating the worst-case "linear" response of groundwater head to sea-level rise. Modeled groundwater heads were then subtracted from high-resolution topographic digital elevation model (DEM) data to obtain the water table depths. Additional details about the groundwater model and data sources are outlined in Befus and others (2020) and in Groundwater_model_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5b8ef008e4b0702d0e7ec72b?name=Groundwater_model_methods.pdf). Methods specific to groundwater head and water table depth products are outlined in Groundwater_head_and_water_table_depth_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bda1563e4b0b3fc5cec39b4?name=Groundwater_head _and_water_table_depth_methods.pdf). Please read the model details, data sources and methods summaries, and inspect model output carefully. Data are complete for the information presented. Users should note that while the metadata Spatial Reference Information/UTM Zone Number in this document is 10, some files in southern California are in UTM Zone 11, as noted in the Format Specification for individual downloadable files. As a result users may need to modify the metadata for automated import and display of Zone 11 datafiles.

Seamless unconfined groundwater heads for coastal California groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (i.e. land surface less than approximately 10 m above mean sea level) areas. In areas where coastal elevations increase rapidly (e.g., bluff stretches), the model boundary was set approximately 1 kilometer inland of the present-day shoreline. Steady-state MODFLOW groundwater flow models were used to obtain detailed (10-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0 to 2 meters (m) in 0.25 m increments, 2.5 m, 3 m, and 5 m) using a range of horizontal hydraulic conductivity (Kh) scenarios (0.1, 1, and 10 m/day). For each SLR/Kh combination, results are provided for two marine boundary conditions, local mean sea level (LMSL) and mean higher-high water (MHHW), and two model versions. In the first model version, groundwater reaching the land surface is removed from the model, simulating loss via natural drainage. In the second model version, groundwater reaching the land surface is retained, simulating the worst-case "linear" response of groundwater head to sea-level rise. Modeled groundwater heads were then subtracted from high-resolution topographic digital elevation model (DEM) data to obtain the water table depths. Additional details about the groundwater model and data sources are outlined in Befus and others (2020) and in Groundwater_model_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5b8ef008e4b0702d0e7ec72b?name=Groundwater_model_methods.pdf). Methods specific to groundwater head and water table depth products are outlined in Groundwater_head_and_water_table_depth_methods.pdf (available at https://www.sciencebase.gov/catalog/file/get/5bda1563e4b0b3fc5cec39b4?name=Groundwater_head _and_water_table_depth_methods.pdf). Please read the model details, data sources and methods summaries, and inspect model output carefully. Data are complete for the information presented. Users should note that while the metadata Spatial Reference Information/UTM Zone Number in this document is 10, some files in southern California are in UTM Zone 11, as noted in the Format Specification for individual downloadable files. As a result users may need to modify the metadata for automated import and display of Zone 11 datafiles.

Seamless unconfined groundwater heads for U.S. coastal North and South Carolina groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K’s, one with published K’s reduced by a factor of 10 (K*0.1), and one with published K’s increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage.

Seamless unconfined groundwater heads for U.S. coastal North and South Carolina groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K’s, one with published K’s reduced by a factor of 10 (K*0.1), and one with published K’s increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage.

Seamless unconfined groundwater heads for U.S. coastal Virginia, Georgia, and Florida (Atlantic and Gulf coast south of Sarasota) groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K's, one with published K's reduced by a factor of 10 (K*0.1), and one with published K's increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage. Similar modeled data for North Carolina and South Carolina are available from Barnard and others, 2023 at https://doi.org/10.5066/P9W91314.

Seamless unconfined groundwater heads for U.S. coastal Virginia, Georgia, and Florida (Atlantic and Gulf coast south of Sarasota) groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K's, one with published K's reduced by a factor of 10 (K*0.1), and one with published K's increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage. Similar modeled data for North Carolina and South Carolina are available from Barnard and others, 2023 at https://doi.org/10.5066/P9W91314.

Tabular data output from a series of groundwater modeling simulations for five counties along the Central Coast of California, USA. We used a spatially explicit state-and-transition simulation model with stocks and flows that integrates climate, land-use change, human water use, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand at 270-m resolution from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future using output from the Basin Characterization Model, a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows to groundwater from excess irrigation. For each combination of the two climate and two land-use change scenarios, we ran 50 Monte Carlo realizations of the model. Results presented here have been aggregated from the individual cell level and summarized by county.