Você pode jogar com os parâmetros GUC
datestyle e timezone dentro de uma função para obter o que deseja. Aqui está um exemplo (no entanto, ele retorna microssegundos, então provavelmente você precisará ajustá-lo um pouco):create or replace function timestamp_iso8601(ts timestamptz, tz text) returns text as $$
declare
res text;
begin
set datestyle = 'ISO';
perform set_config('timezone', tz, true);
res := ts::timestamptz(3)::text;
reset datestyle;
reset timezone;
return replace(res, ' ', 'T') || ':00';
end;
$$ language plpgsql volatile;
Resultados:
test=# select timestamp_iso8601(now()::timestamptz, 'Europe/Moscow');
timestamp_iso8601
-------------------------------
2017-07-12T08:56:58.692985+03:00
test=# select timestamp_iso8601(now()::timestamptz, 'Pacific/Auckland');
timestamp_iso8601
-------------------------------
2017-07-12T17:59:05.863483+12:00
(1 row)
Atualização:editado. Você pode usar
timestamptz(3) , especificando a precisão (por padrão, ele irá com microssegundos, enquanto 3 manterá apenas milissegundos). Como alternativa, você pode usar res := to_char(ts::timestamptz, 'IYYY-MM-DDT HH24:MI:SS:MSOF'); em vez de ::timestamptz(3)::text cadeia de conversão e, neste caso, (3) não será necessário. 
