Proliferation of microplastics in rivers, lakes, estuaries, coastal waters and oceans is a major global challenge and threat to the environment, livelihoods and human health. Reliable predictive tools can play an essential role in developing an improved understanding of microplastics behaviour, exposure and risk in water bodies, and facilitate identification of sources and accumulation hot spots, thereby enabling informed decision-making for targeted prevention and clean-up activities. This study presents a new numerical framework (CaMPSim-3D) for predicting microplastics fate and transport in different aquatic settings, which consists of a Lagrangian, three-dimensional (3D) particle-tracking model (PTM) coupled with an Eulerian-based hydrodynamic modeling system (TELEMAC). The 3D PTM has several innovative features that enable accurate simulation and efficient coupling with TELEMAC, which utilizes an unstructured computational mesh. The PTM is capable of considering spatio-temporally varying diffusivity, and uses an innovative algorithm to locate particles within the Eulerian mesh. Model accuracy associated with different advection schemes was verified by comparing numerical predictions to known analytical solutions for several test cases. The implications of choosing different advection schemes for modeling microplastics transport was then investigated by applying the PTM to simulate particle transport in the lower Saint John River Estuary in eastern Canada. The sensitivity of the PTM predictions to the advection scheme was investigated using six numerical schemes with different levels of complexity. Predicted particle distributions and residence times based on the fourth-order Runge-Kutta (RK4) scheme differed significantly (residence times by up to 100 %) from those computed using the traditional first-order (Euler) method. The Third Order Total Variation Diminishing (TVD3) Runge-Kutta method was found to be optimal, providing the closest results to RK4 with approximately 27 % lower computational cost.

A numerical study was conducted to characterize the probability and intensity of storm surge hazards in Canada’s western Arctic. The utility of the European Centre for Medium-Range Weather Forecasts Reanalysis 5th Generation (ERA5) dataset to force numerical simulations of storm surges was explored. Fifty historical storm surge events that were captured on a tide gauge near Tuktoyaktuk, Northwest Territories, were simulated using a two-dimensional (depth-averaged) hydrodynamic model accounting for the influence of sea ice on air-sea momentum transfer. The extent of sea ice and the duration of the ice season has been reducing in the Arctic region, which may contribute to increasing risk from storm surge-driven hazards. Comparisons between winter storm events under present-day ice concentrations and future open-water scenarios revealed that the decline in ice cover has potential to result in storm surges that are up to three times higher. The numerical model was also used to hindcast a significant surge event that was not recorded by the tide gauge, but for which driftwood lines along the coast provided insights to the high-water marks. Compared to measurements at proximate meteorological stations, the ERA5 reanalysis dataset provided reasonable estimates of atmospheric pressure but did not accurately capture peak wind speeds during storm surge events. By adjusting the wind drag coefficients to compensate, reasonably accurate predictions of storm surges were attained for most of the simulated events. The extreme value probability distributions (i.e., return periods and values) of the storm surges were significantly altered when events absent from the tide gauge record were included in the frequency analysis, demonstrating the value of non-conventional data sources, such as driftwood line surveys, in supporting coastal hazard assessments in remote regions.

Storm surges and waves are key climate-driven parameters affecting the design and operation of ports and other infrastructure on the coast. Reliable predictions of future storm surges and waves are not yet available for the west coast of Canada, and this data gap hinders effective climate risk assessment, planning and adaptation. This paper presents numerical simulations of storm surges and waves in British Columbia coastal waters under a future climate (Representative Concentration Pathway) scenario (RCP8.5). The numerical models were first forced by wind and surface pressure fields from the ERA-5 global reanalysis, and calibrated and validated using historical wave and water level records. The models were then driven by atmospheric data from four regional climate models (RCMs) to investigate potential changes in the frequency and magnitude of storm surges and extreme waves over the 21st century. The model outputs were analyzed to determine the potential impacts of climate change on storm surges and wave effects at key ports and transportation assets in western Canada. The study is the first of its kind to utilize unstructured, computational models to simulate storm surges and waves for the entire western Canada coastal region, while maintaining the high spatial resolution in coastal sub-basins needed to capture local dynamic responses.