If I use a stateful RNN in Keras for processing a sequence of length N divided into N parts (each time step is processed individually),
The back propagation horizon is limited to the second dimension of the input sequence. i.e. if your data is of type (num_sequences, num_time_steps_per_seq, data_dim)
then back prop is done over a time horizon of value num_time_steps_per_seq
Take a look at