From spectrogram one can deduce that in addition to being man made,
- it is for short range (attenuation of 40 kHz signals in water is 10dB/km)
- there are 15-16 distinct bands of independent signals (no harmonics)
- waveforms indicate constant Time-Bandwidth (and constant Q) transmissions aiming at constant signal processing gain
- individual sweeps show volume reverberation, so it is generated at some distance
- no clear multi-path is visible, so signal was most likely generated close to surface or sensor was close to surface
Without a more comprehensive dataset and recoding context it is difficult to speculate on origin and purpose.
Edit on processing gain: It seems to me that upsweeps at lower frequencies that have narrower bandwidth are also somewhat longer that the signals at higher frequencies that have larger bandwidth. If the receiver is match-filtering then the processing gain is determined by the product of duration and bandwidth (time-bandwidth product), which could be close to constant. The increase of bandwidth with frequencies would also indicate a constant Q filtering. Obviously, a more detailed analysis would be facilitated with access to the timeseries. Definitely, a very interesting signal (and in case it was a joke made to you by some technicians it is a very intelligent one).
Edit on multi-path: There are double signals that could be confused with multi-path. However, as they time separation corresponds to the duration of the sweep, it is very unlikely that they are multiple arrivals. More likely they are part of the code. Only reprocessing of the time series with proper matched filters would reveal the medium transfer function (multi-path structure)