Accurate prediction of tuna distribution is essential for sustainable fisheries management. This study develops a two-stage hybrid model combining Convolutional Neural Networks (CNN), Recurrent Neural Networks (RNN), and Random Forest (RF) to predict tuna distribution around drifting fish aggregating devices (DFAD) in the Western and Central Pacific Ocean (WCPO). Echo-sounder buoy data from DFAD were aggregated into 2° × 2° grid cells and matched with oceanographic variables from the Copernicus Marine Service. Random Forest-based variable importance analysis identified primary productivity (27%), chlorophyll-a (22%), and dissolved oxygen (18%) as the three dominant environmental drivers. The CNN-RNN component extracts spatiotemporal features from multi-layer ocean data, while the RF classifier performs binary classification of tuna aggregation zones (high-yield vs. low-yield). All five models (Decision Tree, RF, CNN, Transformer, and CNN-RNN-RF) were evaluated on 557 samples using 5-fold stratified cross-validation, with each fold further split 80:20 for training and validation. The proposed CNN-RNN-RF model achieved the highest performance with an AUC of 0.830, accuracy of 82.6%, and F1-scores of 86.3% (high-yield) and 76.2% (low-yield), outperforming the best baseline model (RF: AUC 0.761, accuracy 75.4%). Predicted high-yield zones showed strong consistency with fishing log records, demonstrating the potential of integrating echo-sounder data with hybrid deep learning for data-driven tuna fisheries management.