<p>This study investigated the effects of fluctuating boundary water levels and surface recharge on groundwater flow within unconfined aquifers. We aimed to understand how changes in recharge patterns and variable boundary water levels, such as those from rivers or canals, affect groundwater levels over time and space. To achieve this, we solved the linearized Boussinesq equation using the time-marching method alongside the generalized integral transformation method. Our analysis focused on how different types of recharge affect groundwater level variations and flow dynamics. We found that boundary effects on groundwater level change propagate from the edges toward the aquifer's center, becoming more pronounced with increased boundary water levels. Over time, the system stabilizes, leading to a steady water table height and flow rate, which depend on the disparity between the boundary water levels. Our analytical model demonstrated flexibility and practical applicability by allowing for the consideration or omission of various influencing factors, thus facilitating complete knowledge about groundwater variations and offering future strategic insights for sustainable groundwater resource management.</p>