## 3 Simulation Results
After addressing the seasonal-trend decomposition for the ammonia cooling system, it should be mentioned that the decomposition residuals average 44 kW and 47 kW in summer and fall, respectively, with a Coefficient of Variation (CV) close to 6%, compared to 19 kW and 21 kW with a CV near 4% in winter and spring. The increase in residuals during the warmer seasons can be attributed to two key factors: (1) higher condensing temperatures, and (2) more frequent fish-handling activities, which result in short but intense bursts of metabolic heat. These results suggest that the residual variability is primarily driven by predictable exogenous events rather than model mis-specification. A similar analysis was performed for the CO2 cooling system, and the decomposition of its power consumption in spring and summer are provided in Figure 5.  
Figure 4: Seasonal-trend decomposition of the power consumption of the ammonia cooling system for autumn and winter. Figure 5: Seasonal-trend decomposition of the power consumption of the CO2 cooling system for spring and summer. The trend component for spring shows that for most hours of the season the cooling system remains in standby mode with minimal power usage, as the breeding pool's water temperature stays below the threshold; however, in the summer greater fluctuations in water temperature lead to more variation in the trend component. Additionally, the high variability in the residual component suggests the influence of other significant factors such as fish metabolism, unexpected events, or increased machinery heat production, making the system's prediction more challenging. Figure 6 presents the same analysis for the following two seasons. It is apparent from Figure 6 that as winter approaches and water temperatures decrease, the power consumption of this unit declines, following a predictable daily cycle; however, the same challenge observed in the previous season persists. High residual values indicate the presence of other influential factors affecting the cooling system's operation. Additionally, similar to the ammonia cooling system in winter, this cooling system operates in standby mode with minimal power consumption due to the lower saltwater temperature. For the CO\({}_{2}\) cascade loop, the residual standard deviation reaches 86 kW in summer and 100 kW in fall, with corresponding CVs of 35% and 52%, respectively. This elevated variability can be explained from two perspectives: first, irregular bursts of metabolic or machinery heat from warm-season activities may contribute to erratic load spikes; second, the trans-critical CO\({}_{2}\) compressor operates near its pressure limit during these periods, making it more sensitive to brief suction-pressure disturbances. In contrast, during winter and spring the system is largely idle and the mean power consumption drops below 10 kW, meaning that even minor startup or defrost pulses result in a CV exceeding 50%. As reviewed in the decomposition analysis, the most prominent factor in the power consumption of the cooling systems associated with breeding pools is the temperature of the water discharged to these pools. In this regard, Figure 7 categorizes the power consumption of the ammonia cooling systems based on water temperature bins. As expected, higher water temperatures correlate with increased power consumption in the cooling system. However, the presence of outliers in certain bins indicates deviations from the expected pattern, suggesting unusual events that warrant further investigation of the recorded data. A similar analysis was conducted for the other cooling system. Figure 8 shows the power consumption distribution across different temperature bins for the CO\({}_{2}\) cooling system. Figure 6: Seasonal-trend decomposition of the power consumption of the CO2 cooling system in fall and winter.