The present study proposes a novel method that merges satellite ocean color bio-optical products with Argo temperature-salinity profiles to infer the vertical distribution of the particulate backscattering coefficient (bbp). This neural network-based method (SOCA-BBP for Satellite Ocean-Color merged with Argo data to infer the vertical distribution of the Particulate Backscattering coefficient) uses three main input components: (1) satellite-based surface estimates of bbp and chlorophyll a concentration matched up in space and time with (2) depth-resolved physical properties derived from temperature-salinity profiles measured by Argo profiling floats and (3) the day of the year of the considered satellite-Argo matchup. The neural network is trained and validated using a database including 4725 simultaneous profiles of temperature-salinity and bio-optical properties collected by Bio-Argo floats, with concomitant satellite-derived products. The Bio-Argo profiles are representative of the global open-ocean in terms of oceanographic conditions, making the proposed method applicable to most open-ocean environments. SOCA-BBP is validated using 20% of the entire database (global error of 21%). We present additional validation results based on two other independent data sets acquired (1) by four Bio-Argo floats deployed in major oceanic basins, not represented in the database used to train the method; and (2) during an AMT (Atlantic Meridional Transect) field cruise in 2009. These validation tests based on two fully independent data sets indicate the robustness of the predicted vertical distribution of bbp. To illustrate the potential of the method, we merged monthly climatological Argo profiles with ocean color products to produce a depth-resolved climatology of bbp for the global ocean.