7+
Ocean Models
8
Core Equations
4
Languages
50+
Model Runs
40+
Yrs of Data

Scientific Programming for Ocean Science

High-performance languages and toolkits for numerical ocean modeling, data analysis, and visualization workflows.

Python

Primary language for oceanographic data processing, model pre/post-processing, and publication-quality visualization of ocean fields and time series.

NumPy SciPy xarray Pandas Matplotlib Cartopy netCDF4 cmocean gsw (TEOS-10) Dask

Julia

High-performance ocean modeling and biogeochemical simulations. Julia's speed rivals Fortran/C while maintaining Python-like syntax, ideal for large-scale eddy-resolving simulations.

Oceananigans.jl OceanBioME.jl Plots.jl DifferentialEquations.jl Distributed.jl

MATLAB / Octave

Classic ocean data analysis, ADCP/CTD processing, spectral analysis, and legacy oceanographic toolboxes that remain standard in the research community.

m_map SeaWater Toolbox ADCP Analysis Signal Processing Time Series

R for Oceanography

Statistical analysis, ecological modeling, and oceanographic visualization in R's rich scientific ecosystem, particularly for biogeochemical and fisheries-linked ocean studies.

oce marmap cmocean ggOceanMaps ncdf4

Ocean & Coastal Hydrodynamic Models

Operational and research-grade ocean circulation models for coastal, regional, and global applications.

ROMS

Regional Ocean Modeling System

Terrain-following sigma-coordinate model widely used for regional and coastal simulations. Supports baroclinic dynamics, nested domains, tidal forcing, and biogeochemical coupling. Applied for Gulf of Mexico circulation, tidal mixing, and sediment transport studies.

Sigma Coordinates Baroclinic Nested Domains Biogeochemistry
FVCOM

Finite Volume Community Ocean Model

Unstructured triangular-grid finite-volume model ideal for estuaries, tidal inlets, and complex coastal geometry. Handles wetting/drying of intertidal flats and incorporates atmosphere-wave-ocean coupling for realistic coastal simulations.

Unstructured Grid Wetting/Drying Estuaries Tidal Channels
ADCIRC

Advanced Circulation Model

High-resolution unstructured-grid model for storm surge, tidal circulation, and coastal flooding. The gold standard for hurricane impact assessment along the Gulf and Atlantic coasts, supporting FEMA flood maps and emergency preparedness planning.

Storm Surge Hurricane Modeling Coastal Flooding FEMA Applications
SCHISM

Semi-implicit Cross-scale Hydroscience Integrated System

Next-generation hybrid V-grid model (unstructured horizontal + terrain-following/z-level vertical) supporting fully baroclinic dynamics, atmosphere-ocean-wave-sediment coupling, and seamless cross-scale simulations from rivers to ocean basins.

Hybrid V-grid Baroclinic Cross-scale Multi-physics
MOM6 / HYCOM

Global & Basin-Scale Ocean Models

Isopycnal-coordinate models used for global ocean reanalysis and forecasting. HYCOM (Hybrid Coordinate Ocean Model) and MOM6 provide open-boundary conditions for regional models and long-term climate simulations with high fidelity in deep-ocean layer representation.

Isopycnal Coords Ocean Reanalysis Global Scale GODAE
WaveWatch III / SWAN

Spectral Wave Models

Third-generation spectral wave models simulating wind-wave generation, propagation, and dissipation. WAVEWATCH III for oceanic scales; SWAN for nearshore wave transformation, shoaling, refraction, and diffraction in coastal zones.

Spectral Waves Wind-Wave Nearshore Wave Shoaling
SWAT / HEC-RAS

Watershed & Riverine Models

Hydrological and hydraulic models for watershed runoff, riverine floodplain inundation, and land-ocean connectivity. Coupled with ADCIRC for compound flooding scenarios integrating riverine discharge, tidal forcing, and storm surge contributions.

Watershed Floodplain Compound Flooding Land-Ocean Coupling

Governing Equations of Ocean Dynamics

The mathematical foundation of numerical ocean modeling — primitive equations, thermodynamics, and transport.

Momentum

Momentum Equations (Horizontal)

The horizontal momentum equations governing large-scale ocean circulation, incorporating the Coriolis force, pressure gradient force, and turbulent friction terms.

Zonal (x) Momentum
$$\frac{Du}{Dt} - fv = -\frac{1}{\rho_0}\frac{\partial p}{\partial x} + \mathcal{F}_u$$
Meridional (y) Momentum
$$\frac{Dv}{Dt} + fu = -\frac{1}{\rho_0}\frac{\partial p}{\partial y} + \mathcal{F}_v$$

where \( f = 2\Omega\sin\phi \) is the Coriolis parameter, \( \rho_0 \) is a reference density, and \( \mathcal{F} \) represents turbulent friction/mixing terms.

Hydrostatic

Hydrostatic Equation

Fundamental vertical force balance in the ocean under the hydrostatic approximation, valid for large horizontal-to-vertical scale ratios typical in ocean basins and continental shelves.

Vertical Pressure Balance
$$\frac{\partial p}{\partial z} = -\rho g$$

This approximation eliminates solving the vertical momentum equation, dramatically reducing computation cost in large-scale ocean models.

Continuity

Continuity Equation (Mass Conservation)

Mass conservation for an incompressible Boussinesq fluid, ensuring the divergence of the velocity field is zero throughout the ocean volume.

Incompressible Continuity
$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$$

In Boussinesq models, density variations are only retained in the buoyancy term; the full mass conservation uses incompressibility.

T & S Conservation

Temperature & Salinity Conservation

Advection-diffusion equations for potential temperature and salinity, the two fundamental tracers determining seawater density and driving thermohaline circulation.

Temperature Conservation
$$\frac{\partial T}{\partial t} + \mathbf{u}\cdot\nabla T = \kappa_T \nabla^2 T + Q_T$$
Salinity Conservation
$$\frac{\partial S}{\partial t} + \mathbf{u}\cdot\nabla S = \kappa_S \nabla^2 S + Q_S$$

\( \kappa_T \) and \( \kappa_S \) are turbulent diffusivities; \( Q_T \) and \( Q_S \) represent surface heat flux and freshwater forcing respectively.

Equation of State

Seawater Equation of State (TEOS-10)

The nonlinear thermodynamic relationship between seawater density, temperature, salinity, and pressure — the closure equation for all ocean models.

Full Nonlinear Form
$$\rho = \rho(T,\, S,\, p)$$
Linear Approximation
$$\rho \approx \rho_0\left[1 - \alpha(T-T_0) + \beta(S-S_0)\right]$$

\( \alpha \) = thermal expansion coefficient, \( \beta \) = haline contraction coefficient. The Gibbs SeaWater (gsw) Python library implements full TEOS-10.

Advection-Diffusion

Advection-Diffusion Equation

General transport equation for any scalar tracer (dissolved oxygen, nutrients, sediment, dye) capturing both advective transport by currents and turbulent mixing.

1-D Scalar Transport
$$\frac{\partial C}{\partial t} + u\frac{\partial C}{\partial x} = K\frac{\partial^2 C}{\partial x^2} + S_C$$
3-D Conservative Form
$$\frac{\partial C}{\partial t} + \nabla\cdot(\mathbf{u}C) = \nabla\cdot(K\nabla C) + S_C$$

\( K \) is the turbulent eddy diffusivity; \( S_C \) is a source/sink term. This equation governs water quality and biogeochemical modeling in coastal systems.

Laplace / Poisson

Laplace & Poisson Equations

Used for the pressure solver in incompressible flow, stream-function computation, and as the classic teaching example in finite-difference methods for numerical oceanography.

Laplace Equation (Steady-State Diffusion)
$$\nabla^2 T = \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} = 0$$
Poisson Equation (Stream Function – Vorticity)
$$\nabla^2 \psi = \zeta$$

Solved numerically via Gauss-Seidel iteration or direct matrix methods. The 2-D heat equation solver is a foundational code exercise in computational ocean science.

Navier-Stokes

Navier-Stokes Equations (Full)

The complete momentum balance for a viscous fluid. Ocean models derive from the Reynolds-averaged Navier-Stokes (RANS) equations with turbulence closure schemes (k-ε, Mellor-Yamada, KPP).

Vector Form
$$\rho\left(\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u}\right) = -\nabla p + \mu\nabla^2\mathbf{u} + \rho\mathbf{g}$$

In ocean modeling, the Coriolis force \( -2\boldsymbol{\Omega}\times\mathbf{u} \) is added to the right-hand side to account for Earth's rotation.

Numerical Methods & Grid Architectures

Discretization schemes, time-stepping methods, and vertical coordinate systems used in operational ocean models.

Finite Difference

Structured grids, simple implementation. Arakawa C-grid staggering for pressure-velocity. Used in ROMS, MOM6, POM.

Finite Element

Unstructured triangular/quad meshes. High flexibility for complex coastal geometry. Used in ADCIRC, SELFE.

Finite Volume

Conservative flux-based formulation. Unstructured grid capability. Used in FVCOM and SCHISM.

Sigma Coordinates

Terrain-following vertical coordinate system. Ideal for shallow coastal flows. Smooth bottom boundary layer.

Z-Level Coordinates

Fixed horizontal layers. Simple for deep ocean. Staircase bottom approximation. Used in MOM and NEMO.

Isopycnal Coordinates

Density-following layers. No diapycnal diffusion. Ideal for basin-scale deep-ocean models (HYCOM, MOM6).

Method Grid Type Time Scheme Key Advantage Example Model
Finite Difference (FD) Structured Leapfrog / AB3 Simple, computationally fast ROMS, MOM6, POM
Finite Element (FE) Unstructured Runge-Kutta / Implicit Handles complex geometry ADCIRC, SELFE
Finite Volume (FV) Unstructured Semi-implicit Flux conservation guaranteed FVCOM, SCHISM
Spectral (SE) Structured Leapfrog High spectral accuracy WaveWatch III

Oceanographic Data Products & Forcing

Observational datasets and reanalysis products used for model initialization, open boundary conditions, and skill validation.

HYCOM Reanalysis / GOFS 3.1

Global Ocean Forecast System at 1/12° resolution providing reanalysis and forecast fields of SSH, temperature, salinity, and currents. Standard open-boundary forcing for regional ROMS and FVCOM configurations.

SSH Temperature Salinity Currents

CMEMS — Copernicus Marine

Comprehensive European ocean data service providing satellite-derived and ensemble reanalysis products: sea level anomaly, SST, chlorophyll, and multi-model ocean forecasts for scientific and operational use.

SLA SST Chl-a Multi-model

Argo Float Profiles

Global array of autonomous profiling floats measuring T/S from the surface down to 2000 m depth. Essential for deep-ocean model validation, mixed layer depth climatologies, and thermohaline structure analysis.

T/S Profiles MLD Deep Ocean 0–2000 m

NOAA CoastWatch / ERDDAP

Satellite-derived coastal ocean products including AVHRR/MODIS SST, ocean color, and wind stress fields. ERDDAP enables standardized programmatic data access for satellite-model matchup studies and validation.

AVHRR SST MODIS Ocean Color ERDDAP API

World Ocean Database (WOD)

NCEI's global in-situ oceanographic database spanning 1770–present with CTD, XBT, Argo, and mooring data. Used for climatological T/S fields (WOA23), model initialization, and long-term Gulf of Mexico trend analysis.

WOA23 CTD Climatology 1770–present

NOAA Tide Gauges & ADCP

Real-time and historical tidal elevation records from NOAA CO-OPS network, combined with ADCP current profiles for tidal harmonic analysis, model skill assessment (RMSE, bias, correlation), and storm surge hindcast validation.

Tidal Harmonics Sea Level Current Profiles Model Skill

Research Projects

Applied numerical oceanography projects combining model simulations, data analysis, and coastal management applications.

2024 – 2027

Mississippi Sound Hydrodynamics & Water Quality

Coupled hydrodynamic–water quality modeling of the Mississippi Sound using FVCOM and CE-QUAL-ICM. Investigating freshwater discharge influence on harmful algal bloom dynamics, seasonal hypoxia, and salinity stratification variability. NSF-funded doctoral research at Mississippi State University.

Mississippi Sound water quality time series
FVCOM Water Quality Hypoxia Mississippi Sound
2022 – 2023

Floodplain Connectivity — Lower Mississippi Alluvial Valley

HEC-RAS and SWAT-coupled floodplain inundation modeling to quantify hydrological connectivity in the Lower Mississippi Alluvial Valley. Assessed habitat restoration scenarios for bottomland hardwood ecosystems and migration corridors for fish and waterfowl.

Floodplain connectivity analysis LMAV
HEC-RAS SWAT Floodplain Connectivity
Gulf of Mexico

ROMS Simulation: Mesoscale Eddies & Loop Current

Regional ROMS simulation of Gulf of Mexico seasonal circulation, focusing on Loop Current dynamics, mesoscale eddy formation, and sea surface temperature anomalies. Validated against satellite altimetry (CMEMS) and HYCOM reanalysis products.

Gulf of Mexico mesoscale eddies simulation
ROMS Mesoscale Eddies Loop Current Gulf of Mexico
Bay of Bengal

Sea Surface Temperature Trend Analysis

Long-term SST trend analysis for the Bay of Bengal using AVHRR, MODIS, and in-situ CTD observations. Identified warming signatures (>10°C seasonal anomalies) linked to altered monsoon-driven circulation and implications for tropical cyclone intensification.

Bay of Bengal SST anomaly analysis
SST Analysis Bay of Bengal AVHRR/MODIS Trend Detection
Ongoing

Python Computational Oceanography Tutorials

Open-source educational materials for computational oceanography: 2-D Laplace heat equation solver (Gauss-Seidel), advection schemes stability analysis, tidal harmonic analysis with t_tide, HYCOM API data access workflows, and xarray-based ocean climatology processing. Published on GitHub and Medium.

Open Source Education Python GitHub

Selected Publications & Reports

Peer-reviewed journal articles and technical reports in numerical and coastal oceanography.

In Preparation

Coupled Hydrodynamic–Water Quality Modeling of Mississippi Sound

Manuscript describing FVCOM-based simulation of coastal circulation and dissolved oxygen dynamics in the Mississippi Sound under varying freshwater discharge and wind forcing regimes. Target journal: Estuarine, Coastal and Shelf Science.

FVCOM Water Quality Mississippi Sound
2022

Hydrological Connectivity in the Lower Mississippi Alluvial Valley

Floodplain inundation frequency and habitat connectivity analysis using HEC-RAS coupled with multi-temporal remote sensing, informing conservation planning for migratory fish corridors and waterbird breeding habitat in the LMAV.

HEC-RAS Connectivity Remote Sensing

Explore the Code

Ocean model scripts, Python notebooks for data analysis, equation solvers, and HYCOM data-access tools are available on GitHub.

View GitHub Repository