The bulk modulus notebook takes the absolute value of r^2 when computing the scoring metric. But, it should not do that - if r^2 is less than zero, it is worse than predicting the mean.
Probably it won't affect this specific notebook, but better not to have that code in there. Very simple fix - just remove np.abs() around any r2 computation.