Define some surface data

In[]:=
bbox={{438400,362700},{439200,363300}};
In[]:=
lidar=
LiDAR
;

Plot the surface data

Use Automatic scaling in the x, y and z directions
ListPlot3D[lidar["Data"],
BoxRatiosAutomatic
,DataRangelidar["DataRange"],ImageSize->Large]

Request the values used for Automatic BoxRatios

In[]:=
AbsoluteOptions[​​ListPlot3D[lidar["Data"],BoxRatios->Automatic,DataRangelidar["DataRange"],ImageSize->Large],​​BoxRatios​​]

In V12.3, I used to get...

In[]:=
{BoxRatios{19.328336901028575`,14.496252675771432`,1.`}};

In V13.0, I now get...

{BoxRatiosAutomatic};
...which is no use at all.

Prior to V13, I could use the absolute value of BoxRatios to set an exaggeration factor for the vertical scale, as follows.

In[]:=
Manipulate[​​ListPlot3D[lidar["Data"],BoxRatios{19.328336901028575`,14.496252675771432`,factor},DataRangelidar["DataRange"],ImageSize->Large],​​{{factor,1},1,10}​​]