Este artículo analiza numéricamente la distribución de los ceros de la función zeta de Riemann a lo largo de la línea crítica (CL). Los ceros se distribuyen de acuerdo con un modelo jerárquico de dos capas, uno determinista y el otro estocástico. Siguiendo una anamorfosis del plano complejo que involucra la función de Lambert, la distribución de ceros a lo largo de la CL transformada sigue la realización de un proceso estocástico de variables aleatorias gaussianas independientes espaciadas regularmente, cada una vinculada a un cero. El valor de la desviación estándar permite la posible superposición de realizaciones adyacentes de las variables aleatorias, sobre un intervalo de confianza estrecho. El modelo jerárquico divide la función en clases de equivalencia secuenciales, con el rango de densidades de probabilidad de las realizaciones coincidiendo con el espectro de estilos de comportamiento de las clases. El modelo tiene como objetivo expresar, en la CL, las coordenadas de las cancelaciones alternantes de las partes real e imaginaria de la función, diseccionar la fórmula para el número de ceros por debajo de un umbral, estimar las leyes estadísticas de dos ceros consecutivos, de máximos de la función y momentos. Esto también ayuda a explicar la ausencia de raíces múltiples.