Combining both processing techniques of horizontal-to-vertical spectral ratio (HVSR) and surface-to-borehole spectral ratio (SBR), using the KiK-net and K-NET database in Japan, could be used in the present study to establish relationships of Vs, Vp, and resonant frequency (f) versus depth to bedrock half-space (h) based on the site-dependent variability in different lithologies. Remarkable correlations of the average Vs and Vp of layers overlying the bedrock half-space (i.e. Vs and Vp) versus h are inversely resembling the relationships of f versus h and their maxima of ≤ 1000 m/s and ≤ 3000 m/s, respectively. Moreover, their ranges are decreasing gradually through site class of B, C, D and E resembling the ranges of the Vs30 according to NEHRP Provision (2000). Underestimation in the correlation coefficients resulted from HVSR at K-NET sites due to shallow and limited ranges of h whenever compared with those resulted from HVSR and SBR at KiK-net sites. Rebuilding the f versus h based on site-dependent variability in different lithologies could strongly decrease scattering and deviations of data points. Although such relationships are site-specific and highly dependent on each region’s geologic conditions, fair comparisons based on site information and site-dependent variability of different bedrock lithologies between previous relationships of f versus h in the literature and the present study relationships are showing remarkable and reasonable similarities. This indicates significant importance of introducing the seismic site classification of different lithologies as a crucial controlling factor in establishing the previous and the present nonlinear regression relationships.