На основе материалов донных тралений с 1977 по 2021 г., проведённых до глубины 400 м, показана высокая статистическая значимость придонной температуры воды T и глубины места D для моделирования распределения биомассы трески в Западно-Беринговоморской зоне во всех проверенных моделях. Наилучшей обобщающей способностью обладали векторные авторегрессионные пространственно-временные (VAST) модели, включающие нелинейные зависимости уловов трески от T и D. Корреляция плотностей трески в тестовом наборе данных с оценками плотностей в моделях VAST была выше, чем таковая с оценками биомасс из более простых моделей, настроенных на полном наборе данных. Использование моделей VAST позволяет получить непрерывные временные ряды биомассы трески с оценкой их неопределенности и статистических весов самих моделей относительно тестовых данных. Полученный усредненный ряд динамики биомассы ансамблевым методом с учетом статистических весов моделей совместно с ранее опубликованными оценками биомасс позволяют установить в обобщенной модели прибавочной продукции в пространстве состояний с Байесовым подходом динамику отклонений биологических процессов от стационарных допущений и приблизительно оценить объем трески, не учитываемый этими процессами. Доля такой трески резко выросла в 2016 г. до 40 % и к 2018 г. достигла максимума в 49 %, что значительно отклоняется от стационарных допущений, но затем эта доля начала снижаться. Анализ годовых тенденций из эмпирических ортогональных функций T выявил резкие изменения основных мод T в эти годы. Таким образом, неоднократно высказанная гипотеза о перераспределении трески Берингова моря из-за изменения площади акватории с низкой температурой воды у дна здесь впервые проверена статистическими методами в пространстве. В связи с найденной высокой ошибкой биологических процессов сделан вывод о невозможности точного прогнозирования динамики биомассы трески без прогнозирования распределения придонной температуры воды в пространстве.