The Coastal Plain in the southern Chesapeake Bay area is becoming increasingly susceptible to nuisance flooding as a result of the combination of sea-level rise and land subsidence associated with aquifer compaction from excessive groundwater pumping. Detailed time-series of cumulative compaction data (land subsidence) from the three U.S. Geological Survey deployed extensometers in the regions, along with cyclical piezometer data, reflect the nature of the complex multi-aquifer/aquitard system in the Coastal Plain. Franklin, Virginia and Suffolk, Virginia extensometers were deactivated in 1995 and were reactivated in 2016 along with the addition of a high-sensitivity borehole extensometer in Nansemond, Virginia in collaboration with the Hampton Roads Sanitation District as a part of the Sustainable Water Initiative for Tomorrow (SWIFT). Yearly compaction rates estimated from the reactivated extensometers are -3.3 mm/year, 15.6 mm/year, and -20.7 mm/year in Franklin, Suffolk, and Nansemond, Virginia respectively. One-dimensional vertical compaction modeling is utilized to estimate the total compaction and differentiate which fine-grained confining units or aquifer interbeds are contributing most to total compaction historically and presently. Additionally, properties of the system can be estimated including the elastic specific storage of the aquitards and aquifers and the inelastic storage of the aquitards. The total cumulative change in aquifer system thickness estimated by the MODFLOW subsidence package can be compared to the observed total cumulative change in aquifer system thickness at each site for validation of hypothesis about the dynamics of the aquifer system to known changes in stress. Subsidence rates and aquifer/aquitard properties can be useful for managing and modeling the groundwater in the Coastal Plain of Virginia. / Master of Science / The Coastal Plain in the southern Chesapeake Bay area is becoming increasingly susceptible to flooding at high tides in low lying areas as a result of the combination of sea-level rise and sinking of the land surface (land subsidence) associated with aquifer compaction from excessive groundwater pumping from buried aquifers. Detailed time-series of land subsidence data from the three U.S. Geological Survey deployed extensometers in the region, along with water level data from nearby wells, reflect the nature of the complex multi-aquifer/aquitard system in the Coastal Plain. Franklin, Virginia and Suffolk, Virginia extensometers were deactivated in 1995 and were reactivated in 2016 along with the addition of a high-sensitivity borehole extensometer in Nansemond, Virginia in collaboration with the Hampton Roads Sanitation District as a part of the Sustainable Water Initiative for Tomorrow (SWIFT). Yearly land subsidence rates estimated from the reactivated extensometers are -3.3 mm/year, 15.6 mm/year, and -20.7 mm/year in Franklin, Suffolk, and Nansemond, Virginia respectively. One-dimensional vertical compaction modeling is utilized to estimate the total sinking of the land surface as well as to differentiate which fine-grained confining units or aquifer interbeds are contributing most to total subsidence historically and presently. Additionally, properties of the system can be estimated including the elastic specific storage of the aquitards and aquifers and the inelastic storage of the aquitards. The total cumulative change in aquifer system thickness estimated by the MODFLOW subsidence package can be compared to the observed total cumulative change in aquifer system thickness at each site for validation of hypothesis about the dynamic changes of the aquifer system with known changes in stress. Subsidence rates, understanding the dynamics of the aquifer system, and aquifer/aquitard properties can be useful for managing groundwater and modeling the aquifer system in the Coastal Plain of Virginia.

The thermal response of a borehole field is often described by non‐dimensional response factors called gfunctions.The g‐function was firstly generated as a numerical solution based on SBM (Superposition BoreholeModel). An analytical approach, the FLS (Finite Line Source), is also accepted for generating the g‐function. In thiswork the potential to numerically produce g‐functions is studied for circular borehole fields using the commercialsoftware COMSOL. The numerical method is flexible and allows the generation of g‐functions for any boreholefield geometry. The approach is partially validated by comparing the solution for a square borehole field containing36 boreholes (6x6) with g‐functions generated with the FLS approach and with the program EED (Earth EnergyDesigner). The latter is based on Eskilsons SBM, one of the first documents where the concept of g‐functions wasintroduced. Once the approach is validated, the square COMSOL model is compared with a circular geometryborehole field developed by the same method, consisting of 3 concentric rings having 6, 12, and 18 boreholes.Finally the influence on the circular geometry g‐function is studied when connecting the boreholes in radial zoneswith different thermal loads. / Den termiska responsen för ett borrhålsfält beskrivs ofta med den dimensionslösa responsfunktionen kallad gfunktion.Responsfunktionen togs först fram som en numerisk lösning med SBM (Superposition Borehole Model).En analytisk metod, FLS (Finite Line Source) är också accepterad för framtagandet av g‐funktioner. I det här arbetetundersöks förutsättningarna att numeriskt ta fram g‐funktioner för cirkulära borrhålsfält genom att använda detkommersiella simuleringsprogrammet COMSOL Multiphysics. Den numeriska metoden är flexibel och kananvändas för alla typer av borrhålsgeometrier. Metoden att använda COMSOL valideras delvis genom att jämföraresultatet för ett kvadratiskt borrhålsfält innehållande 36 borrhål (6x6) med lösningar framtagna med FLS och meddimensioneringsprogrammet EED (Earth Energy Designer). Det senare har sin grund i Eskilsons SBM, ett av deförsta arbeten där begreppet g‐funktion introducerades. När metoden att använda COMSOL verifierats, jämförsden kvadratiska borrhålsmodellen med en cirkulär borrhålskonfiguration, upprättad med samma metod,innehållande 3 koncentriska ringar om vardera 6, 12, 18 borrhål. Slutligen undersöks hur den termiska responsenpåverkas då borrhålen i ett cirkulärt borrhålsfält kopplas samman och grupperas i radiella zoner med olika termiskalaster. / SEEC Scandinavian Energy Efficiency Co.