Measurements of tides are relatively sparse in the Arctic. This paper studies GNSS buoy tracks to complement existing data. Existing methods to perform tidal harmonic analysis of the buoy data are inadequate in the Arctic region because these methods for tidal analysis combine data from multiple buoy tracks, which is often infeasible in the Arctic. Moreover, we find that there are significant spatial and temporal variations in amplitudes and phases in baroclinic zones. To address these complexities, we introduce a new approach–Model-derived Fitting Method–to estimate the tidal current constituents (TCC) from a single buoy trajectory. Our study assesses the proposed method by analyzing GNSS buoy data from three Arctic regions characterized by barotropic or baroclinic tidal currents. Through detailed case studies in the Barents Sea, Chukchi Sea, and Baffin Bay, our approach demonstrates accuracy, robustness, and operational capabilities. In the Barents Sea, TCC estimates from two buoys were compared at a common location within their trajectories and compared against model estimates. In the Chukchi Sea's barotropic dominant zone, our method's estimates were evaluated against nearby ADCP mooring data. In Baffin Bay, known for baroclinic currents, a synthetic evaluation confirmed the method's effectiveness. Our study also highlights that phase variations along buoy trajectories can lead to frequency shifts in the spectrum, similar to the Doppler shift effect, particularly notable in regions with baroclinic tides.